1
|
Insana G, Martin MJ, Pearson WR. Improved selection of canonical proteins for reference proteomes. NAR Genom Bioinform 2024; 6:lqae066. [PMID: 38863529 PMCID: PMC11165316 DOI: 10.1093/nargab/lqae066] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/08/2024] [Revised: 05/04/2024] [Accepted: 05/23/2024] [Indexed: 06/13/2024] Open
Abstract
The 'canonical' protein sets distributed by UniProt are widely used for similarity searching, and functional and structural annotation. For many investigators, canonical sequences are the only version of a protein examined. However, higher eukaryotes often encode multiple isoforms of a protein from a single gene. For unreviewed (UniProtKB/TrEMBL) protein sequences, the longest sequence in a Gene-Centric group is chosen as canonical. This choice can create inconsistencies, selecting >95% identical orthologs with dramatically different lengths, which is biologically unlikely. We describe the ortho2tree pipeline, which examines Reference Proteome canonical and isoform sequences from sets of orthologous proteins, builds multiple alignments, constructs gap-distance trees, and identifies low-cost clades of isoforms with similar lengths. After examining 140 000 proteins from eight mammals in UniProtKB release 2022_05, ortho2tree proposed 7804 canonical changes for release 2023_01, while confirming 53 434 canonicals. Gap distributions for isoforms selected by ortho2tree are similar to those in bacterial and yeast alignments, organisms unaffected by isoform selection, suggesting ortho2tree canonicals more accurately reflect genuine biological variation. 82% of ortho2tree proposed-changes agreed with MANE; for confirmed canonicals, 92% agreed with MANE. Ortho2tree can improve canonical assignment among orthologous sequences that are >60% identical, a group that includes vertebrates and higher plants.
Collapse
Affiliation(s)
- Giuseppe Insana
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton CB10 1SD, UK
| | - Maria J Martin
- European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Hinxton CB10 1SD, UK
| | - William R Pearson
- Department of Biochemistry and Molecular Genetics, University of Virginia School of Medicine, Charlottesville, VA 22908, USA
| |
Collapse
|
2
|
Villanueva-Cañas JL, Fernandez-Fuentes N, Saul D, Kosinsky RL, Teyssier C, Rogalska ME, Pérez FP, Oliva B, Notredame C, Beato M, Sharma P. Evolutionary analysis reveals the role of a non-catalytic domain of peptidyl arginine deiminase 2 in transcriptional regulation. iScience 2024; 27:109584. [PMID: 38623337 PMCID: PMC11016909 DOI: 10.1016/j.isci.2024.109584] [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: 12/21/2023] [Revised: 02/13/2024] [Accepted: 03/25/2024] [Indexed: 04/17/2024] Open
Abstract
Peptidyl arginine deiminases (PADIs) catalyze protein citrullination, a post-translational conversion of arginine to citrulline. The most widely expressed member of this family, PADI2, regulates cellular processes that impact several diseases. We hypothesized that we could gain new insights into PADI2 function through a systematic evolutionary and structural analysis. Here, we identify 20 positively selected PADI2 residues, 16 of which are structurally exposed and maintain PADI2 interactions with cognate proteins. Many of these selected residues reside in non-catalytic regions of PADI2. We validate the importance of a prominent loop in the middle domain that encompasses PADI2 L162, a residue under positive selection. This site is essential for interaction with the transcription elongation factor (P-TEFb) and mediates the active transcription of the oncogenes c-MYC, and CCNB1, as well as impacting cellular proliferation. These insights could be key to understanding and addressing the role of the PADI2 c-MYC axis in cancer progression.
Collapse
Affiliation(s)
- José Luis Villanueva-Cañas
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Narcis Fernandez-Fuentes
- Institute of Biological, Environmental and Rural Sciences, Aberystwyth University, Aberystwyth, Ceredigion, United Kingdom
| | - Dominik Saul
- Division of Endocrinology, Mayo Clinic, Rochester, MN 55905, USA; Robert and Arlene Kogod Center on Aging, Mayo Clinic, Rochester, MN 55905, USA
- Department of Trauma and Reconstructive Surgery, BG Clinic, University of Tübingen, Tübingen, Germany
| | | | - Catherine Teyssier
- Institut de Recherche en Cancérologie de Montpellier (IRCM), INSERM U1194, Université de Montpellier, Institut Du Cancer de Montpellier (ICM), F-34298 Montpellier, France
| | - Malgorzata Ewa Rogalska
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
| | - Ferran Pegenaute Pérez
- Live-Cell Structural Biology Laboratory, Department of Medicine and Life Sciences, E-08005 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Baldomero Oliva
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
- Structural Bioinformatics Laboratory (GRIB-IMIM), Department of Medicine and Life Sciences, E-08003 Barcelona, Spain
| | - Cedric Notredame
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Miguel Beato
- Centre for Genomic Regulation (CRG), The Barcelona Institute of Science and Technology (BIST), Dr. Aiguader 88, 08003 Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | - Priyanka Sharma
- Institut de Pharmacologie et de Biologie Structurale, IPBS, Université de Toulouse, CNRS, UPS, Toulouse, France
| |
Collapse
|
3
|
Feldmeyer B, Bornberg-Bauer E, Dohmen E, Fouks B, Heckenhauer J, Huylmans AK, Jones ARC, Stolle E, Harrison MC. Comparative Evolutionary Genomics in Insects. Methods Mol Biol 2024; 2802:473-514. [PMID: 38819569 DOI: 10.1007/978-1-0716-3838-5_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 06/01/2024]
Abstract
Genome sequencing quality, in terms of both read length and accuracy, is constantly improving. By combining long-read sequencing technologies with various scaffolding techniques, chromosome-level genome assemblies are now achievable at an affordable price for non-model organisms. Insects represent an exciting taxon for studying the genomic underpinnings of evolutionary innovations, due to ancient origins, immense species-richness, and broad phenotypic diversity. Here we summarize some of the most important methods for carrying out a comparative genomics study on insects. We describe available tools and offer concrete tips on all stages of such an endeavor from DNA extraction through genome sequencing, annotation, and several evolutionary analyses. Along the way we describe important insect-specific aspects, such as DNA extraction difficulties or gene families that are particularly difficult to annotate, and offer solutions. We describe results from several examples of comparative genomics analyses on insects to illustrate the fascinating questions that can now be addressed in this new age of genomics research.
Collapse
Affiliation(s)
- Barbara Feldmeyer
- Senckenberg Biodiversity and Climate Research Centre (SBiK-F), Molecular Ecology, Frankfurt, Germany
| | - Erich Bornberg-Bauer
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany
- Department of Protein Evolution, Max Planck Institute for Developmental Biology, Tübingen, Germany
| | - Elias Dohmen
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany
| | - Bertrand Fouks
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany
| | - Jacqueline Heckenhauer
- LOEWE Centre for Translational Biodiversity Genomics (LOEWE-TBG), Frankfurt, Germany
- Department of Terrestrial Zoology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Germany
| | - Ann Kathrin Huylmans
- Institute of Organismic and Molecular Evolution, Johannes Gutenberg University, Mainz, Germany
| | - Alun R C Jones
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany
| | - Eckart Stolle
- Museum Koenig, Leibniz Institute for the Analysis of Biodiversity Change (LIB), Bonn, Germany
| | - Mark C Harrison
- Institute for Evolution and Biodiversity, University of Münster, Münster, Germany.
| |
Collapse
|
4
|
Cui L, Cheng H, Yang Z, Xia C, Zhang L, Kong X. Comparative Analysis Reveals Different Evolutionary Fates and Biological Functions in Wheat Duplicated Genes ( Triticum aestivum L.). PLANTS (BASEL, SWITZERLAND) 2023; 12:3021. [PMID: 37687268 PMCID: PMC10489728 DOI: 10.3390/plants12173021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/31/2023] [Revised: 08/20/2023] [Accepted: 08/21/2023] [Indexed: 09/10/2023]
Abstract
Wheat (Triticum aestivum L.) is a staple food crop that provides 20% of total human calorie consumption. Gene duplication has been considered to play an important role in evolution by providing new genetic resources. However, the evolutionary fates and biological functions of the duplicated genes in wheat remain to be elucidated. In this study, the resulting data showed that the duplicated genes evolved faster with shorter gene lengths, higher codon usage bias, lower expression levels, and higher tissue specificity when compared to non-duplicated genes. Our analysis further revealed functions of duplicated genes in various biological processes with significant enrichment to environmental stresses. In addition, duplicated genes derived from dispersed, proximal, tandem, transposed, and whole-genome duplication differed in abundance, evolutionary rate, gene compactness, expression pattern, and genetic diversity. Tandem and proximal duplicates experienced stronger selective pressure and showed a more compact gene structure with diverse expression profiles than other duplication modes. Moreover, genes derived from different duplication modes showed an asymmetrical evolutionary pattern for wheat A, B, and D subgenomes. Several candidate duplication hotspots associated with wheat domestication or polyploidization were characterized as potential targets for wheat molecular breeding. Our comprehensive analysis revealed the evolutionary trajectory of duplicated genes and laid the foundation for future functional studies on wheat.
Collapse
Affiliation(s)
- Licao Cui
- Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (L.C.); (H.C.); (Z.Y.); (C.X.); (L.Z.)
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang 330045, China
| | - Hao Cheng
- Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (L.C.); (H.C.); (Z.Y.); (C.X.); (L.Z.)
- State Key Laboratory of Crop Stress Biology for Arid Areas, College of Life Sciences, Northwest A&F University, Yangling 712100, China
| | - Zhe Yang
- Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (L.C.); (H.C.); (Z.Y.); (C.X.); (L.Z.)
| | - Chuan Xia
- Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (L.C.); (H.C.); (Z.Y.); (C.X.); (L.Z.)
| | - Lichao Zhang
- Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (L.C.); (H.C.); (Z.Y.); (C.X.); (L.Z.)
| | - Xiuying Kong
- Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing 100081, China; (L.C.); (H.C.); (Z.Y.); (C.X.); (L.Z.)
| |
Collapse
|
5
|
Lian ZJ, Phan MD, Hancock SJ, Nhu NTK, Paterson DL, Schembri MA. Genetic basis of I-complex plasmid stability and conjugation. PLoS Genet 2023; 19:e1010773. [PMID: 37347771 DOI: 10.1371/journal.pgen.1010773] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2023] [Accepted: 05/05/2023] [Indexed: 06/24/2023] Open
Abstract
Plasmids are major drivers of increasing antibiotic resistance, necessitating an urgent need to understand their biology. Here we describe a detailed dissection of the molecular components controlling the genetics of I-complex plasmids, a group of antibiotic resistance plasmids found frequently in pathogenic Escherichia coli and other Enterobacteriaceae that cause significant human disease. We show these plasmids cluster into four distinct subgroups, with the prototype IncI1 plasmid R64 subgroup displaying low nucleotide sequence conservation to other I-complex plasmids. Using pMS7163B, an I-complex plasmid distantly related to R64, we performed a high-resolution transposon-based genetic screen and defined genes involved in replication, stability, and conjugative transfer. We identified the replicon and a partitioning system as essential for replication/stability. Genes required for conjugation included the type IV secretion system, relaxosome, and several uncharacterised genes located in the pMS7163B leading transfer region that exhibited an upstream strand-specific transposon insertion bias. The overexpression of these genes severely impacted host cell growth or reduced fitness during mixed competitive growth, demonstrating that their expression must be controlled to avoid deleterious impacts. These genes were present in >80% of all I-complex plasmids and broadly conserved across multiple plasmid incompatibility groups, implicating an important role in plasmid dissemination.
Collapse
Affiliation(s)
- Zheng Jie Lian
- Institute for Molecular Bioscience (IMB), The University of Queensland, Brisbane, Queensland, Australia
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia
| | - Minh-Duy Phan
- Institute for Molecular Bioscience (IMB), The University of Queensland, Brisbane, Queensland, Australia
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia
| | - Steven J Hancock
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia
| | - Nguyen Thi Khanh Nhu
- Institute for Molecular Bioscience (IMB), The University of Queensland, Brisbane, Queensland, Australia
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia
| | - David L Paterson
- The University of Queensland Centre for Clinical Research, Brisbane, Australia
| | - Mark A Schembri
- Institute for Molecular Bioscience (IMB), The University of Queensland, Brisbane, Queensland, Australia
- School of Chemistry and Molecular Biosciences, The University of Queensland, Brisbane, Queensland, Australia
- Australian Infectious Diseases Research Centre, The University of Queensland, Brisbane, Queensland, Australia
| |
Collapse
|
6
|
Tao W, Li R, Li T, Li Z, Li Y, Cui L. The evolutionary patterns, expression profiles, and genetic diversity of expanded genes in barley. FRONTIERS IN PLANT SCIENCE 2023; 14:1168124. [PMID: 37180392 PMCID: PMC10171312 DOI: 10.3389/fpls.2023.1168124] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/17/2023] [Accepted: 03/28/2023] [Indexed: 05/16/2023]
Abstract
Gene duplication resulting from whole-genome duplication (WGD), small-scale duplication (SSD), or unequal hybridization plays an important role in the expansion of gene families. Gene family expansion can also mediate species formation and adaptive evolution. Barley (Hordeum vulgare) is the world's fourth largest cereal crop, and it contains valuable genetic resources due to its ability to tolerate various types of environmental stress. In this study, 27,438 orthogroups in the genomes of seven Poaceae were identified, and 214 of them were significantly expanded in barley. The evolutionary rates, gene properties, expression profiles, and nucleotide diversity between expanded and non-expanded genes were compared. Expanded genes evolved more rapidly and experienced lower negative selection. Expanded genes, including their exons and introns, were shorter, they had fewer exons, their GC content was lower, and their first exons were longer compared with non-expanded genes. Codon usage bias was also lower for expanded genes than for non-expanded genes; the expression levels of expanded genes were lower than those of non-expanded genes, and the expression of expanded genes showed higher tissue specificity than that of non-expanded genes. Several stress-response-related genes/gene families were identified, and these genes could be used to breed barley plants with greater resistance to environmental stress. Overall, our analysis revealed evolutionary, structural, and functional differences between expanded and non-expanded genes in barley. Additional studies are needed to clarify the functions of the candidate genes identified in our study and evaluate their utility for breeding barley plants with greater stress resistance.
Collapse
Affiliation(s)
- Wenjing Tao
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, China
- State Key Laboratory of Cellular Stress Biology, School of Life Sciences, Faculty of Medicine and Life Sciences, Xiamen University, Xiamen, Fujian, China
| | - Ruiying Li
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, China
| | - Tingting Li
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, China
| | - Zhimin Li
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, China
| | - Yihan Li
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, China
- *Correspondence: Yihan Li, ; Licao Cui,
| | - Licao Cui
- College of Bioscience and Engineering, Jiangxi Agricultural University, Nanchang, Jiangxi, China
- *Correspondence: Yihan Li, ; Licao Cui,
| |
Collapse
|
7
|
Kaczanowska J, Ganglberger F, Chernomor O, Kargl D, Galik B, Hess A, Moodley Y, von Haeseler A, Bühler K, Haubensak W. Molecular archaeology of human cognitive traits. Cell Rep 2022; 40:111287. [PMID: 36044840 DOI: 10.1016/j.celrep.2022.111287] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/24/2019] [Revised: 05/20/2022] [Accepted: 08/05/2022] [Indexed: 01/06/2023] Open
Abstract
The brains and minds of our human ancestors remain inaccessible for experimental exploration. Therefore, we reconstructed human cognitive evolution by projecting nonsynonymous/synonymous rate ratios (ω values) in mammalian phylogeny onto the anatomically modern human (AMH) brain. This atlas retraces human neurogenetic selection and allows imputation of ancestral evolution in task-related functional networks (FNs). Adaptive evolution (high ω values) is associated with excitatory neurons and synaptic function. It shifted from FNs for motor control in anthropoid ancestry (60-41 mya) to attention in ancient hominoids (26-19 mya) and hominids (19-7.4 mya). Selection in FNs for language emerged with an early hominin ancestor (7.4-1.7 mya) and was later accompanied by adaptive evolution in FNs for strategic thinking during recent (0.8 mya-present) speciation of AMHs. This pattern mirrors increasingly complex cognitive demands and suggests that co-selection for language alongside strategic thinking may have separated AMHs from their archaic Denisovan and Neanderthal relatives.
Collapse
Affiliation(s)
- Joanna Kaczanowska
- Research Institute of Molecular Pathology (IMP), Vienna Biocenter (VBC), Campus-Vienna-Biocenter 1, 1030 Vienna, Austria
| | | | - Olga Chernomor
- Center for Integrative Bioinformatics Vienna (CIBIV), Max Perutz Labs, University of Vienna, Medical University of Vienna, Dr. Bohr Gasse 9, 1030 Vienna, Austria
| | - Dominic Kargl
- Research Institute of Molecular Pathology (IMP), Vienna Biocenter (VBC), Campus-Vienna-Biocenter 1, 1030 Vienna, Austria; Department of Neuronal Cell Biology, Center for Brain Research, Medical University of Vienna, Vienna, Austria
| | - Bence Galik
- Bioinformatics and Scientific Computing, Vienna Biocenter Core Facilities (VBCF), Dr. Bohr Gasse 3, 1030 Vienna, Austria
| | - Andreas Hess
- Institute of Experimental and Clinical Pharmacology and Toxicology, Friedrich-Alexander University Erlangen-Nuremberg, Fahrstrasse 17, 91054 Erlangen, Germany
| | - Yoshan Moodley
- Department of Zoology, University of Venda, Private Bag X5050, Thohoyandou, Republic of South Africa
| | - Arndt von Haeseler
- Center for Integrative Bioinformatics Vienna (CIBIV), Max Perutz Labs, University of Vienna, Medical University of Vienna, Dr. Bohr Gasse 9, 1030 Vienna, Austria; Faculty of Computer Science, University of Vienna, Währinger Str. 29, 1090 Vienna, Austria
| | - Katja Bühler
- VRVis Research Center, Donau-City Strasse 11, 1220 Vienna, Austria
| | - Wulf Haubensak
- Research Institute of Molecular Pathology (IMP), Vienna Biocenter (VBC), Campus-Vienna-Biocenter 1, 1030 Vienna, Austria; Department of Neuronal Cell Biology, Center for Brain Research, Medical University of Vienna, Vienna, Austria.
| |
Collapse
|
8
|
Garrido-Gala J, Higuera JJ, Rodríguez-Franco A, Muñoz-Blanco J, Amil-Ruiz F, Caballero JL. A Comprehensive Study of the WRKY Transcription Factor Family in Strawberry. PLANTS 2022; 11:plants11121585. [PMID: 35736736 PMCID: PMC9229891 DOI: 10.3390/plants11121585] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 05/06/2022] [Revised: 06/10/2022] [Accepted: 06/11/2022] [Indexed: 11/16/2022]
Abstract
WRKY transcription factors play critical roles in plant growth and development or stress responses. Using up-to-date genomic data, a total of 64 and 257 WRKY genes have been identified in the diploid woodland strawberry, Fragaria vesca, and the more complex allo-octoploid commercial strawberry, Fragaria × ananassa cv. Camarosa, respectively. The completeness of the new genomes and annotations has enabled us to perform a more detailed evolutionary and functional study of the strawberry WRKY family members, particularly in the case of the cultivated hybrid, in which homoeologous and paralogous FaWRKY genes have been characterized. Analysis of the available expression profiles has revealed that many strawberry WRKY genes show preferential or tissue-specific expression. Furthermore, significant differential expression of several FaWRKY genes has been clearly detected in fruit receptacles and achenes during the ripening process and pathogen challenged, supporting a precise functional role of these strawberry genes in such processes. Further, an extensive analysis of predicted development, stress and hormone-responsive cis-acting elements in the strawberry WRKY family is shown. Our results provide a deeper and more comprehensive knowledge of the WRKY gene family in strawberry.
Collapse
Affiliation(s)
| | - José-Javier Higuera
- Departamento de Bioquímica y Biología Molecular, Campus Universitario de Rabanales y Campus de Excelencia Internacional Agroalimentario ceiA3, Edificio Severo Ochoa-C6, Universidad de Córdoba, 14071 Córdoba, Spain; (J.-J.H.); (A.R.-F.); (J.M.-B.)
| | - Antonio Rodríguez-Franco
- Departamento de Bioquímica y Biología Molecular, Campus Universitario de Rabanales y Campus de Excelencia Internacional Agroalimentario ceiA3, Edificio Severo Ochoa-C6, Universidad de Córdoba, 14071 Córdoba, Spain; (J.-J.H.); (A.R.-F.); (J.M.-B.)
| | - Juan Muñoz-Blanco
- Departamento de Bioquímica y Biología Molecular, Campus Universitario de Rabanales y Campus de Excelencia Internacional Agroalimentario ceiA3, Edificio Severo Ochoa-C6, Universidad de Córdoba, 14071 Córdoba, Spain; (J.-J.H.); (A.R.-F.); (J.M.-B.)
| | - Francisco Amil-Ruiz
- Unidad de Bioinformática, Servicio Central de Apoyo a la Investigación (SCAI), Universidad de Córdoba, 14071 Córdoba, Spain;
| | - José L. Caballero
- Departamento de Bioquímica y Biología Molecular, Campus Universitario de Rabanales y Campus de Excelencia Internacional Agroalimentario ceiA3, Edificio Severo Ochoa-C6, Universidad de Córdoba, 14071 Córdoba, Spain; (J.-J.H.); (A.R.-F.); (J.M.-B.)
- Correspondence:
| |
Collapse
|
9
|
Yoshioka Y, Suzuki G, Zayasu Y, Yamashita H, Shinzato C. Comparative genomics highlight the importance of lineage-specific gene families in evolutionary divergence of the coral genus, Montipora. BMC Ecol Evol 2022; 22:71. [PMID: 35624412 PMCID: PMC9145168 DOI: 10.1186/s12862-022-02023-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Accepted: 05/13/2022] [Indexed: 12/23/2022] Open
Abstract
Background Scleractinian corals of the genus Montipora (Anthozoa, Cnidaria) possess some unusual biological traits, such as vertical transmission of algal symbionts; however, the genetic bases for those traits remain unknown. We performed extensive comparative genomic analyses among members of the family Acroporidae (Montipora, Acropora, and Astreopora) to explore genomic novelties that might explain unique biological traits of Montipora using improved genome assemblies and gene predictions for M. cactus, M. efflorescens and Astreopora myriophthalma. Results We obtained genomic data for the three species of comparable high quality to other published coral genomes. Comparative genomic analyses revealed that the gene families restricted to Montipora are significantly more numerous than those of Acropora and Astreopora, but their functions are largely unknown. The number of gene families specifically expanded in Montipora was much lower than the number specifically expanded in Acropora. In addition, we found that evolutionary rates of the Montipora-specific gene families were significantly higher than other gene families shared with Acropora and/or Astreopora. Of 40 gene families under positive selection (Ka/Ks ratio > 1) in Montipora, 30 were specifically detected in Montipora-specific gene families. Comparative transcriptomic analysis of early life stages of Montipora, which possesses maternally inherited symbionts, and Acropora, which lacks them, revealed that most gene families continuously expressed in Montipora, but not expressed in Acropora do not have orthologs in Acropora. Among the 30 Montipora-specific gene families under positive selection, 27 are expressed in early life stages. Conclusions Lineage-specific gene families were important to establish the genus Montipora, particularly genes expressed throughout early life stages, which under positive selection, gave rise to biological traits unique to Montipora. Our findings highlight evolutionarily acquired genomic bases that may support symbiosis in these stony corals and provide novel insights into mechanisms of coral-algal symbiosis, the physiological foundation of coral reefs. Supplementary Information The online version contains supplementary material available at 10.1186/s12862-022-02023-8.
Collapse
Affiliation(s)
- Yuki Yoshioka
- Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba, Japan.,Graduate School of Frontier Science, The University of Tokyo, Kashiwa, Chiba, Japan
| | - Go Suzuki
- Fisheries Technology Institute, Japan Fisheries Research and Education Agency, Ishigaki, Okinawa, Japan
| | - Yuna Zayasu
- Marine Genomics Unit, Okinawa Institute of Science and Technology Graduate University, Onna, Okinawa, Japan
| | - Hiroshi Yamashita
- Fisheries Technology Institute, Japan Fisheries Research and Education Agency, Ishigaki, Okinawa, Japan
| | - Chuya Shinzato
- Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba, Japan.
| |
Collapse
|
10
|
Genomic signatures of high-altitude adaptation and chromosomal polymorphism in geladas. Nat Ecol Evol 2022; 6:630-643. [PMID: 35332281 PMCID: PMC9090980 DOI: 10.1038/s41559-022-01703-4] [Citation(s) in RCA: 11] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2021] [Accepted: 02/15/2022] [Indexed: 01/31/2023]
Abstract
Primates have adapted to numerous environments and lifestyles, but very few species are native to high elevations. Here, we investigated high-altitude adaptations in the gelada (Theropithecus gelada), a monkey endemic to the Ethiopian Plateau. We examined genome-wide variation in conjunction with measurements of hematological and morphological traits. Our new gelada reference genome is highly intact and assembled at chromosome-length levels. Unexpectedly, we identified a chromosomal polymorphism in geladas that could potentially contribute to reproductive barriers between populations. Compared to baboons at low altitude, we found that high-altitude geladas exhibit significantly expanded chest circumferences, potentially allowing for greater lung surface area for increased oxygen diffusion. We identified gelada-specific amino acid substitutions in the alpha-chain subunit of adult hemoglobin but found that gelada hemoglobin does not exhibit markedly altered oxygenation properties compared to lowland primates. We also found that geladas at high altitude do not exhibit elevated blood hemoglobin concentrations, in contrast to the normal acclimatization response to hypoxia in lowland primates. The absence of altitude-related polycythemia suggests that geladas are able to sustain adequate tissue-oxygen delivery despite environmental hypoxia. Finally, we identified numerous genes and genomic regions exhibiting accelerated rates of evolution, as well as gene families exhibiting expansions in the gelada lineage, potentially reflecting altitude-related selection. Our findings lend insight into putative mechanisms of high-altitude adaptation while suggesting promising avenues for functional hypoxia research.
Collapse
|
11
|
Stevison LS, Bailey NP, Szpiech ZA, Novak TE, Melnick DJ, Evans BJ, Wall JD. Evolution of genes involved in the unusual genitals of the bear macaque, Macaca arctoides. Ecol Evol 2022; 12:e8897. [PMID: 35646310 PMCID: PMC9130562 DOI: 10.1002/ece3.8897] [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: 04/01/2022] [Accepted: 04/05/2022] [Indexed: 11/30/2022] Open
Abstract
Genital divergence is thought to contribute to reproductive barriers by establishing a “lock‐and‐key" mechanism for reproductive compatibility. One such example, Macaca arctoides, the bear macaque, has compensatory changes in both male and female genital morphology as compared to close relatives. M. arctoides also has a complex evolutionary history, having extensive introgression between the fascicularis and sinica macaque species groups. Here, phylogenetic relationships were analyzed via whole‐genome sequences from five species, including M. arctoides, and two species each from the putative parental species groups. This analysis revealed ~3x more genomic regions supported placement in the sinica species group as compared to the fascicularis species group. Additionally, introgression analysis of the M. arctoides genome revealed it is a mosaic of recent polymorphisms shared with both species groups. To examine the evolution of their unique genital morphology further, the prevalence of candidate genes involved in genital morphology was compared against genome‐wide outliers in various population genetic metrics of diversity, divergence, introgression, and selection, while accounting for background variation in recombination rate. This analysis identified 67 outlier genes, including several genes that influence baculum morphology in mice, which were of interest since the bear macaque has the longest primate baculum. The mean of four of the seven population genetic metrics was statistically different in the candidate genes as compared to the rest of the genome, suggesting that genes involved in genital morphology have increased divergence and decreased diversity beyond expectations. These results highlight specific genes that may have played a role in shaping the unique genital morphology in the bear macaque.
Collapse
Affiliation(s)
- Laurie S Stevison
- Department of Biological Sciences Auburn University Auburn Alabama USA
| | - Nick P Bailey
- Department of Biological Sciences Auburn University Auburn Alabama USA
| | - Zachary A Szpiech
- Department of Biological Sciences Auburn University Auburn Alabama USA.,Department of Biology Pennsylvania State University University Park Pennsylvania USA.,Institute for Computational and Data Sciences Pennsylvania State University University Park Pennsylvania USA
| | - Taylor E Novak
- Department of Biological Sciences Auburn University Auburn Alabama USA
| | - Don J Melnick
- Department of Ecology, Evolution, and Environmental Biology Columbia University New York New York USA
| | - Ben J Evans
- Biology Department McMaster University Hamilton Ontario Canada
| | - Jeffrey D Wall
- Institute for Human Genetics University of California, San Francisco San Francisco California USA
| |
Collapse
|
12
|
Luchetti A, Forni G, Martelossi J, Savojardo C, Martelli PL, Casadio R, Skaist AM, Wheelan SJ, Mantovani B. Comparative genomics of tadpole shrimps (Crustacea, Branchiopoda, Notostraca): Dynamic genome evolution against the backdrop of morphological stasis. Genomics 2021; 113:4163-4172. [PMID: 34748900 DOI: 10.1016/j.ygeno.2021.11.001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/14/2021] [Revised: 10/26/2021] [Accepted: 11/02/2021] [Indexed: 12/21/2022]
Abstract
This analysis presents five genome assemblies of four Notostraca taxa. Notostraca origin dates to the Permian/Upper Devonian and the extant forms show a striking morphological similarity to fossil taxa. The comparison of sequenced genomes with other Branchiopoda genomes shows that, despite the morphological stasis, Notostraca share a dynamic genome evolution with high turnover for gene families' expansion/contraction and a transposable elements content comparable to other branchiopods. While Notostraca substitutions rate appears similar or lower in comparison to other branchiopods, a subset of genes shows a faster evolutionary pace, highlighting the difficulty of generalizing about genomic stasis versus dynamism. Moreover, we found that the variation of Triops cancriformis transposable elements content appeared linked to reproductive strategies, in line with theoretical expectations. Overall, besides providing new genomic resources for the study of these organisms, which appear relevant for their ecology and evolution, we also confirmed the decoupling of morphological and molecular evolution.
Collapse
Affiliation(s)
- Andrea Luchetti
- Department of Biological, Geological and Environmental Sciences, University of Bologna, via Selmi 3, 40126 Bologna, Italy.
| | - Giobbe Forni
- Department of Biological, Geological and Environmental Sciences, University of Bologna, via Selmi 3, 40126 Bologna, Italy
| | - Jacopo Martelossi
- Department of Biological, Geological and Environmental Sciences, University of Bologna, via Selmi 3, 40126 Bologna, Italy
| | - Castrense Savojardo
- Biocomputing Group, Department of Pharmacy and Biotechnology, University of Bologna, Italy
| | - Pier Luigi Martelli
- Biocomputing Group, Department of Pharmacy and Biotechnology, University of Bologna, Italy
| | - Rita Casadio
- Biocomputing Group, Department of Pharmacy and Biotechnology, University of Bologna, Italy
| | - Alyza M Skaist
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, The Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Sarah J Wheelan
- Department of Oncology, Sidney Kimmel Comprehensive Cancer Center, The Johns Hopkins University School of Medicine, Baltimore, MD 21205, USA
| | - Barbara Mantovani
- Department of Biological, Geological and Environmental Sciences, University of Bologna, via Selmi 3, 40126 Bologna, Italy
| |
Collapse
|
13
|
Bailey NP, Stevison LS. Mitonuclear conflict in a macaque species exhibiting phylogenomic discordance. J Evol Biol 2021; 34:1568-1579. [PMID: 34379829 DOI: 10.1111/jeb.13914] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2021] [Accepted: 07/28/2021] [Indexed: 12/19/2022]
Abstract
Speciation and hybridization are intertwined processes in the study of evolution. Hybridization between sufficiently diverged populations can result in genomic conflict within offspring, causing reduced viability and fertility, thus increasing divergence between populations. Conflicts between mitochondrial and nuclear genes are increasingly found to play a role in this process in various systems. We examine the possibility of this conflict in the bear macaque, Macaca arctoides (Primates: Cercopithecidae), a primate species exhibiting mitonuclear discordance due to extensive hybridization with species in the sinica and fascicularis groups. Here, divergence, introgression and natural selection of mitonuclear genes (N = 160) relative to nuclear control genes (N = 144) were analysed to determine whether there are evolutionary processes involved in resolving the potential conflict caused by mitonuclear discordance. Nucleotide divergence of mitonuclear genes is increased relative to control nuclear genes between M. arctoides and the species sharing its nuclear ancestry (p = 0.007), consistent with genetic conflict. However, measures of introgression and selection do not identify large-scale co-introgression or co-evolution as means to resolve mitonuclear conflict. Nonetheless, mitochondrial tRNA synthetases stand out in analyses using dN/dS and extended branch lengths as potential targets of selection. The methodology implemented provides a framework that can be used to examine the effects of mitonuclear co-introgression and co-evolution on a genomic scale in a variety of systems.
Collapse
Affiliation(s)
- Nick P Bailey
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| | - Laurie S Stevison
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| |
Collapse
|
14
|
New Approaches for Inferring Phylogenies in the Presence of Paralogs. Trends Genet 2021; 37:174-187. [DOI: 10.1016/j.tig.2020.08.012] [Citation(s) in RCA: 28] [Impact Index Per Article: 9.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/08/2020] [Revised: 08/13/2020] [Accepted: 08/19/2020] [Indexed: 12/18/2022]
|
15
|
Phylomitogenomics provides new perspectives on the Euphasmatodea radiation (Insecta: Phasmatodea). Mol Phylogenet Evol 2020; 155:106983. [PMID: 33059069 DOI: 10.1016/j.ympev.2020.106983] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2019] [Revised: 09/11/2020] [Accepted: 10/06/2020] [Indexed: 11/20/2022]
Abstract
Phasmatodea species diversity lies almost entirely within its suborder Euphasmatodea, which exhibits a pantropical distribution and is considered to derive from a recent and rapid evolutionary radiation. To shed light on Euphasmatodea origins and diversification, we assembled the mitogenomes of 17 species from transcriptomic sequencing data and analysed them along with 22 already available Phasmatodea mitogenomes and 33 mitogenomes representing most of the Polyneoptera lineages. Maximum Likelihood and Bayesian Inference approaches retrieved consistent topologies, both showing the widespread conflict between phylogenetic approaches and traditional systematics. We performed a divergence time analysis leveraging ten fossil specimens representative of most polyneopteran lineages: the time tree obtained supports an older radiation of the clade with respect to previous hypotheses. Euphasmatodea diversification is inferred to have started ~ 187 million years ago, suggesting that the Triassic-Jurassic mass extinction and the breakup of Pangea could have contributed to the process. We also investigated Euphasmatodea mitogenomes patterns of dN, dS and dN/dS ratio throughout our time-tree, trying to characterize the selective regime which may have shaped the clade evolution.
Collapse
|
16
|
Cardoso GA, Deszo MS, Torres TT. Evolution of coding sequence and gene expression of blowflies and botflies with contrasting feeding habits. Genomics 2020; 113:699-706. [PMID: 33022358 DOI: 10.1016/j.ygeno.2020.09.066] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/15/2020] [Revised: 09/05/2020] [Accepted: 09/28/2020] [Indexed: 10/23/2022]
Abstract
The Oestroidea superfamily is characterized by the diversity of feeding preferences among closely-related species; these flies are saprophagous, obligate parasites, or facultative parasites. We used gene expression and coding sequence data from five species (Cochliomyia hominivorax, Chrysomya megacephala, Lucilia cuprina, Dermatobia hominis, and Oestrus ovis) to identify underlying genetic differences involved in the diverse lifestyles. We tested whether 1287 orthologs have different expression and evolutionary constraints under different scenarios. We found two up-regulated genes; one in species causing cutaneous myiasis that is involved in iron transportation/metabolization (ferritin), and another in species causing traumatic myiasis that responds to reduced oxygen levels (anoxia up-regulated-like). Our evolutionary analysis showed a similar result. In the Co. hominivorax branch, we found one gene with the same function as ferritin that may be evolving under positive selection, spook. This is the first step towards understanding origins and evolution of parasitic strategy diversity in Oestroidea.
Collapse
Affiliation(s)
- Gisele Antoniazzi Cardoso
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo (USP), São Paulo, SP 05508-090, Brazil
| | - Marina Santos Deszo
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo (USP), São Paulo, SP 05508-090, Brazil
| | - Tatiana Teixeira Torres
- Departamento de Genética e Biologia Evolutiva, Instituto de Biociências, Universidade de São Paulo (USP), São Paulo, SP 05508-090, Brazil.
| |
Collapse
|
17
|
David KT, Oaks JR, Halanych KM. Patterns of gene evolution following duplications and speciations in vertebrates. PeerJ 2020; 8:e8813. [PMID: 32266119 PMCID: PMC7120047 DOI: 10.7717/peerj.8813] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Accepted: 02/27/2020] [Indexed: 11/24/2022] Open
Abstract
BACKGROUND Eukaryotic genes typically form independent evolutionary lineages through either speciation or gene duplication events. Generally, gene copies resulting from speciation events (orthologs) are expected to maintain similarity over time with regard to sequence, structure and function. After a duplication event, however, resulting gene copies (paralogs) may experience a broader set of possible fates, including partial (subfunctionalization) or complete loss of function, as well as gain of new function (neofunctionalization). This assumption, known as the Ortholog Conjecture, is prevalent throughout molecular biology and notably plays an important role in many functional annotation methods. Unfortunately, studies that explicitly compare evolutionary processes between speciation and duplication events are rare and conflicting. METHODS To provide an empirical assessment of ortholog/paralog evolution, we estimated ratios of nonsynonymous to synonymous substitutions (ω = dN/dS) for 251,044 lineages in 6,244 gene trees across 77 vertebrate taxa. RESULTS Overall, we found ω to be more similar between lineages descended from speciation events (p < 0.001) than lineages descended from duplication events, providing strong support for the Ortholog Conjecture. The asymmetry in ω following duplication events appears to be largely driven by an increase along one of the paralogous lineages, while the other remains similar to the parent. This trend is commonly associated with neofunctionalization, suggesting that gene duplication is a significant mechanism for generating novel gene functions.
Collapse
Affiliation(s)
- Kyle T. David
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| | - Jamie R. Oaks
- Department of Biological Sciences, Auburn University, Auburn, AL, USA
| | | |
Collapse
|
18
|
Okamura Y, Sato A, Tsuzuki N, Murakami M, Heidel‐Fischer H, Vogel H. Molecular signatures of selection associated with host plant differences in
Pieris
butterflies. Mol Ecol 2019; 28:4958-4970. [DOI: 10.1111/mec.15268] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/07/2019] [Revised: 09/24/2019] [Accepted: 09/30/2019] [Indexed: 11/29/2022]
Affiliation(s)
- Yu Okamura
- Department of Entomology Max Planck Institute for Chemical Ecology Jena Germany
- Community Ecology Lab Faculty of Science Chiba University Chiba Japan
| | - Ai Sato
- Community Ecology Lab Faculty of Science Chiba University Chiba Japan
| | - Natsumi Tsuzuki
- Community Ecology Lab Faculty of Science Chiba University Chiba Japan
| | - Masashi Murakami
- Community Ecology Lab Faculty of Science Chiba University Chiba Japan
| | - Hanna Heidel‐Fischer
- Department of Entomology Max Planck Institute for Chemical Ecology Jena Germany
- Leibniz Institute for Natural Product Research and Infection Biology Hans Knöll Institute (HKI) Jena Germany
| | - Heiko Vogel
- Department of Entomology Max Planck Institute for Chemical Ecology Jena Germany
| |
Collapse
|
19
|
Ely ZA, Moon JM, Sliwoski GR, Sangha AK, Shen XX, Labella AL, Meiler J, Capra JA, Rokas A. The Impact of Natural Selection on the Evolution and Function of Placentally Expressed Galectins. Genome Biol Evol 2019; 11:2574-2592. [PMID: 31504490 PMCID: PMC6751361 DOI: 10.1093/gbe/evz183] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 08/19/2019] [Indexed: 01/03/2023] Open
Abstract
Immunity genes have repeatedly experienced natural selection during mammalian evolution. Galectins are carbohydrate-binding proteins that regulate diverse immune responses, including maternal-fetal immune tolerance in placental pregnancy. Seven human galectins, four conserved across vertebrates and three specific to primates, are involved in placental development. To comprehensively study the molecular evolution of these galectins, both across mammals and within humans, we conducted a series of between- and within-species evolutionary analyses. By examining patterns of sequence evolution between species, we found that primate-specific galectins showed uniformly high substitution rates, whereas two of the four other galectins experienced accelerated evolution in primates. By examining human population genomic variation, we found that galectin genes and variants, including variants previously linked to immune diseases, showed signatures of recent positive selection in specific human populations. By examining one nonsynonymous variant in Galectin-8 previously associated with autoimmune diseases, we further discovered that it is tightly linked to three other nonsynonymous variants; surprisingly, the global frequency of this four-variant haplotype is ∼50%. To begin understanding the impact of this major haplotype on Galectin-8 protein structure, we modeled its 3D protein structure and found that it differed substantially from the reference protein structure. These results suggest that placentally expressed galectins experienced both ancient and more recent selection in a lineage- and population-specific manner. Furthermore, our discovery that the major Galectin-8 haplotype is structurally distinct from and more commonly found than the reference haplotype illustrates the significance of understanding the evolutionary processes that sculpted variants associated with human genetic disease.
Collapse
Affiliation(s)
- Zackery A Ely
- Department of Biological Sciences, Vanderbilt University
| | - Jiyun M Moon
- Department of Biological Sciences, Vanderbilt University
| | | | - Amandeep K Sangha
- Department of Chemistry, Vanderbilt University
- Center for Structural Biology, Vanderbilt University
| | - Xing-Xing Shen
- Department of Biological Sciences, Vanderbilt University
| | | | - Jens Meiler
- Department of Chemistry, Vanderbilt University
- Center for Structural Biology, Vanderbilt University
| | - John A Capra
- Department of Biological Sciences, Vanderbilt University
- Department of Biomedical Informatics, Vanderbilt University School of Medicine
| | - Antonis Rokas
- Department of Biological Sciences, Vanderbilt University
- Department of Biomedical Informatics, Vanderbilt University School of Medicine
| |
Collapse
|
20
|
Marks RA, Smith JJ, Cronk Q, Grassa CJ, McLetchie DN. Genome of the tropical plant Marchantia inflexa: implications for sex chromosome evolution and dehydration tolerance. Sci Rep 2019; 9:8722. [PMID: 31217536 PMCID: PMC6584576 DOI: 10.1038/s41598-019-45039-9] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Accepted: 05/29/2019] [Indexed: 01/29/2023] Open
Abstract
We present a draft genome assembly for the tropical liverwort, Marchantia inflexa, which adds to a growing body of genomic resources for bryophytes and provides an important perspective on the evolution and diversification of land plants. We specifically address questions related to sex chromosome evolution, sexual dimorphisms, and the genomic underpinnings of dehydration tolerance. This assembly leveraged the recently published genome of related liverwort, M. polymorpha, to improve scaffolding and annotation, aid in the identification of sex-linked sequences, and quantify patterns of sequence differentiation within Marchantia. We find that genes on sex chromosomes are under greater diversifying selection than autosomal and organellar genes. Interestingly, this is driven primarily by divergence of male-specific genes, while divergence of other sex-linked genes is similar to autosomal genes. Through analysis of sex-specific read coverage, we identify and validate genetic sex markers for M. inflexa, which will enable diagnosis of sex for non-reproductive individuals. To investigate dehydration tolerance, we capitalized on a difference between genetic lines, which allowed us to identify multiple dehydration associated genes two of which were sex-linked, suggesting that dehydration tolerance may be impacted by sex-specific genes.
Collapse
Affiliation(s)
- Rose A Marks
- Department of Biology, University of Kentucky, 101 Thomas Hunt Morgan Building, Lexington, KY, 40506, USA.
| | - Jeramiah J Smith
- Department of Biology, University of Kentucky, 101 Thomas Hunt Morgan Building, Lexington, KY, 40506, USA
| | - Quentin Cronk
- Department of Botany, University of British Columbia, 6270 University Boulevard, Vancouver, BC, V6T 1Z4, Canada
| | - Christopher J Grassa
- Department of Botany, University of British Columbia, 6270 University Boulevard, Vancouver, BC, V6T 1Z4, Canada
- Department of Organismic and Evolutionary Biology, Harvard University Herbaria, 22 Divinity Avenue, Cambridge, MA, 02138, USA
| | - D Nicholas McLetchie
- Department of Biology, University of Kentucky, 101 Thomas Hunt Morgan Building, Lexington, KY, 40506, USA
| |
Collapse
|
21
|
Little Evidence of Antagonistic Selection in the Evolutionary Strata of Fungal Mating-Type Chromosomes ( Microbotryum lychnidis-dioicae). G3-GENES GENOMES GENETICS 2019; 9:1987-1998. [PMID: 31015196 PMCID: PMC6553529 DOI: 10.1534/g3.119.400242] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 11/18/2022]
Abstract
Recombination suppression on sex chromosomes often extends in a stepwise manner, generating evolutionary strata of differentiation between sex chromosomes. Sexual antagonism is a widely accepted explanation for evolutionary strata, postulating that sets of genes beneficial in only one sex are successively linked to the sex-determining locus. The anther-smut fungus Microbotryum lychnidis-dioicae has mating-type chromosomes with evolutionary strata, only some of which link mating-type genes. Male and female roles are non-existent in this fungus, but mating-type antagonistic selection can also generate evolutionary strata, although the life cycle of the fungus suggests it should be restricted to few traits. Here, we tested the hypothesis that mating-type antagonism may have triggered recombination suppression beyond mating-type genes in M. lychnidis-dioicae by searching for footprints of antagonistic selection in evolutionary strata not linking mating-type loci. We found that these evolutionary strata (i) were not enriched in genes upregulated in the haploid phase, where cells are of alternative mating types, (ii) carried no gene differentially expressed between mating types, and (iii) carried no genes displaying footprints of specialization in terms of protein sequences (dN/dS) between mating types after recommended filtering. Without filtering, eleven genes showed signs of positive selection in the strata not linking mating-type genes, which constituted an enrichment compared to autosomes, but their functions were not obviously involved in antagonistic selection. Thus, we found no strong evidence that antagonistic selection has contributed to extending recombination suppression beyond mating-type genes. Alternative hypotheses should therefore be explored to improve our understanding of the sex-related chromosome evolution.
Collapse
|
22
|
Abstract
An attractive and long-standing hypothesis regarding the evolution of genes after duplication posits that the duplication event creates new evolutionary possibilities by releasing a copy of the gene from constraint. Apparent support was found in numerous analyses, particularly, the observation of higher rates of evolution in duplicated as compared with singleton genes. Could it, instead, be that more duplicable genes (owing to mutation, fixation, or retention biases) are intrinsically faster evolving? To uncouple the measurement of rates of evolution from the determination of duplicate or singleton status, we measure the rates of evolution in singleton genes in outgroup primate lineages but classify these genes as to whether they have duplicated or not in a crown group of great apes. We find that rates of evolution are higher in duplicable genes prior to the duplication event. In part this is owing to a negative correlation between coding sequence length and rate of evolution, coupled with a bias toward smaller genes being more duplicable. The effect is masked by difference in expression rate between duplicable genes and singletons. Additionally, in contradiction to the classical assumption, we find no convincing evidence for an increase in dN/dS after duplication, nor for rate asymmetry between duplicates. We conclude that high rates of evolution of duplicated genes are not solely a consequence of the duplication event, but are rather a predictor of duplicability. These results are consistent with a model in which successful gene duplication events in mammals are skewed toward events of minimal phenotypic impact.
Collapse
Affiliation(s)
- Áine N O'Toole
- Department of Genetics, Smurfit Institute of Genetics, Trinity College Dublin, Dublin, Ireland
| | - Laurence D Hurst
- The Milner Centre for Evolution, Department of Biology and Biochemistry, University of Bath, Bath, Somerset, United Kingdom
| | - Aoife McLysaght
- Department of Genetics, Smurfit Institute of Genetics, Trinity College Dublin, Dublin, Ireland
| |
Collapse
|
23
|
Marcionetti A, Rossier V, Roux N, Salis P, Laudet V, Salamin N. Insights into the Genomics of Clownfish Adaptive Radiation: Genetic Basis of the Mutualism with Sea Anemones. Genome Biol Evol 2019; 11:869-882. [PMID: 30830203 PMCID: PMC6430985 DOI: 10.1093/gbe/evz042] [Citation(s) in RCA: 21] [Impact Index Per Article: 4.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 02/28/2019] [Indexed: 02/06/2023] Open
Abstract
Clownfishes are an iconic group of coral reef fishes, especially known for their mutualism with sea anemones. This mutualism is particularly interesting as it likely acted as the key innovation that triggered clownfish adaptive radiation. Indeed, after the acquisition of the mutualism, clownfishes diversified into multiple ecological niches linked with host and habitat use. However, despite the importance of this mutualism, the genetic mechanisms allowing clownfishes to interact with sea anemones are still unclear. Here, we used a comparative genomics and molecular evolutionary analyses to investigate the genetic basis of clownfish mutualism with sea anemones. We assembled and annotated the genome of nine clownfish species and one closely related outgroup. Orthologous genes inferred between these species and additional publicly available teleost genomes resulted in almost 16,000 genes that were tested for positively selected substitutions potentially involved in the adaptation of clownfishes to live in sea anemones. We identified 17 genes with a signal of positive selection at the origin of clownfish radiation. Two of them (Versican core protein and Protein O-GlcNAse) show particularly interesting functions associated with N-acetylated sugars, which are known to be involved in sea anemone discharge of toxins. This study provides the first insights into the genetic mechanisms of clownfish mutualism with sea anemones. Indeed, we identified the first candidate genes likely to be associated with clownfish protection form sea anemones, and thus the evolution of their mutualism. Additionally, the genomic resources acquired represent a valuable resource for further investigation of the genomic basis of clownfish adaptive radiation.
Collapse
Affiliation(s)
- Anna Marcionetti
- Department of Computational Biology, Génopode, University of Lausanne, Switzerland
| | - Victor Rossier
- Department of Computational Biology, Génopode, University of Lausanne, Switzerland
| | - Natacha Roux
- Observatoire Océanologique de Banyuls-sur-Mer, UMR CNRS 7232 BIOM, Sorbonne University, Banyuls-sur-Mer, France
| | - Pauline Salis
- Observatoire Océanologique de Banyuls-sur-Mer, UMR CNRS 7232 BIOM, Sorbonne University, Banyuls-sur-Mer, France
| | - Vincent Laudet
- Observatoire Océanologique de Banyuls-sur-Mer, UMR CNRS 7232 BIOM, Sorbonne University, Banyuls-sur-Mer, France
| | - Nicolas Salamin
- Department of Computational Biology, Génopode, University of Lausanne, Switzerland
| |
Collapse
|
24
|
Mori S, Matsunami M. Signature of positive selection in mitochondrial DNA in Cetartiodactyla. Genes Genet Syst 2018; 93:65-73. [PMID: 29643269 DOI: 10.1266/ggs.17-00015] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/23/2022] Open
Abstract
Acceleration of the amino acid substitution rate is a good indicator of positive selection in adaptive evolutionary changes of functional genes. Genomic information about mammals has become readily available in recent years, as many researchers have attempted to clarify the adaptive evolution of mammals by examining evolutionary rate change based on multiple loci. The order Cetartiodactyla (Artiodactyla and Cetacea) is one of the most diverse orders of mammals. Species in this order are found throughout all continents and seas, except Antarctica, and they exhibit wide variation in morphology and habitat. Here, we focused on the metabolism-related genes of mitochondrial DNA (mtDNA) in species of the order Cetartiodactyla using 191 mtDNA sequences available in databases. Based on comparisons of the dN/dS ratio (ω) in 12 protein-coding genes, ATP8 was shown to have a higher ω value (ω = 0.247) throughout Cetartiodactyla than the other 11 genes (ω < 0.05). In a branch-site analysis of ATP8 sequences, a markedly higher ω value of 0.801 was observed in the ancestral lineage of the clade of Cetacea, which is indicative of adaptive evolution. Through efforts to detect positively selected amino acids, codon positions 52 and 54 of ATP8 were shown to have experienced positive selective pressure during the course of evolution; multiple substitutions have occurred at these sites throughout the cetacean lineage. At position 52, glutamic acid was replaced with asparagine, and, at position 54, lysine was replaced with non-charged amino acids. These sites are conserved in most Artiodactyla. These results imply that the ancestor of cetaceans underwent accelerated amino acid changes in ATP8 and replacements at codons 52 and 54, which adjusted metabolism to adapt to the marine environment.
Collapse
Affiliation(s)
- Satoko Mori
- Laboratory of Ecology and Genetics, Graduate School of Environmental Science, Hokkaido University
| | - Masatoshi Matsunami
- Laboratory of Ecology and Genetics, Graduate School of Environmental Science, Hokkaido University.,Graduate School of Medicine, University of the Ryukyus
| |
Collapse
|
25
|
Sahm A, Bens M, Platzer M, Szafranski K. PosiGene: automated and easy-to-use pipeline for genome-wide detection of positively selected genes. Nucleic Acids Res 2017; 45:e100. [PMID: 28334822 PMCID: PMC5499814 DOI: 10.1093/nar/gkx179] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/29/2016] [Accepted: 03/09/2017] [Indexed: 11/12/2022] Open
Abstract
Many comparative genomics studies aim to find the genetic basis of species-specific phenotypic traits. A prevailing strategy is to search genome-wide for genes that evolved under positive selection based on the non-synonymous to synonymous substitution ratio. However, incongruent results largely due to high false positive rates indicate the need for standardization of quality criteria and software tools. Main challenges are the ortholog and isoform assignment, the high sensitivity of the statistical models to alignment errors and the imperative to parallelize large parts of the software. We developed the software tool PosiGene that (i) detects positively selected genes (PSGs) on genome-scale, (ii) allows analysis of specific evolutionary branches, (iii) can be used in arbitrary species contexts and (iv) offers visualization of the results for further manual validation and biological interpretation. We exemplify PosiGene's performance using simulated and real data. In the simulated data approach, we determined a false positive rate <1%. With real data, we found that 68.4% of the PSGs detected by PosiGene, were shared by at least one previous study that used the same set of species. PosiGene is a user-friendly, reliable tool for reproducible genome-wide identification of PSGs and freely available at https://github.com/gengit/PosiGene.
Collapse
Affiliation(s)
- Arne Sahm
- Leibniz Institute on Aging, Fritz Lipmann Institute, 07745 Jena, Germany
| | - Martin Bens
- Leibniz Institute on Aging, Fritz Lipmann Institute, 07745 Jena, Germany
| | - Matthias Platzer
- Leibniz Institute on Aging, Fritz Lipmann Institute, 07745 Jena, Germany
| | - Karol Szafranski
- Leibniz Institute on Aging, Fritz Lipmann Institute, 07745 Jena, Germany
| |
Collapse
|
26
|
De La Torre AR, Li Z, Van de Peer Y, Ingvarsson PK. Contrasting Rates of Molecular Evolution and Patterns of Selection among Gymnosperms and Flowering Plants. Mol Biol Evol 2017; 34:1363-1377. [PMID: 28333233 PMCID: PMC5435085 DOI: 10.1093/molbev/msx069] [Citation(s) in RCA: 111] [Impact Index Per Article: 15.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
The majority of variation in rates of molecular evolution among seed plants remains both unexplored and unexplained. Although some attention has been given to flowering plants, reports of molecular evolutionary rates for their sister plant clade (gymnosperms) are scarce, and to our knowledge differences in molecular evolution among seed plant clades have never been tested in a phylogenetic framework. Angiosperms and gymnosperms differ in a number of features, of which contrasting reproductive biology, life spans, and population sizes are the most prominent. The highly conserved morphology of gymnosperms evidenced by similarity of extant species to fossil records and the high levels of macrosynteny at the genomic level have led scientists to believe that gymnosperms are slow-evolving plants, although some studies have offered contradictory results. Here, we used 31,968 nucleotide sites obtained from orthologous genes across a wide taxonomic sampling that includes representatives of most conifers, cycads, ginkgo, and many angiosperms with a sequenced genome. Our results suggest that angiosperms and gymnosperms differ considerably in their rates of molecular evolution per unit time, with gymnosperm rates being, on average, seven times lower than angiosperm species. Longer generation times and larger genome sizes are some of the factors explaining the slow rates of molecular evolution found in gymnosperms. In contrast to their slow rates of molecular evolution, gymnosperms possess higher substitution rate ratios than angiosperm taxa. Finally, our study suggests stronger and more efficient purifying and diversifying selection in gymnosperm than in angiosperm species, probably in relation to larger effective population sizes.
Collapse
Affiliation(s)
- Amanda R De La Torre
- Department of Plant Sciences, University of California-Davis, Davis, CA.,Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden
| | - Zhen Li
- Department of Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Yves Van de Peer
- Department of Plant Systems Biology, VIB, Ghent, Belgium.,Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium.,Genomics Research Institute, University of Pretoria, Hatfield Campus, Pretoria, South Africa
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Umeå, Sweden.,Department of Plant Biology, Uppsala Biocenter, Swedish University of Agricultural Sciences, Uppsala, Sweden
| |
Collapse
|
27
|
Ishengoma E, Agaba M, Cavener DR. Evolutionary analysis of vision genes identifies potential drivers of visual differences between giraffe and okapi. PeerJ 2017; 5:e3145. [PMID: 28396824 PMCID: PMC5385128 DOI: 10.7717/peerj.3145] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/13/2016] [Accepted: 03/04/2017] [Indexed: 12/20/2022] Open
Abstract
Background The capacity of visually oriented species to perceive and respond to visual signal is integral to their evolutionary success. Giraffes are closely related to okapi, but the two species have broad range of phenotypic differences including their visual capacities. Vision studies rank giraffe’s visual acuity higher than all other artiodactyls despite sharing similar vision ecological determinants with many of them. The extent to which the giraffe’s unique visual capacity and its difference with okapi is reflected by changes in their vision genes is not understood. Methods The recent availability of giraffe and okapi genomes provided opportunity to identify giraffe and okapi vision genes. Multiple strategies were employed to identify thirty-six candidate mammalian vision genes in giraffe and okapi genomes. Quantification of selection pressure was performed by a combination of branch-site tests of positive selection and clade models of selection divergence through comparing giraffe and okapi vision genes and orthologous sequences from other mammals. Results Signatures of selection were identified in key genes that could potentially underlie giraffe and okapi visual adaptations. Importantly, some genes that contribute to optical transparency of the eye and those that are critical in light signaling pathway were found to show signatures of adaptive evolution or selection divergence. Comparison between giraffe and other ruminants identifies significant selection divergence in CRYAA and OPN1LW. Significant selection divergence was identified in SAG while positive selection was detected in LUM when okapi is compared with ruminants and other mammals. Sequence analysis of OPN1LW showed that at least one of the sites known to affect spectral sensitivity of the red pigment is uniquely divergent between giraffe and other ruminants. Discussion By taking a systemic approach to gene function in vision, the results provide the first molecular clues associated with giraffe and okapi vision adaptations. At least some of the genes that exhibit signature of selection may reflect adaptive response to differences in giraffe and okapi habitat. We hypothesize that requirement for long distance vision associated with predation and communication with conspecifics likely played an important role in the adaptive pressure on giraffe vision genes.
Collapse
Affiliation(s)
- Edson Ishengoma
- The School of Life Sciences and Bio-Engineering, Nelson Mandela African Institution of Science and Technology, Arusha, Tanzania; Mkwawa University of College of Education, University of Dar-es-Salaam, Iringa, Tanzania
| | - Morris Agaba
- The School of Life Sciences and Bio-Engineering, Nelson Mandela African Institution of Science and Technology , Arusha , Tanzania
| | - Douglas R Cavener
- Department of Biology and the Huck Institute of Life Sciences, Pennsylvania State University , University Park , PA , United States
| |
Collapse
|
28
|
Philippon H, Souvane A, Brochier-Armanet C, Perrière G. IsoSel: Protein Isoform Selector for phylogenetic reconstructions. PLoS One 2017; 12:e0174250. [PMID: 28323858 PMCID: PMC5360266 DOI: 10.1371/journal.pone.0174250] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2016] [Accepted: 03/06/2017] [Indexed: 12/05/2022] Open
Abstract
The reliability of molecular phylogenies is strongly dependent on the quality of the assembled datasets. In the case of eukaryotes, the selection of only one protein isoform per genomic locus is mandatory to avoid biases linked to redundancy. Here, we present IsoSel, a tool devoted to the selection of alternative isoforms in the context of phylogenetic reconstruction. It provides a better alternative to the widely used approach consisting in the selection of the longest isoforms and it performs better than Guidance, its only available counterpart. IsoSel is publicly available at http://doua.prabi.fr/software/isosel.
Collapse
Affiliation(s)
- Héloïse Philippon
- Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive, UMR 5558, F-69622, Villeurbanne, France
| | - Alexia Souvane
- Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive, UMR 5558, F-69622, Villeurbanne, France
| | - Céline Brochier-Armanet
- Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive, UMR 5558, F-69622, Villeurbanne, France
| | - Guy Perrière
- Univ. Lyon, Université Claude Bernard Lyon 1, CNRS, Laboratoire de Biométrie et Biologie Evolutive, UMR 5558, F-69622, Villeurbanne, France
| |
Collapse
|
29
|
Mu H, Sun J, Heras H, Chu KH, Qiu JW. An integrated proteomic and transcriptomic analysis of perivitelline fluid proteins in a freshwater gastropod laying aerial eggs. J Proteomics 2017; 155:22-30. [PMID: 28095328 DOI: 10.1016/j.jprot.2017.01.006] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2016] [Revised: 01/04/2017] [Accepted: 01/08/2017] [Indexed: 11/27/2022]
Abstract
Proteins of the egg perivitelline fluid (PVF) that surrounds the embryo are critical for embryonic development in many animals, but little is known about their identities. Using an integrated proteomic and transcriptomic approach, we identified 64 proteins from the PVF of Pomacea maculata, a freshwater snail adopting aerial oviposition. Proteins were classified into eight functional groups: major multifunctional perivitellin subunits, immune response, energy metabolism, protein degradation, oxidation-reduction, signaling and binding, transcription and translation, and others. Comparison of gene expression levels between tissues showed that 22 PVF genes were exclusively expressed in albumen gland, the female organ that secretes PVF. Base substitution analysis of PVF and housekeeping genes between P. maculata and its closely related species Pomacea canaliculata showed that the reproductive proteins had a higher mean evolutionary rate. Predicted 3D structures of selected PVF proteins showed that some nonsynonymous substitutions are located at or near the binding regions that may affect protein function. The proteome and sequence divergence analysis revealed a substantial amount of maternal investment in embryonic nutrition and defense, and higher adaptive selective pressure on PVF protein-coding genes when compared with housekeeping genes, providing insight into the adaptations associated with the unusual reproductive strategy in these mollusks. SIGNIFICANCE There has been great interest in studying reproduction-related proteins as such studies may not only answer fundamental questions about speciation and evolution, but also solve practical problems of animal infertility and pest outbreak. Our study has demonstrated the effectiveness of an integrated proteomic and transcriptomic approach in understanding the heavy maternal investment of proteins in the eggs of a non-model snail, and how the reproductive proteins may have evolved during the transition from laying underwater eggs to aerial eggs.
Collapse
Affiliation(s)
- Huawei Mu
- Department of Biology, Hong Kong Baptist University, Hong Kong, China
| | - Jin Sun
- Division of Life Science, Hong Kong University of Science and Technology, Hong Kong, China
| | - Horacio Heras
- Instituto de Investigaciones Bioquímicas de La Plata (INIBIOLP), Universidad Nacional de La Plata (UNLP)-CONICET CCT-La Plata, La Plata, Argentina; Cátedra de Química Biológica, Facultad de Ciencias Naturales y Museo, UNLP, Argentina
| | - Ka Hou Chu
- Simon F.S. Li Marine Science Laboratory, School of Life Sciences, The Chinese University of Hong Kong, Hong Kong, China
| | - Jian-Wen Qiu
- Department of Biology, Hong Kong Baptist University, Hong Kong, China.
| |
Collapse
|
30
|
Abascal F, Corvelo A, Cruz F, Villanueva-Cañas JL, Vlasova A, Marcet-Houben M, Martínez-Cruz B, Cheng JY, Prieto P, Quesada V, Quilez J, Li G, García F, Rubio-Camarillo M, Frias L, Ribeca P, Capella-Gutiérrez S, Rodríguez JM, Câmara F, Lowy E, Cozzuto L, Erb I, Tress ML, Rodriguez-Ales JL, Ruiz-Orera J, Reverter F, Casas-Marce M, Soriano L, Arango JR, Derdak S, Galán B, Blanc J, Gut M, Lorente-Galdos B, Andrés-Nieto M, López-Otín C, Valencia A, Gut I, García JL, Guigó R, Murphy WJ, Ruiz-Herrera A, Marques-Bonet T, Roma G, Notredame C, Mailund T, Albà MM, Gabaldón T, Alioto T, Godoy JA. Extreme genomic erosion after recurrent demographic bottlenecks in the highly endangered Iberian lynx. Genome Biol 2016; 17:251. [PMID: 27964752 PMCID: PMC5155386 DOI: 10.1186/s13059-016-1090-1] [Citation(s) in RCA: 85] [Impact Index Per Article: 10.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/26/2016] [Accepted: 10/25/2016] [Indexed: 11/21/2022] Open
Abstract
BACKGROUND Genomic studies of endangered species provide insights into their evolution and demographic history, reveal patterns of genomic erosion that might limit their viability, and offer tools for their effective conservation. The Iberian lynx (Lynx pardinus) is the most endangered felid and a unique example of a species on the brink of extinction. RESULTS We generate the first annotated draft of the Iberian lynx genome and carry out genome-based analyses of lynx demography, evolution, and population genetics. We identify a series of severe population bottlenecks in the history of the Iberian lynx that predate its known demographic decline during the 20th century and have greatly impacted its genome evolution. We observe drastically reduced rates of weak-to-strong substitutions associated with GC-biased gene conversion and increased rates of fixation of transposable elements. We also find multiple signatures of genetic erosion in the two remnant Iberian lynx populations, including a high frequency of potentially deleterious variants and substitutions, as well as the lowest genome-wide genetic diversity reported so far in any species. CONCLUSIONS The genomic features observed in the Iberian lynx genome may hamper short- and long-term viability through reduced fitness and adaptive potential. The knowledge and resources developed in this study will boost the research on felid evolution and conservation genomics and will benefit the ongoing conservation and management of this emblematic species.
Collapse
Affiliation(s)
- Federico Abascal
- Structural Biology and Biocomputing Programme, Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
| | - André Corvelo
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - Fernando Cruz
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
- Department of Integrative Ecology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), C/ Americo Vespucio, s/n, 41092, Sevilla, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - José L Villanueva-Cañas
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Anna Vlasova
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Marina Marcet-Houben
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Begoña Martínez-Cruz
- Department of Integrative Ecology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), C/ Americo Vespucio, s/n, 41092, Sevilla, Spain
| | - Jade Yu Cheng
- Bioinformatics Research Centre, Aarhus University, C.F. Møllers Allé 8, 8000, Aarhus, Denmark
| | - Pablo Prieto
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Víctor Quesada
- Departamento de Bioquímica y Biología Molecular, Instituto Universitario de Oncología (IUOPA), Universidad de Oviedo, 33006, Oviedo, Spain
| | - Javier Quilez
- Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, PRBB, Doctor Aiguader, 88, 08003, Barcelona, Spain
| | - Gang Li
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine, Texas A&M University, College Station, TX, 77843, USA
| | - Francisca García
- Servei de Cultius Cel.lulars (SCC, SCAC), Universitat Autònoma de Barcelona, Barcelona, Spain
| | - Miriam Rubio-Camarillo
- Structural Biology and Biocomputing Programme, Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
| | - Leonor Frias
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - Paolo Ribeca
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - Salvador Capella-Gutiérrez
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - José M Rodríguez
- Structural Biology and Biocomputing Programme, Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
- National Bioinformatics Institute (INB), Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
| | - Francisco Câmara
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Ernesto Lowy
- Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Luca Cozzuto
- Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Ionas Erb
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Michael L Tress
- Structural Biology and Biocomputing Programme, Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
| | - Jose L Rodriguez-Ales
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Jorge Ruiz-Orera
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Ferran Reverter
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Mireia Casas-Marce
- Department of Integrative Ecology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), C/ Americo Vespucio, s/n, 41092, Sevilla, Spain
| | - Laura Soriano
- Department of Integrative Ecology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), C/ Americo Vespucio, s/n, 41092, Sevilla, Spain
| | - Javier R Arango
- Departamento de Bioquímica y Biología Molecular, Instituto Universitario de Oncología (IUOPA), Universidad de Oviedo, 33006, Oviedo, Spain
| | - Sophia Derdak
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - Beatriz Galán
- Department of Environmental Biology, Center for Biological Research (CIB), Spanish National Research Council (CSIC), Ramiro de Maeztu 9, 28040, Madrid, Spain
| | - Julie Blanc
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - Marta Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - Belen Lorente-Galdos
- Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, PRBB, Doctor Aiguader, 88, 08003, Barcelona, Spain
| | - Marta Andrés-Nieto
- Institut de Biotecnologia i de Biomedicina, Universitat Autònoma de Barcelona, 08193, Cerdanyola del Vallès, Spain
| | - Carlos López-Otín
- Departamento de Bioquímica y Biología Molecular, Instituto Universitario de Oncología (IUOPA), Universidad de Oviedo, 33006, Oviedo, Spain
| | - Alfonso Valencia
- Structural Biology and Biocomputing Programme, Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
- National Bioinformatics Institute (INB), Spanish National Cancer Research Centre (CNIO), Madrid, 28029, Spain
| | - Ivo Gut
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
| | - José L García
- Department of Environmental Biology, Center for Biological Research (CIB), Spanish National Research Council (CSIC), Ramiro de Maeztu 9, 28040, Madrid, Spain
| | - Roderic Guigó
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
- Computational Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - William J Murphy
- Department of Veterinary Integrative Biosciences, College of Veterinary Medicine, Texas A&M University, College Station, TX, 77843, USA
| | - Aurora Ruiz-Herrera
- Institut de Biotecnologia i de Biomedicina, Universitat Autònoma de Barcelona, 08193, Cerdanyola del Vallès, Spain
- Departament de Biologia Cel.lular, Fisiologia i Immunologia, Universitat Autònoma de Barcelona, 08193, Cerdanyola del Vallès, Spain
| | - Tomas Marques-Bonet
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
- Institut de Biologia Evolutiva (UPF-CSIC), Universitat Pompeu Fabra, PRBB, Doctor Aiguader, 88, 08003, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, 08010, Barcelona, Spain
| | - Guglielmo Roma
- Bioinformatics Core Facility, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Cedric Notredame
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - Thomas Mailund
- Bioinformatics Research Centre, Aarhus University, C.F. Møllers Allé 8, 8000, Aarhus, Denmark
| | - M Mar Albà
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, 08010, Barcelona, Spain
| | - Toni Gabaldón
- Bioinformatics and Genomics Programme, Centre for Genomic Regulation (CRG), Dr. Aiguader 88, 08003, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
- Institució Catalana de Recerca i Estudis Avançats (ICREA), Pg. Lluís Companys 23, 08010, Barcelona, Spain
| | - Tyler Alioto
- CNAG-CRG, Centre for Genomic Regulation (CRG), Barcelona Institute of Science and Technology (BIST), Baldiri i Reixac 4, 08028, Barcelona, Spain
- Universitat Pompeu Fabra (UPF), Dr. Aiguader 88, 08003, Barcelona, Spain
| | - José A Godoy
- Department of Integrative Ecology, Doñana Biological Station (EBD), Spanish National Research Council (CSIC), C/ Americo Vespucio, s/n, 41092, Sevilla, Spain.
| |
Collapse
|
31
|
Faherty SL, Villanueva-Cañas JL, Klopfer PH, Albà MM, Yoder AD. Gene Expression Profiling in the Hibernating Primate, Cheirogaleus Medius. Genome Biol Evol 2016; 8:2413-26. [PMID: 27412611 PMCID: PMC5010898 DOI: 10.1093/gbe/evw163] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 07/01/2016] [Indexed: 12/24/2022] Open
Abstract
Hibernation is a complex physiological response that some mammalian species employ to evade energetic demands. Previous work in mammalian hibernators suggests that hibernation is activated not by a set of genes unique to hibernators, but by differential expression of genes that are present in all mammals. This question of universal genetic mechanisms requires further investigation and can only be tested through additional investigations of phylogenetically dispersed species. To explore this question, we use RNA-Seq to investigate gene expression dynamics as they relate to the varying physiological states experienced throughout the year in a group of primate hibernators-Madagascar's dwarf lemurs (genus Cheirogaleus). In a novel experimental approach, we use longitudinal sampling of biological tissues as a method for capturing gene expression profiles from the same individuals throughout their annual hibernation cycle. We identify 90 candidate genes that have variable expression patterns when comparing two active states (Active 1 and Active 2) with a torpor state. These include genes that are involved in metabolic pathways, feeding behavior, and circadian rhythms, as might be expected to correlate with seasonal physiological state changes. The identified genes appear to be critical for maintaining the health of an animal that undergoes prolonged periods of metabolic depression concurrent with the hibernation phenotype. By focusing on these differentially expressed genes in dwarf lemurs, we compare gene expression patterns in previously studied mammalian hibernators. Additionally, by employing evolutionary rate analysis, we find that hibernation-related genes do not evolve under positive selection in hibernating species relative to nonhibernators.
Collapse
Affiliation(s)
| | - José Luis Villanueva-Cañas
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona, Spain
| | | | - M Mar Albà
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona, Spain Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| | | |
Collapse
|
32
|
Duguet TB, Charvet CL, Forrester SG, Wever CM, Dent JA, Neveu C, Beech RN. Recent Duplication and Functional Divergence in Parasitic Nematode Levamisole-Sensitive Acetylcholine Receptors. PLoS Negl Trop Dis 2016; 10:e0004826. [PMID: 27415016 PMCID: PMC4945070 DOI: 10.1371/journal.pntd.0004826] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/18/2016] [Accepted: 06/16/2016] [Indexed: 01/07/2023] Open
Abstract
Helminth parasites rely on fast-synaptic transmission in their neuromusculature to experience the outside world and respond to it. Acetylcholine plays a pivotal role in this and its receptors are targeted by a wide variety of both natural and synthetic compounds used in human health and for the control of parasitic disease. The model, Caenorhabditis elegans is characterized by a large number of acetylcholine receptor subunit genes, a feature shared across the nematodes. This dynamic family is characterized by both gene duplication and loss between species. The pentameric levamisole-sensitive acetylcholine receptor has been characterized from C. elegans, comprised of five different subunits. More recently, cognate receptors have been reconstituted from multiple parasitic nematodes that are found to vary in subunit composition. In order to understand the implications of receptor composition change and the origins of potentially novel drug targets, we investigated a specific example of subunit duplication based on analysis of genome data for 25 species from the 50 helminth genome initiative. We found multiple independent duplications of the unc-29, acetylcholine receptor subunit, where codon substitution rate analysis identified positive, directional selection acting on amino acid positions associated with subunit assembly. Characterization of four gene copies from a model parasitic nematode, Haemonchus contortus, demonstrated that each copy has acquired unique functional characteristics based on phenotype rescue of transgenic C. elegans and electrophysiology of receptors reconstituted in Xenopus oocytes. We found evidence that a specific incompatibility has evolved for two subunits co-expressed in muscle. We demonstrated that functional divergence of acetylcholine receptors, driven by directional selection, can occur more rapidly than previously thought and may be mediated by alteration of receptor assembly. This phenomenon is common among the clade V parasitic nematodes and this work provides a foundation for understanding the broader context of changing anthelmintic drug targets across the parasitic nematodes.
Collapse
Affiliation(s)
- Thomas B. Duguet
- Institute of Parasitology, McGill University, Sainte-Anne-de-Bellevue, Quebec, Canada
| | - Claude L. Charvet
- INRA, UMR1282 Infectiologie Animale et Santé Publique, Nouzilly, France
- Université François Rabelais de Tours, UMR1282, Infectiologie Santé Publique, Tours, France
| | - Sean G. Forrester
- Faculty of Science, University of Ontario Institute of Technology, Oshawa, Ontario, Canada
| | - Claudia M. Wever
- Department of Biology, McGill University, Montreal, Quebec, Canada
| | - Joseph A. Dent
- Department of Biology, McGill University, Montreal, Quebec, Canada
- Centre for Host-Parasite Interactions, Ste-Anne-de-Bellevue, Quebec, Canada
| | - Cedric Neveu
- INRA, UMR1282 Infectiologie Animale et Santé Publique, Nouzilly, France
- Université François Rabelais de Tours, UMR1282, Infectiologie Santé Publique, Tours, France
| | - Robin N. Beech
- Institute of Parasitology, McGill University, Sainte-Anne-de-Bellevue, Quebec, Canada
- Centre for Host-Parasite Interactions, Ste-Anne-de-Bellevue, Quebec, Canada
- * E-mail:
| |
Collapse
|
33
|
Calvo-Martín JM, Librado P, Aguadé M, Papaceit M, Segarra C. Adaptive selection and coevolution at the proteins of the Polycomb repressive complexes in Drosophila. Heredity (Edinb) 2016; 116:213-23. [PMID: 26486609 PMCID: PMC4806890 DOI: 10.1038/hdy.2015.91] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2015] [Revised: 07/23/2015] [Accepted: 08/10/2015] [Indexed: 11/08/2022] Open
Abstract
Polycomb group (PcG) proteins are important epigenetic regulatory proteins that modulate the chromatin state through posttranslational histone modifications. These interacting proteins form multimeric complexes that repress gene expression. Thus, PcG proteins are expected to evolve coordinately, which might be reflected in their phylogenetic trees by concordant episodes of positive selection and by a correlation in evolutionary rates. In order to detect these signals of coevolution, the molecular evolution of 17 genes encoding the subunits of five Polycomb repressive complexes has been analyzed in the Drosophila genus. The observed distribution of divergence differs substantially among and along proteins. Indeed, CAF1 is uniformly conserved, whereas only the established protein domains are conserved in other proteins, such as PHO, PHOL, PSC, PH-P and ASX. Moreover, regions with a low divergence not yet described as protein domains are present, for instance, in SFMBT and SU(Z)12. Maximum likelihood methods indicate an acceleration in the nonsynonymous substitution rate at the lineage ancestral to the obscura group species in most genes encoding subunits of the Pcl-PRC2 complex and in genes Sfmbt, Psc and Kdm2. These methods also allow inferring the action of positive selection in this lineage at genes E(z) and Sfmbt. Finally, the protein interaction network predicted from the complete proteomes of 12 Drosophila species using a coevolutionary approach shows two tight PcG clusters. These clusters include well-established binary interactions among PcG proteins as well as new putative interactions.
Collapse
Affiliation(s)
- J M Calvo-Martín
- Facultat de Biologia, Departament de Genètica, and Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - P Librado
- Facultat de Biologia, Departament de Genètica, and Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - M Aguadé
- Facultat de Biologia, Departament de Genètica, and Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - M Papaceit
- Facultat de Biologia, Departament de Genètica, and Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - C Segarra
- Facultat de Biologia, Departament de Genètica, and Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| |
Collapse
|
34
|
Radó-Trilla N, Arató K, Pegueroles C, Raya A, de la Luna S, Albà MM. Key Role of Amino Acid Repeat Expansions in the Functional Diversification of Duplicated Transcription Factors. Mol Biol Evol 2015; 32:2263-72. [PMID: 25931513 PMCID: PMC4540963 DOI: 10.1093/molbev/msv103] [Citation(s) in RCA: 20] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
The high regulatory complexity of vertebrates has been related to two rounds of whole genome duplication (2R-WGD) that occurred before the divergence of the major vertebrate groups. Following these events, many developmental transcription factors (TFs) were retained in multiple copies and subsequently specialized in diverse functions, whereas others reverted to their singleton state. TFs are known to be generally rich in amino acid repeats or low-complexity regions (LCRs), such as polyalanine or polyglutamine runs, which can evolve rapidly and potentially influence the transcriptional activity of the protein. Here we test the hypothesis that LCRs have played a major role in the diversification of TF gene duplicates. We find that nearly half of the TF gene families originated during the 2R-WGD contains LCRs. The number of gene duplicates with LCRs is 155 out of 550 analyzed (28%), about twice as many as the number of single copy genes with LCRs (15 out of 115, 13%). In addition, duplicated TFs preferentially accumulate certain LCR types, the most prominent of which are alanine repeats. We experimentally test the role of alanine-rich LCRs in two different TF gene families, PHOX2A/PHOX2B and LHX2/LHX9. In both cases, the presence of the alanine-rich LCR in one of the copies (PHOX2B and LHX2) significantly increases the capacity of the TF to activate transcription. Taken together, the results provide strong evidence that LCRs are important driving forces of evolutionary change in duplicated genes.
Collapse
Affiliation(s)
- Núria Radó-Trilla
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Barcelona, Spain
| | - Krisztina Arató
- Department of Experimental and Health Sciences, Universitat Pompeu Fabra (UPF), Barcelona, Spain Centre for Genomic Regulation (CRG), Barcelona, Spain Centro de Investigación Biomèdica en Red en Enfermedades Raras (CIBERER), Barcelona, Spain
| | - Cinta Pegueroles
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Barcelona, Spain Centre for Genomic Regulation (CRG), Barcelona, Spain
| | - Alicia Raya
- Department of Experimental and Health Sciences, Universitat Pompeu Fabra (UPF), Barcelona, Spain Centre for Genomic Regulation (CRG), Barcelona, Spain Centro de Investigación Biomèdica en Red en Enfermedades Raras (CIBERER), Barcelona, Spain
| | - Susana de la Luna
- Department of Experimental and Health Sciences, Universitat Pompeu Fabra (UPF), Barcelona, Spain Centre for Genomic Regulation (CRG), Barcelona, Spain Centro de Investigación Biomèdica en Red en Enfermedades Raras (CIBERER), Barcelona, Spain Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
| | - M Mar Albà
- Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Barcelona, Spain Department of Experimental and Health Sciences, Universitat Pompeu Fabra (UPF), Barcelona, Spain Institució Catalana de Recerca i Estudis Avançats (ICREA), Barcelona, Spain
| |
Collapse
|
35
|
Meng Y, Zhang W, Zhou J, Liu M, Chen J, Tian S, Zhuo M, Zhang Y, Zhong Y, Du H, Wang X. Genome-wide analysis of positively selected genes in seasonal and non-seasonal breeding species. PLoS One 2015; 10:e0126736. [PMID: 26000771 PMCID: PMC4441472 DOI: 10.1371/journal.pone.0126736] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/18/2014] [Accepted: 04/07/2015] [Indexed: 01/04/2023] Open
Abstract
Some mammals breed throughout the year, while others breed only at certain times of year. These differences in reproductive behavior can be explained by evolution. We identified positively-selected genes in two sets of species with different degrees of relatedness including seasonal and non-seasonal breeding species, using branch-site models. After stringent filtering by sum of pairs scoring, we revealed that more genes underwent positive selection in seasonal compared with non-seasonal breeding species. Positively-selected genes were verified by cDNA mapping of the positive sites with the corresponding cDNA sequences. The design of the evolutionary analysis can effectively lower the false-positive rate and thus identify valid positive genes. Validated, positively-selected genes, including CGA, DNAH1, INVS, and CD151, were related to reproductive behaviors such as spermatogenesis and cell proliferation in non-seasonal breeding species. Genes in seasonal breeding species, including THRAP3, TH1L, and CMTM6, may be related to the evolution of sperm and the circadian rhythm system. Identification of these positively-selected genes might help to identify the molecular mechanisms underlying seasonal and non-seasonal reproductive behaviors.
Collapse
Affiliation(s)
- Yuhuan Meng
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Wenlu Zhang
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Jinghui Zhou
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Mingyu Liu
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Junhui Chen
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Shuai Tian
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Min Zhuo
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Yu Zhang
- Guangdong Key Laboratory of Laboratory Animals/Guangdong laboratory animals monitoring institution, Guangzhou, China
| | - Yang Zhong
- School of Life Sciences, Fudan University, Shanghai, China
- Institute of Biodiversity Science, Tibet University, Lhasa, China
| | - Hongli Du
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
| | - Xiaoning Wang
- School of Bioscience and Bioengineering, Guangdong Provincial Key Laboratory of Fermentation and Enzyme Engineering, South China University of Technology, Guangzhou, China
- Chinese PLA General Hospital, Beijing, China
| |
Collapse
|
36
|
De La Torre AR, Lin YC, Van de Peer Y, Ingvarsson PK. Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families. Genome Biol Evol 2015; 7:1002-15. [PMID: 25747252 PMCID: PMC4419791 DOI: 10.1093/gbe/evv044] [Citation(s) in RCA: 50] [Impact Index Per Article: 5.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/13/2023] Open
Abstract
The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.
Collapse
Affiliation(s)
| | - Yao-Cheng Lin
- Department of Plant Systems Biology, VIB, and Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium
| | - Yves Van de Peer
- Department of Plant Systems Biology, VIB, and Department of Plant Biotechnology and Bioinformatics, Ghent University, Ghent, Belgium Genomics Research Institute, University of Pretoria, South Africa
| | - Pär K Ingvarsson
- Department of Ecology and Environmental Science, Umeå University, Sweden Umeå Plant Science Centre, Umeå, Sweden
| |
Collapse
|
37
|
Gayà-Vidal M, Albà MM. Uncovering adaptive evolution in the human lineage. BMC Genomics 2014; 15:599. [PMID: 25030307 PMCID: PMC4124166 DOI: 10.1186/1471-2164-15-599] [Citation(s) in RCA: 31] [Impact Index Per Article: 3.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/03/2014] [Accepted: 07/10/2014] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The recent increase in human polymorphism data, together with the availability of genome sequences from several primate species, provides an unprecedented opportunity to investigate how natural selection has shaped human evolution. RESULTS We compared human branch-specific substitutions with variation data in the current human population to measure the impact of adaptive evolution on human protein coding genes. The use of single nucleotide polymorphisms (SNPs) with high derived allele frequencies (DAFs) minimized the influence of segregating slightly deleterious mutations and improved the estimation of the number of adaptive sites. Using DAF ≥ 60% we showed that the proportion of adaptive substitutions is 0.2% in the complete gene set. However, the percentage rose to 40% when we focused on genes that are specifically accelerated in the human branch with respect to the chimpanzee branch, or on genes that show signatures of adaptive selection at the codon level by the maximum likelihood based branch-site test. In general, neural genes are enriched in positive selection signatures. Genes with multiple lines of evidence of positive selection include taxilin beta, which is involved in motor nerve regeneration and syntabulin, and is required for the formation of new presynaptic boutons. CONCLUSIONS We combined several methods to detect adaptive evolution in human coding sequences at a genome-wide level. The use of variation data, in addition to sequence divergence information, uncovered previously undetected positive selection signatures in neural genes.
Collapse
Affiliation(s)
| | - M Mar Albà
- Evolutionary Genomics Group IMIM-UPF Research Programme on Biomedical Informatics, Barcelona Biomedical Research Park (PRBB), Aiguader 88, 08003 Barcelona, Spain.
| |
Collapse
|
38
|
Villanueva-Cañas JL, Faherty SL, Yoder AD, Albà MM. Comparative genomics of mammalian hibernators using gene networks. Integr Comp Biol 2014; 54:452-62. [PMID: 24881044 DOI: 10.1093/icb/icu048] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
In recent years, the study of the molecular processes involved in mammalian hibernation has shifted from investigating a few carefully selected candidate genes to large-scale analysis of differential gene expression. The availability of high-throughput data provides an unprecedented opportunity to ask whether phylogenetically distant species show similar mechanisms of genetic control, and how these relate to particular genes and pathways involved in the hibernation phenotype. In order to address these questions, we compare 11 datasets of differentially expressed (DE) genes from two ground squirrel species, one bat species, and the American black bear, as well as a list of genes extracted from the literature that previously have been correlated with the drastic physiological changes associated with hibernation. We identify several genes that are DE in different species, indicating either ancestral adaptations or evolutionary convergence. When we use a network approach to expand the original datasets of DE genes to large gene networks using available interactome data, a higher agreement between datasets is achieved. This indicates that the same key pathways are important for activating and maintaining the hibernation phenotype. Functional-term-enrichment analysis identifies several important metabolic and mitochondrial processes that are critical for hibernation, such as fatty acid beta-oxidation and mitochondrial transport. We do not detect any enrichment of positive selection signatures in the coding sequences of genes from the networks of hibernation-associated genes, supporting the hypothesis that the genetic processes shaping the hibernation phenotype are driven primarily by changes in gene regulation.
Collapse
Affiliation(s)
- José Luis Villanueva-Cañas
- *Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona 08003, Spain; Department of Biology, Duke University, Durham, NC 27708, USA; Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| | - Sheena L Faherty
- *Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona 08003, Spain; Department of Biology, Duke University, Durham, NC 27708, USA; Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| | - Anne D Yoder
- *Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona 08003, Spain; Department of Biology, Duke University, Durham, NC 27708, USA; Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| | - M Mar Albà
- *Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona 08003, Spain; Department of Biology, Duke University, Durham, NC 27708, USA; Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain*Evolutionary Genomics Group, Research Programme on Biomedical Informatics (GRIB), Hospital del Mar Research Institute (IMIM), Universitat Pompeu Fabra (UPF), Barcelona 08003, Spain; Department of Biology, Duke University, Durham, NC 27708, USA; Catalan Institution for Research and Advanced Studies (ICREA), Barcelona, Spain
| |
Collapse
|
39
|
Redelings B. Erasing errors due to alignment ambiguity when estimating positive selection. Mol Biol Evol 2014; 31:1979-93. [PMID: 24866534 DOI: 10.1093/molbev/msu174] [Citation(s) in RCA: 39] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Current estimates of diversifying positive selection rely on first having an accurate multiple sequence alignment. Simulation studies have shown that under biologically plausible conditions, relying on a single estimate of the alignment from commonly used alignment software can lead to unacceptably high false-positive rates in detecting diversifying positive selection. We present a novel statistical method that eliminates excess false positives resulting from alignment error by jointly estimating the degree of positive selection and the alignment under an evolutionary model. Our model treats both substitutions and insertions/deletions as sequence changes on a tree and allows site heterogeneity in the substitution process. We conduct inference starting from unaligned sequence data by integrating over all alignments. This approach naturally accounts for ambiguous alignments without requiring ambiguously aligned sites to be identified and removed prior to analysis. We take a Bayesian approach and conduct inference using Markov chain Monte Carlo to integrate over all alignments on a fixed evolutionary tree topology. We introduce a Bayesian version of the branch-site test and assess the evidence for positive selection using Bayes factors. We compare two models of differing dimensionality using a simple alternative to reversible-jump methods. We also describe a more accurate method of estimating the Bayes factor using Rao-Blackwellization. We then show using simulated data that jointly estimating the alignment and the presence of positive selection solves the problem with excessive false positives from erroneous alignments and has nearly the same power to detect positive selection as when the true alignment is known. We also show that samples taken from the posterior alignment distribution using the software BAli-Phy have substantially lower alignment error compared with MUSCLE, MAFFT, PRANK, and FSA alignments.
Collapse
Affiliation(s)
- Benjamin Redelings
- Biology Department, Duke UniversityThe National Evolutionary Synthesis Center, Durham, NC
| |
Collapse
|