1
|
Dong Z, Wang C, Qu Q. WGCCRR: a web-based tool for genome-wide screening of convergent indels and substitutions of amino acids. BIOINFORMATICS ADVANCES 2024; 4:vbae070. [PMID: 38808070 PMCID: PMC11132816 DOI: 10.1093/bioadv/vbae070] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 08/14/2023] [Revised: 04/05/2024] [Accepted: 05/23/2024] [Indexed: 05/30/2024]
Abstract
Summary Genome-wide analyses of proteincoding gene sequences are being employed to examine the genetic basis of adaptive evolution in many organismal groups. Previous studies have revealed that convergent/parallel adaptive evolution may be caused by convergent/parallel amino acid changes. Similarly, detailed analysis of lineage-specific amino acid changes has shown correlations with certain lineage-specific traits. However, experimental validation remains the ultimate measure of causality. With the increasing availability of genomic data, a streamlined tool for such analyses would facilitate and expedite the screening of genetic loci that hold potential for adaptive evolution, while alleviating the bioinformatic burden for experimental biologists. In this study, we present a user-friendly web-based tool called WGCCRR (Whole Genome Comparative Coding Region Read) designed to screen both convergent/parallel and lineage-specific amino acid changes on a genome-wide scale. Our tool allows users to replicate previous analyses with just a few clicks, and the exported results are straightforward to interpret. In addition, we have also included amino acid indels that are usually neglected in previous work. Our website provides an efficient platform for screening candidate loci for downstream experimental tests. Availability and Implementation The tool is available at: https://fishevo.xmu.edu.cn/.
Collapse
Affiliation(s)
- Zheng Dong
- State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Xiamen University, Xià-Mén, Fú-Jiàn 361102, China
| | - Chen Wang
- State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Xiamen University, Xià-Mén, Fú-Jiàn 361102, China
| | - Qingming Qu
- State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Xiamen University, Xià-Mén, Fú-Jiàn 361102, China
| |
Collapse
|
2
|
Roh YR, Yim HS, Park K, Lee JH. Molecular characterization of positively selected genes contributing aquatic adaptation in marine mammals. Genes Genomics 2024:10.1007/s13258-023-01487-2. [PMID: 38733518 DOI: 10.1007/s13258-023-01487-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2023] [Accepted: 12/20/2023] [Indexed: 05/13/2024]
Abstract
BACKGROUND Marine mammals, which have evolved independently into three distinct lineages, share common physiological features that contribute to their adaptation to the marine environment. OBJECTIVE To identify positively selected genes (PSGs) for adaptation to the marine environment using available genomic data from three taxonomic orders: cetaceans, pinnipeds, and sirenians. METHODS Based on the genomes within each group of Artiodactyla, Carnivora and Afrotheria, we performed selection analysis using the branch-site model in CODEML. RESULTS Based on the branch-site model, 460, 614, and 359 PSGs were predicted for the cetaceans, pinnipeds, and sirenians, respectively. Functional enrichment analysis indicated that genes associated with hemostasis were positively selected across all lineages of marine mammals. We observed positive selection signals for the hemostasis and coagulation-related genes plasminogen activator, urokinase (PLAU), multimerin 1 (MMRN1), gamma-glutamyl carboxylase (GGCX), and platelet endothelial aggregation receptor 1 (PEAR1). Additionally, we found out that the sodium voltage-gated channel alpha subunit 9 (SCN9A), serine/arginine repetitive matrix 4 (SRRM4), and Ki-ras-induced actin-interacting protein (KRAP) are under positive selection pressure and are associated with cognition, neurite outgrowth, and IP3-mediated Ca2 + release, respectively. CONCLUSION This study will contribute to our understanding of the adaptive evolution of marine mammals by providing information on a group of candidate genes that are predicted to influence adaptation to aquatic environments, as well as their functional characteristics.
Collapse
Affiliation(s)
- Yoo-Rim Roh
- Marine Biotechnology Research Center, Korea Institute of Ocean Science and Technology, 385 Haeyang-ro, Yeongdo-gu, Busan, 49111, Republic of Korea
- Department of Marine Biotechnology, Korea National University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon, 34113, Republic of Korea
| | - Hyung-Soon Yim
- Marine Biotechnology Research Center, Korea Institute of Ocean Science and Technology, 385 Haeyang-ro, Yeongdo-gu, Busan, 49111, Republic of Korea
- Department of Marine Biotechnology, Korea National University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon, 34113, Republic of Korea
| | - Kiejung Park
- Cheonan Industry-Academic Collaboration Foundation, Sangmyung University, 31 Sangmyeongdae-gil, Dongnam-gu, Cheonan, 31066, Republic of Korea.
| | - Jung-Hyun Lee
- Marine Biotechnology Research Center, Korea Institute of Ocean Science and Technology, 385 Haeyang-ro, Yeongdo-gu, Busan, 49111, Republic of Korea.
- Department of Marine Biotechnology, Korea National University of Science and Technology, 217 Gajeong-ro, Yuseong-gu, Daejeon, 34113, Republic of Korea.
| |
Collapse
|
3
|
Morel M, Zhukova A, Lemoine F, Gascuel O. Accurate Detection of Convergent Mutations in Large Protein Alignments With ConDor. Genome Biol Evol 2024; 16:evae040. [PMID: 38451738 PMCID: PMC10986858 DOI: 10.1093/gbe/evae040] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/09/2023] [Revised: 01/30/2024] [Accepted: 02/22/2024] [Indexed: 03/09/2024] Open
Abstract
Evolutionary convergences are observed at all levels, from phenotype to DNA and protein sequences, and changes at these different levels tend to be correlated. Notably, convergent mutations can lead to convergent changes in phenotype, such as changes in metabolism, drug resistance, and other adaptations to changing environments. We propose a two-component approach to detect mutations subject to convergent evolution in protein alignments. The "Emergence" component selects mutations that emerge more often than expected, while the "Correlation" component selects mutations that correlate with the convergent phenotype under study. With regard to Emergence, a phylogeny deduced from the alignment is provided by the user and is used to simulate the evolution of each alignment position. These simulations allow us to estimate the expected number of mutations in a neutral model, which is compared to the observed number of mutations in the data studied. In Correlation, a comparative phylogenetic approach, is used to measure whether the presence of each of the observed mutations is correlated with the convergent phenotype. Each component can be used on its own, for example Emergence when no phenotype is available. Our method is implemented in a standalone workflow and a webserver, called ConDor. We evaluate the properties of ConDor using simulated data, and we apply it to three real datasets: sedge PEPC proteins, HIV reverse transcriptase, and fish rhodopsin. The results show that the two components of ConDor complement each other, with an overall accuracy that compares favorably to other available tools, especially on large datasets.
Collapse
Affiliation(s)
- Marie Morel
- Institut Pasteur, Université Paris Cité, Unité Bioinformatique Evolutive, Paris, France
- Université Claude Bernard Lyon 1, LBBE, UMR 5558, CNRS, VAS, Villeurbanne, 69100, France
| | - Anna Zhukova
- Institut Pasteur, Université Paris Cité, Unité Bioinformatique Evolutive, Paris, France
- Institut Pasteur, Université Paris Cité, Bioinformatics and Biostatistics Hub, Paris, France
| | - Frédéric Lemoine
- Institut Pasteur, Université Paris Cité, Unité Bioinformatique Evolutive, Paris, France
- Institut Pasteur, Université Paris Cité, Bioinformatics and Biostatistics Hub, Paris, France
- Institut Pasteur, Université Paris Cité, CNR Virus Des Infections Respiratoires, Paris, France
| | - Olivier Gascuel
- Institut Pasteur, Université Paris Cité, Unité Bioinformatique Evolutive, Paris, France
- Institut de Systématique, Evolution, Biodiversité (UMR 7205—CNRS, Muséum National d’Histoire Naturelle, SU, EPHE, UA), Paris, France
| |
Collapse
|
4
|
Redlich R, Kowalczyk A, Tene M, Sestili HH, Foley K, Saputra E, Clark N, Chikina M, Meyer WK, Pfenning A. RERconverge Expansion: Using Relative Evolutionary Rates to Study Complex Categorical Trait Evolution. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.12.06.570425. [PMID: 38106136 PMCID: PMC10723433 DOI: 10.1101/2023.12.06.570425] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 12/19/2023]
Abstract
Comparative genomics approaches seek to associate evolutionary genetic changes with the evolution of phenotypes across a phylogeny. Many of these methods, including our evolutionary rates based method, RERconverge, lack the capability of analyzing non-ordinal, multicategorical traits. To address this limitation, we introduce an expansion to RERconverge that associates shifts in evolutionary rates with the convergent evolution of multi-categorical traits. The categorical RERconverge expansion includes methods for performing categorical ancestral state reconstruction, statistical tests for associating relative evolutionary rates with categorical variables, and a new method for performing phylogenetic permulations on multi-categorical traits. In addition to demonstrating our new method on a three-category diet phenotype, we compare its performance to naive pairwise binary RERconverge analyses and two existing methods for comparative genomic analyses of categorical traits: phylogenetic simulations and a phylogenetic signal based method. We also present a diagnostic analysis of the new permulations approach demonstrating how the method scales with the number of species and the number of categories included in the analysis. Our results show that our new categorical method outperforms phylogenetic simulations at identifying genes and enriched pathways significantly associated with the diet phenotype and that the new ancestral reconstruction drives an improvement in our ability to capture diet-related enriched pathways. Our categorical permulations were able to account for non-uniform null distributions and correct for non-independence in gene rank during pathway enrichment analysis. The categorical expansion to RERconverge will provide a strong foundation for applying the comparative method to categorical traits on larger data sets with more species and more complex trait evolution.
Collapse
|
5
|
Rahi ML, Mather PB, de Bello Cioffi M, Ezaz T, Hurwood DA. Genomic Basis of Freshwater Adaptation in the Palaemonid Prawn Genus Macrobrachium: Convergent Evolution Following Multiple Independent Colonization Events. J Mol Evol 2023; 91:976-989. [PMID: 38010517 DOI: 10.1007/s00239-023-10149-6] [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: 10/14/2022] [Accepted: 11/14/2023] [Indexed: 11/29/2023]
Abstract
Adaptation to different salinity environments can enhance morphological and genomic divergence between related aquatic taxa. Species of prawns in the genus Macrobrachium naturally inhabit different osmotic niches and possess distinctive lifecycle traits associated with salinity tolerance. This study was conducted to investigate the patterns of adaptive genomic divergence during freshwater colonization in 34 Macrobrachium species collected from four continents; Australia, Asia, North and South America. Genotyping-by-sequencing (GBS) technique identified 5018 loci containing 82,636 single nucleotide polymorphisms (SNPs) that were used to reconstruct a phylogenomic tree. An additional phylogeny was reconstructed based on 43 candidate genes, previously identified as being potentially associated with freshwater adaptation. Comparison of the two phylogenetic trees revealed contrasting topologies. The GBS tree indicated multiple independent continent-specific invasions into freshwater by Macrobrachium lineages following common marine ancestry, as species with abbreviated larval development (ALD), i.e., species having a full freshwater life history, appeared reciprocally monophyletic within each continent. In contrast, the candidate gene tree showed convergent evolution for all ALD species worldwide, forming a single, well-supported clade. This latter pattern is likely the result of common evolutionary pressures selecting key mutations favored in continental freshwater habitats Results suggest that following multiple independent invasions into continental freshwaters at different evolutionary timescales, Macrobrachium taxa experienced adaptive genomic divergence, and in particular, convergence in the same genomic regions with parallel shifts in specific conserved phenotypic traits, such as evolution of larger eggs with abbreviated larval developmental.
Collapse
Affiliation(s)
- Md Lifat Rahi
- Fisheries and Marine Resource Technology Discipline, Khulna University, Khulna, Bangladesh.
| | - Peter B Mather
- Faculty of Science, Queensland University of Technology (QUT), Brisbane, QLD, 4001, Australia
| | - Marcelo de Bello Cioffi
- Department of Genetics and Evolution, Federal University of Sao Carlos, São Carlos, SP, Brazil
| | - Tariq Ezaz
- Institute for Applied Ecology (IAE), University of Canberra (UC), Canberra, ACT, 2617, Australia
| | - David A Hurwood
- Faculty of Science, Queensland University of Technology (QUT), Brisbane, QLD, 4001, Australia
| |
Collapse
|
6
|
Wu B, Gao X, Hu M, Hu J, Lan T, Xue T, Xu W, Zhu C, Yuan Y, Zheng J, Qin T, Xin P, Li Y, Gong L, Feng C, He S, Liu H, Li H, Wang Q, Ma Z, Qiu Q, Wang K. Distinct and shared endothermic strategies in the heat producing tissues of tuna and other teleosts. SCIENCE CHINA. LIFE SCIENCES 2023; 66:2629-2645. [PMID: 37273070 DOI: 10.1007/s11427-022-2312-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/06/2022] [Accepted: 02/28/2023] [Indexed: 06/06/2023]
Abstract
Although most fishes are ectothermic, some, including tuna and billfish, achieve endothermy through specialized heat producing tissues that are modified muscles. How these heat producing tissues evolved, and whether they share convergent molecular mechanisms, remain unresolved. Here, we generated a high-quality genome from the mackerel tuna (Euthynnus affinis) and investigated the heat producing tissues of this fish by single-nucleus and bulk RNA sequencing. Compared with other teleosts, tuna-specific genetic variation is strongly associated with muscle differentiation. Single-nucleus RNA-seq revealed a high proportion of specific slow skeletal muscle cell subtypes in the heat producing tissues of tuna. Marker genes of this cell subtype are associated with the relative sliding of actin and myosin, suggesting that tuna endothermy is mainly based on shivering thermogenesis. In contrast, cross-species transcriptome analysis indicated that endothermy in billfish relies mainly on non-shivering thermogenesis. Nevertheless, the heat producing tissues of the different species do share some tissue-specific genes, including vascular-related and mitochondrial genes. Overall, although tunas and billfishes differ in their thermogenic strategies, they share similar expression patterns in some respects, highlighting the complexity of convergent evolution.
Collapse
Affiliation(s)
- Baosheng Wu
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Xueli Gao
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Mingling Hu
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Jing Hu
- South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou, 510300, China
| | - Tianming Lan
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, 518083, China
- BGI Life Science Joint Research Center, Northeast Forestry University, Harbin, 150006, China
| | - Tingfeng Xue
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Wenjie Xu
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Chenglong Zhu
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Yuan Yuan
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Jiangmin Zheng
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Tao Qin
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Peidong Xin
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Ye Li
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Li Gong
- National Engineering Laboratory of Marine Germplasm Resources Exploration and Utilization, Marine Science and Technology College, Zhejiang Ocean University, Zhoushan, 316022, China
| | - Chenguang Feng
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China
| | - Shunping He
- Institute of Deep-Sea Science and Engineering, Chinese Academy of Sciences, Sanya, 572000, China
- The Key Laboratory of Aquatic Biodiversity and Conservation of Chinese Academy of Sciences, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, 430072, China
| | - Huan Liu
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, 518083, China
- BGI Life Science Joint Research Center, Northeast Forestry University, Harbin, 150006, China
| | - Haimeng Li
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, 518083, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Qing Wang
- State Key Laboratory of Agricultural Genomics, BGI-Shenzhen, Shenzhen, 518083, China
- College of Life Sciences, University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Zhenhua Ma
- South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences, Guangzhou, 510300, China.
| | - Qiang Qiu
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China.
| | - Kun Wang
- School of Ecology and Environment, Northwestern Polytechnical University, Xi'an, 710072, China.
| |
Collapse
|
7
|
Herbert AL, Allard CAH, McCoy MJ, Wucherpfennig JI, Krueger SP, Chen HI, Gourlay AN, Jackson KD, Abbo LA, Bennett SH, Sears JD, Rhyne AL, Bellono NW, Kingsley DM. The genetic basis of novel trait gain in walking fish. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.10.14.562356. [PMID: 37873105 PMCID: PMC10592820 DOI: 10.1101/2023.10.14.562356] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 10/25/2023]
Abstract
A major goal in biology is to understand how organisms evolve novel traits. Multiple studies have identified genes contributing to regressive evolution, the loss of structures that existed in a recent ancestor. However, fewer examples exist for genes underlying constructive evolution, the gain of novel structures and capabilities in lineages that previously lacked them. Sea robins are fish that have evolved enlarged pectoral fins, six mobile locomotory fin rays (legs) and six novel macroscopic lobes in the central nervous system (CNS) that innervate the corresponding legs. Here, we establish successful husbandry and use a combination of transcriptomics, CRISPR-Cas9 editing, and behavioral assays to identify key transcription factors that are required for leg formation and function in sea robins. We also generate hybrids between two sea robin species with distinct leg morphologies and use allele-specific expression analysis and gene editing to explore the genetic basis of species-specific trait diversity, including a novel sensory gain of function. Collectively, our study establishes sea robins as a new model for studying the genetic basis of novel organ formation, and demonstrates a crucial role for the conserved limb gene tbx3a in the evolution of chemosensory legs in walking fish.
Collapse
Affiliation(s)
- Amy L Herbert
- Department of Developmental Biology, Stanford University School of Medicine, Stanford CA 94305 USA
| | - Corey AH Allard
- Department of Molecular and Cellular Biology, Harvard University, Cambridge MA 02138 USA
| | - Matthew J McCoy
- Department of Pathology, Stanford University School of Medicine, Stanford CA 94305 USA
| | - Julia I Wucherpfennig
- Department of Developmental Biology, Stanford University School of Medicine, Stanford CA 94305 USA
| | - Stephanie P Krueger
- Department of Molecular and Cellular Biology, Harvard University, Cambridge MA 02138 USA
| | - Heidi I Chen
- Department of Developmental Biology, Stanford University School of Medicine, Stanford CA 94305 USA
| | | | - Kohle D Jackson
- Department of Developmental Biology, Stanford University School of Medicine, Stanford CA 94305 USA
| | - Lisa A Abbo
- Marine Biological Laboratory, Woods Hole, MA, 02543 USA
| | | | | | | | - Nicholas W Bellono
- Department of Molecular and Cellular Biology, Harvard University, Cambridge MA 02138 USA
| | - David M Kingsley
- Howard Hughes Medical Institute, Stanford University School of Medicine, Stanford CA 94305 USA
| |
Collapse
|
8
|
Chen HI, Turakhia Y, Bejerano G, Kingsley DM. Whole-genome Comparisons Identify Repeated Regulatory Changes Underlying Convergent Appendage Evolution in Diverse Fish Lineages. Mol Biol Evol 2023; 40:msad188. [PMID: 37739926 PMCID: PMC10516590 DOI: 10.1093/molbev/msad188] [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] [Indexed: 09/24/2023] Open
Abstract
Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 percomorph fish species that span over 100 million years of evolution and either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA, a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, syngnathid, and synbranchid species with caudal fin reduction, and alters caudal fin development when targeted by genome editing. Our study illustrates a broadly applicable strategy for mapping phenotypes to genotypes across a tree of vertebrate species and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes across 100 million years of fish evolution.
Collapse
Affiliation(s)
- Heidi I Chen
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
| | - Yatish Turakhia
- Department of Electrical and Computer Engineering, University of California, San Diego, CA, USA
| | - Gill Bejerano
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
- Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, CA, USA
- Department of Computer Science, Stanford University School of Engineering, Stanford, CA, USA
- Department of Pediatrics, Stanford University School of Medicine, Stanford, CA, USA
| | - David M Kingsley
- Department of Developmental Biology, Stanford University School of Medicine, Stanford, CA, USA
- Howard Hughes Medical Institute, Stanford University, Stanford, CA, USA
| |
Collapse
|
9
|
Tong H, Shao G, Wang L, Li J, Wang T, Zhang L, Lv Y, Ye F, Fu C, Jin Y. Association of a single amino acid replacement with dorsal pigmentation in a lizard from the Qinghai-Tibetan Plateau. Int J Biol Macromol 2023; 242:124907. [PMID: 37230451 DOI: 10.1016/j.ijbiomac.2023.124907] [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: 02/15/2023] [Revised: 05/08/2023] [Accepted: 05/13/2023] [Indexed: 05/27/2023]
Abstract
Reptiles can evolve adaptive colors in different environments, but relatively little is known about the genetic mechanisms. Here, we identified the MC1R gene and its association with intraspecific color variation in the lizard Phrynocephalus erythrurus. Analysis of the MC1R sequence in 143 individuals from dark South Qiangtang Plateau (SQP) and light North Qiangtang plateau (NQP) populations, revealed two amino acid sites that showed significant differences in frequency between two areas. One SNP, corresponding to Glu183Lys residue, was found to be a highly significant outlier and differentially fixed for SQP and NQP populations. This residue is located in an extracellular area in the second small extracellular loop within the secondary structure of MC1R, which represents an "attachment pocket" part of the 3D structure. Cytological expression of MC1R alleles with the Glu183Lys replacement showed a 39 % increase in intracellular agonist-induced cyclic AMP levels and a 23.18 % greater cell surface expression of MC1R protein in the SQP relative to the NQP allele. Further in silico 3D modeling and in vitro binding experiments indicated a higher MC1R-α-MSH binding for the SQP allele, and elevated melanin synthesis. We provide an overview of how a single amino acid replacement leads to fundamental changes in MC1R function, and hence shapes variation in dorsal pigmentation in lizards from different environments.
Collapse
Affiliation(s)
- Haojie Tong
- College of Life Sciences, China Jiliang University, Hangzhou 310018, PR China
| | - Gang Shao
- College of Life Sciences, China Jiliang University, Hangzhou 310018, PR China
| | - Leijie Wang
- College of Life Sciences, China Jiliang University, Hangzhou 310018, PR China
| | - Jiasheng Li
- College of Life Sciences, China Jiliang University, Hangzhou 310018, PR China
| | - Tao Wang
- Zhejiang Provincial Key Laboratory of Silkworm Bioreactor and Biomedicine, College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China
| | - Lun Zhang
- Zhejiang Provincial Key Laboratory of Silkworm Bioreactor and Biomedicine, College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China
| | - Yudie Lv
- College of Life Sciences, China Jiliang University, Hangzhou 310018, PR China
| | - Fei Ye
- Zhejiang Provincial Key Laboratory of Silkworm Bioreactor and Biomedicine, College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China
| | - Caiyun Fu
- Zhejiang Provincial Key Laboratory of Silkworm Bioreactor and Biomedicine, College of Life Sciences and Medicine, Zhejiang Sci-Tech University, Hangzhou 310018, China
| | - Yuanting Jin
- College of Life Sciences, China Jiliang University, Hangzhou 310018, PR China.
| |
Collapse
|
10
|
Christmas MJ, Kaplow IM, Genereux DP, Dong MX, Hughes GM, Li X, Sullivan PF, Hindle AG, Andrews G, Armstrong JC, Bianchi M, Breit AM, Diekhans M, Fanter C, Foley NM, Goodman DB, Goodman L, Keough KC, Kirilenko B, Kowalczyk A, Lawless C, Lind AL, Meadows JRS, Moreira LR, Redlich RW, Ryan L, Swofford R, Valenzuela A, Wagner F, Wallerman O, Brown AR, Damas J, Fan K, Gatesy J, Grimshaw J, Johnson J, Kozyrev SV, Lawler AJ, Marinescu VD, Morrill KM, Osmanski A, Paulat NS, Phan BN, Reilly SK, Schäffer DE, Steiner C, Supple MA, Wilder AP, Wirthlin ME, Xue JR, Birren BW, Gazal S, Hubley RM, Koepfli KP, Marques-Bonet T, Meyer WK, Nweeia M, Sabeti PC, Shapiro B, Smit AFA, Springer MS, Teeling EC, Weng Z, Hiller M, Levesque DL, Lewin HA, Murphy WJ, Navarro A, Paten B, Pollard KS, Ray DA, Ruf I, Ryder OA, Pfenning AR, Lindblad-Toh K, Karlsson EK. Evolutionary constraint and innovation across hundreds of placental mammals. Science 2023; 380:eabn3943. [PMID: 37104599 PMCID: PMC10250106 DOI: 10.1126/science.abn3943] [Citation(s) in RCA: 45] [Impact Index Per Article: 45.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/23/2021] [Accepted: 12/16/2022] [Indexed: 04/29/2023]
Abstract
Zoonomia is the largest comparative genomics resource for mammals produced to date. By aligning genomes for 240 species, we identify bases that, when mutated, are likely to affect fitness and alter disease risk. At least 332 million bases (~10.7%) in the human genome are unusually conserved across species (evolutionarily constrained) relative to neutrally evolving repeats, and 4552 ultraconserved elements are nearly perfectly conserved. Of 101 million significantly constrained single bases, 80% are outside protein-coding exons and half have no functional annotations in the Encyclopedia of DNA Elements (ENCODE) resource. Changes in genes and regulatory elements are associated with exceptional mammalian traits, such as hibernation, that could inform therapeutic development. Earth's vast and imperiled biodiversity offers distinctive power for identifying genetic variants that affect genome function and organismal phenotypes.
Collapse
Affiliation(s)
- Matthew J. Christmas
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Irene M. Kaplow
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | | | - Michael X. Dong
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Graham M. Hughes
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Xue Li
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Morningside Graduate School of Biomedical Sciences, UMass Chan Medical School, Worcester, MA 01605, USA
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
| | - Patrick F. Sullivan
- Department of Genetics, University of North Carolina Medical School, Chapel Hill, NC 27599, USA
- Department of Medical Epidemiology and Biostatistics, Karolinska Institutet, Stockholm, Sweden
| | - Allyson G. Hindle
- School of Life Sciences, University of Nevada Las Vegas, Las Vegas, NV 89154, USA
| | - Gregory Andrews
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
| | - Joel C. Armstrong
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Matteo Bianchi
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Ana M. Breit
- School of Biology and Ecology, University of Maine, Orono, ME 04469, USA
| | - Mark Diekhans
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Cornelia Fanter
- School of Life Sciences, University of Nevada Las Vegas, Las Vegas, NV 89154, USA
| | - Nicole M. Foley
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX 77843, USA
| | - Daniel B. Goodman
- Department of Microbiology and Immunology, University of California San Francisco, San Francisco, CA 94143, USA
| | | | - Kathleen C. Keough
- Fauna Bio, Inc., Emeryville, CA 94608, USA
- Department of Epidemiology and Biostatistics, University of California San Francisco, San Francisco, CA 94158, USA
- Gladstone Institutes, San Francisco, CA 94158, USA
| | - Bogdan Kirilenko
- Faculty of Biosciences, Goethe-University, 60438 Frankfurt, Germany
- LOEWE Centre for Translational Biodiversity Genomics, 60325 Frankfurt, Germany
- Senckenberg Research Institute, 60325 Frankfurt, Germany
| | - Amanda Kowalczyk
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Colleen Lawless
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Abigail L. Lind
- Department of Epidemiology and Biostatistics, University of California San Francisco, San Francisco, CA 94158, USA
- Gladstone Institutes, San Francisco, CA 94158, USA
| | - Jennifer R. S. Meadows
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Lucas R. Moreira
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
| | - Ruby W. Redlich
- Department of Biological Sciences, Mellon College of Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Louise Ryan
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Ross Swofford
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
| | - Alejandro Valenzuela
- Department of Experimental and Health Sciences, Institute of Evolutionary Biology (UPF-CSIC), Universitat Pompeu Fabra, 08003 Barcelona, Spain
| | - Franziska Wagner
- Museum of Zoology, Senckenberg Natural History Collections Dresden, 01109 Dresden, Germany
| | - Ola Wallerman
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Ashley R. Brown
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Joana Damas
- The Genome Center, University of California Davis, Davis, CA 95616, USA
| | - Kaili Fan
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
| | - John Gatesy
- Division of Vertebrate Zoology, American Museum of Natural History, New York, NY 10024, USA
| | - Jenna Grimshaw
- Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409, USA
| | - Jeremy Johnson
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
| | - Sergey V. Kozyrev
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Alyssa J. Lawler
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Department of Biological Sciences, Mellon College of Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Voichita D. Marinescu
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
| | - Kathleen M. Morrill
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Morningside Graduate School of Biomedical Sciences, UMass Chan Medical School, Worcester, MA 01605, USA
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
| | - Austin Osmanski
- Medical Scientist Training Program, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Nicole S. Paulat
- Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409, USA
| | - BaDoi N. Phan
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Medical Scientist Training Program, University of Pittsburgh School of Medicine, Pittsburgh, PA 15261, USA
| | - Steven K. Reilly
- Department of Genetics, Yale School of Medicine, New Haven, CT 06510, USA
| | - Daniel E. Schäffer
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Cynthia Steiner
- Conservation Genetics, San Diego Zoo Wildlife Alliance, Escondido, CA 92027, USA
| | - Megan A. Supple
- Department of Ecology and Evolutionary Biology, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Aryn P. Wilder
- Conservation Genetics, San Diego Zoo Wildlife Alliance, Escondido, CA 92027, USA
| | - Morgan E. Wirthlin
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Allen Institute for Brain Science, Seattle, WA 98109, USA
| | - James R. Xue
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | | | - Bruce W. Birren
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
| | - Steven Gazal
- Keck School of Medicine, University of Southern California, Los Angeles, CA 90033, USA
| | | | - Klaus-Peter Koepfli
- Center for Species Survival, Smithsonian’s National Zoo and Conservation Biology Institute, Washington, DC 20008, USA
- Computer Technologies Laboratory, ITMO University, St. Petersburg 197101, Russia
- Smithsonian-Mason School of Conservation, George Mason University, Front Royal, VA 22630, USA
| | - Tomas Marques-Bonet
- Catalan Institution of Research and Advanced Studies (ICREA), 08010 Barcelona, Spain
- CNAG-CRG, Centre for Genomic Regulation, Barcelona Institute of Science and Technology (BIST), 08036 Barcelona, Spain
- Department of Medicine and Life Sciences, Institute of Evolutionary Biology (UPF-CSIC), Universitat Pompeu Fabra, 08003 Barcelona, Spain
- Institut Català de Paleontologia Miquel Crusafont, Universitat Autònoma de Barcelona, 08193 Cerdanyola del Vallès, Barcelona, Spain
| | - Wynn K. Meyer
- Department of Biological Sciences, Lehigh University, Bethlehem, PA 18015, USA
| | - Martin Nweeia
- Department of Comprehensive Care, School of Dental Medicine, Case Western Reserve University, Cleveland, OH 44106, USA
- Department of Vertebrate Zoology, Canadian Museum of Nature, Ottawa, Ontario K2P 2R1, Canada
- Department of Vertebrate Zoology, Smithsonian Institution, Washington, DC 20002, USA
- Narwhal Genome Initiative, Department of Restorative Dentistry and Biomaterials Sciences, Harvard School of Dental Medicine, Boston, MA 02115, USA
| | - Pardis C. Sabeti
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
- Howard Hughes Medical Institute, Harvard University, Cambridge, MA 02138, USA
| | - Beth Shapiro
- Department of Ecology and Evolutionary Biology, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Howard Hughes Medical Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | | | - Mark S. Springer
- Department of Evolution, Ecology and Organismal Biology, University of California Riverside, Riverside, CA 92521, USA
| | - Emma C. Teeling
- School of Biology and Environmental Science, University College Dublin, Belfield, Dublin 4, Ireland
| | - Zhiping Weng
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
| | - Michael Hiller
- Faculty of Biosciences, Goethe-University, 60438 Frankfurt, Germany
- LOEWE Centre for Translational Biodiversity Genomics, 60325 Frankfurt, Germany
- Senckenberg Research Institute, 60325 Frankfurt, Germany
| | | | - Harris A. Lewin
- The Genome Center, University of California Davis, Davis, CA 95616, USA
- Department of Evolution and Ecology, University of California Davis, Davis, CA 95616, USA
- John Muir Institute for the Environment, University of California Davis, Davis, CA 95616, USA
| | - William J. Murphy
- Veterinary Integrative Biosciences, Texas A&M University, College Station, TX 77843, USA
| | - Arcadi Navarro
- Catalan Institution of Research and Advanced Studies (ICREA), 08010 Barcelona, Spain
- Department of Medicine and Life Sciences, Institute of Evolutionary Biology (UPF-CSIC), Universitat Pompeu Fabra, 08003 Barcelona, Spain
- BarcelonaBeta Brain Research Center, Pasqual Maragall Foundation, 08005 Barcelona, Spain
- CRG, Centre for Genomic Regulation, Barcelona Institute of Science and Technology (BIST), 08003 Barcelona, Spain
| | - Benedict Paten
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Katherine S. Pollard
- Department of Epidemiology and Biostatistics, University of California San Francisco, San Francisco, CA 94158, USA
- Gladstone Institutes, San Francisco, CA 94158, USA
- Chan Zuckerberg Biohub, San Francisco, CA 94158, USA
| | - David A. Ray
- Department of Biological Sciences, Texas Tech University, Lubbock, TX 79409, USA
| | - Irina Ruf
- Division of Messel Research and Mammalogy, Senckenberg Research Institute and Natural History Museum Frankfurt, 60325 Frankfurt am Main, Germany
| | - Oliver A. Ryder
- Conservation Genetics, San Diego Zoo Wildlife Alliance, Escondido, CA 92027, USA
- Department of Evolution, Behavior and Ecology, School of Biological Sciences, University of California San Diego, La Jolla, CA 92039, USA
| | - Andreas R. Pfenning
- Department of Computational Biology, School of Computer Science, Carnegie Mellon University, Pittsburgh, PA 15213, USA
- Neuroscience Institute, Carnegie Mellon University, Pittsburgh, PA 15213, USA
| | - Kerstin Lindblad-Toh
- Department of Medical Biochemistry and Microbiology, Science for Life Laboratory, Uppsala University, 751 32 Uppsala, Sweden
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
| | - Elinor K. Karlsson
- Broad Institute of MIT and Harvard, Cambridge, MA 02139, USA
- Program in Bioinformatics and Integrative Biology, UMass Chan Medical School, Worcester, MA 01605, USA
- Program in Molecular Medicine, UMass Chan Medical School, Worcester, MA 01605, USA
| |
Collapse
|
11
|
Lipovsek M, Elgoyhen AB. The evolutionary tuning of hearing. Trends Neurosci 2023; 46:110-123. [PMID: 36621369 DOI: 10.1016/j.tins.2022.12.002] [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: 09/12/2022] [Revised: 11/29/2022] [Accepted: 12/06/2022] [Indexed: 01/08/2023]
Abstract
After the transition to life on land, tympanic middle ears emerged separately in different groups of tetrapods, facilitating the efficient detection of airborne sounds and paving the way for high frequency sensitivity. The processes that brought about high-frequency hearing in mammals are tightly linked to the accumulation of coding sequence changes in inner ear genes; many of which were selected during evolution. These include proteins involved in hair bundle morphology, mechanotransduction and high endolymphatic potential, somatic electromotility for sound amplification, ribbon synapses for high-fidelity transmission of sound stimuli, and efferent synapses for the modulation of sound amplification. Here, we review the molecular evolutionary processes behind auditory functional innovation. Overall, the evidence to date supports the hypothesis that changes in inner ear proteins were central to the fine tuning of mammalian hearing.
Collapse
Affiliation(s)
- Marcela Lipovsek
- Ear Institute, Faculty of Brain Sciences, University College London, London, UK.
| | - Ana Belén Elgoyhen
- Instituto de Investigaciones en Ingeniería Genética y Biología Molecular "Dr. Héctor N. Torres" (INGEBI), Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET), Buenos Aires, Argentina.
| |
Collapse
|
12
|
Chen HI, Turakhia Y, Bejerano G, Kingsley DM. Whole-genome comparisons identify repeated regulatory changes underlying convergent appendage evolution in diverse fish lineages. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.01.30.526059. [PMID: 36778215 PMCID: PMC9915506 DOI: 10.1101/2023.01.30.526059] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 02/04/2023]
Abstract
Fins are major functional appendages of fish that have been repeatedly modified in different lineages. To search for genomic changes underlying natural fin diversity, we compared the genomes of 36 wild fish species that either have complete or reduced pelvic and caudal fins. We identify 1,614 genomic regions that are well-conserved in fin-complete species but missing from multiple fin-reduced lineages. Recurrent deletions of conserved sequences (CONDELs) in wild fin-reduced species are enriched for functions related to appendage development, suggesting that convergent fin reduction at the organismal level is associated with repeated genomic deletions near fin-appendage development genes. We used sequencing and functional enhancer assays to confirm that PelA , a Pitx1 enhancer previously linked to recurrent pelvic loss in sticklebacks, has also been independently deleted and may have contributed to the fin morphology in distantly related pelvic-reduced species. We also identify a novel enhancer that is conserved in the majority of percomorphs, drives caudal fin expression in transgenic stickleback, is missing in tetraodontiform, s yngnathid, and synbranchid species with caudal fin reduction, and which alters caudal fin development when targeted by genome editing. Our study illustrates a general strategy for mapping phenotypes to genotypes across a tree of vertebrate species, and highlights notable new examples of regulatory genomic hotspots that have been used to evolve recurrent phenotypes during 100 million years of fish evolution.
Collapse
Affiliation(s)
- Heidi I. Chen
- Department of Developmental Biology, Stanford University School of Medicine, CA
| | - Yatish Turakhia
- Department of Electrical and Computer Engineering, University of California, San Diego, San Diego, CA
| | - Gill Bejerano
- Department of Developmental Biology, Stanford University School of Medicine, CA
- Department of Biomedical Data Science, Stanford University School of Medicine, CA
- Department of Computer Science, Stanford University School of Engineering, CA
- Department of Pediatrics, Stanford University School of Medicine, CA
| | - David M. Kingsley
- Department of Developmental Biology, Stanford University School of Medicine, CA
- Howard Hughes Medical Institute, Stanford University, CA
| |
Collapse
|
13
|
Duchemin L, Lanore V, Veber P, Boussau B. Evaluation of Methods to Detect Shifts in Directional Selection at the Genome Scale. Mol Biol Evol 2022; 40:6889995. [PMID: 36510704 PMCID: PMC9940701 DOI: 10.1093/molbev/msac247] [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] [Received: 03/16/2022] [Revised: 10/24/2022] [Accepted: 10/26/2022] [Indexed: 12/15/2022] Open
Abstract
Identifying the footprints of selection in coding sequences can inform about the importance and function of individual sites. Analyses of the ratio of nonsynonymous to synonymous substitutions (dN/dS) have been widely used to pinpoint changes in the intensity of selection, but cannot distinguish them from changes in the direction of selection, that is, changes in the fitness of specific amino acids at a given position. A few methods that rely on amino-acid profiles to detect changes in directional selection have been designed, but their performances have not been well characterized. In this paper, we investigate the performance of six of these methods. We evaluate them on simulations along empirical phylogenies in which transition events have been annotated and compare their ability to detect sites that have undergone changes in the direction or intensity of selection to that of a widely used dN/dS approach, codeml's branch-site model A. We show that all methods have reduced performance in the presence of biased gene conversion but not CpG hypermutability. The best profile method, Pelican, a new implementation of Tamuri AU, Hay AJ, Goldstein RA. (2009. Identifying changes in selective constraints: host shifts in influenza. PLoS Comput Biol. 5(11):e1000564), performs as well as codeml in a range of conditions except for detecting relaxations of selection, and performs better when tree length increases, or in the presence of persistent positive selection. It is fast, enabling genome-scale searches for site-wise changes in the direction of selection associated with phenotypic changes.
Collapse
Affiliation(s)
| | - Vincent Lanore
- Laboratoire de Biométrie et Biologie Evolutive, Univ Lyon, Univ Lyon 1, CNRS, VetAgro Sup, UMR5558, Villeurbanne, France
| | - Philippe Veber
- Laboratoire de Biométrie et Biologie Evolutive, Univ Lyon, Univ Lyon 1, CNRS, VetAgro Sup, UMR5558, Villeurbanne, France
| | | |
Collapse
|
14
|
Rosvall KA. Evolutionary endocrinology and the problem of Darwin's tangled bank. Horm Behav 2022; 146:105246. [PMID: 36029721 DOI: 10.1016/j.yhbeh.2022.105246] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/12/2021] [Revised: 06/20/2022] [Accepted: 08/10/2022] [Indexed: 11/04/2022]
Abstract
Like Darwin's tangled bank of biodiversity, the endocrine mechanisms that give rise to phenotypic diversity also exhibit nearly endless forms. This tangled bank of mechanistic diversity can prove problematic as we seek general principles on the role of endocrine mechanisms in phenotypic evolution. A key unresolved question is therefore: to what degree are specific endocrine mechanisms re-used to bring about replicated phenotypic evolution? Related areas of inquiry are booming in molecular ecology, but behavioral traits are underrepresented in this literature. Here, I leverage the rich comparative tradition in evolutionary endocrinology to evaluate whether and how certain mechanisms may be repeated hotspots of behavioral evolutionary change. At one extreme, mechanisms may be parallel, such that evolution repeatedly uses the same gene or pathway to arrive at multiple independent (or, convergent) origins of a particular behavioral trait. At the other extreme, the building blocks of behavior may be unique, such that outwardly similar phenotypes are generated via lineage-specific mechanisms. This review synthesizes existing case studies, phylogenetic analyses, and experimental evolutionary research on mechanistic parallelism in animal behavior. These examples show that the endocrine building blocks of behavior have some elements of parallelism across replicated evolutionary events. However, support for parallelism is variable among studies, at least some of which relates to the level of complexity at which we consider sameness (i.e. pathway vs. gene level). Moving forward, we need continued experimentation and better testing of neutral models to understand whether, how - and critically, why - mechanism A is used in one lineage and mechanism B is used in another. We also need continued growth of large-scale comparative analyses, especially those that can evaluate which endocrine parameters are more or less likely to undergo parallel evolution alongside specific behavioral traits. These efforts will ultimately deepen understanding of how and why hormone-mediated behaviors are constructed the way that they are.
Collapse
Affiliation(s)
- Kimberly A Rosvall
- Indiana University, Bloomington, USA; Department of Biology, USA; Center for the Integrative Study of Animal Behavior, USA.
| |
Collapse
|
15
|
Bittner NKJ, Mack KL, Nachman MW. Shared Patterns of Gene Expression and Protein Evolution Associated with Adaptation to Desert Environments in Rodents. Genome Biol Evol 2022; 14:6765154. [PMID: 36268582 PMCID: PMC9648513 DOI: 10.1093/gbe/evac155] [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: 10/17/2022] [Indexed: 01/18/2023] Open
Abstract
Desert specialization has arisen multiple times across rodents and is often associated with a suite of convergent phenotypes, including modification of the kidneys to mitigate water loss. However, the extent to which phenotypic convergence in desert rodents is mirrored at the molecular level is unknown. Here, we sequenced kidney mRNA and assembled transcriptomes for three pairs of rodent species to search for shared differences in gene expression and amino acid sequence associated with adaptation to deserts. We conducted phylogenetically independent comparisons between a desert specialist and a non-desert relative in three families representing ∼70 million years of evolution. Overall, patterns of gene expression faithfully recapitulated the phylogeny of these six taxa providing a strong evolutionary signal in levels of mRNA abundance. We also found that 8.6% of all genes showed shared patterns of expression divergence between desert and non-desert taxa, much of which likely reflects convergent evolution, and representing more than expected by chance under a model of independent gene evolution. In addition to these shared changes, we observed many species-pair-specific changes in gene expression indicating that instances of adaptation to deserts include a combination of unique and shared changes. Patterns of protein evolution revealed a small number of genes showing evidence of positive selection, the majority of which did not show shared changes in gene expression. Overall, our results suggest that convergent changes in gene regulation play an important role in the complex trait of desert adaptation in rodents.
Collapse
Affiliation(s)
- Noëlle K J Bittner
- Department of Integrative Biology and Museum of Vertebrate Zoology, 3101 Valley Life Sciences Building, University of California Berkeley, California 94720
| | - Katya L Mack
- Present address: Department of Biology, Stanford University, CA 94305
| | | |
Collapse
|
16
|
A high-quality genome of the dobsonfly Neoneuromus ignobilis reveals molecular convergences in aquatic insects. Genomics 2022; 114:110437. [PMID: 35902070 DOI: 10.1016/j.ygeno.2022.110437] [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: 05/08/2022] [Revised: 07/03/2022] [Accepted: 07/21/2022] [Indexed: 11/22/2022]
Abstract
Neoneuromus ignobilis is an archaic holometabolous aquatic predatory insect. However, a lack of genomic resources hinders the use of whole genome sequencing to explore their genetic basis and molecular mechanisms for adaptive evolution. Here, we provided a high-contiguity, chromosome-level genome assembly of N. ignobilis using high coverage Nanopore and PacBio reads with the Hi-C technique. The final assembly is 480.67 MB in size, containing 12 telomere-ended pseudochromosomes with only 17 gaps. We compared 42 hexapod species genomes including six independent lineages comprising 11 aquatic insects, and found convergent expansions of long wavelength-sensitive and blue-sensitive opsins, thermal stress response TRP channels, and sulfotransferases in aquatic insects, which may be related to their aquatic adaptation. We also detected strong nonrandom signals of convergent amino acid substitutions in aquatic insects. Collectively, our comparative genomic analysis revealed the evidence of molecular convergences in aquatic insects during both gene family evolution and convergent amino acid substitutions.
Collapse
|
17
|
Font-Porterias N, McNelis MG, Comas D, Hlusko LJ. Evidence of selection in the ectodysplasin pathway among endangered aquatic mammals. Integr Org Biol 2022; 4:obac018. [PMID: 35874492 PMCID: PMC9299678 DOI: 10.1093/iob/obac018] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 05/06/2022] [Accepted: 05/21/2022] [Indexed: 11/13/2022] Open
Abstract
Synopsis The ectodysplasin pathway has been a target of evolution repeatedly. Genetic variation in the key genes of this pathway (EDA, EDAR, and EDARADD) results in a rich source of pleiotropic effects across ectodermally-derived structures, including teeth, hair, sweat glands, and mammary glands. In addition, a non-canonical Wnt pathway has a very similar functional role, making variation in the WNT10A gene also of evolutionary significance. The adaptation of mammals to aquatic environments has occurred independently in at least 4 orders, whose species occupy a wide geographic range (from equatorial to polar regions) and exhibit great phenotypic variation in ectodermally-derived structures, including the presence or absence of fur and extreme lactational strategies. The role of the ectodysplasin pathway in the adaptation to aquatic environments has been never explored in mammalian species. In the present study, we analyze the genetic variation in orthologous coding sequences from EDA, EDAR, EDARADD, and WNT10A genes together with ectodermally-derived phenotypic variation from 34 aquatic and non-aquatic mammalian species to assess signals of positive selection, gene-trait coevolution, and genetic convergence. Our study reveals strong evidence of positive selection in a proportion of coding sites in EDA and EDAR genes in 3 endangered aquatic mammals (the Hawaiian monk seal, the Yangtze finless porpoise, and the sea otter). We hypothesize functional implications potentially related to the adaptation to the low-latitude aquatic environment in the Hawaiian monk seal and the freshwater in the Yangtze finless porpoise. The signal in the sea otter is likely the result of an increased genetic drift after an intense bottleneck and reduction of genetic diversity. Besides positive selection, we have not detected robust signals of gene-trait coevolution or convergent amino acid shifts in the ectodysplasin pathway associated with shared phenotypic traits among aquatic mammals. This study provides new evidence of the evolutionary role of the ectodysplasin pathway and encourages further investigation, including functional studies, to fully resolve its relationship with mammalian aquatic adaptation. Spanish La vía de la ectodisplasina ha sido objeto de la evolución repetidamente. La variación genética en los principales genes de esta vía (EDA, EDAR y EDARADD) da como resultado una gran diversidad de efectos pleiotrópicos en las estructuras derivadas del ectodermo, incluidos los dientes, el cabello, las glándulas sudoríparas y las glándulas mamarias. Además, una vía wnt no canónica tiene un papel funcional muy similar, por lo que la variación en el gen WNT10A también tiene importancia evolutiva. La adaptación de los mamíferos a los entornes acuáticos se ha producido de forma independiente en al menos cuatro órdenes, cuyas especies ocupan un amplio rango geográfico (desde regiones ecuatoriales a polares) y presentan una gran variación fenotípica en las estructuras derivadas del ectodermo, incluyendo la presencia o ausencia de pelaje y estrategias de lactancia muy diferentes. El papel de la vía de la ectodisplasina en la adaptación a entornos acuáticos no se ha explorado nunca en especies de mamíferos. En este estudio, analizamos la variación genética en las secuencias codificantes ortólogas de los genes EDA, EDAR, EDARADD y WNT10A junto con la variación fenotípica derivada del ectodermo de 34 especies de mamíferos acuáticos y no acuáticos para evaluar señales de selección positiva, coevolución gen-rasgo y convergencia genética. Nuestro estudio revela señales de selección positiva en regiones de las secuencias codificantes de los genes EDA y EDAR en tres mamíferos acuáticos en peligro de extinción (la foca monje de Hawái, la marsopa lisa y la nutria marina). Estas señales podrían tener implicaciones funcionales potencialmente relacionadas con la adaptación al entorno acuático de baja latitud en la foca monje de Hawái y el agua dulce en la marsopa lisa. La señal en la nutria marina es probablemente el resultado de una mayor deriva genética tras un intenso un cuello de botella y una reducción de la diversidad genética. A parte de selección positiva, no hemos detectado señales sólidas de coevolución gen-rasgo o cambios convergentes de aminoácidos en la vía de la ectodisplasina asociados a rasgos fenotípicos compartidos entre mamíferos acuáticos. Este estudio proporciona nuevas evidencias del papel evolutivo de la vía de la ectodisplasina y quiere promover futuras investigaciones con estudios funcionales para acabar de resolver la relación de esta vía con la adaptación acuática de los mamíferos.
Collapse
Affiliation(s)
- Neus Font-Porterias
- Departament de Ciències Experimentals i de la Salut, Universitat Pompeu Fabra, Institut de Biologia Evolutiva (UPF-CSIC) , Barcelona , Spain
| | - Madeline G McNelis
- Department of Integrative Biology, University of California Berkeley , California , USA
| | - David Comas
- Departament de Ciències Experimentals i de la Salut, Universitat Pompeu Fabra, Institut de Biologia Evolutiva (UPF-CSIC) , Barcelona , Spain
| | - Leslea J Hlusko
- Department of Integrative Biology, University of California Berkeley , California , USA
- National Research Center on Human Evolution (CENIEH) , Burgos , Spain
| |
Collapse
|
18
|
Prestin-Mediated Frequency Selectivity Does not Cover Ultrahigh Frequencies in Mice. Neurosci Bull 2022; 38:769-784. [DOI: 10.1007/s12264-022-00839-4] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/30/2021] [Accepted: 12/24/2021] [Indexed: 02/08/2023] Open
|
19
|
Schull JK, Turakhia Y, Hemker JA, Dally WJ, Bejerano G. OUP accepted manuscript. Genome Biol Evol 2022; 14:6529394. [PMID: 35171243 PMCID: PMC8920512 DOI: 10.1093/gbe/evac013] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 01/10/2022] [Indexed: 11/14/2022] Open
Abstract
We present Champagne, a whole-genome method for generating character matrices for phylogenomic analysis using large genomic indel events. By rigorously picking orthologous genes and locating large insertion and deletion events, Champagne delivers a character matrix that considerably reduces homoplasy compared with morphological and nucleotide-based matrices, on both established phylogenies and difficult-to-resolve nodes in the mammalian tree. Champagne provides ample evidence in the form of genomic structural variation to support incomplete lineage sorting and possible introgression in Paenungulata and human–chimp–gorilla which were previously inferred primarily through matrices composed of aligned single-nucleotide characters. Champagne also offers further evidence for Myomorpha as sister to Sciuridae and Hystricomorpha in the rodent tree. Champagne harbors distinct theoretical advantages as an automated method that produces nearly homoplasy-free character matrices on the whole-genome scale.
Collapse
Affiliation(s)
- James K Schull
- Department of Computer Science, Stanford University, USA
| | - Yatish Turakhia
- Department of Electrical and Computer Engineering, University of California San Diego, USA
| | - James A Hemker
- Department of Computer Science, Stanford University, USA
| | - William J Dally
- Department of Computer Science, Stanford University, USA
- NVIDIA, Santa Clara, California, USA
- Department of Electrical Engineering, Stanford University, USA
| | - Gill Bejerano
- Department of Computer Science, Stanford University, USA
- Department of Developmental Biology, Stanford University, USA
- Department of Biomedical Data Science, Stanford University, USA
- Department of Pediatrics, Stanford University, USA
- Corresponding author: E-mail:
| |
Collapse
|
20
|
Upadhyay M, Kunz E, Sandoval-Castellanos E, Hauser A, Krebs S, Graf A, Blum H, Dotsev A, Okhlopkov I, Shakhin A, Bagirov V, Brem G, Fries R, Zinovieva N, Medugorac I. Whole genome sequencing reveals a complex introgression history and the basis of adaptation to subarctic climate in wild sheep. Mol Ecol 2021; 30:6701-6717. [PMID: 34534381 DOI: 10.1111/mec.16184] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2021] [Revised: 09/01/2021] [Accepted: 09/07/2021] [Indexed: 11/29/2022]
Abstract
To predict species responses to anthropogenic disturbances and climate change, it is reasonable to use species with high sensitivity to such factors. Snow sheep (Ovis nivicola) could represent a good candidate for this; as the only large herbivore species adapted to the cold and alpine habitats of northeastern Siberia, it plays a crucial role in its ecosystem. Despite having an extensive geographical distribution among all ovine species, it is one of the least studied. In this study, we sequenced and analysed six genomes of snow sheep in combination with all other wild sheep species to infer key aspects of their evolutionary history and unveil the genetic basis of their adaptation to subarctic environments. Despite their large census population size, snow sheep genomes showed remarkably low heterozygosity, which could reflect the effect of isolation and historical bottlenecks that we inferred using the pairwise sequential Markovian coalescent and runs of homozygosity. F4 -statistics indicated instances of introgression involving snow sheep with argali (Ovis ammon) and Dall (Ovis dalli) sheep, suggesting that these species might have been more widespread during the Pleistocene. Furthermore, the introgressed segments, which were identified using mainly minimum relative node depth, covered genes associated with immunity, adipogenesis and morphology-related traits, representing potential targets of adaptive introgression. Genes related to mitochondrial functions and thermogenesis associated with adipose tissue were identified to be under selection. Overall, our data suggest introgression as a mechanism facilitating adaptation in wild sheep species and provide insights into the genetic mechanisms underlying cold adaptation in snow sheep.
Collapse
Affiliation(s)
- Maulik Upadhyay
- Population Genomics Group, Department of Veterinary Sciences, LMU Munich, Munich, Germany
| | - Elisabeth Kunz
- Population Genomics Group, Department of Veterinary Sciences, LMU Munich, Munich, Germany
| | | | - Andreas Hauser
- Laboratory for Functional Genome Analysis, Gene Center, LMU Munich, Munich, Germany
| | - Stefan Krebs
- Laboratory for Functional Genome Analysis, Gene Center, LMU Munich, Munich, Germany
| | - Alexander Graf
- Laboratory for Functional Genome Analysis, Gene Center, LMU Munich, Munich, Germany
| | - Helmut Blum
- Laboratory for Functional Genome Analysis, Gene Center, LMU Munich, Munich, Germany
| | - Arsen Dotsev
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | | | - Alexey Shakhin
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Vugar Bagirov
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Gottfried Brem
- Institute of Animal Breeding and Genetics, VMU, Vienna, Austria
| | - Ruedi Fries
- Lehrstuhl für Tierzucht, Technische Universität München, Freising, Germany
| | - Natalia Zinovieva
- L.K. Ernst Federal Science Center for Animal Husbandry, Podolsk, Russia
| | - Ivica Medugorac
- Population Genomics Group, Department of Veterinary Sciences, LMU Munich, Munich, Germany
| |
Collapse
|
21
|
Comparative genomics provides insights into the aquatic adaptations of mammals. Proc Natl Acad Sci U S A 2021; 118:2106080118. [PMID: 34503999 PMCID: PMC8449357 DOI: 10.1073/pnas.2106080118] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/26/2021] [Indexed: 12/30/2022] Open
Abstract
Divergent lineages can respond to common environmental factors through convergent processes involving shared genomic components or pathways, but the molecular mechanisms are poorly understood. Here, we provide genomic resources and insights into the evolution of mammalian lineages adapting to aquatic life. Our data suggest convergent evolution, for example, in association with thermoregulation through genes associated with a surface heat barrier (NFIA) and internal heat exchange (SEMA3E). Combined with the support of previous reports showing that the UCP1 locus has been lost in many marine mammals independently, our results suggest that the thermostatic strategy of marine mammals shifted from enhancing heat production to limiting heat loss. The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.
Collapse
|
22
|
Treaster S, Daane JM, Harris MP. Refining Convergent Rate Analysis with Topology in Mammalian Longevity and Marine Transitions. Mol Biol Evol 2021; 38:5190-5203. [PMID: 34324001 PMCID: PMC8557430 DOI: 10.1093/molbev/msab226] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
The quest to map the genetic foundations of phenotypes has been empowered by the modern diversity, quality, and availability of genomic resources. Despite these expanding resources, the abundance of variation within lineages makes it challenging to associate genetic change to specific phenotypes, without an a priori means of isolating the changes from background genomic variation. Evolution provides this means through convergence-i.e., the shared variation that may result from replicate evolutionary experiments across independent trait occurrences. To leverage these opportunities, we developed TRACCER: Topologically Ranked Analysis of Convergence via Comparative Evolutionary Rates. Compared to current methods, this software empowers rate convergence analysis by factoring in topological relationships, because genetic variation between phylogenetically proximate trait changes is more likely to be facilitating the trait. Comparisons are performed not with singular branches, but with the complete paths to the most recent common ancestor for each pair of lineages. This ensures that comparisons represent a single context diverging over the same timeframe while obviating the problematic requirement of assigning ancestral states. We applied TRACCER to two case studies: mammalian transitions to marine environments, an unambiguous collection of traits which have independently evolved three times; and the evolution of mammalian longevity, a less delineated trait but with more instances to compare. By factoring in topology, TRACCER identifies highly significant, convergent genetic signals, with important incongruities and statistical resolution when compared to existing approaches. These improvements in sensitivity and specificity of convergence analysis generates refined targets for downstream validation and identification of genotype-phenotype relationships.
Collapse
Affiliation(s)
- Stephen Treaster
- Department of Orthopaedic Research, Boston Children's Hospital, Boston, MA, 02124, USA.,Department of Genetics, Harvard Medical School, Boston, MA, 02124, USA
| | - Jacob M Daane
- Department of Orthopaedic Research, Boston Children's Hospital, Boston, MA, 02124, USA.,Department of Genetics, Harvard Medical School, Boston, MA, 02124, USA.,Department of Marine and Environmental Sciences, Northeastern University Marine Science Center, Nahant, MA, 01908, USA
| | - Matthew P Harris
- Department of Orthopaedic Research, Boston Children's Hospital, Boston, MA, 02124, USA.,Department of Genetics, Harvard Medical School, Boston, MA, 02124, USA
| |
Collapse
|
23
|
Turakhia Y, Chen HI, Marcovitz A, Bejerano G. A fully-automated method discovers loss of mouse-lethal and human-monogenic disease genes in 58 mammals. Nucleic Acids Res 2020; 48:e91. [PMID: 32614390 PMCID: PMC7498332 DOI: 10.1093/nar/gkaa550] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2020] [Revised: 05/23/2020] [Accepted: 06/23/2020] [Indexed: 01/20/2023] Open
Abstract
Gene losses provide an insightful route for studying the morphological and physiological adaptations of species, but their discovery is challenging. Existing genome annotation tools focus on annotating intact genes and do not attempt to distinguish nonfunctional genes from genes missing annotation due to sequencing and assembly artifacts. Previous attempts to annotate gene losses have required significant manual curation, which hampers their scalability for the ever-increasing deluge of newly sequenced genomes. Using extreme sequence erosion (amino acid deletions and substitutions) and sister species support as an unambiguous signature of loss, we developed an automated approach for detecting high-confidence gene loss events across a species tree. Our approach relies solely on gene annotation in a single reference genome, raw assemblies for the remaining species to analyze, and the associated phylogenetic tree for all organisms involved. Using human as reference, we discovered over 400 unique human ortholog erosion events across 58 mammals. This includes dozens of clade-specific losses of genes that result in early mouse lethality or are associated with severe human congenital diseases. Our discoveries yield intriguing potential for translational medical genetics and evolutionary biology, and our approach is readily applicable to large-scale genome sequencing efforts across the tree of life.
Collapse
Affiliation(s)
- Yatish Turakhia
- Department of Electrical Engineering, Stanford University, Stanford, CA 94305, USA
| | - Heidi I Chen
- Department of Developmental Biology, Stanford University, Stanford, CA 94305, USA
| | - Amir Marcovitz
- Department of Developmental Biology, Stanford University, Stanford, CA 94305, USA
| | - Gill Bejerano
- Department of Developmental Biology, Stanford University, Stanford, CA 94305, USA
- Department of Computer Science, Stanford University, Stanford, CA 94305, USA
- Department of Biomedical Data Science, Stanford University, Stanford, CA 94305, USA
- Department of Pediatrics, Stanford University, Stanford, CA 94305, USA
| |
Collapse
|
24
|
Pyott SJ, van Tuinen M, Screven LA, Schrode KM, Bai JP, Barone CM, Price SD, Lysakowski A, Sanderford M, Kumar S, Santos-Sacchi J, Lauer AM, Park TJ. Functional, Morphological, and Evolutionary Characterization of Hearing in Subterranean, Eusocial African Mole-Rats. Curr Biol 2020; 30:4329-4341.e4. [PMID: 32888484 DOI: 10.1016/j.cub.2020.08.035] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/03/2020] [Revised: 06/09/2020] [Accepted: 08/07/2020] [Indexed: 12/26/2022]
Abstract
Naked mole-rats are highly vocal, eusocial, subterranean rodents with, counterintuitively, poor hearing. The causes underlying their altered hearing are unknown. Moreover, whether altered hearing is degenerate or adaptive to their unique lifestyles is controversial. We used various methods to identify the factors contributing to altered hearing in naked and the related Damaraland mole-rats and to examine whether these alterations result from relaxed or adaptive selection. Remarkably, we found that cochlear amplification was absent from both species despite normal prestin function in outer hair cells isolated from naked mole-rats. Instead, loss of cochlear amplification appears to result from abnormal hair bundle morphologies observed in both species. By exploiting a well-curated deafness phenotype-genotype database, we identified amino acid substitutions consistent with abnormal hair bundle morphology and reduced hearing sensitivity. Amino acid substitutions were found in unique groups of six hair bundle link proteins. Molecular evolutionary analyses revealed shifts in selection pressure at both the gene and the codon level for five of these six hair bundle link proteins. Substitutions in three of these proteins are associated exclusively with altered hearing. Altogether, our findings identify the likely mechanism of altered hearing in African mole-rats, making them the only identified mammals naturally lacking cochlear amplification. Moreover, our findings suggest that altered hearing in African mole-rats is adaptive, perhaps tailoring hearing to eusocial and subterranean lifestyles. Finally, our work reveals multiple, unique evolutionary trajectories in African mole-rat hearing and establishes species members as naturally occurring disease models to investigate human hearing loss.
Collapse
Affiliation(s)
- Sonja J Pyott
- University Medical Center Groningen and University of Groningen, Department of Otorhinolaryngology and Head/Neck Surgery, 9713GZ Groningen, the Netherlands.
| | - Marcel van Tuinen
- University Medical Center Groningen and University of Groningen, Department of Otorhinolaryngology and Head/Neck Surgery, 9713GZ Groningen, the Netherlands
| | - Laurel A Screven
- Johns Hopkins School of Medicine, Department of Otolaryngology, Baltimore, MD 21205, USA
| | - Katrina M Schrode
- Johns Hopkins School of Medicine, Department of Otolaryngology, Baltimore, MD 21205, USA
| | - Jun-Ping Bai
- Yale University School of Medicine, Department of Neurology, 333 Cedar Street, New Haven, CT 06510, USA
| | - Catherine M Barone
- University of Illinois at Chicago, Department of Biological Sciences, Chicago, IL 60612, USA
| | - Steven D Price
- University of Illinois at Chicago, Department of Anatomy and Cell Biology, Chicago, IL 60612, USA
| | - Anna Lysakowski
- University of Illinois at Chicago, Department of Anatomy and Cell Biology, Chicago, IL 60612, USA
| | - Maxwell Sanderford
- Temple University, Institute for Genomics and Evolutionary Medicine and Department of Biology, Philadelphia, PA 19122, USA
| | - Sudhir Kumar
- Temple University, Institute for Genomics and Evolutionary Medicine and Department of Biology, Philadelphia, PA 19122, USA; King Abdulaziz University, Center for Excellence in Genome Medicine and Research, Jeddah, Saudi Arabia
| | - Joseph Santos-Sacchi
- Yale University School of Medicine, Department of Surgery (Otolaryngology) and Department of Neuroscience and Cellular and Molecular Physiology, 333 Cedar Street, New Haven, CT 06510, USA
| | - Amanda M Lauer
- Johns Hopkins School of Medicine, Department of Otolaryngology, Baltimore, MD 21205, USA
| | - Thomas J Park
- University of Illinois at Chicago, Department of Biological Sciences, Chicago, IL 60612, USA
| |
Collapse
|
25
|
Burskaia V, Naumenko S, Schelkunov M, Bedulina D, Neretina T, Kondrashov A, Yampolsky L, Bazykin GA. Excessive Parallelism in Protein Evolution of Lake Baikal Amphipod Species Flock. Genome Biol Evol 2020; 12:1493-1503. [PMID: 32653919 DOI: 10.1093/gbe/evaa138] [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] [Accepted: 07/03/2020] [Indexed: 11/12/2022] Open
Abstract
Repeated emergence of similar adaptations is often explained by parallel evolution of underlying genes. However, evidence of parallel evolution at amino acid level is limited. When the analyzed species are highly divergent, this can be due to epistatic interactions underlying the dynamic nature of the amino acid preferences: The same amino acid substitution may have different phenotypic effects on different genetic backgrounds. Distantly related species also often inhabit radically different environments, which makes the emergence of parallel adaptations less likely. Here, we hypothesize that parallel molecular adaptations are more prevalent between closely related species. We analyze the rate of parallel evolution in genome-size sets of orthologous genes in three groups of species with widely ranging levels of divergence: 46 species of the relatively recent lake Baikal amphipod radiation, a species flock of very closely related cichlids, and a set of significantly more divergent vertebrates. Strikingly, in genes of amphipods, the rate of parallel substitutions at nonsynonymous sites exceeded that at synonymous sites, suggesting rampant selection driving parallel adaptation. At sites of parallel substitutions, the intraspecies polymorphism is low, suggesting that parallelism has been driven by positive selection and is therefore adaptive. By contrast, in cichlids, the rate of nonsynonymous parallel evolution was similar to that at synonymous sites, whereas in vertebrates, this rate was lower than that at synonymous sites, indicating that in these groups of species, parallel substitutions are mainly fixed by drift.
Collapse
Affiliation(s)
- Valentina Burskaia
- Center of Life Sciences, Skolkovo Institute of Science and Technology, Moscow, Moscow Oblast, Russia
| | - Sergey Naumenko
- Institute for Information Transmission Problems of the Russian Academy of Sciences (Kharkevitch Institute), Moscow, Russia
- Department of Biostatistics, Harvard T.H. Chan School of Public Health, Boston, Massachusetts
| | - Mikhail Schelkunov
- Center of Life Sciences, Skolkovo Institute of Science and Technology, Moscow, Moscow Oblast, Russia
- Institute for Information Transmission Problems of the Russian Academy of Sciences (Kharkevitch Institute), Moscow, Russia
| | - Daria Bedulina
- Institute of Biology, Irkutsk State University, Russia
- Baikal Research Centre, Irkutsk, Russia
| | - Tatyana Neretina
- Institute for Information Transmission Problems of the Russian Academy of Sciences (Kharkevitch Institute), Moscow, Russia
- N.A. Pertsov White Sea Biological Station, Lomonosov Moscow State University, Primorskiy, Russia
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Russia
| | - Alexey Kondrashov
- Belozersky Institute of Physico-Chemical Biology, Lomonosov Moscow State University, Russia
- Department of Ecology and Evolutionary Biology, University of Michigan
| | - Lev Yampolsky
- Department of Biological Sciences, East Tennessee State University
| | - Georgii A Bazykin
- Center of Life Sciences, Skolkovo Institute of Science and Technology, Moscow, Moscow Oblast, Russia
- Institute for Information Transmission Problems of the Russian Academy of Sciences (Kharkevitch Institute), Moscow, Russia
| |
Collapse
|
26
|
Smith SD, Pennell MW, Dunn CW, Edwards SV. Phylogenetics is the New Genetics (for Most of Biodiversity). Trends Ecol Evol 2020; 35:415-425. [DOI: 10.1016/j.tree.2020.01.005] [Citation(s) in RCA: 24] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Revised: 01/15/2020] [Accepted: 01/20/2020] [Indexed: 12/15/2022]
|