1
|
Calamari ZT, Song A, Cohen E, Akter M, Roy RD, Hallikas O, Christensen MM, Li P, Marangoni P, Jernvall J, Klein OD. Vole genomics links determinate and indeterminate growth of teeth. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2023.12.18.572015. [PMID: 38187646 PMCID: PMC10769287 DOI: 10.1101/2023.12.18.572015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/09/2024]
Abstract
Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells responsible for tooth crown growth are lost at the onset of tooth root formation. Genetic signaling that initiates this loss is difficult to study with the ever-growing incisor and rooted molars of mice, the most common mammalian dental model species, because signals for root formation overlap with signals that pattern tooth size and shape (i.e., cusp patterns). Different species of voles (Cricetidae, Rodentia, Glires) have evolved rooted and unrooted molars that have similar size and shape, providing alternative models for studying roots. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes. We identified 15 dental genes with changing synteny relationships and six dental genes undergoing positive selection across Glires, two of which were undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Decreased expression of both genes in prairie voles with unrooted molars compared to bank voles supports the presence of positive selection and may underlie differences in root formation. Bulk transcriptomics analyses of embryonic molar development in bank voles also demonstrated conserved patterns of dental gene expression compared to mice, with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars. Our results support ongoing evolution of dental genes across Glires, revealing the complex evolutionary background of convergent evolution for ever-growing molars.
Collapse
Affiliation(s)
- Zachary T. Calamari
- Baruch College, City University of New York, One Bernard Baruch Way, New York, NY 10010, USA
- The Graduate Center, City University of New York, 365 Fifth Ave, New York, NY 10016, USA
- Program in Craniofacial Biology and Department of Orofacial Sciences, University of California, San Francisco, San Francisco, CA 94158, USA
- Division of Paleontology, American Museum of Natural History, Central Park West at 79th Street, New York, NY, 10024, USA
| | - Andrew Song
- Baruch College, City University of New York, One Bernard Baruch Way, New York, NY 10010, USA
- Cornell University, 616 Thurston Ave, Ithaca, NY 14853, USA
| | - Emily Cohen
- Baruch College, City University of New York, One Bernard Baruch Way, New York, NY 10010, USA
- New York University College of Dentistry, 345 E 34th St, New York, NY 10010
| | - Muspika Akter
- Baruch College, City University of New York, One Bernard Baruch Way, New York, NY 10010, USA
| | - Rishi Das Roy
- Institute of Biotechnology, University of Helsinki, FI-00014 Helsinki, Finland
| | - Outi Hallikas
- Institute of Biotechnology, University of Helsinki, FI-00014 Helsinki, Finland
| | - Mona M. Christensen
- Institute of Biotechnology, University of Helsinki, FI-00014 Helsinki, Finland
| | - Pengyang Li
- Program in Craniofacial Biology and Department of Orofacial Sciences, University of California, San Francisco, San Francisco, CA 94158, USA
- Department of Pediatrics, Cedars-Sinai Guerin Children’s, 8700 Beverly Blvd., Suite 2416, Los Angeles, CA 90048, USA
| | - Pauline Marangoni
- Program in Craniofacial Biology and Department of Orofacial Sciences, University of California, San Francisco, San Francisco, CA 94158, USA
- Department of Pediatrics, Cedars-Sinai Guerin Children’s, 8700 Beverly Blvd., Suite 2416, Los Angeles, CA 90048, USA
| | - Jukka Jernvall
- Institute of Biotechnology, University of Helsinki, FI-00014 Helsinki, Finland
- Department of Geosciences and Geography, University of Helsinki, FI-00014 Helsinki, Finland
| | - Ophir D. Klein
- Program in Craniofacial Biology and Department of Orofacial Sciences, University of California, San Francisco, San Francisco, CA 94158, USA
- Department of Pediatrics, Cedars-Sinai Guerin Children’s, 8700 Beverly Blvd., Suite 2416, Los Angeles, CA 90048, USA
| |
Collapse
|
2
|
Vasquez YM, Li Z, Xue AZ, Bennett GM. Chromosome-level genome assembly of the aster leafhopper (Macrosteles quadrilineatus) reveals the role of environment and microbial symbiosis in shaping pest insect genome evolution. Mol Ecol Resour 2024; 24:e13919. [PMID: 38146900 DOI: 10.1111/1755-0998.13919] [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: 08/22/2023] [Revised: 11/12/2023] [Accepted: 12/13/2023] [Indexed: 12/27/2023]
Abstract
Leafhoppers comprise over 20,000 plant-sap feeding species, many of which are important agricultural pests. Most species rely on two ancestral bacterial symbionts, Sulcia and Nasuia, for essential nutrition lacking in their phloem and xylem plant sap diets. To understand how pest leafhopper genomes evolve and are shaped by microbial symbioses, we completed a chromosomal-level assembly of the aster leafhopper's genome (ALF; Macrosteles quadrilineatus). We compared ALF's genome to three other pest leafhoppers, Nephotettix cincticeps, Homalodisca vitripennis, and Empoasca onukii, which have distinct ecologies and symbiotic relationships. Despite diverging ~155 million years ago, leafhoppers have high levels of chromosomal synteny and gene family conservation. Conserved genes include those involved in plant chemical detoxification, resistance to various insecticides, and defence against environmental stress. Positive selection acting upon these genes further points to ongoing adaptive evolution in response to agricultural environments. In relation to leafhoppers' general dependence on symbionts, species that retain the ancestral symbiont, Sulcia, displayed gene enrichment of metabolic processes in their genomes. Leafhoppers with both Sulcia and its ancient partner, Nasuia, showed genomic enrichment in genes related to microbial population regulation and immune responses. Finally, horizontally transferred genes (HTGs) associated with symbiont support of Sulcia and Nasuia are only observed in leafhoppers that maintain symbionts. In contrast, HTGs involved in non-symbiotic functions are conserved across all species. The high-quality ALF genome provides deep insights into how host ecology and symbioses shape genome evolution and a wealth of genetic resources for pest control targets.
Collapse
Affiliation(s)
- Yumary M Vasquez
- Department of Life and Environmental Sciences, University of California, Merced, Merced, California, USA
| | - Zheng Li
- Department of Integrative Biology, University of Texas at Austin, Austin, Texas, USA
| | - Allen Z Xue
- Department of Integrative Biology, University of Texas at Austin, Austin, Texas, USA
| | - Gordon M Bennett
- Department of Life and Environmental Sciences, University of California, Merced, Merced, California, USA
| |
Collapse
|
3
|
Liu Q, Zhao RM, Wang DY, Li P, Qu YF, Ji X. Genome-wide characterization of the TGF-β gene family and their expression in different tissues during tail regeneration in the Schlegel's Japanese gecko Gekko japonicus. Int J Biol Macromol 2024; 255:128127. [PMID: 37984573 DOI: 10.1016/j.ijbiomac.2023.128127] [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] [Received: 05/12/2023] [Revised: 10/19/2023] [Accepted: 11/08/2023] [Indexed: 11/22/2023]
Abstract
The transforming growth factor-β (TGF-β) gene family is unique to animals and is involved in various important processes including tissue regeneration. Here, we identified 52 TGF-β family genes based on genome sequences of the gecko (Gekko japonicus), compared TGF-β genes between G. japonicus and other four reptilian species, and evaluated the expression of 14 randomly selected genes in muscle, kidney, liver, heart, and brain during tail regeneration to investigate whether their expression was tissue-dependent. We detected 23 conserved domains, 13 in the TGF-β ligand subfamily, and 10 in the receptor subfamily. The pattern of higher genetic variation in the ligand subfamily than in the receptor subfamily in vertebrates might result from the precise localization of agonists and antagonists in the cell surface and intracellular compartment. TGF-β genes were unevenly distributed across 15 chromosomes in G. japonicus, presumably resulting from gene losses and gains during evolution. Genes in the TGF-β receptor subfamily (ACVR2A, ACVR2B, ACVR1, BMPR1A, ACVRL1, BMPR2 and TGFBR1) played a vital role in the TGF-β signal pathway. The expression of all 14 randomly selected TGF-β genes was tissue-specific. Our study supports the speculation that some TGF-β family genes are involved in the early stages of tail regeneration.
Collapse
Affiliation(s)
- Qian Liu
- College of Life Sciences, Nanjing Normal University, Nanjing 210023, China
| | - Ru-Meng Zhao
- College of Life Sciences, Nanjing Normal University, Nanjing 210023, China
| | - Dan-Yan Wang
- College of Life Sciences, Nanjing Normal University, Nanjing 210023, China
| | - Peng Li
- College of Life Sciences, Nanjing Normal University, Nanjing 210023, China
| | - Yan-Fu Qu
- College of Life Sciences, Nanjing Normal University, Nanjing 210023, China.
| | - Xiang Ji
- Zhejiang Provincial Key Laboratory for Water Environment and Marine Biological Resources Protection, College of Life and Environmental Sciences, Wenzhou University, Wenzhou 325035, China.
| |
Collapse
|
4
|
Habib MT, Rahman S, Afrad MH, Howlader AM, Khan MH, Khanam F, Alam AN, Chowdhury EK, Rahman Z, Rahman M, Shirin T, Qadri F. Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh. Front Genet 2023; 14:1220906. [PMID: 37621704 PMCID: PMC10446972 DOI: 10.3389/fgene.2023.1220906] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/11/2023] [Accepted: 07/20/2023] [Indexed: 08/26/2023] Open
Abstract
Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.
Collapse
Affiliation(s)
| | - Saikt Rahman
- Institute for Developing Science and Health Initiatives, Dhaka, Bangladesh
| | | | | | | | - Farhana Khanam
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | - Ahmed Nawsher Alam
- Institute of Epidemiology, Disease Control and Research, Dhaka, Bangladesh
| | - Emran Kabir Chowdhury
- Department of Biochemistry and Molecular Biochemistry, University of Dhaka, Dhaka, Bangladesh
| | - Ziaur Rahman
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | - Mustafizur Rahman
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| | - Tahmina Shirin
- Institute of Epidemiology, Disease Control and Research, Dhaka, Bangladesh
| | - Firdausi Qadri
- Institute for Developing Science and Health Initiatives, Dhaka, Bangladesh
- International Centre for Diarrhoeal Disease Research, Dhaka, Bangladesh
| |
Collapse
|
5
|
Wang Y, Yang Y, Kong L, Sasaki T, Li Q. Phylogenomic resolution of Imparidentia (Mollusca: Bivalvia) diversification through mitochondrial genomes. MARINE LIFE SCIENCE & TECHNOLOGY 2023; 5:326-336. [PMID: 37637250 PMCID: PMC10449738 DOI: 10.1007/s42995-023-00178-x] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 09/23/2021] [Accepted: 04/25/2023] [Indexed: 08/29/2023]
Abstract
Despite significant advances in the phylogenomics of bivalves over the past decade, the higher-level phylogeny of Imparidentia (a superorder of Heterodonta) remains elusive. Here, a total of five new mitochondrial sequences (Chama asperella, Chama limbula, Chama dunkeri, Barnea manilensis and Ctena divergens) was added to provide resolution in nodes that required additional study. Although the monophyly of Lucinida remains less clear, the results revealed the overall backbone of the Imparidentia tree and the monophyly of Imparidentia. Likewise, most relationships among the five major Imparidentia lineages-Lucinida, Cardiida, Adapedonta, Myida and Venerida-were addressed with a well-supported topology. Basal relationships of Imparidentia recovered Lucinidae as the sister group to all remaining imparidentian taxa. Thyasiridae is a sister group to other imparidentian bivalves (except Lucinidae species) which is split into Cardiida, Adapedonta and the divergent clade of Neoheterodontei. Neoheterodontei was comprised of Venerida and Myida, the former of which now also contains Chamidae as the sister group to all the remaining venerid taxa. Moreover, molecular divergence times were inferred by calibrating nine nodes in the Imparidentia tree of life by extinct taxa. The origin of these major clades ranged from Ordovician to Permian with the diversification through the Palaeozoic to Mesozoic. Overall, the results obtained in this study demonstrate a better-resolved Imparidentia phylogeny based on mitochondrial genomes. Supplementary Information The online version contains supplementary material available at 10.1007/s42995-023-00178-x.
Collapse
Affiliation(s)
- Yu Wang
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, 266003 China
| | - Yi Yang
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, 266003 China
| | - Lingfeng Kong
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, 266003 China
| | - Takenori Sasaki
- The University Museum, The University of Tokyo, Tokyo, 113-0033 Japan
| | - Qi Li
- Key Laboratory of Mariculture, Ministry of Education, Ocean University of China, Qingdao, 266003 China
- Laboratory for Marine Fisheries Science and Food Production Processes, Laoshan Laboratory, Qingdao, 266237 China
| |
Collapse
|
6
|
Nanni AV, Martinez N, Graze R, Morse A, Newman JRB, Jain V, Vlaho S, Signor S, Nuzhdin SV, Renne R, McIntyre LM. Sex-Biased Expression Is Associated With Chromatin State in Drosophila melanogaster and Drosophila simulans. Mol Biol Evol 2023; 40:msad078. [PMID: 37116218 PMCID: PMC10162771 DOI: 10.1093/molbev/msad078] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Revised: 02/24/2023] [Accepted: 03/13/2023] [Indexed: 04/30/2023] Open
Abstract
In Drosophila melanogaster and D. simulans head tissue, 60% of orthologous genes show evidence of sex-biased expression in at least one species. Of these, ∼39% (2,192) are conserved in direction. We hypothesize enrichment of open chromatin in the sex where we see expression bias and closed chromatin in the opposite sex. Male-biased orthologs are significantly enriched for H3K4me3 marks in males of both species (∼89% of male-biased orthologs vs. ∼76% of unbiased orthologs). Similarly, female-biased orthologs are significantly enriched for H3K4me3 marks in females of both species (∼90% of female-biased orthologs vs. ∼73% of unbiased orthologs). The sex-bias ratio in female-biased orthologs was similar in magnitude between the two species, regardless of the closed chromatin (H3K27me2me3) marks in males. However, in male-biased orthologs, the presence of H3K27me2me3 in both species significantly reduced the correlation between D. melanogaster sex-bias ratio and the D. simulans sex-bias ratio. Male-biased orthologs are enriched for evidence of positive selection in the D. melanogaster group. There are more male-biased genes than female-biased genes in both species. For orthologs with gains/losses of sex-bias between the two species, there is an excess of male-bias compared to female-bias, but there is no consistent pattern in the relationship between H3K4me3 or H3K27me2me3 chromatin marks and expression. These data suggest chromatin state is a component of the maintenance of sex-biased expression and divergence of sex-bias between species is reflected in the complexity of the chromatin status.
Collapse
Affiliation(s)
- Adalena V Nanni
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL
| | - Natalie Martinez
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
| | - Rita Graze
- Department of Biological Sciences, Auburn University, Auburn, AL
| | - Alison Morse
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL
| | - Jeremy R B Newman
- University of Florida Genetics Institute, University of Florida, Gainesville, FL
| | - Vaibhav Jain
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
| | - Srna Vlaho
- Department of Biological Sciences, University of Southern California, Los Angeles, CA
| | - Sarah Signor
- Department of Biological Sciences, North Dakota State University, Fargo, ND
| | - Sergey V Nuzhdin
- Department of Biological Sciences, University of Southern California, Los Angeles, CA
| | - Rolf Renne
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL
| | - Lauren M McIntyre
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL
| |
Collapse
|
7
|
Ramos NI, DeLeo DM, Horowitz J, McFadden CS, Quattrini AM. Selection in coral mitogenomes, with insights into adaptations in the deep sea. Sci Rep 2023; 13:6016. [PMID: 37045882 PMCID: PMC10097804 DOI: 10.1038/s41598-023-31243-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/19/2022] [Accepted: 03/08/2023] [Indexed: 04/14/2023] Open
Abstract
Corals are a dominant benthic fauna that occur across a vast range of depths from just below the ocean's surface to the abyssopelagic zone. However, little is known about the evolutionary mechanisms that enable them to inhabit such a wide range of environments. The mitochondrial (mt) genome, which is involved in energetic pathways, may be subject to selection pressures at greater depths to meet the metabolic demands of that environment. Here, we use a phylogenomic framework combined with codon-based models to evaluate whether mt protein-coding genes (PCGs) associated with cellular energy functions are under positive selection across depth in three groups of corals: Octocorallia, Scleractinia, and Antipatharia. The results demonstrated that mt PCGs of deep- and shallow-water species of all three groups were primarily under strong purifying selection (0.0474 < ω < 0.3123), with the exception of positive selection in atp6 (ω = 1.3263) of deep-sea antipatharians. We also found evidence for positive selection at fifteen sites across cox1, mtMutS, and nad1 in deep-sea octocorals and nad3 of deep-sea antipatharians. These results contribute to our limited understanding of mt adaptations as a function of depth and provide insight into the molecular response of corals to the extreme deep-sea environment.
Collapse
Affiliation(s)
- Nina I Ramos
- Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20560, USA
| | - Danielle M DeLeo
- Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20560, USA
| | - Jeremy Horowitz
- Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20560, USA
| | | | - Andrea M Quattrini
- Department of Invertebrate Zoology, National Museum of Natural History, Smithsonian Institution, Washington, DC, 20560, USA.
| |
Collapse
|
8
|
Álvarez-Carretero S, Kapli P, Yang Z. Beginner's Guide on the Use of PAML to Detect Positive Selection. Mol Biol Evol 2023; 40:7140562. [PMID: 37096789 PMCID: PMC10127084 DOI: 10.1093/molbev/msad041] [Citation(s) in RCA: 34] [Impact Index Per Article: 34.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 04/26/2023] Open
Abstract
The CODEML program in the PAML package has been widely used to analyze protein-coding gene sequences to estimate the synonymous and nonsynonymous rates (dS and dN) and to detect positive Darwinian selection driving protein evolution. For users not familiar with molecular evolutionary analysis, the program is known to have a steep learning curve. Here, we provide a step-by-step protocol to illustrate the commonly used tests available in the program, including the branch models, the site models, and the branch-site models, which can be used to detect positive selection driving adaptive protein evolution affecting particular lineages of the species phylogeny, affecting a subset of amino acid residues in the protein, and affecting a subset of sites along prespecified lineages, respectively. A data set of the myxovirus (Mx) genes from ten mammal and two bird species is used as an example. We discuss a new feature in CODEML that allows users to perform positive selection tests for multiple genes for the same set of taxa, as is common in modern genome-sequencing projects. The PAML package is distributed at https://github.com/abacus-gene/paml under the GNU license, with support provided at its discussion site (https://groups.google.com/g/pamlsoftware). Data files used in this protocol are available at https://github.com/abacus-gene/paml-tutorial.
Collapse
Affiliation(s)
- Sandra Álvarez-Carretero
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| | - Paschalia Kapli
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| | - Ziheng Yang
- Department of Genetics, Evolution and Environment, University College London, London, United Kingdom
| |
Collapse
|
9
|
Latrille T, Rodrigue N, Lartillot N. Genes and sites under adaptation at the phylogenetic scale also exhibit adaptation at the population-genetic scale. Proc Natl Acad Sci U S A 2023; 120:e2214977120. [PMID: 36897968 PMCID: PMC10089192 DOI: 10.1073/pnas.2214977120] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2022] [Accepted: 02/11/2023] [Indexed: 03/12/2023] Open
Abstract
Adaptation in protein-coding sequences can be detected from multiple sequence alignments across species or alternatively by leveraging polymorphism data within a population. Across species, quantification of the adaptive rate relies on phylogenetic codon models, classically formulated in terms of the ratio of nonsynonymous over synonymous substitution rates. Evidence of an accelerated nonsynonymous substitution rate is considered a signature of pervasive adaptation. However, because of the background of purifying selection, these models are potentially limited in their sensitivity. Recent developments have led to more sophisticated mutation-selection codon models aimed at making a more detailed quantitative assessment of the interplay between mutation, purifying, and positive selection. In this study, we conducted a large-scale exome-wide analysis of placental mammals with mutation-selection models, assessing their performance at detecting proteins and sites under adaptation. Importantly, mutation-selection codon models are based on a population-genetic formalism and thus are directly comparable to the McDonald and Kreitman test at the population level to quantify adaptation. Taking advantage of this relationship between phylogenetic and population genetics analyses, we integrated divergence and polymorphism data across the entire exome for 29 populations across 7 genera and showed that proteins and sites detected to be under adaptation at the phylogenetic scale are also under adaptation at the population-genetic scale. Altogether, our exome-wide analysis shows that phylogenetic mutation-selection codon models and the population-genetic test of adaptation can be reconciled and are congruent, paving the way for integrative models and analyses across individuals and populations.
Collapse
Affiliation(s)
- Thibault Latrille
- Université de Lyon, Université Lyon 1, CNRS, VetAgro Sup, Laboratoire de Biométrie et Biologie Evolutive, UMR5558, 69100Villeurbanne, France
- École Normale Supérieure de Lyon, Université de Lyon, 69342Lyon, France
- Department of Computational Biology, Université de Lausanne, 1015Lausanne, Switzerland
| | - Nicolas Rodrigue
- Department of Biology, Institute of Biochemistry, and School of Mathematics and Statistics, Carleton University, K1S 5B6Ottawa, Canada
| | - Nicolas Lartillot
- Université de Lyon, Université Lyon 1, CNRS, VetAgro Sup, Laboratoire de Biométrie et Biologie Evolutive, UMR5558, 69100Villeurbanne, France
| |
Collapse
|
10
|
Genome-wide identification, expression profile and evolutionary relationships of TPS genes in the neotropical fruit tree species Psidium cattleyanum. Sci Rep 2023; 13:3930. [PMID: 36894661 PMCID: PMC9998390 DOI: 10.1038/s41598-023-31061-5] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/22/2022] [Accepted: 03/06/2023] [Indexed: 03/11/2023] Open
Abstract
Terpenoids are essential for plant growth, development, defense, and adaptation mechanisms. Psidium cattleyanum (Myrtaceae) is a fleshy fruit tree species endemics from Atlantic Forest, known for its pleasant fragrance and sweet taste, attributed to terpenoids in its leaves and fruits. In this study, we conducted genome-wide identification, evolutionary and expression analyses of the terpene synthase gene (TPS) family in P. cattleyanum red guava (var. cattleyanum), and yellow guava (var. lucidum Hort.) morphotypes. We identified 32 full-length TPS in red guava (RedTPS) and 30 in yellow guava (YlwTPS). We showed different expression patterns of TPS paralogous in the two morphotypes, suggesting the existence of distinct gene regulation mechanisms and their influence on the final essential oil content in both morphotypes. Moreover, the oil profile of red guava was dominated by 1,8-cineole and linalool and yellow guava was enriched in α-pinene, coincident in proportion to TPS-b1 genes, which encode enzymes that produce cyclic monoterpenes, suggesting a lineage-specific subfamily expansion of this family. Finally, we identified amino acid residues near the catalytic center and functional areas under positive selection. Our findings provide valuable insights into the terpene biosynthesis in a Neotropical Myrtaceae species and their potential involvement in adaptation mechanisms.
Collapse
|
11
|
Divergence time of mites of the family Laelapidae based on mitochondrial barcoding region. PLoS One 2023; 18:e0279598. [PMID: 36787294 PMCID: PMC9928082 DOI: 10.1371/journal.pone.0279598] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/20/2022] [Accepted: 12/12/2022] [Indexed: 02/15/2023] Open
Abstract
Using the mitochondrial barcoding region to correlate research with 58 species in 19 genera of the family Laelapidae with the aim of determining the origin, phylogenetic relationships, and biogeographic historical distribution characteristics of mites in the family Laelapidae. Phylogenetic trees were obtained using Bayesian inference (BI) and Maximum-likelihood (ML) methods, based on three fossil records calibrated as molecular clock nodes, to estimate the divergence time of mites in the family Laelapidae as well as to apply Dispersal-Extinction-Cladogenesis (DEC) analyses to obtain biogeographic history inferences. The result showed species of the genera Hyperlaelaps and Haemolaelaps and some species of the genus Androlaelaps in the family Laelapidae were divided into clades of the genus Laelaps in both the BI and ML trees. Divergence time estimates and biogeographic history analysis revealed that the family Laelapidae likely diverged from other taxa during the Middle Jurassic (ca. 156.73 Mya), with Asia considered the most likely ancestral region for the family Laelapidae. Species of various genera began to undergo massive diversification events during the Cenozoic Tertiary. The results suggest that some genera in the family Laelapidae need to be re-defined or new genera need to be established; the Late Cretaceous to Late Neogene warm period would have promoted the divergence and expansion of species in the family Laelapidae. The divergence and dispersal of the family Laelapidae species is most likely a joint response to the continued northward drift of the Indian plate away from the Gondwana paleo-continent and gradually closer to Asia during the Late Cretaceous and the geological activity of the Tibetan Plateau during the Cenozoic Tertiary. The results strengthen our understanding of the origin and evolution of species in the family Laelapidae.
Collapse
|
12
|
Nanni AV, Martinez N, Graze R, Morse A, Newman JRB, Jain V, Vlaho S, Signor S, Nuzhdin SV, Renne R, McIntyre LM. Sex-biased expression is associated with chromatin state in D. melanogaster and D. simulans. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.13.523946. [PMID: 36711631 PMCID: PMC9882225 DOI: 10.1101/2023.01.13.523946] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/15/2023]
Abstract
We propose a new model for the association of chromatin state and sex-bias in expression. We hypothesize enrichment of open chromatin in the sex where we see expression bias (OS) and closed chromatin in the opposite sex (CO). In this study of D. melanogaster and D. simulans head tissue, sex-bias in expression is associated with H3K4me3 (open mark) in males for male-biased genes and in females for female-biased genes in both species. Sex-bias in expression is also largely conserved in direction and magnitude between the two species on the X and autosomes. In male-biased orthologs, the sex-bias ratio is more divergent between species if both species have H3K27me2me3 marks in females compared to when either or neither species has H3K27me2me3 in females. H3K27me2me3 marks in females are associated with male-bias in expression on the autosomes in both species, but on the X only in D. melanogaster . In female-biased orthologs the relationship between the species for the sex-bias ratio is similar regardless of the H3K27me2me3 marks in males. Female-biased orthologs are more similar in the ratio of sex-bias than male-biased orthologs and there is an excess of male-bias in expression in orthologs that gain/lose sex-bias. There is an excess of male-bias in sex-limited expression in both species suggesting excess male-bias is due to rapid evolution between the species. The X chromosome has an enrichment in male-limited H3K4me3 in both species and an enrichment of sex-bias in expression compared to the autosomes.
Collapse
Affiliation(s)
- Adalena V Nanni
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Natalie Martinez
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
| | - Rita Graze
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| | - Alison Morse
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Jeremy R B Newman
- University of Florida Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Vaibhav Jain
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
| | - Srna Vlaho
- Department of Biological Sciences, University of Southern California, Los Angeles, CA, USA
| | - Sarah Signor
- Department of Biological Sciences, North Dakota State University, Fargo, ND, USA
| | - Sergey V Nuzhdin
- Department of Biological Sciences, University of Southern California, Los Angeles, CA, USA
| | - Rolf Renne
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL, USA
| | - Lauren M McIntyre
- Department of Molecular Genetics and Microbiology, University of Florida, Gainesville, FL
- University of Florida Genetics Institute, University of Florida, Gainesville, FL, USA
| |
Collapse
|
13
|
Ceron-Noriega A, Almeida MV, Levin M, Butter F. Nematode gene annotation by machine-learning-assisted proteotranscriptomics enables proteome-wide evolutionary analysis. Genome Res 2023; 33:112-128. [PMID: 36653121 PMCID: PMC9977148 DOI: 10.1101/gr.277070.122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2022] [Accepted: 11/18/2022] [Indexed: 01/19/2023]
Abstract
Nematodes encompass more than 24,000 described species, which were discovered in almost every ecological habitat, and make up >80% of metazoan taxonomic diversity in soils. The last common ancestor of nematodes is believed to date back to ∼650-750 million years, generating a large and phylogenetically diverse group to be explored. However, for most species high-quality gene annotations are incomprehensive or missing. Combining short-read RNA sequencing with mass spectrometry-based proteomics and machine-learning quality control in an approach called proteotranscriptomics, we improve gene annotations for nine genome-sequenced nematode species and provide new gene annotations for three additional species without genome assemblies. Emphasizing the sensitivity of our methodology, we provide evidence for two hitherto undescribed genes in the model organism Caenorhabditis elegans Extensive phylogenetic systems analysis using this comprehensive proteome annotation provides new insights into evolutionary processes of this metazoan group.
Collapse
Affiliation(s)
| | | | - Michal Levin
- Institute of Molecular Biology (IMB), 55128 Mainz, Germany
| | - Falk Butter
- Institute of Molecular Biology (IMB), 55128 Mainz, Germany
| |
Collapse
|
14
|
Kiemel K, Gurke M, Paraskevopoulou S, Havenstein K, Weithoff G, Tiedemann R. Variation in heat shock protein 40 kDa relates to divergence in thermotolerance among cryptic rotifer species. Sci Rep 2022; 12:22626. [PMID: 36587065 PMCID: PMC9805463 DOI: 10.1038/s41598-022-27137-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/16/2022] [Accepted: 12/26/2022] [Indexed: 01/01/2023] Open
Abstract
Genetic divergence and the frequency of hybridization are central for defining species delimitations, especially among cryptic species where morphological differences are merely absent. Rotifers are known for their high cryptic diversity and therefore are ideal model organisms to investigate such patterns. Here, we used the recently resolved Brachionus calyciflorus species complex to investigate whether previously observed between species differences in thermotolerance and gene expression are also reflected in their genomic footprint. We identified a Heat Shock Protein gene (HSP 40 kDa) which exhibits cross species pronounced sequence variation. This gene exhibits species-specific fixed sites, alleles, and sites putatively under positive selection. These sites are located in protein binding regions involved in chaperoning and may therefore reflect adaptive diversification. By comparing three genetic markers (ITS, COI, HSP 40 kDa), we revealed hybridization events between the cryptic species. The low frequency of introgressive haplotypes/alleles suggest a tight, but not fully impermeable boundary between the cryptic species.
Collapse
Affiliation(s)
- K. Kiemel
- grid.11348.3f0000 0001 0942 1117Unit of Evolutionary Biology/Systematic Zoology, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht Straße 24-25, 14476 Potsdam, Germany
| | - M. Gurke
- grid.422371.10000 0001 2293 9957Museum für Naturkunde – Leibniz Institute for Evolution and Biodiversity Science, Invalidenstraße 43, 10115 Berlin, Germany ,grid.7468.d0000 0001 2248 7639Department of Biology, Humboldt-University, Invalidenstraße 42, 10115 Berlin, Germany
| | - S. Paraskevopoulou
- grid.4514.40000 0001 0930 2361Department of Biology, Lund University, Microbiology Group, Sölvegatan 35, 223 62 Lund, Sweden
| | - K. Havenstein
- grid.11348.3f0000 0001 0942 1117Unit of Evolutionary Biology/Systematic Zoology, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht Straße 24-25, 14476 Potsdam, Germany
| | - G. Weithoff
- grid.11348.3f0000 0001 0942 1117Unit of Ecology and Ecosystem Modelling, Institute for Biochemistry and Biology, University of Potsdam, Am Neuen Palais 10, 14469 Potsdam, Germany
| | - R. Tiedemann
- grid.11348.3f0000 0001 0942 1117Unit of Evolutionary Biology/Systematic Zoology, Institute for Biochemistry and Biology, University of Potsdam, Karl-Liebknecht Straße 24-25, 14476 Potsdam, Germany
| |
Collapse
|
15
|
Thompson CE, Brisolara-Corrêa L, Thompson HN, Stassen H, de Freitas LB. Evolutionary and structural aspects of Solanaceae RNases T2. Genet Mol Biol 2022; 46:e20220115. [PMID: 36534953 PMCID: PMC9762611 DOI: 10.1590/1678-4685-gmb-2022-0115] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/24/2022] [Accepted: 10/20/2022] [Indexed: 12/23/2022] Open
Abstract
Plant RNases T2 are involved in several physiological and developmental processes, including inorganic phosphate starvation, senescence, wounding, defense against pathogens, and the self-incompatibility system. Solanaceae RNases form three main clades, one composed exclusively of S-RNases and two that include S-like RNases. We identified several positively selected amino acids located in highly flexible regions of these molecules, mainly close to the B1 and B2 substrate-binding sites in S-like RNases and the hypervariable regions of S-RNases. These differences between S- and S-like RNases in the flexibility of amino acids in substrate-binding regions are essential to understand the RNA-binding process. For example, in the S-like RNase NT, two positively selected amino acid residues (Tyr156 and Asn134) are located at the most flexible sites on the molecular surface. RNase NT is induced in response to tobacco mosaic virus infection; these sites may thus be regions of interaction with pathogen proteins or viral RNA. Differential selective pressures acting on plant ribonucleases have increased amino acid variability and, consequently, structural differences within and among S-like RNases and S-RNases that seem to be essential for these proteins play different functions.
Collapse
Affiliation(s)
- Claudia Elizabeth Thompson
- Universidade Federal de Ciências da Saúde de Porto Alegre,
Departamento de Farmacociências, Porto Alegre, RS, Brazil
| | - Lauís Brisolara-Corrêa
- Universidade Federal do Rio Grande do Sul, Departamento de Genética,
Porto Alegre, RS, Brazil
| | - Helen Nathalia Thompson
- Universidade Federal do Rio Grande do Sul, Instituto de Química,
Departamento de Fisico-Química, Laboratório de Química Teórica e Computacional,
Porto Alegre, RS, Brazil
| | - Hubert Stassen
- Universidade Federal do Rio Grande do Sul, Instituto de Química,
Departamento de Fisico-Química, Laboratório de Química Teórica e Computacional,
Porto Alegre, RS, Brazil
| | | |
Collapse
|
16
|
Zhao D, Guo Y, Gao Y. Natural selection drives the evolution of mitogenomes in Acrossocheilus. PLoS One 2022; 17:e0276056. [PMID: 36227932 PMCID: PMC9560497 DOI: 10.1371/journal.pone.0276056] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/26/2022] [Accepted: 09/27/2022] [Indexed: 11/10/2022] Open
Abstract
The mitochondrial genome plays a crucial role in the balance of energy and heat production in organisms and, thus, may be under natural selection due to its potential role in adaptive divergence and speciation. Here, we studied natural selection on the mitogenome of Acrossocheilus (Cypriniformes: Cyprinidae), a genus of fish that inhabits a broad latitudinal distribution ranging from the tropics and subtropics through temperate regions. Specifically, we used 25 published mitogenome sequences of Acrossocheilus species to investigate phylogenetic relationships in this genus and detected signals of positive selection on 13 protein-coding, mitochondrial genes. We found that relaxed purifying selection and genetic drift were the predominant evolutionary forces acting on the analyzed mitogenomes. However, we also found evidence of diversifying selection on some codons, indicating episodes of positive selection. Additionally, we analyzed the mitogenomic data within an environmental modeling framework and found that the Ka/Ks ratio of ATP6 may correlated with a mean diurnal temperature range (p = 0.0449), while the Ka/Ks ratio of COX2 may correlated with precipitation during the driest month (p = 0.00761). These results suggest that the mitogenomes of Acrossocheilus species may be involved in evolutionary adaptations to different habitats. Based on this, we believe that our study provides a new insight into the role of the mitochondrial genome of Acrossocheilus species in adaptation to different environments. During our study, we also discovered several cases of paraphyly and polyphyly among accessions of species and their putative synonyms. Thus, our study suggests that a careful reassessment of the taxonomy of Acrossocheilus is using high-quality molecular data merited.
Collapse
Affiliation(s)
- Dan Zhao
- Fisheries College, Zhejiang Ocean University, Zhoushan, Zhejiang, China
| | - Yudong Guo
- Fisheries College, Zhejiang Ocean University, Zhoushan, Zhejiang, China
| | - Yang Gao
- Fisheries College, Zhejiang Ocean University, Zhoushan, Zhejiang, China
- * E-mail:
| |
Collapse
|
17
|
Heyduk K, McAssey EV, Leebens‐Mack J. Differential timing of gene expression and recruitment in independent origins of CAM in the Agavoideae (Asparagaceae). THE NEW PHYTOLOGIST 2022; 235:2111-2126. [PMID: 35596719 PMCID: PMC9796715 DOI: 10.1111/nph.18267] [Citation(s) in RCA: 8] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/01/2021] [Accepted: 05/09/2022] [Indexed: 06/15/2023]
Abstract
Crassulacean acid metabolism (CAM) photosynthesis has evolved repeatedly across the plant tree of life, however our understanding of the genetic convergence across independent origins remains hampered by the lack of comparative studies. Here, we explore gene expression profiles in eight species from the Agavoideae (Asparagaceae) encompassing three independent origins of CAM. Using comparative physiology and transcriptomics, we examined the variable modes of CAM in this subfamily and the changes in gene expression across time of day and between well watered and drought-stressed treatments. We further assessed gene expression and the molecular evolution of genes encoding phosphoenolpyruvate carboxylase (PPC), an enzyme required for primary carbon fixation in CAM. Most time-of-day expression profiles are largely conserved across all eight species and suggest that large perturbations to the central clock are not required for CAM evolution. By contrast, transcriptional response to drought is highly lineage specific. Yucca and Beschorneria have CAM-like expression of PPC2, a copy of PPC that has never been shown to be recruited for CAM in angiosperms. Together the physiological and transcriptomic comparison of closely related C3 and CAM species reveals similar gene expression profiles, with the notable exception of differential recruitment of carboxylase enzymes for CAM function.
Collapse
Affiliation(s)
- Karolina Heyduk
- School of Life SciencesUniversity of Hawaiʻi at MānoaHonoluluHI96822USA
- Department of Plant BiologyUniversity of GeorgiaAthensGA30602USA
- Department of Ecology and Evolutionary BiologyYale UniversityNew HavenCT06520USA
| | - Edward V. McAssey
- School of Life SciencesUniversity of Hawaiʻi at MānoaHonoluluHI96822USA
| | - Jim Leebens‐Mack
- Department of Plant BiologyUniversity of GeorgiaAthensGA30602USA
| |
Collapse
|
18
|
Johnson BD, Anderson AP, Small CM, Rose E, Flanagan SP, Hendrickson-Rose C, Jones AG. The evolution of the testis transcriptome in pregnant male pipefishes and seahorses. Evolution 2022; 76:2162-2180. [PMID: 35863060 DOI: 10.1111/evo.14579] [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: 01/10/2022] [Revised: 06/17/2022] [Accepted: 06/22/2022] [Indexed: 01/22/2023]
Abstract
In many animals, sperm competition and sexual conflict are thought to drive the rapid evolution of male-specific genes, especially those expressed in the testes. A potential exception occurs in the male pregnant pipefishes, where females transfer eggs to the males, eliminating testes from participating in these processes. Here, we show that testis-related genes differ dramatically in their rates of molecular evolution and expression patterns in pipefishes and seahorses (Syngnathidae) compared to other fish. Genes involved in testis or sperm function within syngnathids experience weaker selection in comparison to their orthologs in spawning and livebearing fishes. An assessment of gene turnover and expression in the testis transcriptome suggests that syngnathids have lost (or significantly reduced expression of) important classes of genes from their testis transcriptomes compared to other fish. Our results indicate that more than 50 million years of male pregnancy have removed syngnathid testes from the molecular arms race that drives the rapid evolution of male reproductive genes in other taxa.
Collapse
Affiliation(s)
| | | | - Clayton M Small
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, 97403
| | - Emily Rose
- Department of Biology, Valdosta State University, Valdosta, Georgia, 31698
| | - Sarah P Flanagan
- School of Biological Sciences, University of Canterbury, Christchurch, 8041, New Zealand
| | | | - Adam G Jones
- Department of Biological Sciences, University of Idaho, Moscow, Idaho, 83844
| |
Collapse
|
19
|
Meslin C, Mainet P, Montagné N, Robin S, Legeai F, Bretaudeau A, Johnston JS, Koutroumpa F, Persyn E, Monsempès C, François MC, Jacquin-Joly E. Spodoptera littoralis genome mining brings insights on the dynamic of expansion of gustatory receptors in polyphagous noctuidae. G3 (BETHESDA, MD.) 2022; 12:6598846. [PMID: 35652787 PMCID: PMC9339325 DOI: 10.1093/g3journal/jkac131] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/01/2022] [Accepted: 05/17/2022] [Indexed: 11/13/2022]
Abstract
The bitter taste, triggered via gustatory receptors, serves as an important natural defense against the ingestion of poisonous foods in animals, and the increased host breadth is usually linked to an increase in the number of gustatory receptor genes. This has been especially observed in polyphagous insect species, such as noctuid species from the Spodoptera genus. However, the dynamic and physical mechanisms leading to these gene expansions and the evolutionary pressures behind them remain elusive. Among major drivers of genome dynamics are the transposable elements but, surprisingly, their potential role in insect gustatory receptor expansion has not been considered yet. In this work, we hypothesized that transposable elements and possibly positive selection would be involved in the highly dynamic evolution of gustatory receptor in Spodoptera spp. We first sequenced de novo the full 465 Mb genome of S. littoralis, and manually annotated the main chemosensory genes, including a large repertoire of 373 gustatory receptor genes (including 19 pseudogenes). We also improved the completeness of S. frugiperda and S. litura gustatory receptor gene repertoires. Then, we annotated transposable elements and revealed that a particular category of class I retrotransposons, the SINE transposons, was significantly enriched in the vicinity of gustatory receptor gene clusters, suggesting a transposon-mediated mechanism for the formation of these clusters. Selection pressure analyses indicated that positive selection within the gustatory receptor gene family is cryptic, only 7 receptors being identified as positively selected. Altogether, our data provide a new good quality Spodoptera genome, pinpoint interesting gustatory receptor candidates for further functional studies and bring valuable genomic information on the mechanisms of gustatory receptor expansions in polyphagous insect species.
Collapse
Affiliation(s)
- Camille Meslin
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France
| | - Pauline Mainet
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France
| | - Nicolas Montagné
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France
| | - Stéphanie Robin
- INRAE, UMR Institut de Génétique, Environnement et Protection des Plantes (IGEPP), BioInformatics Platform for Agroecosystems Arthropods (BIPAA), Campus Beaulieu, 35042 Rennes, France.,INRIA, IRISA, GenOuest Core Facility, Campus de Beaulieu, Rennes 5042, France
| | - Fabrice Legeai
- INRAE, UMR Institut de Génétique, Environnement et Protection des Plantes (IGEPP), BioInformatics Platform for Agroecosystems Arthropods (BIPAA), Campus Beaulieu, 35042 Rennes, France.,INRIA, IRISA, GenOuest Core Facility, Campus de Beaulieu, Rennes 5042, France
| | - Anthony Bretaudeau
- INRAE, UMR Institut de Génétique, Environnement et Protection des Plantes (IGEPP), BioInformatics Platform for Agroecosystems Arthropods (BIPAA), Campus Beaulieu, 35042 Rennes, France.,INRIA, IRISA, GenOuest Core Facility, Campus de Beaulieu, Rennes 5042, France
| | - J Spencer Johnston
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA
| | - Fotini Koutroumpa
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France.,Present address: INRAE, Université Tours, Infectiologie et Santé Publique (ISP), 37380 Nouzilly, France
| | - Emma Persyn
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France.,CIRAD, UMR PVBMT, Réunion, France
| | - Christelle Monsempès
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France
| | - Marie-Christine François
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France
| | - Emmanuelle Jacquin-Joly
- INRAE, Sorbonne Université, CNRS, IRD, UPEC, Université de Paris, Institut d'Ecologie et des Sciences de l'Environnement de Paris (iEES-Paris), 78026 Versailles, France
| |
Collapse
|
20
|
Nunney L. Cancer suppression and the evolution of multiple retrogene copies of TP53 in elephants: a re‐evaluation. Evol Appl 2022; 15:891-901. [PMID: 35603034 PMCID: PMC9108310 DOI: 10.1111/eva.13383] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/04/2021] [Revised: 03/27/2022] [Accepted: 04/02/2022] [Indexed: 11/28/2022] Open
Abstract
Evolving to become bigger and/or longer lived should increase cancer susceptibility, but this predicted increase is not observed, a contradiction named Peto's paradox. A solution is that cancer suppression evolves to minimize cancer susceptibility, and the discovery of 19 retrogene (RTG) copies of the tumor suppressor gene TP53 in the African elephant (Loxodonta africana) is increasingly cited as a classic example of such adaptive suppression. However, classic examples need rigorous evaluation and an alternative hypothesis is that the RTGs spread by genetic drift. This study shows that before its duplication, the ancestral elephant RTG was already truncated from 390 amino acids to 157 by a frameshift mutation, and that 14 of the 19 copies are now truncated to ≤88 amino acids. There was no compelling evidence of either positive or negative selection acting on these 88 codons, and the pattern of RTG accumulation fits a neutral model with a duplication rate of ~10−6 per generation. It is concluded that there is no evidence supporting the hypothesis that the 19 elephant RTGs spread to fixation by selection; instead, the evidence indicates that these RTGs accumulated primarily by segmental duplication and drift. It is shown that the evolutionary multistage model of carcinogenesis (EMMC) predicts the recruitment of 1–2 independently acting tumor suppressor genes to suppress the increased cancer risk in elephants, so it is possible that one or a few RTGs may have been favored by selection resulting in the known enhanced sensitivity of elephant cells to DNA damage. However, the analysis does not provide any support for either a direct (via conserved TP53 activity) or indirect (via supporting canonical TP53 function) role of the RTGs sequences, so that the presence of multiple copies of TP53 retrogenes in elephants needs to be further justified before being used as a classic example of tumor suppression in large‐bodied animals.
Collapse
Affiliation(s)
- Leonard Nunney
- Department of Evolution, Ecology, and Organismal Biology University of California Riverside 900 University Avenue Riverside CA 92521 USA
| |
Collapse
|
21
|
Wang F, Tekle YI. Variation of natural selection in the Amoebozoa reveals heterogeneity across the phylogeny and adaptive evolution in diverse lineages. Front Ecol Evol 2022; 10:851816. [PMID: 36874909 PMCID: PMC9980437 DOI: 10.3389/fevo.2022.851816] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The evolution and diversity of the supergroup Amoebozoa is complex and poorly understood. The supergroup encompasses predominantly amoeboid lineages characterized by extreme diversity in phenotype, behavior and genetics. The study of natural selection, a driving force of diversification, within and among species of Amoebozoa will play a crucial role in understanding the evolution of the supergroup. In this study, we searched for traces of natural selection based on a set of highly conserved protein-coding genes in a phylogenetic framework from a broad sampling of amoebozoans. Using these genes, we estimated substitution rates and inferred patterns of selective pressure in lineages and sites with various models. We also examined the effect of selective pressure on codon usage bias and potential correlations with observed biological traits and habitat. Results showed large heterogeneity of selection across lineages of Amoebozoa, indicating potential species-specific optimization of adaptation to their diverse ecological environment. Overall, lineages in Tubulinea had undergone stronger purifying selection with higher average substitution rates compared to Discosea and Evosea. Evidence of adaptive evolution was observed in some representative lineages and in a gene (Rpl7a) within Evosea, suggesting potential innovation and beneficial mutations in these lineages. Our results revealed that members of the fast-evolving lineages, Entamoeba and Cutosea, all underwent strong purifying selection but had distinct patterns of codon usage bias. For the first time, this study revealed an overall pattern of natural selection across the phylogeny of Amoebozoa and provided significant implications on their distinctive evolutionary processes.
Collapse
Affiliation(s)
- Fang Wang
- Department of Biology, Spelman College, Atlanta, GA, United States
| | - Yonas I Tekle
- Department of Biology, Spelman College, Atlanta, GA, United States
| |
Collapse
|
22
|
Chen Q, Yang H, Feng X, Chen Q, Shi S, Wu CI, He Z. Two decades of suspect evidence for adaptive molecular evolution – Negative selection confounding positive selection signals. Natl Sci Rev 2021; 9:nwab217. [PMID: 35663241 PMCID: PMC9154339 DOI: 10.1093/nsr/nwab217] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/09/2021] [Accepted: 11/21/2021] [Indexed: 11/21/2022] Open
Abstract
There has been a large literature in the last two decades affirming adaptive DNA sequence evolution between species. The main lines of evidence are from (i) the McDonald-Kreitman (MK) test, which compares divergence and polymorphism data, and (ii) the phylogenetic analysis by maximum likelihood (PAML) test, which analyzes multispecies divergence data. Here, we apply these two tests concurrently to genomic data of Drosophila and Arabidopsis. To our surprise, the >100 genes identified by the two tests do not overlap beyond random expectation. Because the non-concordance could be due to low powers leading to high false negatives, we merge every 20–30 genes into a ‘supergene’. At the supergene level, the power of detection is large but the calls still do not overlap. We rule out methodological reasons for the non-concordance. In particular, extensive simulations fail to find scenarios whereby positive selection can only be detected by either MK or PAML, but not both. Since molecular evolution is governed by positive and negative selection concurrently, a fundamental assumption for estimating one of these (say, positive selection) is that the other is constant. However, in a broad survey of primates, birds, Drosophila and Arabidopsis, we found that negative selection rarely stays constant for long in evolution. As a consequence, the variation in negative selection is often misconstrued as a signal of positive selection. In conclusion, MK, PAML and any method that examines genomic sequence evolution has to explicitly address the variation in negative selection before estimating positive selection. In a companion study, we propose a possible path forward in two stages—first, by mapping out the changes in negative selection and then using this map to estimate positive selection. For now, the large literature on positive selection between species has to await reassessment.
Collapse
Affiliation(s)
- Qipian Chen
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| | - Hao Yang
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| | - Xiao Feng
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| | - Qingjian Chen
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| | - Suhua Shi
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| | - Chung-I Wu
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| | - Ziwen He
- State Key Laboratory of Biocontrol, Guangdong Key Lab of Plant Resources, School of Life Sciences, Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai), Sun Yat-sen University, Guangzhou, China
| |
Collapse
|
23
|
Duarte MA, Fernandes CR, Heckel G, da Luz Mathias M, Bastos-Silveira C. Variation and Selection in the Putative Sperm-Binding Region of ZP3 in Muroid Rodents: A Comparison between Cricetids and Murines. Genes (Basel) 2021; 12:genes12091450. [PMID: 34573431 PMCID: PMC8469249 DOI: 10.3390/genes12091450] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/19/2021] [Revised: 09/15/2021] [Accepted: 09/16/2021] [Indexed: 11/16/2022] Open
Abstract
In mammals, the zona pellucida glycoprotein 3 (ZP3) is considered a primary sperm receptor of the oocyte and is hypothesized to be involved in reproductive isolation. We investigated patterns of diversity and selection in the putative sperm-binding region (pSBR) of mouse ZP3 across Cricetidae and Murinae, two hyperdiverse taxonomic groups within muroid rodents. In murines, the pSBR is fairly conserved, in particular the serine-rich stretch containing the glycosylation sites proposed as essential for sperm binding. In contrast, cricetid amino acid sequences of the pSBR were much more variable and the serine-rich motif, typical of murines, was generally substantially modified. Overall, our results suggest a general lack of species specificity of the pSBR across the two muroid families. We document statistical evidence of positive selection acting on exons 6 and 7 of ZP3 and identified several amino acid sites that are likely targets of selection, with most positively selected sites falling within or adjacent to the pSBR.
Collapse
Affiliation(s)
- Margarida Alexandra Duarte
- Champalimaud Centre for the Uknown, Champalimaud Research, Champalimaud Foundation, Avenida Brasília, 1400-038 Lisboa, Portugal
- Museu Nacional de História Natural e da Ciência, Departamento de Zoologia e Antropologia, Universidade de Lisboa, Rua da Escola Politécnica, 58, Lisboa, 1250-102 Lisboa, Portugal
- Departamento de Biologia Animal, Faculdade de Ciências da Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal;
- Centro de Estudos de Ambiente e Mar, Faculdade de Ciências da Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal
- Correspondence:
| | - Carlos Rodríguez Fernandes
- cE3c-Centre for Ecology, Evolution and Environmental Changes, Departamento de Biologia Animal, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal; (C.R.F.); (C.B.-S.)
- Faculdade de Psicologia, Universidade de Lisboa, Alameda da Universidade, 1649-013 Lisboa, Portugal
| | - Gerald Heckel
- Institute of Ecology and Evolution, University of Bern, Baltzerstrasse 6, CH-3012 Bern, Switzerland;
- SIB Swiss Institute of Bioinformatics, Quartier Sorge-Batiment Amphipole, CH-1015 Lausanne, Switzerland
| | - Maria da Luz Mathias
- Departamento de Biologia Animal, Faculdade de Ciências da Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal;
- Centro de Estudos de Ambiente e Mar, Faculdade de Ciências da Universidade de Lisboa, Campo Grande, 1749-016 Lisboa, Portugal
| | - Cristiane Bastos-Silveira
- cE3c-Centre for Ecology, Evolution and Environmental Changes, Departamento de Biologia Animal, Faculdade de Ciências, Universidade de Lisboa, 1749-016 Lisboa, Portugal; (C.R.F.); (C.B.-S.)
| |
Collapse
|
24
|
Zhang X, Chi H, Li G, Irwin DM, Zhang S, Rossiter SJ, Liu Y. Parallel Independent Losses of G-Type Lysozyme Genes in Hairless Aquatic Mammals. Genome Biol Evol 2021; 13:6358722. [PMID: 34450623 PMCID: PMC8449827 DOI: 10.1093/gbe/evab201] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/23/2021] [Indexed: 12/03/2022] Open
Abstract
Lysozyme enzymes provide classic examples of molecular adaptation and parallel evolution, however, nearly all insights to date come from chicken-type (c-type) lysozymes. Goose-type (g-type) lysozymes occur in diverse vertebrates, with multiple independent duplications reported. Most mammals possess two g-type lysozyme genes (Lyg1 and Lyg2), the result of an early duplication, although some lineages are known to have subsequently lost one copy. Here we examine g-type lysozyme evolution across >250 mammals and reveal widespread losses of either Lyg1 or Lyg2 in several divergent taxa across the mammal tree of life. At the same time, we report strong evidence of extensive losses of both gene copies in cetaceans and sirenians, with an additional putative case of parallel loss in the tarsier. To validate these findings, we inspected published short-read data and confirmed the presence of loss of function mutations. Despite these losses, comparisons of selection pressures between intact g- and c-type lysozyme genes showed stronger purifying selection in the former, indicative of conserved function. Although the reasons for the evolutionary loss of g-type lysozymes in fully aquatic mammals are not known, we suggest that this is likely to at least partially relate to their hairlessness. Indeed, although Lyg1 does not show tissue-specific expression, recent studies have linked Lyg2 expression to anagen hair follicle development and hair loss. Such a role for g-type lysozyme would explain why the Lyg2 gene became obsolete when these taxa lost their body hair.
Collapse
Affiliation(s)
- Xiaoqing Zhang
- College of Life Sciences, Shaanxi Normal University, Xi'an, China.,College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, China
| | - Hai Chi
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Gang Li
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - David M Irwin
- Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Canada
| | - Shuyi Zhang
- College of Animal Science and Veterinary Medicine, Shenyang Agricultural University, Shenyang, China
| | - Stephen J Rossiter
- School of Biological and Chemical Sciences, Queen Mary University of London, London, United Kingdom
| | - Yang Liu
- College of Life Sciences, Shaanxi Normal University, Xi'an, China.,Key Laboratory of Zoonosis of Liaoning Province, Shenyang Agricultural University, Shenyang, China
| |
Collapse
|
25
|
Toda Y, Ko MC, Liang Q, Miller ET, Rico-Guevara A, Nakagita T, Sakakibara A, Uemura K, Sackton T, Hayakawa T, Sin SYW, Ishimaru Y, Misaka T, Oteiza P, Crall J, Edwards SV, Buttemer W, Matsumura S, Baldwin MW. Early origin of sweet perception in the songbird radiation. Science 2021; 373:226-231. [PMID: 34244416 DOI: 10.1126/science.abf6505] [Citation(s) in RCA: 24] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2020] [Accepted: 05/19/2021] [Indexed: 12/24/2022]
Abstract
Early events in the evolutionary history of a clade can shape the sensory systems of descendant lineages. Although the avian ancestor may not have had a sweet receptor, the widespread incidence of nectar-feeding birds suggests multiple acquisitions of sugar detection. In this study, we identify a single early sensory shift of the umami receptor (the T1R1-T1R3 heterodimer) that conferred sweet-sensing abilities in songbirds, a large evolutionary radiation containing nearly half of all living birds. We demonstrate sugar responses across species with diverse diets, uncover critical sites underlying carbohydrate detection, and identify the molecular basis of sensory convergence between songbirds and nectar-specialist hummingbirds. This early shift shaped the sensory biology of an entire radiation, emphasizing the role of contingency and providing an example of the genetic basis of convergence in avian evolution.
Collapse
Affiliation(s)
- Yasuka Toda
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, University of Tokyo, Tokyo 113-8657, Japan.,Department of Agricultural Chemistry, School of Agriculture, Meiji University, Kawasaki, Kanagawa 214-8571, Japan.,Japan Society for the Promotion of Science, Tokyo 102-0083, Japan
| | - Meng-Ching Ko
- Evolution of Sensory Systems Research Group, Max Planck Institute for Ornithology, Seewiesen, Germany
| | - Qiaoyi Liang
- Evolution of Sensory Systems Research Group, Max Planck Institute for Ornithology, Seewiesen, Germany
| | - Eliot T Miller
- Macaulay Library, Cornell Lab of Ornithology, Ithaca, NY, USA
| | - Alejandro Rico-Guevara
- Department of Biology, University of Washington, Seattle, WA 98105, USA.,Burke Museum of Natural History and Culture, University of Washington, Seattle, WA 98105, USA
| | - Tomoya Nakagita
- Proteo-Science Center, Ehime University, Matsuyama, Ehime 790-8577, Japan
| | - Ayano Sakakibara
- Faculty of Applied Biological Sciences, Gifu University, Gifu, 501-1193, Japan
| | - Kana Uemura
- Faculty of Applied Biological Sciences, Gifu University, Gifu, 501-1193, Japan
| | | | - Takashi Hayakawa
- Faculty of Environmental Earth Science, Hokkaido University, Sapporo, Hokkaido 060-0810, Japan.,Japan Monkey Centre, Inuyama, Aichi 484-0081, Japan
| | - Simon Yung Wa Sin
- School of Biological Sciences, The University of Hong Kong, Pok Fu Lam Road, Hong Kong.,Department of Organismic and Evolutionary Biology and the Museum of Comparative Zoology, Harvard University, 26 Oxford Street, Cambridge, MA, USA
| | - Yoshiro Ishimaru
- Department of Agricultural Chemistry, School of Agriculture, Meiji University, Kawasaki, Kanagawa 214-8571, Japan
| | - Takumi Misaka
- Department of Applied Biological Chemistry, Graduate School of Agricultural and Life Sciences, University of Tokyo, Tokyo 113-8657, Japan
| | - Pablo Oteiza
- Flow Sensing Research Group, Max Planck Institute for Ornithology, Seewiesen Germany
| | - James Crall
- Department of Organismic and Evolutionary Biology and the Museum of Comparative Zoology, Harvard University, 26 Oxford Street, Cambridge, MA, USA.,Department of Entomology, University of Wisconsin-Madison, WI, USA
| | - Scott V Edwards
- Department of Organismic and Evolutionary Biology and the Museum of Comparative Zoology, Harvard University, 26 Oxford Street, Cambridge, MA, USA
| | - William Buttemer
- Centre for Integrative Ecology, Deakin University, Geelong, Victoria, Australia.,School of Earth, Atmospheric and Life Sciences, University of Wollongong, Wollongong, New South Wales, Australia
| | - Shuichi Matsumura
- Faculty of Applied Biological Sciences, Gifu University, Gifu, 501-1193, Japan
| | - Maude W Baldwin
- Evolution of Sensory Systems Research Group, Max Planck Institute for Ornithology, Seewiesen, Germany. .,Department of Organismic and Evolutionary Biology and the Museum of Comparative Zoology, Harvard University, 26 Oxford Street, Cambridge, MA, USA
| |
Collapse
|
26
|
Junker N, Gossmann TI. Adaptation-Driven Evolution of Sirtuin 1 (SIRT1), a Key Regulator of Metabolism and Aging, in Marmot Species. Front Ecol Evol 2021. [DOI: 10.3389/fevo.2021.666564] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The sirtuin protein family plays a role in the lifespan of various species and is involved in numerous key metabolic processes. To understand the evolutionary role of sirtuins in marmots, a long-living rodent species group with remarkable metabolic shutdown during hibernation, we conducted a phylogeny-based substitution rate analysis of coding genes based on genetic information of seven marmot species. We show that sirtuin 1 (SIRT1) has evolved under positive selection in the marmot lineage. We pinpoint three amino acid changes in four different marmot species that underlie the signal of positive selection and that may favor increased longevity in marmots. Based on a computational structural analysis we can show that all three substitutions affect the secondary structure of the same region in human SIRT1. We propose that the identified region is close to the catalytic domain and that the potential structural changes may impact the catalytic activity of the enzyme and therefore might be playing a functional role in marmot's extended lifespan and metabolic shutdown.
Collapse
|
27
|
Pervaiz N, Kang H, Bao Y, Abbasi AA. Molecular evolutionary analysis of human primary microcephaly genes. BMC Ecol Evol 2021; 21:76. [PMID: 33941077 PMCID: PMC8091745 DOI: 10.1186/s12862-021-01801-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/13/2020] [Accepted: 04/22/2021] [Indexed: 12/05/2022] Open
Abstract
Background There has been a rapid increase in the brain size relative to body size during mammalian evolutionary history. In particular, the enlarged and globular brain is the most distinctive anatomical feature of modern humans that set us apart from other extinct and extant primate species. Genetic basis of large brain size in modern humans has largely remained enigmatic. Genes associated with the pathological reduction of brain size (primary microcephaly-MCPH) have the characteristics and functions to be considered ideal candidates to unravel the genetic basis of evolutionary enlargement of human brain size. For instance, the brain size of microcephaly patients is similar to the brain size of Pan troglodyte and the very early hominids like the Sahelanthropus tchadensis and Australopithecus afarensis. Results The present study investigates the molecular evolutionary history of subset of autosomal recessive primary microcephaly (MCPH) genes; CEP135, ZNF335, PHC1, SASS6, CDK6, MFSD2A, CIT, and KIF14 across 48 mammalian species. Codon based substitutions site analysis indicated that ZNF335, SASS6, CIT, and KIF14 have experienced positive selection in eutherian evolutionary history. Estimation of divergent selection pressure revealed that almost all of the MCPH genes analyzed in the present study have maintained their functions throughout the history of placental mammals. Contrary to our expectations, human-specific adoptive evolution was not detected for any of the MCPH genes analyzed in the present study. Conclusion Based on these data it can be inferred that protein-coding sequence of MCPH genes might not be the sole determinant of increase in relative brain size during primate evolutionary history. Supplementary Information The online version contains supplementary material available at 10.1186/s12862-021-01801-0.
Collapse
Affiliation(s)
- Nashaiman Pervaiz
- National Center for Bioinformatics, Program of Comparative and Evolutionary Genomics, Faculty of Biological Sciences, Quaid-i-Azam University, Islamabad, 45320, Pakistan
| | - Hongen Kang
- China National Center for Bioinformation and National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, 100101, China
| | - Yiming Bao
- China National Center for Bioinformation and National Genomics Data Center, Beijing Institute of Genomics, Chinese Academy of Sciences, Beijing, 100101, China.
| | - Amir Ali Abbasi
- National Center for Bioinformatics, Program of Comparative and Evolutionary Genomics, Faculty of Biological Sciences, Quaid-i-Azam University, Islamabad, 45320, Pakistan.
| |
Collapse
|
28
|
Salzman S, Crook D, Calonje M, Stevenson DW, Pierce NE, Hopkins R. Cycad-Weevil Pollination Symbiosis Is Characterized by Rapidly Evolving and Highly Specific Plant-Insect Chemical Communication. FRONTIERS IN PLANT SCIENCE 2021; 12:639368. [PMID: 33995438 PMCID: PMC8121082 DOI: 10.3389/fpls.2021.639368] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/08/2020] [Accepted: 04/12/2021] [Indexed: 06/12/2023]
Abstract
Coevolution between plants and insects is thought to be responsible for generating biodiversity. Extensive research has focused largely on antagonistic herbivorous relationships, but mutualistic pollination systems also likely contribute to diversification. Here we describe an example of chemically-mediated mutualistic species interactions affecting trait evolution and lineage diversification. We show that volatile compounds produced by closely related species of Zamia cycads are more strikingly different from each other than are other phenotypic characters, and that two distantly related pollinating weevil species have specialized responses only to volatiles from their specific host Zamia species. Plant transcriptomes show that approximately a fifth of genes related to volatile production are evolving under positive selection, but we find no differences in the relative proportion of genes under positive selection in different categories. The importance of phenotypic divergence coupled with chemical communication for the maintenance of this obligate mutualism highlights chemical signaling as a key mechanism of coevolution between cycads and their weevil pollinators.
Collapse
Affiliation(s)
- Shayla Salzman
- Plant Sciences, Cornell University, Ithaca, NY, United States
- Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, United States
| | - Damon Crook
- Otis Laboratory, USDA-APHIS-PPQ CPHST, Otis ANGB, MA, United States
| | - Michael Calonje
- Montgomery Botanical Center, Coral Gables, FL, United States
| | | | - Naomi E. Pierce
- Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, United States
| | - Robin Hopkins
- Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, United States
| |
Collapse
|
29
|
Razban RM, Dasmeh P, Serohijos AWR, Shakhnovich EI. Avoidance of protein unfolding constrains protein stability in long-term evolution. Biophys J 2021; 120:2413-2424. [PMID: 33932438 PMCID: PMC8390877 DOI: 10.1016/j.bpj.2021.03.042] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/16/2020] [Revised: 02/24/2021] [Accepted: 03/17/2021] [Indexed: 11/28/2022] Open
Abstract
Every amino acid residue can influence a protein's overall stability, making stability highly susceptible to change throughout evolution. We consider the distribution of protein stabilities evolutionarily permittable under two previously reported protein fitness functions: flux dynamics and misfolding avoidance. We develop an evolutionary dynamics theory and find that it agrees better with an extensive protein stability data set for dihydrofolate reductase orthologs under the misfolding avoidance fitness function rather than the flux dynamics fitness function. Further investigation with ribonuclease H data demonstrates that not any misfolded state is avoided; rather, it is only the unfolded state. At the end, we discuss how our work pertains to the universal protein abundance-evolutionary rate correlation seen across organisms' proteomes. We derive a closed-form expression relating protein abundance to evolutionary rate that captures Escherichia coli, Saccharomyces cerevisiae, and Homo sapiens experimental trends without fitted parameters.
Collapse
Affiliation(s)
- Rostam M Razban
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts
| | - Pouria Dasmeh
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts; Departement de Biochimie, Université de Montréal, Montreal, Quebec, Canada
| | | | - Eugene I Shakhnovich
- Department of Chemistry and Chemical Biology, Harvard University, Cambridge, Massachusetts.
| |
Collapse
|
30
|
Genome-Wide Analysis of Terpene Synthase Gene Family in Mentha longifolia and Catalytic Activity Analysis of a Single Terpene Synthase. Genes (Basel) 2021; 12:genes12040518. [PMID: 33918244 PMCID: PMC8066702 DOI: 10.3390/genes12040518] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Revised: 03/30/2021] [Accepted: 03/31/2021] [Indexed: 11/17/2022] Open
Abstract
Terpenoids are a wide variety of natural products and terpene synthase (TPS) plays a key role in the biosynthesis of terpenoids. Mentha plants are rich in essential oils, whose main components are terpenoids, and their biosynthetic pathways have been basically elucidated. However, there is a lack of systematic identification and study of TPS in Mentha plants. In this work, we genome-widely identified and analyzed the TPS gene family in Mentha longifolia, a model plant for functional genomic research in the genus Mentha. A total of 63 TPS genes were identified in the M. longifolia genome sequence assembly, which could be divided into six subfamilies. The TPS-b subfamily had the largest number of genes, which might be related to the abundant monoterpenoids in Mentha plants. The TPS-e subfamily had 18 members and showed a significant species-specific expansion compared with other sequenced Lamiaceae plant species. The 63 TPS genes could be mapped to nine scaffolds of the M. longifolia genome sequence assembly and the distribution of these genes is uneven. Tandem duplicates and fragment duplicates contributed greatly to the increase in the number of TPS genes in M. longifolia. The conserved motifs (RR(X)8W, NSE/DTE, RXR, and DDXXD) were analyzed in M. longifolia TPSs, and significant differentiation was found between different subfamilies. Adaptive evolution analysis showed that M. longifolia TPSs were subjected to purifying selection after the species-specific expansion, and some amino acid residues under positive selection were identified. Furthermore, we also cloned and analyzed the catalytic activity of a single terpene synthase, MlongTPS29, which belongs to the TPS-b subfamily. MlongTPS29 could encode a limonene synthase and catalyze the biosynthesis of limonene, an important precursor of essential oils from the genus Mentha. This study provides useful information for the biosynthesis of terpenoids in the genus Mentha.
Collapse
|
31
|
Wisotsky SR, Kosakovsky Pond SL, Shank SD, Muse SV. Synonymous Site-to-Site Substitution Rate Variation Dramatically Inflates False Positive Rates of Selection Analyses: Ignore at Your Own Peril. Mol Biol Evol 2021; 37:2430-2439. [PMID: 32068869 PMCID: PMC7403620 DOI: 10.1093/molbev/msaa037] [Citation(s) in RCA: 31] [Impact Index Per Article: 10.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/11/2022] Open
Abstract
Most molecular evolutionary studies of natural selection maintain the decades-old assumption that synonymous substitution rate variation (SRV) across sites within genes occurs at levels that are either nonexistent or negligible. However, numerous studies challenge this assumption from a biological perspective and show that SRV is comparable in magnitude to that of nonsynonymous substitution rate variation. We evaluated the impact of this assumption on methods for inferring selection at the molecular level by incorporating SRV into an existing method (BUSTED) for detecting signatures of episodic diversifying selection in genes. Using simulated data we found that failing to account for even moderate levels of SRV in selection testing is likely to produce intolerably high false positive rates. To evaluate the effect of the SRV assumption on actual inferences we compared results of tests with and without the assumption in an empirical analysis of over 13,000 Euteleostomi (bony vertebrate) gene alignments from the Selectome database. This exercise reveals that close to 50% of positive results (i.e., evidence for selection) in empirical analyses disappear when SRV is modeled as part of the statistical analysis and are thus candidates for being false positives. The results from this work add to a growing literature establishing that tests of selection are much more sensitive to certain model assumptions than previously believed.
Collapse
Affiliation(s)
- Sadie R Wisotsky
- Bioinformatics Research Center, North Carolina State University, Raleigh, NC.,Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, PA
| | | | - Stephen D Shank
- Institute for Genomics and Evolutionary Medicine, Temple University, Philadelphia, PA
| | - Spencer V Muse
- Bioinformatics Research Center, North Carolina State University, Raleigh, NC.,Department of Statistics, North Carolina State University, Raleigh, NC
| |
Collapse
|
32
|
Nowak MD, Birkeland S, Mandáková T, Roy Choudhury R, Guo X, Gustafsson ALS, Gizaw A, Schrøder‐Nielsen A, Fracassetti M, Brysting AK, Rieseberg L, Slotte T, Parisod C, Lysak MA, Brochmann C. The genome of Draba nivalis shows signatures of adaptation to the extreme environmental stresses of the Arctic. Mol Ecol Resour 2021; 21:661-676. [PMID: 33058468 PMCID: PMC7983928 DOI: 10.1111/1755-0998.13280] [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/01/2020] [Revised: 08/26/2020] [Accepted: 10/09/2020] [Indexed: 01/04/2023]
Abstract
The Arctic is one of the most extreme terrestrial environments on the planet. Here, we present the first chromosome-scale genome assembly of a plant adapted to the high Arctic, Draba nivalis (Brassicaceae), an attractive model species for studying plant adaptation to the stresses imposed by this harsh environment. We used an iterative scaffolding strategy with data from short-reads, single-molecule long reads, proximity ligation data, and a genetic map to produce a 302 Mb assembly that is highly contiguous with 91.6% assembled into eight chromosomes (the base chromosome number). To identify candidate genes and gene families that may have facilitated adaptation to Arctic environmental stresses, we performed comparative genomic analyses with nine non-Arctic Brassicaceae species. We show that the D. nivalis genome contains expanded suites of genes associated with drought and cold stress (e.g., related to the maintenance of oxidation-reduction homeostasis, meiosis, and signaling pathways). The expansions of gene families associated with these functions appear to be driven in part by the activity of transposable elements. Tests of positive selection identify suites of candidate genes associated with meiosis and photoperiodism, as well as cold, drought, and oxidative stress responses. Our results reveal a multifaceted landscape of stress adaptation in the D. nivalis genome, offering avenues for the continued development of this species as an Arctic model plant.
Collapse
Affiliation(s)
| | | | | | | | - Xinyi Guo
- CEITECMasaryk UniversityBrnoCzech Republic
| | | | - Abel Gizaw
- Natural History MuseumUniversity of OsloOsloNorway
| | | | - Marco Fracassetti
- Science for Life Laboratory and Department of EcologyEnvironment and Plant ScienceStockholm UniversityStockholmSweden
| | - Anne K. Brysting
- Centre for Ecological and Evolutionary SynthesisDepartment of BiosciencesUniversity of OsloOsloNorway
| | - Loren Rieseberg
- Department of BotanyThe University of British ColumbiaVancouverBCCanada
| | - Tanja Slotte
- Science for Life Laboratory and Department of EcologyEnvironment and Plant ScienceStockholm UniversityStockholmSweden
| | | | | | | |
Collapse
|
33
|
Birkeland S, Gustafsson ALS, Brysting AK, Brochmann C, Nowak MD. Multiple Genetic Trajectories to Extreme Abiotic Stress Adaptation in Arctic Brassicaceae. Mol Biol Evol 2021; 37:2052-2068. [PMID: 32167553 PMCID: PMC7306683 DOI: 10.1093/molbev/msaa068] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 02/23/2020] [Accepted: 03/10/2020] [Indexed: 12/11/2022] Open
Abstract
Extreme environments offer powerful opportunities to study how different organisms have adapted to similar selection pressures at the molecular level. Arctic plants have adapted to some of the coldest and driest biomes on Earth and typically possess suites of similar morphological and physiological adaptations to extremes in light and temperature. Here, we compare patterns of molecular evolution in three Brassicaceae species that have independently colonized the Arctic and present some of the first genetic evidence for plant adaptations to the Arctic environment. By testing for positive selection and identifying convergent substitutions in orthologous gene alignments for a total of 15 Brassicaceae species, we find that positive selection has been acting on different genes, but similar functional pathways in the three Arctic lineages. The positively selected gene sets identified in the three Arctic species showed convergent functional profiles associated with extreme abiotic stress characteristic of the Arctic. However, there was little evidence for independently fixed mutations at the same sites and for positive selection acting on the same genes. The three species appear to have evolved similar suites of adaptations by modifying different components in similar stress response pathways, implying that there could be many genetic trajectories for adaptation to the Arctic environment. By identifying candidate genes and functional pathways potentially involved in Arctic adaptation, our results provide a framework for future studies aimed at testing for the existence of a functional syndrome of Arctic adaptation in the Brassicaceae and perhaps flowering plants in general.
Collapse
Affiliation(s)
- Siri Birkeland
- Natural History Museum, University of Oslo, Oslo, Norway
| | | | - Anne K Brysting
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | | | | |
Collapse
|
34
|
Andersen MJ, McCullough JM, Gyllenhaal EF, Mapel XM, Haryoko T, Jønsson KA, Joseph L. Complex histories of gene flow and a mitochondrial capture event in a nonsister pair of birds. Mol Ecol 2021; 30:2087-2103. [PMID: 33615597 PMCID: PMC8252742 DOI: 10.1111/mec.15856] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2019] [Revised: 02/02/2021] [Accepted: 02/12/2021] [Indexed: 01/03/2023]
Abstract
Hybridization, introgression, and reciprocal gene flow during speciation, specifically the generation of mitonuclear discordance, are increasingly observed as parts of the speciation process. Genomic approaches provide insight into where, when, and how adaptation operates during and after speciation and can measure historical and modern introgression. Whether adaptive or neutral in origin, hybridization can cause mitonuclear discordance by placing the mitochondrial genome of one species (or population) in the nuclear background of another species. The latter, introgressed species may eventually have its own mtDNA replaced or “captured” by other species across its entire geographical range. Intermediate stages in the capture process should be observable. Two nonsister species of Australasian monarch‐flycatchers, Spectacled Monarch (Symposiachrus trivirgatus) mostly of Australia and Indonesia and Spot‐winged Monarch (S. guttula) of New Guinea, present an opportunity to observe this process. We analysed thousands of single nucleotide polymorphisms (SNPs) derived from ultraconserved elements of all subspecies of both species. Mitochondrial DNA sequences of Australian populations of S. trivirgatus form two paraphyletic clades, one being sister to and presumably introgressed by S. guttula despite little nuclear signal of introgression. Population genetic analyses (e.g., tests for modern and historical gene flow and selection) support at least one historical gene flow event between S. guttula and Australian S. trivirgatus. We also uncovered introgression from the Maluku Islands subspecies of S. trivirgatus into an island population of S. guttula, resulting in apparent nuclear paraphyly. We find that neutral demographic processes, not adaptive introgression, are the most likely cause of these complex population histories. We suggest that a Pleistocene extinction of S. guttula from mainland Australia resulted from range expansion by S. trivirgatus.
Collapse
Affiliation(s)
- Michael J Andersen
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Jenna M McCullough
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Ethan F Gyllenhaal
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA
| | - Xena M Mapel
- Department of Biology and Museum of Southwestern Biology, University of New Mexico, Albuquerque, New Mexico, USA.,Animal Genomics, ETH Zürich, Lindau, Switzerland
| | - Tri Haryoko
- Museum Zoologicum Bogoriense, Research Centre for Biology, Indonesian Institute of Sciences (LIPI), Cibinong, Indonesia
| | - Knud A Jønsson
- Natural History Museum of Denmark, University of Copenhagen, Copenhagen Ø, Denmark
| | - Leo Joseph
- Australian National Wildlife Collection, CSIRO National Research Collections, Canberra, Australian Capital Territory, Australia
| |
Collapse
|
35
|
Resequencing and SNP discovery of Amur ide (Leuciscus waleckii) provides insights into local adaptations to extreme environments. Sci Rep 2021; 11:5064. [PMID: 33658614 PMCID: PMC7930030 DOI: 10.1038/s41598-021-84652-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2020] [Accepted: 02/18/2021] [Indexed: 01/31/2023] Open
Abstract
Amur ide (Leuciscus waleckii), a Cyprinid species, is broadly distributed in Northeast Asia. Different from its freshwater counterparts, the population in Lake Dali Nor has a strong alkalinity tolerance and can adapt to extremely alkali-saline water with bicarbonate over 50 mmol/L. To uncover the genetic basis of its alkaline adaptation, three populations, including one alkali form from Lake Dali Nor (DL), one freshwater form from its adjacent sister Lake Ganggeng Nor (GG), and one freshwater form from its historical origin, namely, the Songhua River (SH), were analyzed using genome resequencing technology. A total of 679.82 Gb clean data and 38,091,163 high-quality single-nucleotide polymorphism (SNP) loci were detected in the three populations. Nucleotide diversity and population structure analysis revealed that the DL and GG populations have lower nucleotide diversities and different genetic structures than those of the SH population. Selective sweeping showed 21 genes involved in osmoregulatory regulation (DLG1, VIPR1, AKT1, and GNAI1), inflammation and immune responses (DLG1, BRINP1, CTSL, TRAF6, AKT1, STAT3, GNAI1, SEC22b, and PSME4b), and cardiorespiratory development (TRAF6, PSME4b, STAT3, AKT1, and COL9A1) to be associated with alkaline adaption of the DL population. Interestingly, selective pressure (CodeML, MEME, and FEL) methods identified two functional codon sites of VIPR1 to be under positive selection in the DL population. The subsequent 3D protein modeling confirmed that these selected sites will incur changes in protein structure and function in the DL population. In brief, this study provides molecular evidence of population divergence and alkaline adaptation, which will be very useful for revealing the genetic basis of alkaline adaptation in Amur ide.
Collapse
|
36
|
Jones CT, Youssef N, Susko E, Bielawski JP. A Phenotype-Genotype Codon Model for Detecting Adaptive Evolution. Syst Biol 2021; 69:722-738. [PMID: 31730199 DOI: 10.1093/sysbio/syz075] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/07/2019] [Revised: 11/09/2019] [Accepted: 11/11/2019] [Indexed: 01/03/2023] Open
Abstract
A central objective in biology is to link adaptive evolution in a gene to structural and/or functional phenotypic novelties. Yet most analytic methods make inferences mainly from either phenotypic data or genetic data alone. A small number of models have been developed to infer correlations between the rate of molecular evolution and changes in a discrete or continuous life history trait. But such correlations are not necessarily evidence of adaptation. Here, we present a novel approach called the phenotype-genotype branch-site model (PG-BSM) designed to detect evidence of adaptive codon evolution associated with discrete-state phenotype evolution. An episode of adaptation is inferred under standard codon substitution models when there is evidence of positive selection in the form of an elevation in the nonsynonymous-to-synonymous rate ratio $\omega$ to a value $\omega > 1$. As it is becoming increasingly clear that $\omega > 1$ can occur without adaptation, the PG-BSM was formulated to infer an instance of adaptive evolution without appealing to evidence of positive selection. The null model makes use of a covarion-like component to account for general heterotachy (i.e., random changes in the evolutionary rate at a site over time). The alternative model employs samples of the phenotypic evolutionary history to test for phenomenological patterns of heterotachy consistent with specific mechanisms of molecular adaptation. These include 1) a persistent increase/decrease in $\omega$ at a site following a change in phenotype (the pattern) consistent with an increase/decrease in the functional importance of the site (the mechanism); and 2) a transient increase in $\omega$ at a site along a branch over which the phenotype changed (the pattern) consistent with a change in the site's optimal amino acid (the mechanism). Rejection of the null is followed by post hoc analyses to identify sites with strongest evidence for adaptation in association with changes in the phenotype as well as the most likely evolutionary history of the phenotype. Simulation studies based on a novel method for generating mechanistically realistic signatures of molecular adaptation show that the PG-BSM has good statistical properties. Analyses of real alignments show that site patterns identified post hoc are consistent with the specific mechanisms of adaptation included in the alternate model. Further simulation studies show that the covarion-like component of the PG-BSM plays a crucial role in mitigating recently discovered statistical pathologies associated with confounding by accounting for heterotachy-by-any-cause. [Adaptive evolution; branch-site model; confounding; mutation-selection; phenotype-genotype.].
Collapse
Affiliation(s)
- Christopher T Jones
- Department of Mathematics and Statistics, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada
| | - Noor Youssef
- Department of Biology, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada
| | - Edward Susko
- Department of Mathematics and Statistics, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada.,Centre for Comparative Genomics and Evolutionary Bioinformatics, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada
| | - Joseph P Bielawski
- Department of Mathematics and Statistics, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada.,Department of Biology, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada.,Centre for Comparative Genomics and Evolutionary Bioinformatics, Dalhousie University, 1233 LeMarchant Street, B3H 4R2, Halifax, Nova Scotia, Canada
| |
Collapse
|
37
|
Tian X, Guo J, Zhou X, Ma K, Ma Y, Shi T, Shi Y. Comparative and Evolutionary Analyses on the Complete Plastomes of Five Kalanchoe Horticultural Plants. FRONTIERS IN PLANT SCIENCE 2021; 12:705874. [PMID: 34512691 PMCID: PMC8429837 DOI: 10.3389/fpls.2021.705874] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/06/2021] [Accepted: 08/02/2021] [Indexed: 05/11/2023]
Abstract
Many species of the genus Kalanchoe are important horticultural plants. They have evolved the Crassulacean acid metabolism (CAM) photosynthetic pathway to allow them to be better adapted to dry environments. Despite their importance, it is still debating whether Kalanchoe is monophyletic, and understanding the past diversification of this genus requires a tremendous amount of effort and work being devoted to the studies of morphological and molecular characters of this genus. However, molecular information, plastic sequence data, in particular, reported on Kalanchoe species is scarce, and this has posed a great challenge in trying to interpret the evolutionary history of this genus. In this study, plastomes of the five Kalanchoe species, including Kalanchoe daigremontiana, Kalanchoe delagoensis, Kalanchoe fedtschenkoi, Kalanchoe longiflora, and Kalanchoe pinnata, were sequenced and analyzed. The results indicate that the five plastomes are comparable in size, guanine-cytosine (GC) contents and the number of genes, which also demonstrate an insignificant difference in comparison with other species from the family Crassulaceae. About 224 simple sequence repeats (SSRs) and 144 long repeats were identified in the five plastomes, and most of these are distributed in the inverted repeat regions. In addition, highly divergent regions containing either single nucleotide polymorphism (SNP) or insertion or deletion (InDel) mutations are discovered, which could be potentially used for establishing phylogenetic relationships among members of the Kalanchoe genus in future studies. Furthermore, phylogenetic analyses suggest that Bryophyllum should be placed into one single genus as Kalanchoe. Further genomic analyses also reveal that several genes are undergone positive selection. Among them, 11 genes are involved in important cellular processes, such as cell survival, electron transfer, and may have played indispensable roles in the adaptive evolution of Kalanchoe to dry environments.
Collapse
Affiliation(s)
- Xiangyu Tian
- School of Life Sciences, Zhengzhou University, Zhengzhou, China
| | - Jia Guo
- School of Life Sciences, Zhengzhou University, Zhengzhou, China
| | - Xiaojiao Zhou
- School of Life Sciences, Zhengzhou University, Zhengzhou, China
| | - Ke Ma
- School of Life Sciences, Zhengzhou University, Zhengzhou, China
| | - Yonggui Ma
- Key Laboratory of Medicinal Animal and Plant Resources of Qinghai-Tibetan Plateau, Qinghai Normal University, Xining, China
| | - Tuansheng Shi
- School of Life Sciences, Zhengzhou University, Zhengzhou, China
| | - Yuhua Shi
- School of Life Sciences, Zhengzhou University, Zhengzhou, China
- *Correspondence: Yuhua Shi
| |
Collapse
|
38
|
Campbell SK, Cortés-Ortiz L. Oxytocin amino acid variation within Neotropical primates: new genetic variants in hormone and receptor sequences and evidence for evolutionary forces driving this unexpected diversity. Biol J Linn Soc Lond 2020. [DOI: 10.1093/biolinnean/blaa173] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022]
Abstract
Abstract
Oxytocin is a mammalian neuropeptide hormone that mediates behaviours important to reproduction. Despite almost universal amino acid sequence conservation across most groups of mammals, several unique forms have been reported across Neotropical primates. To explore sequence diversity, we investigated the genes encoding oxytocin and its receptor across the Atelidae, which was known to contain at least three unique oxytocin sequences. Additionally, we included the genus Cebus, within the Cebidae, to further explore the ubiquity of the Pro8 variant in this family. We found a novel amino acid variant (Val3) within the Atelidae radiation, bringing the total number of oxytocin sequences within Neotropical primates to seven. Analyses of physicochemical properties revealed conservative substitutions that are likely tolerated within the selective constraints imposed by receptor binding. Furthermore, we report radical substitutions at the eighth codon and evidence for co-evolution between Pro8 and a ligand-binding region of the oxytocin receptor in the Atelidae, supporting the notion that this variant may affect binding specificity. Overall, we suggest that selective constraint on binding specificity may maintain proper oxytocin function and that the diversification of amino acid sequence is likely due to a variety of processes such as relaxed constraint, neutral mutation, positive selection and coevolution.
Collapse
Affiliation(s)
- Susanna K Campbell
- Department of Ecology and Evolutionary Biology, The University of Michigan, Ann Arbor, MI, USA
| | - Liliana Cortés-Ortiz
- Department of Ecology and Evolutionary Biology, The University of Michigan, Ann Arbor, MI, USA
| |
Collapse
|
39
|
Ajayi OO, Showalter AM. Systems identification and characterization of β-glucuronosyltransferase genes involved in arabinogalactan-protein biosynthesis in plant genomes. Sci Rep 2020; 10:20562. [PMID: 33239665 PMCID: PMC7689455 DOI: 10.1038/s41598-020-72658-4] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2020] [Accepted: 08/24/2020] [Indexed: 11/10/2022] Open
Abstract
Utilizing plant biomass for bioethanol production requires an understanding of the molecular mechanisms involved in plant cell wall assembly. Arabinogalactan-proteins (AGPs) are glycoproteins that interact with other cell wall polymers to influence plant growth and developmental processes. Glucuronic acid, which is transferred to the AGP glycan by β-glucuronosyltransferases (GLCATs), is the only acidic sugar in AGPs with the ability to bind calcium. We carried out a comprehensive genome-wide analysis of a putative GLCAT gene family involved in AGP biosynthesis by examining its sequence diversity, genetic architecture, phylogenetic and motif characteristics, selection pressure and gene expression in plants. We report the identification of 161 putative GLCAT genes distributed across 14 plant genomes and a widely conserved GLCAT catalytic domain. We discovered a phylogenetic clade shared between bryophytes and higher land plants of monocot grass and dicot lineages and identified positively selected sites that do not result in functional divergence of GLCATs. RNA-seq and microarray data analyses of the putative GLCAT genes revealed gene expression signatures that likely influence the assembly of plant cell wall polymers which is critical to the overall growth and development of edible and bioenergy crops.
Collapse
Affiliation(s)
- Oyeyemi Olugbenga Ajayi
- Department of Environmental and Plant Biology, Ohio University, Athens, 45701 USA
- Molecular and Cellular Biology Program, Ohio University, Athens, 45701 USA
| | - Allan M. Showalter
- Department of Environmental and Plant Biology, Ohio University, Athens, 45701 USA
- Molecular and Cellular Biology Program, Ohio University, Athens, 45701 USA
| |
Collapse
|
40
|
Rowe M, Whittington E, Borziak K, Ravinet M, Eroukhmanoff F, Sætre GP, Dorus S. Molecular Diversification of the Seminal Fluid Proteome in a Recently Diverged Passerine Species Pair. Mol Biol Evol 2020; 37:488-506. [PMID: 31665510 PMCID: PMC6993853 DOI: 10.1093/molbev/msz235] [Citation(s) in RCA: 28] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Seminal fluid proteins (SFPs) mediate an array of postmating reproductive processes that influence fertilization and fertility. As such, it is widely held that SFPs may contribute to postmating, prezygotic reproductive barriers between closely related taxa. We investigated seminal fluid (SF) diversification in a recently diverged passerine species pair (Passer domesticus and Passer hispaniolensis) using a combination of proteomic and comparative evolutionary genomic approaches. First, we characterized and compared the SF proteome of the two species, revealing consistencies with known aspects of SFP biology and function in other taxa, including the presence and diversification of proteins involved in immunity and sperm maturation. Second, using whole-genome resequencing data, we assessed patterns of genomic differentiation between house and Spanish sparrows. These analyses detected divergent selection on immunity-related SF genes and positive selective sweeps in regions containing a number of SF genes that also exhibited protein abundance diversification between species. Finally, we analyzed the molecular evolution of SFPs across 11 passerine species and found a significantly higher rate of positive selection in SFPs compared with the rest of the genome, as well as significant enrichments for functional pathways related to immunity in the set of positively selected SF genes. Our results suggest that selection on immunity pathways is an important determinant of passerine SF composition and evolution. Assessing the role of immunity genes in speciation in other recently diverged taxa should be prioritized given the potential role for immunity-related proteins in reproductive incompatibilities in Passer sparrows.
Collapse
Affiliation(s)
- Melissah Rowe
- Natural History Museum, University of Oslo, Oslo, Norway.,Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway.,Department of Animal Ecology, Netherlands Institute of Ecology (NIOO-KNAW), Wageningen, The Netherlands
| | - Emma Whittington
- Center for Reproductive Evolution, Department of Biology, Syracuse University, Syracuse, NY
| | - Kirill Borziak
- Center for Reproductive Evolution, Department of Biology, Syracuse University, Syracuse, NY
| | - Mark Ravinet
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Fabrice Eroukhmanoff
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Glenn-Peter Sætre
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Steve Dorus
- Center for Reproductive Evolution, Department of Biology, Syracuse University, Syracuse, NY
| |
Collapse
|
41
|
Ma S, Liu H, Sun W, Mustafa A, Xi Y, Pu F, Li Y, Han C, Bai L, Hua H. Molecular evolution of the ATP-binding cassette subfamily G member 2 gene subfamily and its paralogs in birds. BMC Evol Biol 2020; 20:85. [PMID: 32664916 PMCID: PMC7362505 DOI: 10.1186/s12862-020-01654-z] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/15/2019] [Accepted: 07/07/2020] [Indexed: 12/11/2022] Open
Abstract
Background ATP-binding cassette (ABC) transporters are involved in the active transportation of various endogenous or exogenous substances. Two ABCG2 gene subfamily members have been identified in birds. A detailed comparative study of the ABCG2 and ABCG2-like genes aid our understanding of their evolutionary history at the molecular level and provide a theoretical reference for studying the specific functions of ABCG2 and ABCG2-like genes in birds. Results We first identified 77 ABCG2/ABCG2-like gene sequences in the genomes of 41 birds. Further analysis showed that both the nucleic acid and amino acid sequences of ABCG2 and ABCG2-like genes were highly conserved and exhibited high homology in birds. However, significant differences in the N-terminal structure were found between the ABCG2 and ABCG2-like amino acid sequences. A selective pressure analysis showed that the ABCG2 and ABCG2-like genes were affected by purifying selection during the process of bird evolution. Conclusions We believe that multiple members of the ABCG2 gene subfamily exist on chromosome 4 in the ancestors of birds. Over the long course of evolution, only the ABCG2 gene was retained on chromosome 4 in birds. The ABCG2-like gene on chromosome 6 might have originated from chromosome replication or fusion. The structural differences between the N terminus of ABCG2 protein and those of ABCG2-like proteins might lead to functional differences between the corresponding genes.
Collapse
Affiliation(s)
- Shengchao Ma
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - Hehe Liu
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China.
| | - Wenqiang Sun
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - Ahsan Mustafa
- Institute of Animal Nutrition, Key Laboratory for Animal Disease-Resistance Nutrition of China, Ministry of Education, Sichuan Agricultural University, Chengdu, P.R. China
| | - Yang Xi
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - Fajun Pu
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - Yanying Li
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - Chunchun Han
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - Lili Bai
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| | - He Hua
- Farm Animal Genetic Resources Exploration and Innovation Key Laboratory of Sichuan Province, Sichuan Agricultural University, Chengdu, Sichuan, 611130, P.R. China
| |
Collapse
|
42
|
Geyer LB, Zigler KS, Tiozzo S, Lessios HA. Slow evolution under purifying selection in the gamete recognition protein bindin of the sea urchin Diadema. Sci Rep 2020; 10:9834. [PMID: 32555217 PMCID: PMC7299941 DOI: 10.1038/s41598-020-66390-2] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2020] [Accepted: 05/06/2020] [Indexed: 01/22/2023] Open
Abstract
Bindin is a sperm protein that mediates attachment and membrane fusion of gametes. The mode of bindin evolution varies across sea urchin genera studied to date. In three genera it evolves under positive selection, in four under mostly purifying selection, and in one, results have been mixed. We studied bindin evolution in the pantropical sea urchin Diadema, which split from other studied genera 250 million years ago. We found that Diadema bindin is structurally similar to that of other genera, but much longer (418 amino acids). In seven species of Diadema, bindin evolves under purifying selection, more slowly than in any other sea urchin genus. Only bindin of the recently rediscovered D. clarki shows evidence of positive selection. As D. clarki is sympatric with D. setosum and D. savignyi, positive selection could arise from avoidance of maladaptive hybridization. However, D. setosum and D. savignyi overlap in the Indo-West Pacific, yet their bindins show no evidence of positive selection, possibly because the two species spawn at different times. Bindin in the East Pacific D. mexicanum, the West Atlantic D. antillarum, the East Atlantic D. africanum, and the Indo-Pacific D. paucispinum also evolves slowly under purifying selection.
Collapse
Affiliation(s)
- L B Geyer
- Smithsonian Tropical Research Institute, Apartado Postal 0843-03092, Balboa, Ancon, Panama.
| | - K S Zigler
- Department of Biology, Sewanee: University of the South, 735 University Ave., Sewanee, TN, 37383, United States
| | - S Tiozzo
- Sorbonne Universite, CNRS, Laboratoire de Biologie du Developpement de Villefranche-sur-mer (LBDV), 06230, Paris, France
| | - H A Lessios
- Smithsonian Tropical Research Institute, Apartado Postal 0843-03092, Balboa, Ancon, Panama
| |
Collapse
|
43
|
Mingrone J, Susko E, Bielawski JP. ModL: exploring and restoring regularity when testing for positive selection. Bioinformatics 2020; 35:2545-2554. [PMID: 30541063 DOI: 10.1093/bioinformatics/bty1019] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2017] [Revised: 11/02/2018] [Accepted: 12/11/2018] [Indexed: 11/14/2022] Open
Abstract
MOTIVATION Likelihood ratio tests are commonly used to test for positive selection acting on proteins. They are usually applied with thresholds for declaring a protein under positive selection determined from a chi-square or mixture of chi-square distributions. Although it is known that such distributions are not strictly justified due to the statistical irregularity of the problem, the hope has been that the resulting tests are conservative and do not lose much power in comparison with the same test using the unknown, correct threshold. We show that commonly used thresholds need not yield conservative tests, but instead give larger than expected Type I error rates. Statistical regularity can be restored by using a modified likelihood ratio test. RESULTS We give theoretical results to prove that, if the number of sites is not too small, the modified likelihood ratio test gives approximately correct Type I error probabilities regardless of the parameter settings of the underlying null hypothesis. Simulations show that modification gives Type I error rates closer to those stated without a loss of power. The simulations also show that parameter estimation for mixture models of codon evolution can be challenging in certain data-generation settings with very different mixing distributions giving nearly identical site pattern distributions unless the number of taxa and tree length are large. Because mixture models are widely used for a variety of problems in molecular evolution, the challenges and general approaches to solving them presented here are applicable in a broader context. AVAILABILITY AND IMPLEMENTATION https://github.com/jehops/codeml_modl. SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Joseph Mingrone
- Department of Mathematics and Statistics, Dalhousie University, Halifax, NS, Canada
- Centre for Comparative Genomics and Evolutionary Bioinformatics, Dalhousie University, Halifax, NS, Canada
| | - Edward Susko
- Department of Mathematics and Statistics, Dalhousie University, Halifax, NS, Canada
- Centre for Comparative Genomics and Evolutionary Bioinformatics, Dalhousie University, Halifax, NS, Canada
| | - Joseph P Bielawski
- Department of Mathematics and Statistics, Dalhousie University, Halifax, NS, Canada
- Centre for Comparative Genomics and Evolutionary Bioinformatics, Dalhousie University, Halifax, NS, Canada
- Department of Biology, Dalhousie University, Halifax, NS, Canada
| |
Collapse
|
44
|
Zhou J, Ren H, Hu M, Zhou J, Li B, Kong N, Zhang Q, Jin Y, Liang L, Yue J. Characterization of Burkholderia cepacia Complex Core Genome and the Underlying Recombination and Positive Selection. Front Genet 2020; 11:506. [PMID: 32528528 PMCID: PMC7253759 DOI: 10.3389/fgene.2020.00506] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2019] [Accepted: 04/24/2020] [Indexed: 11/13/2022] Open
Abstract
Recombination and positive selection are two key factors that play a vital role in pathogenic microorganisms’ population adaptation and diversification. The Burkholderia cepacia complex (Bcc) represents bacterial species with high similarity, which can cause severe infections among cases suffering from the chronic granulomatous disorder and cystic fibrosis (CF). At present, no genome-wide study has been carried out focusing on investigating the core genome of Bcc associated with the two evolutionary forces. The general characteristics of the core genome of Bcc species remain scarce as well. In this study, we explored the core orthologous genes of 116 Bcc strains using comparative genomic analysis and studied the two adaptive evolutionary forces: recombination and positive selection. We estimated 1005 orthogroups consisting entirely of single copy genes. These single copy orthologous genes in some Cluster of Orthologous Groups (COG) categories showed significant differences in the comparison of several evolutionary properties, and the encoding proteins were relatively simple and compact. Our findings showed that 5.8% of the core orthologous genes strongly supported recombination; in the meantime, 1.1% supported positive selection. We found that genes involved in protein synthesis as well as material transport and metabolism are favored by selection pressure. More importantly, homologous recombination contributed more genetic variation to a large number of genes and largely maintained the genetic cohesion in Bcc. This high level of recombination between Bcc species blurs their taxonomic boundaries, which leads Bcc species to be difficult or impossible to distinguish phenotypically and genotypically.
Collapse
Affiliation(s)
- Jianglin Zhou
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Hongguang Ren
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Mingda Hu
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Jing Zhou
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Beiping Li
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Na Kong
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China.,Institutes of Physical Science and Information Technology, Anhui University, Hefei, China
| | - Qi Zhang
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Yuan Jin
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Long Liang
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| | - Junjie Yue
- State Key Laboratory of Pathogen and Biosecurity, Beijing Institute of Biotechnology, Beijing, China
| |
Collapse
|
45
|
Zhang Y, Feng X, Wang L, Su Y, Chu Z, Sun Y. The structure, functional evolution, and evolutionary trajectories of the H +-PPase gene family in plants. BMC Genomics 2020; 21:195. [PMID: 32122295 PMCID: PMC7053079 DOI: 10.1186/s12864-020-6604-2] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Accepted: 02/21/2020] [Indexed: 08/16/2023] Open
Abstract
BACKGROUND The H+-PPase (pyrophosphatase) gene family is an important class of proton transporters that play key roles in plant development and stress resistance. Although the physiological and biochemical functions of H+-PPases are well characterized, the structural evolution and functional differentiation of this gene family remain unclear. RESULTS We identified 124 H+-PPase members from 27 plant species using complete genomic data obtained from algae to angiosperms. We found that all analyzed plants carried H+-PPase genes, and members were not limited to the two main types (type I and II). Differentiation of this gene family occurred early in evolutionary history, probably prior to the emergence of algae. The type I and II H+-PPase genes were retained during the subsequent evolution of higher plants, and their copy numbers increased rapidly in some angiosperms following whole-genome duplication (WGD) events, with obvious expression pattern differentiation among the new copies. We found significant functional divergence between type I and II H+-PPase genes, with both showing evidence for positive selection pressure. We classified angiosperm type I H+-PPases into subtypes Ia and non-Ia, which probably differentiated at an early stage of angiosperm evolution. Compared with non-Ia subtype, the Ia subtype appears to confer some advantage in angiosperms, as it is highly conserved and abundantly expressed, but shows no evidence for positive selection. CONCLUSIONS We hypothesized that there were many types of H+-PPase genes in the plant ancestral genome, and that different plant groups retained different types of these genes. In the early stages of angiosperm evolution, the type I H+-PPase genes differentiated into various subtypes. In addition, the expression pattern varied not only among genes of different types or subtypes, but also among copies of the same subtype. Based on the expression patterns and copy numbers of H+-PPase genes in higher plants, we propose two possible evolutionary trajectories for this gene family.
Collapse
Affiliation(s)
- Yiming Zhang
- College of Life Sciences, Langfang Normal University, Langfang, 065000, China
| | - Xue Feng
- College of Life Sciences, Langfang Normal University, Langfang, 065000, China
| | - Lihui Wang
- College of Plant Protection, Fujian Agricultural and Forestry University, Fuzhou, 350000, China
| | - Yanping Su
- College of Life Sciences, Langfang Normal University, Langfang, 065000, China
| | - Zhuodong Chu
- College of Life Sciences, Langfang Normal University, Langfang, 065000, China
| | - Yanxiang Sun
- College of Life Sciences, Langfang Normal University, Langfang, 065000, China.
| |
Collapse
|
46
|
Chong RA, Park H, Moran NA. Genome Evolution of the Obligate Endosymbiont Buchnera aphidicola. Mol Biol Evol 2020; 36:1481-1489. [PMID: 30989224 DOI: 10.1093/molbev/msz082] [Citation(s) in RCA: 56] [Impact Index Per Article: 14.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/25/2022] Open
Abstract
An evolutionary consequence of uniparentally transmitted symbiosis is degradation of symbiont genomes. We use the system of aphids and their maternally inherited obligate endosymbiont, Buchnera aphidicola, to explore the evolutionary process of genome degradation. We compared complete genome sequences for 39 Buchnera strains, including 23 newly sequenced symbiont genomes from diverse aphid hosts. We reconstructed the genome of the most recent shared Buchnera ancestor, which contained 616 protein-coding genes, and 39 RNA genes. The extent of subsequent gene loss varied across lineages, resulting in modern genomes ranging from 412 to 646 kb and containing 354-587 protein-coding genes. Loss events were highly nonrandom across loci. Genes involved in replication, transcription, translation, and amino acid biosynthesis are largely retained, whereas genes underlying ornithine biosynthesis, stress responses, and transcriptional regulation were lost repeatedly. Aside from losses, gene order is almost completely stable. The main exceptions involve movement between plasmid and chromosome locations of genes underlying tryptophan and leucine biosynthesis and supporting nutrition of aphid hosts. This set of complete genomes enabled tests for signatures of positive diversifying selection. Of 371 Buchnera genes tested, 29 genes show strong support for ongoing positive selection. These include genes encoding outer membrane porins that are expected to be involved in direct interactions with hosts. Collectively, these results indicate that extensive genome reduction occurred in the ancestral Buchnera prior to aphid diversification and that reduction has continued since, with losses greater in some lineages and for some loci.
Collapse
Affiliation(s)
- Rebecca A Chong
- Department of Integrative Biology, University of Texas at Austin, Austin, TX
- Department of Biology, University of Hawaii at Mānoa, Honolulu, HI
| | - Hyunjin Park
- Department of Integrative Biology, University of Texas at Austin, Austin, TX
| | - Nancy A Moran
- Department of Integrative Biology, University of Texas at Austin, Austin, TX
| |
Collapse
|
47
|
Hermida-Carrera C, Fares MA, Font-Carrascosa M, Kapralov MV, Koch MA, Mir A, Molins A, Ribas-Carbó M, Rocha J, Galmés J. Exploring molecular evolution of Rubisco in C 3 and CAM Orchidaceae and Bromeliaceae. BMC Evol Biol 2020; 20:11. [PMID: 31969115 PMCID: PMC6977233 DOI: 10.1186/s12862-019-1551-8] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/03/2019] [Accepted: 11/29/2019] [Indexed: 01/02/2023] Open
Abstract
BACKGROUND The CO2-concentrating mechanism associated to Crassulacean acid metabolism (CAM) alters the catalytic context for Rubisco by increasing CO2 availability and provides an advantage in particular ecological conditions. We hypothesized about the existence of molecular changes linked to these particular adaptations in CAM Rubisco. We investigated molecular evolution of the Rubisco large (L-) subunit in 78 orchids and 144 bromeliads with C3 and CAM photosynthetic pathways. The sequence analyses were complemented with measurements of Rubisco kinetics in some species with contrasting photosynthetic mechanism and differing in the L-subunit sequence. RESULTS We identified potential positively selected sites and residues with signatures of co-adaptation. The implementation of a decision tree model related Rubisco specific variable sites to the leaf carbon isotopic composition of the species. Differences in the Rubisco catalytic traits found among C3 orchids and between strong CAM and C3 bromeliads suggested Rubisco had evolved in response to differing CO2 concentration. CONCLUSIONS The results revealed that the variability in the Rubisco L-subunit sequence in orchids and bromeliads is composed of coevolving sites under potential positive adaptive signal. The sequence variability was related to δ13C in orchids and bromeliads, however it could not be linked to the variability found in the kinetic properties of the studied species.
Collapse
Affiliation(s)
- Carmen Hermida-Carrera
- Research Group on Plant Biology under Mediterranean Conditions, Universitat de les Illes Balears-INAGEA, Ctra. Valldemossa km. 7.5, 07122 Palma, Illes Balears Spain
| | - Mario A. Fares
- Integrative and Systems Biology Group, Department of Abiotic Stress, Instituto de Biología Molecular y Celular de Plantas (CSIC–UPV), 46022 Valencia, Spain
- Department of Genetics, Trinity College Dublin, University of Dublin, Dublin 2, Ireland
| | - Marcel Font-Carrascosa
- Research Group on Plant Biology under Mediterranean Conditions, Universitat de les Illes Balears-INAGEA, Ctra. Valldemossa km. 7.5, 07122 Palma, Illes Balears Spain
| | - Maxim V. Kapralov
- School of Natural and Environmental Sciences, Newcastle University, Newcastle upon Tyne, NE1 7RU United Kingdom
| | - Marcus A. Koch
- Department of Biodiversity and Plant Systematics, Centre for Organismal Studies (COS) Heidelberg, Heidelberg University, Im Neuenheimer Feld 345, 9120 Heidelberg, Germany
| | - Arnau Mir
- Computational Biology and Bioinformatics Research Group, Department of Mathematics and Computer Science, Universitat de les Illes Balears, 07122 Palma, Illes Balears Spain
| | - Arántzazu Molins
- Research Group on Plant Biology under Mediterranean Conditions, Universitat de les Illes Balears-INAGEA, Ctra. Valldemossa km. 7.5, 07122 Palma, Illes Balears Spain
| | - Miquel Ribas-Carbó
- Research Group on Plant Biology under Mediterranean Conditions, Universitat de les Illes Balears-INAGEA, Ctra. Valldemossa km. 7.5, 07122 Palma, Illes Balears Spain
| | - Jairo Rocha
- Computational Biology and Bioinformatics Research Group, Department of Mathematics and Computer Science, Universitat de les Illes Balears, 07122 Palma, Illes Balears Spain
| | - Jeroni Galmés
- Research Group on Plant Biology under Mediterranean Conditions, Universitat de les Illes Balears-INAGEA, Ctra. Valldemossa km. 7.5, 07122 Palma, Illes Balears Spain
| |
Collapse
|
48
|
Chaintoutis SC, Papa A, Pervanidou D, Dovas CI. Evolutionary dynamics of lineage 2 West Nile virus in Europe, 2004–2018: Phylogeny, selection pressure and phylogeography. Mol Phylogenet Evol 2019; 141:106617. [DOI: 10.1016/j.ympev.2019.106617] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/30/2019] [Revised: 09/10/2019] [Accepted: 09/11/2019] [Indexed: 12/29/2022]
|
49
|
Wang Y, Jiao H, Jiang P, Zhao H. Functional divergence of bitter taste receptors in a nectar-feeding bird. Biol Lett 2019; 15:20190461. [PMID: 31551065 DOI: 10.1098/rsbl.2019.0461] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/02/2023] Open
Abstract
Nectar may contain many secondary metabolites that are commonly toxic and bitter-tasting. It has been hypothesized that such bitter-tasting secondary metabolites might keep the nectar exclusive to only a few pollinators. To test this hypothesis, we examined functional changes of bitter taste receptor genes (Tas2rs) in a species of nectar-feeding bird (Anna's hummingbird) by comparing these genes with those from two closely related insect-feeding species (chimney swift and chuck-will's widow). We previously identified a larger number of Tas2rs in the hummingbird than in its close insectivorous relatives. In the present study, we demonstrate higher sensitivity and new functions in the hummingbird Tas2r gene copies generated by a lineage-specific duplication, which has been shaped by positive selection. These results suggest that the bitter taste may lead to increased sensitivities and specialized abilities of the hummingbird to detect bitter-tasting nectar. Moreover, this study potentially supports the hypothesis that bitter-tasting nectar may have been specialized for some pollinators, thus enforcing plant-pollinator mutualism.
Collapse
Affiliation(s)
- Yi Wang
- Department of Ecology, Hubei Key Laboratory of Cell Homeostasis, College of Life Sciences, Wuhan University, Wuhan 430072, People's Republic of China.,Monell Chemical Senses Center, Philadelphia, PA 19104, USA
| | - Hengwu Jiao
- Department of Ecology, Hubei Key Laboratory of Cell Homeostasis, College of Life Sciences, Wuhan University, Wuhan 430072, People's Republic of China
| | - Peihua Jiang
- Monell Chemical Senses Center, Philadelphia, PA 19104, USA
| | - Huabin Zhao
- Department of Ecology, Hubei Key Laboratory of Cell Homeostasis, College of Life Sciences, Wuhan University, Wuhan 430072, People's Republic of China
| |
Collapse
|
50
|
Jacquot M, Rao PP, Yadav S, Nomikou K, Maan S, Jyothi YK, Reddy N, Putty K, Hemadri D, Singh KP, Maan NS, Hegde NR, Mertens P, Biek R. Contrasting selective patterns across the segmented genome of bluetongue virus in a global reassortment hotspot. Virus Evol 2019; 5:vez027. [PMID: 31392031 PMCID: PMC6680063 DOI: 10.1093/ve/vez027] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022] Open
Abstract
For segmented viruses, rapid genomic and phenotypic changes can occur through the process of reassortment, whereby co-infecting strains exchange entire segments creating novel progeny virus genotypes. However, for many viruses with segmented genomes, this process and its effect on transmission dynamics remain poorly understood. Here, we assessed the consequences of reassortment for selection on viral diversity through time using bluetongue virus (BTV), a segmented arbovirus that is the causative agent of a major disease of ruminants. We analysed ninety-two BTV genomes isolated across four decades from India, where BTV diversity, and thus opportunities for reassortment, are among the highest in the world. Our results point to frequent reassortment and segment turnover, some of which appear to be driven by selective sweeps and serial hitchhiking. Particularly, we found evidence for a recent selective sweep affecting segment 5 and its encoded NS1 protein that has allowed a single variant to essentially invade the full range of BTV genomic backgrounds and serotypes currently circulating in India. In contrast, diversifying selection was found to play an important role in maintaining genetic diversity in genes encoding outer surface proteins involved in virus interactions (VP2 and VP5, encoded by segments 2 and 6, respectively). Our results support the role of reassortment in driving rapid phenotypic change in segmented viruses and generate testable hypotheses for in vitro experiments aiming at understanding the specific mechanisms underlying differences in fitness and selection across viral genomes.
Collapse
Affiliation(s)
- Maude Jacquot
- Institute of Biodiversity, Animal Health and Comparative Medicine, Boyd Orr Centre for Population and Ecosystem Health, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, UK
| | - Pavuluri P Rao
- Ella Foundation, Genome Valley Hyderabad, Hyderabad, Telangana, India
| | - Sarita Yadav
- The Pirbright Institute, Pirbright, Woking, Surrey, UK
| | - Kyriaki Nomikou
- MRC-University of Glasgow Centre for Virus Research, Glasgow, UK
| | - Sushila Maan
- College of Veterinary Sciences, LLR University of Veterinary and Animal Sciences, Hisar, Haryana, India
| | - Y Krishna Jyothi
- Veterinary Biological and Research Institute, Vijayawada, Andhra Pradesh, India
| | - Narasimha Reddy
- PVNR Telangana Veterinary University, Hyderabad, Telangana, India
| | - Kalyani Putty
- PVNR Telangana Veterinary University, Hyderabad, Telangana, India
| | - Divakar Hemadri
- ICAR-National Institute of Veterinary Epidemiology and Disease Informatics, Bengaluru, Karnataka, India
| | - Karam P Singh
- Centre for Animal Disease Research and Diagnosis, Division of Pathology, ICAR-Indian Veterinary Research Institute, Izatnagar, Bareilly, Uttar Pradesh, India
| | - Narender Singh Maan
- College of Veterinary Sciences, LLR University of Veterinary and Animal Sciences, Hisar, Haryana, India
| | - Nagendra R Hegde
- Ella Foundation, Genome Valley Hyderabad, Hyderabad, Telangana, India
| | - Peter Mertens
- The Pirbright Institute, Pirbright, Woking, Surrey, UK.,The School of Veterinary Medicine and Science, University of Nottingham, Sutton Bonington, Leicestershire, UK
| | - Roman Biek
- Institute of Biodiversity, Animal Health and Comparative Medicine, Boyd Orr Centre for Population and Ecosystem Health, College of Medical, Veterinary and Life Sciences, University of Glasgow, Glasgow, UK
| |
Collapse
|