1
|
Zeng ZH, Zhong L, Sun HY, Wu ZK, Wang X, Wang H, Li DZ, Barrett SCH, Zhou W. Parallel evolution of morphological and genomic selfing syndromes accompany the breakdown of heterostyly. THE NEW PHYTOLOGIST 2024; 242:302-316. [PMID: 38214455 DOI: 10.1111/nph.19522] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/07/2023] [Accepted: 12/18/2023] [Indexed: 01/13/2024]
Abstract
Evolutionary transitions from outcrossing to selfing in flowering plants have convergent morphological and genomic signatures and can involve parallel evolution within related lineages. Adaptive evolution of morphological traits is often assumed to evolve faster than nonadaptive features of the genomic selfing syndrome. We investigated phenotypic and genomic changes associated with transitions from distyly to homostyly in the Primula oreodoxa complex. We determined whether the transition to selfing occurred more than once and investigated stages in the evolution of morphological and genomic selfing syndromes using 22 floral traits and both nuclear and plastid genomic data from 25 populations. Two independent transitions were detected representing an earlier and a more recently derived selfing lineage. The older lineage exhibited classic features of the morphological and genomic selfing syndrome. Although features of both selfing syndromes were less developed in the younger selfing lineage, they exhibited parallel development with the older selfing lineage. This finding contrasts with the prediction that some genomic changes should lag behind adaptive changes to morphological traits. Our findings highlight the value of comparative studies on the timing and extent of transitions from outcrossing to selfing between related lineages for investigating the tempo of morphological and molecular evolution.
Collapse
Affiliation(s)
- Zhi-Hua Zeng
- Germplasm Bank of Wild Species, Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Li Zhong
- Germplasm Bank of Wild Species, Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
| | - Hua-Ying Sun
- School of Chinese Materia Medica, Yunnan University of Chinese Medicine, Kunming, Yunnan, 650500, China
| | - Zhi-Kun Wu
- Department of Pharmacy, Guizhou University of Traditional Chinese Medicine, Guiyang, Guizhou, 550002, China
| | - Xin Wang
- Germplasm Bank of Wild Species, Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
- University of Chinese Academy of Sciences, Beijing, 100049, China
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
| | - Hong Wang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
| | - De-Zhu Li
- Germplasm Bank of Wild Species, Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
| | - Spencer C H Barrett
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, ON, M5S 3B2, Canada
| | - Wei Zhou
- Germplasm Bank of Wild Species, Yunnan Key Laboratory of Crop Wild Relatives Omics, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, 650201, China
- Lijiang Forest Biodiversity National Observation and Research Station, Kunming Institute of Botany, Chinese Academy of Sciences, Lijiang, Yunnan, 674100, China
| |
Collapse
|
2
|
Fujimori C, Sugimoto K, Ishida M, Yang C, Kayo D, Tomihara S, Sano K, Akazome Y, Oka Y, Kanda S. Long-lasting redundant gnrh1/3 expression in GnRH neurons enabled apparent switching of paralog usage during evolution. iScience 2024; 27:109304. [PMID: 38464591 PMCID: PMC10924128 DOI: 10.1016/j.isci.2024.109304] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Revised: 12/08/2023] [Accepted: 02/16/2024] [Indexed: 03/12/2024] Open
Abstract
Expressed subtype of paralogous genes in functionally homologous cells sometimes show differences across species, the reasons for which have not been explained. The present study examined hypophysiotropic gonadotropin-releasing hormone (GnRH) neurons in vertebrates to investigate this mechanism. These neurons express either gnrh1 or gnrh3 paralogs, depending on the species, and apparent switching of the expressed paralogs in them occurred at least four times in vertebrate evolution. First, we found redundant expression of gnrh1 and gnrh3 in a single neuron in piranha and hypothesized that it may represent an ancestral GnRH system. Moreover, the gnrh1/gnrh3 enhancer of piranha induced reporter RFP/GFP co-expression in a single hypophysiotropic GnRH neuron in both zebrafish and medaka, whose GnRH neurons only express either gnrh3 or gnrh1. Thus, we propose that redundant expression of gnrh1/3 of relatively recent common ancestors may be the key to apparent switching of the paralog usage among present-day species.
Collapse
Affiliation(s)
- Chika Fujimori
- Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba, Japan
- Optics and Imaging Facility, National Institute for Basic Biology, Okazaki, Aichi, Japan
| | - Kohei Sugimoto
- Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba, Japan
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Mio Ishida
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Christopher Yang
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
- Department of Neuroscience, Johns Hopkins University, Baltimore, MD, USA
| | - Daichi Kayo
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
- Graduate School of Life Sciences, Tohoku University, Sendai, Miyagi, Japan
| | - Soma Tomihara
- Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba, Japan
- Nagahama Institute of Bio-Science and Technology, Nagahama, Shiga, Japan
| | - Kaori Sano
- Department of Chemistry, Faculty of Science, Josai University, Sakado, Saitama, Japan
| | - Yasuhisa Akazome
- Department of Anatomy, St. Marianna University School of Medicine, Kawasaki, Kanagawa, Japan
| | - Yoshitaka Oka
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| | - Shinji Kanda
- Atmosphere and Ocean Research Institute, The University of Tokyo, Kashiwa, Chiba, Japan
- Department of Biological Sciences, Graduate School of Science, The University of Tokyo, Bunkyo-ku, Tokyo, Japan
| |
Collapse
|
3
|
de Jong MJ, van Oosterhout C, Hoelzel AR, Janke A. Moderating the neutralist-selectionist debate: exactly which propositions are we debating, and which arguments are valid? Biol Rev Camb Philos Soc 2024; 99:23-55. [PMID: 37621151 DOI: 10.1111/brv.13010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2022] [Revised: 08/04/2023] [Accepted: 08/07/2023] [Indexed: 08/26/2023]
Abstract
Half a century after its foundation, the neutral theory of molecular evolution continues to attract controversy. The debate has been hampered by the coexistence of different interpretations of the core proposition of the neutral theory, the 'neutral mutation-random drift' hypothesis. In this review, we trace the origins of these ambiguities and suggest potential solutions. We highlight the difference between the original, the revised and the nearly neutral hypothesis, and re-emphasise that none of them equates to the null hypothesis of strict neutrality. We distinguish the neutral hypothesis of protein evolution, the main focus of the ongoing debate, from the neutral hypotheses of genomic and functional DNA evolution, which for many species are generally accepted. We advocate a further distinction between a narrow and an extended neutral hypothesis (of which the latter posits that random non-conservative amino acid substitutions can cause non-ecological phenotypic divergence), and we discuss the implications for evolutionary biology beyond the domain of molecular evolution. We furthermore point out that the debate has widened from its initial focus on point mutations, and also concerns the fitness effects of large-scale mutations, which can alter the dosage of genes and regulatory sequences. We evaluate the validity of neutralist and selectionist arguments and find that the tested predictions, apart from being sensitive to violation of underlying assumptions, are often derived from the null hypothesis of strict neutrality, or equally consistent with the opposing selectionist hypothesis, except when assuming molecular panselectionism. Our review aims to facilitate a constructive neutralist-selectionist debate, and thereby to contribute to answering a key question of evolutionary biology: what proportions of amino acid and nucleotide substitutions and polymorphisms are adaptive?
Collapse
Affiliation(s)
- Menno J de Jong
- Senckenberg Biodiversity and Climate Research Institute (SBiK-F), Georg-Voigt-Strasse 14-16, Frankfurt am Main, 60325, Germany
| | - Cock van Oosterhout
- Centre for Ecology, Evolution and Conservation, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ, UK
| | - A Rus Hoelzel
- Department of Biosciences, Durham University, South Road, Durham, DH1 3LE, UK
| | - Axel Janke
- Senckenberg Biodiversity and Climate Research Institute (SBiK-F), Georg-Voigt-Strasse 14-16, Frankfurt am Main, 60325, Germany
- Institute for Ecology, Evolution and Diversity, Goethe University, Max-von-Laue-Strasse 9, Frankfurt am Main, 60438, Germany
- LOEWE-Centre for Translational Biodiversity Genomics (TBG), Senckenberg Nature Research Society, Georg-Voigt-Straße 14-16, Frankfurt am Main, 60325, Germany
| |
Collapse
|
4
|
Xu X, Wang C, Xu C, Yuan J, Wang G, Wu Y, Huang C, Jing H, Yang P, Xu L, Peng S, Shan F, Xia X, Jin F, Hou F, Wang J, Mi D, Ren Y, Liu Y, Irwin DM, Li X, Chen W, Li G. Genomic evolution of island birds from the view of the Swinhoe's pheasant (Lophura swinhoii). Mol Ecol Resour 2024; 24:e13896. [PMID: 37955396 DOI: 10.1111/1755-0998.13896] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/05/2023] [Revised: 10/26/2023] [Accepted: 10/31/2023] [Indexed: 11/14/2023]
Abstract
Island endemic birds account for the majority of extinct vertebrates in the past few centuries. To date, the evolutionary characteristics of island endemic bird's is poorly known. In this research, we de novo assembled a high-quality chromosome-level reference genome for the Swinhoe's pheasant, which is a typical endemic island bird. Results of collinearity tests suggest rapid ancient chromosome rearrangement that may have contributed to the initial species radiation within Phasianidae, and a role for the insertions of CR1 transposable elements in rearranging chromosomes in Phasianidae. During the evolution of the Swinhoe's pheasant, natural selection positively selected genes involved in fecundity and body size functions, at both the species and population levels, which reflect genetic variation associated with island adaptation. We further tested for variation in population genomic traits between the Swinhoe's pheasant and its phylogenetically closely related mainland relative the silver pheasant, and found higher levels of genetic drift and inbreeding in the Swinhoe's pheasant genome. Divergent demographic histories of insular and mainland bird species during the last glacial period may reflect the differing impact of insular and continental climates on the evolution of species. Our research interprets the natural history and population genetic characteristics of the insular endemic bird the Swinhoe's pheasant, at a genome-wide scale, provides a broader perspective on insular speciation, and adaptive evolution and contributes to the genetic conservation of island endemic birds.
Collapse
Affiliation(s)
- Xiao Xu
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Chen Wang
- Guangzhou Zoo, Guangzhou, China
- Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape, Guangzhou, China
| | - Chunzhong Xu
- Shanghai Wild Animal Park Development Co., Ltd, Shanghai, China
| | - Jiaqing Yuan
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Guiqiang Wang
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Yajiang Wu
- Guangzhou Zoo, Guangzhou, China
- Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape, Guangzhou, China
| | - Chen Huang
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Haohao Jing
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Peng Yang
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Lulu Xu
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Shiming Peng
- Guangzhou Zoo, Guangzhou, China
- Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape, Guangzhou, China
| | - Fen Shan
- Guangzhou Zoo, Guangzhou, China
- Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape, Guangzhou, China
| | - Xiaochao Xia
- Guangdong Wildlife Monitoring, Rescue and Conservation Center, Guangzhou, China
| | - Fuyuan Jin
- Guangdong Maoming Forest Park Administrative Office, Maoming, China
| | - Fanghui Hou
- Shanghai Wild Animal Park Development Co., Ltd, Shanghai, China
| | - Jinhong Wang
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Da Mi
- Xi'an Haorui Genomics Technology Co., Ltd, Xi'an, China
- Key Laboratory of Biomedical Information Engineering of Ministry of Education, School of Life Science and Technology, Xi'an Jiaotong University, Xi'an, China
| | - Yandong Ren
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Yang Liu
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - David M Irwin
- Department of Laboratory Medicine and Pathobiology, University of Toronto, Toronto, Ontario, Canada
| | - Xuejuan Li
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
| | - Wu Chen
- Guangzhou Zoo, Guangzhou, China
- Guangzhou Collaborative Innovation Center on Science-Tech of Ecology and Landscape, Guangzhou, China
| | - Gang Li
- College of Life Sciences, Shaanxi Normal University, Xi'an, China
- Guangzhou Zoo, Guangzhou, China
| |
Collapse
|
5
|
Lei H, Li J, Zhao B, Kou SH, Xiao F, Chen T, Wang SM. Evolutionary origin of germline pathogenic variants in human DNA mismatch repair genes. Hum Genomics 2024; 18:5. [PMID: 38287404 PMCID: PMC10823654 DOI: 10.1186/s40246-024-00573-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/27/2023] [Accepted: 01/17/2024] [Indexed: 01/31/2024] Open
Abstract
BACKGROUND Mismatch repair (MMR) system is evolutionarily conserved for genome stability maintenance. Germline pathogenic variants (PVs) in MMR genes that lead to MMR functional deficiency are associated with high cancer risk. Knowing the evolutionary origin of germline PVs in human MMR genes will facilitate understanding the biological base of MMR deficiency in cancer. However, systematic knowledge is lacking to address the issue. In this study, we performed a comprehensive analysis to know the evolutionary origin of human MMR PVs. METHODS We retrieved MMR gene variants from the ClinVar database. The genomes of 100 vertebrates were collected from the UCSC genome browser and ancient human sequencing data were obtained through comprehensive data mining. Cross-species conservation analysis was performed based on the phylogenetic relationship among 100 vertebrates. Rescaled ancient sequencing data were used to perform variant calling for archeological analysis. RESULTS Using the phylogenetic approach, we traced the 3369 MMR PVs identified in modern humans in 99 non-human vertebrate genomes but found no evidence for cross-species conservation as the source for human MMR PVs. Using the archeological approach, we searched the human MMR PVs in over 5000 ancient human genomes dated from 45,045 to 100 years before present and identified a group of MMR PVs shared between modern and ancient humans mostly within 10,000 years with similar quantitative patterns. CONCLUSION Our study reveals that MMR PVs in modern humans were arisen within the recent human evolutionary history.
Collapse
Affiliation(s)
- Huijun Lei
- Ministry of Education Frontiers Science Center for Precision Oncology, Cancer Centre and Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Taipa, Macau SAR, 999078, China
- Hangzhou Institute of Medicine (HIM), Chinese Academy of Sciences, Hangzhou, 310018, Zhejiang, China
- Department of Cancer Prevention, Zhejiang Cancer Hospital, Hangzhou, 310022, Zhejiang, China
| | - Jiaheng Li
- Ministry of Education Frontiers Science Center for Precision Oncology, Cancer Centre and Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Taipa, Macau SAR, 999078, China
| | - Bojin Zhao
- Ministry of Education Frontiers Science Center for Precision Oncology, Cancer Centre and Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Taipa, Macau SAR, 999078, China
| | - Si Hoi Kou
- Ministry of Education Frontiers Science Center for Precision Oncology, Cancer Centre and Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Taipa, Macau SAR, 999078, China
| | - Fengxia Xiao
- Ministry of Education Frontiers Science Center for Precision Oncology, Cancer Centre and Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Taipa, Macau SAR, 999078, China
| | - Tianhui Chen
- Hangzhou Institute of Medicine (HIM), Chinese Academy of Sciences, Hangzhou, 310018, Zhejiang, China.
- Department of Cancer Prevention, Zhejiang Cancer Hospital, Hangzhou, 310022, Zhejiang, China.
| | - San Ming Wang
- Ministry of Education Frontiers Science Center for Precision Oncology, Cancer Centre and Institute of Translational Medicine, Faculty of Health Sciences, University of Macau, Taipa, Macau SAR, 999078, China.
| |
Collapse
|
6
|
Tsunenari K, Ito T, Yokota M, Shibabayashi M, Endo C, Chung KF, Suyama Y, Matsuo A, Abe A, Naiki A, Setoguchi H, Makino T, Isagi Y. Double migration of the endangered Tricyrtis formosana (Liliaceae) in Japan. Sci Rep 2024; 14:957. [PMID: 38200076 PMCID: PMC10781951 DOI: 10.1038/s41598-024-51431-x] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/10/2023] [Accepted: 01/04/2024] [Indexed: 01/12/2024] Open
Abstract
The Ryukyu Islands of Japan are a biodiversity hotspot due to geographical and historical factors. Tricyrtis formosana is a perennial herbaceous plant that commonly found in Taiwan. But only a few populations have been identified in a limited habitat on Iriomote Island, while populations of unknown origin occur near human settlements in an area on the main island of Okinawa. To better understand these populations of the phylogenetic uniqueness and intrinsic vulnerability, we conducted comparative analyses including (1) phylogeny and population structure with MIG-seq data, (2) photosynthesis-related traits of plants grown under common conditions and (3) transcriptome analysis to detect deleterious variations. Results revealed that T. formosana was split into two clades by the congeners and that Iriomote and Okinawa populations independently derived from ancestral Taiwanese populations in each clade. Photosynthetic efficiency was lowest in the Iriomote population, followed by Okinawa and Taiwan. Transcriptome analysis showed that the Iriomote population accumulated more deleterious variations, suggesting intrinsic vulnerability. These results indicate that each T. formosana population in Japan is phylogenetically unique and has been independently dispersed from Taiwan, and that the Iriomote population presents a high conservation difficulty with a unique photosynthesis-related characteristic and a larger amount of deleterious variations.
Collapse
Affiliation(s)
- Kaori Tsunenari
- Graduate School of Agriculture, Kyoto University, Kyoto, Japan.
- Japan Broadcasting Cooporation, Tokyo, Japan.
| | - Takuro Ito
- The Center for Academic Resources and Archives, Tohoku University, Sendai, Japan
| | | | | | - Chiharu Endo
- Graduate School of Agriculture, Kyoto University, Kyoto, Japan
| | - Kuo-Fang Chung
- School of Forestry and Resource Conservation, National Taiwan University, Taipei, Taiwan
- Biodiversity Research Center, Academia Sinica, Taipei, Taiwan
| | - Yoshihisa Suyama
- Graduate School of Agricultural Science, Tohoku University, Sendai, Japan
| | - Ayumi Matsuo
- Graduate School of Agricultural Science, Tohoku University, Sendai, Japan
| | - Atsushi Abe
- Okinawa Churashima Foundation Research Institute, Botanical Laboratory, Okinawa, Japan
| | - Akiyo Naiki
- Tropical Biosphere Research Center, University of the Ryukyus, Taketomi, Okinawa, Japan
| | - Hiroaki Setoguchi
- Graduate School of Human and Environmental Studies, Kyoto University, Kyoto, Japan
| | - Takashi Makino
- Graduate School of Life Science, Tohoku University, Sendai, Japan
| | - Yuji Isagi
- Graduate School of Agriculture, Kyoto University, Kyoto, Japan.
| |
Collapse
|
7
|
Rivas-Santisteban J, Yubero P, Robaina-Estévez S, González JM, Tamames J, Pedrós-Alió C. Quantifying microbial guilds. ISME COMMUNICATIONS 2024; 4:ycae042. [PMID: 38707845 PMCID: PMC11069341 DOI: 10.1093/ismeco/ycae042] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 02/05/2024] [Revised: 03/22/2024] [Accepted: 03/22/2024] [Indexed: 05/07/2024]
Abstract
The ecological role of microorganisms is of utmost importance due to their multiple interactions with the environment. However, assessing the contribution of individual taxonomic groups has proven difficult despite the availability of high throughput data, hindering our understanding of such complex systems. Here, we propose a quantitative definition of guild that is readily applicable to metagenomic data. Our framework focuses on the functional character of protein sequences, as well as their diversifying nature. First, we discriminate functional sequences from the whole sequence space corresponding to a gene annotation to then quantify their contribution to the guild composition across environments. In addition, we identify and distinguish functional implementations, which are sequence spaces that have different ways of carrying out the function. In contrast, we found that orthology delineation did not consistently align with ecologically (or functionally) distinct implementations of the function. We demonstrate the value of our approach with two case studies: the ammonia oxidation and polyamine uptake guilds from the Malaspina circumnavigation cruise, revealing novel ecological dynamics of the latter in marine ecosystems. Thus, the quantification of guilds helps us to assess the functional role of different taxonomic groups with profound implications on the study of microbial communities.
Collapse
Affiliation(s)
- Juan Rivas-Santisteban
- Microbiome Analysis Laboratory, Centro Nacional de Biotecnología (CNB), CSIC, Calle Darwin no. 3, Madrid, 28049, Spain
| | - Pablo Yubero
- Logic of Genomic Systems Laboratory, Centro Nacional de Biotecnología (CNB), CSIC, Spain
| | | | | | - Javier Tamames
- Microbiome Analysis Laboratory, Centro Nacional de Biotecnología (CNB), CSIC, Calle Darwin no. 3, Madrid, 28049, Spain
| | - Carlos Pedrós-Alió
- Microbiome Analysis Laboratory, Centro Nacional de Biotecnología (CNB), CSIC, Calle Darwin no. 3, Madrid, 28049, Spain
| |
Collapse
|
8
|
Li Z, Huang Z, Wan X, Yu J, Dong H, Zhang J, Zhang C, Wang S. Complete chloroplast genome sequence of Rhododendronmariesii and comparative genomics of related species in the family Ericaeae. COMPARATIVE CYTOGENETICS 2023; 17:163-180. [PMID: 37650109 PMCID: PMC10464601 DOI: 10.3897/compcytogen.17.101427] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 02/03/2023] [Accepted: 07/26/2023] [Indexed: 09/01/2023]
Abstract
Rhododendronmariesii Hemsley et Wilson, 1907, a typical member of the family Ericaeae, possesses valuable medicinal and horticultural properties. In this research, the complete chloroplast (cp) genome of R.mariesii was sequenced and assembled, which proved to be a typical quadripartite structure with the length of 203,480 bp. In particular, the lengths of the large single copy region (LSC), small single copy region (SSC), and inverted repeat regions (IR) were 113,715 bp, 7,953 bp, and 40,918 bp, respectively. Among the 151 unique genes, 98 were protein-coding genes, 8 were tRNA genes, and 45 were rRNA genes. The structural characteristics of the R.mariesiicp genome was similar to other angiosperms. Leucine was the most representative amino acid, while cysteine was the lowest representative. Totally, 30 codons showed obvious codon usage bias, and most were A/U-ending codons. Six highly variable regions were observed, such as trnK-pafI and atpE-rpoB, which could serve as potential markers for future barcoding and phylogenetic research of R.mariesii species. Coding regions were more conserved than non-coding regions. Expansion and contraction in the IR region might be the main length variation in R.mariesii and related Ericaeae species. Maximum-likelihood (ML) phylogenetic analysis revealed that R.mariesii was relatively closed to the R.simsii Planchon, 1853 and R.pulchrum Sweet,1831. This research will supply rich genetic resource for R.mariesii and related species of the Ericaeae.
Collapse
Affiliation(s)
- Zhiliang Li
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| | - Zhiwei Huang
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| | - Xuchun Wan
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| | - Jiaojun Yu
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| | - Hongjin Dong
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| | - Jialiang Zhang
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| | - Chunyu Zhang
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
- College of Plant Science & Technology, Huazhong Agricultural University, Wuhan, 430070, Hubei Province, ChinaHuazhong Agricultural UniversityWuhanChina
| | - Shuzhen Wang
- College of Biology and Agricultural Resources, Huanggang Normal University, Huanggang, 438000, Hubei Province, ChinaHuanggang Normal UniversityHuanggangChina
| |
Collapse
|
9
|
Hajihassan Z, Afsharian NP, Ansari-Pour N. In silico engineering a CD80 variant with increased affinity to CTLA-4 and decreased affinity to CD28 for optimized cancer immunotherapy. J Immunol Methods 2023; 513:113425. [PMID: 36638881 DOI: 10.1016/j.jim.2023.113425] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/25/2022] [Revised: 11/20/2022] [Accepted: 01/08/2023] [Indexed: 01/11/2023]
Abstract
CD80 or cluster of differentiation 80, also known as B7-1, is a member of the immunoglobulin super family, which binds to CTLA-4 and CD28 T cell receptors and induces inhibitory and inductive signals respectively. Although CTLA-4 and CD28 receptors belong to the same protein family, slight differences in their structures leads to CD80 having a higher binding affinity to CTLA-4 (-14.55 kcal/mol) compared with CD28(-12.51 kcal/mol). In this study, we constructed a variant of CD80 protein with increased binding affinity to CTLA-4 and decreased binding affinity to CD28. This variant has no signaling capability, and can act as a cap for these receptors to protect them from natural CD80 proteins existing in the body. The first step was the evolutionary and alanine scanning analysis of CD80 protein to determine conserved regions in this protein. Next, complex alanine scanning technique was employed to determine CD80 protein hotspots in CD80-CTLA-4 and CD80-CD28 protein complexes. This information was fed into a computational model developed in R for in silico mutagenesis and CD80 variant library construction. The 3D structures of variants were modeled using the Swiss model webserver. After modeling the 3D structures, HADDOCK server was employed to build all protein-protein complexes, which contain CTLA-4-CD80 variant complexes, Wild type CD80-CD28 complexes and CD28-CD80 variant complexes. Protein-protein binding free energy was determined using FoldX and the variant number 316 with mutations at 29, 31, 33 positions showed increased binding affinity to CTLA-4 (-21.43 kcal/mol) and decreased binding affinity to CD28 (- 9.54 kcal/mol). Finally, molecular dynamics (MD) simulations confirmed the stability of variant 316. In conclusion, we designed a new CD80 protein variant with potential immunotherapeutic applications.
Collapse
Affiliation(s)
- Zahra Hajihassan
- Department of Life Science Engineering, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran.
| | - Nessa Pesaran Afsharian
- Department of Life Science Engineering, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran
| | - Naser Ansari-Pour
- Department of Life Science Engineering, Faculty of New Sciences & Technologies, University of Tehran, Tehran, Iran; MRC Molecular Haematology Unit, Weatherall Institute of Molecular Medicine, University of Oxford, Oxford, UK.
| |
Collapse
|
10
|
Gu X. d N/d S-H, a New Test to Distinguish Different Selection Modes in Protein Evolution and Cancer Evolution. J Mol Evol 2022; 90:342-351. [PMID: 35920867 DOI: 10.1007/s00239-022-10064-2] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/19/2022] [Accepted: 06/14/2022] [Indexed: 11/25/2022]
Abstract
One of the most popular measures in the analysis of protein sequence evolution is the ratio of nonsynonymous distance (dN) to synonymous distance (dS). Under the assumption that synonymous substitutions in the coding region are selectively neutral, the dN/dS ratio can be used to statistically detect the adaptive evolution (or purifying selection) if dN/dS > 1 (or dN/dS < 1) significantly. However, due to strong structural constraints and/or variable functional constraints imposed on amino acid sites, most encoding genes in most species have demonstrated dN/dS < 1. Consequently, the statistical power for testing dN/dS = 1 may be insufficient to distinguish between different selection modes. In this paper, we propose a more powerful test, called dN/dS-H, in which a new parameter H, a relative measure of rate variation among sites, was introduced. Given the condition of strong purifying selections at some sites, the dN/dS-H model predicts dN/dS = 1-H for neutral evolution, dN/dS < 1-H for nearly neutral selection, and dN/dS > 1-H for adaptive evolution. The potential of this new method for resolving the neutral-adaptive debates is illustrated by the protein sequence evolution in vertebrates, Drosophila and yeasts, as well as somatic cancer evolution (specialized as the CN/CS-H test).
Collapse
Affiliation(s)
- Xun Gu
- The Laurence H. Baker Center in Bioinformatics on Biological Statistics, Iowa State University, Ames, IA, 50011, USA. .,Department of Genetics, Development and Cell Biology, Iowa State University, Ames, IA, 50011, USA. .,Program of Ecological and Evolutionary Biology, Iowa State University, Ames, IA, 50011, USA.
| |
Collapse
|
11
|
Melissa MJ, Good BH, Fisher DS, Desai MM. Population genetics of polymorphism and divergence in rapidly evolving populations. Genetics 2022; 221:6564664. [PMID: 35389471 PMCID: PMC9339298 DOI: 10.1093/genetics/iyac053] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2021] [Accepted: 03/19/2022] [Indexed: 11/14/2022] Open
Abstract
In rapidly evolving populations, numerous beneficial and deleterious mutations can arise and segregate within a population at the same time. In this regime, evolutionary dynamics cannot be analyzed using traditional population genetic approaches that assume that sites evolve independently. Instead, the dynamics of many loci must be analyzed simultaneously. Recent work has made progress by first analyzing the fitness variation within a population, and then studying how individual lineages interact with this traveling fitness wave. However, these "traveling wave" models have previously been restricted to extreme cases where selection on individual mutations is either much faster or much slower than the typical coalescent timescale Tc. In this work, we show how the traveling wave framework can be extended to intermediate regimes in which the scaled fitness effects of mutations (Tcs) are neither large nor small compared to one. This enables us to describe the dynamics of populations subject to a wide range of fitness effects, and in particular, in cases where it is not immediately clear which mutations are most important in shaping the dynamics and statistics of genetic diversity. We use this approach to derive new expressions for the fixation probabilities and site frequency spectra of mutations as a function of their scaled fitness effects, along with related results for the coalescent timescale Tc and the rate of adaptation or Muller's ratchet. We find that competition between linked mutations can have a dramatic impact on the proportions of neutral and selected polymorphisms, which is not simply summarized by the scaled selection coefficient Tcs. We conclude by discussing the implications of these results for population genetic inferences.
Collapse
Affiliation(s)
- Matthew J Melissa
- Department of Organismic and Evolutionary Biology, Department of Physics, Quantitative Biology Initiative, and NSF-Simons Center for Mathematical and Statistical Analysis of Biology, Harvard University, Cambridge MA 02138, USA
| | - Benjamin H Good
- Department of Applied Physics and Department of Bioengineering, Stanford University, Stanford CA 94305, USA
| | - Daniel S Fisher
- Department of Applied Physics and Department of Bioengineering, Stanford University, Stanford CA 94305, USA
| | - Michael M Desai
- Department of Organismic and Evolutionary Biology, Department of Physics, Quantitative Biology Initiative, and NSF-Simons Center for Mathematical and Statistical Analysis of Biology, Harvard University, Cambridge MA 02138, USA
| |
Collapse
|
12
|
DeGiorgio M, Szpiech ZA. A spatially aware likelihood test to detect sweeps from haplotype distributions. PLoS Genet 2022; 18:e1010134. [PMID: 35404934 PMCID: PMC9022890 DOI: 10.1371/journal.pgen.1010134] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/18/2021] [Revised: 04/21/2022] [Accepted: 03/04/2022] [Indexed: 01/13/2023] Open
Abstract
The inference of positive selection in genomes is a problem of great interest in evolutionary genomics. By identifying putative regions of the genome that contain adaptive mutations, we are able to learn about the biology of organisms and their evolutionary history. Here we introduce a composite likelihood method that identifies recently completed or ongoing positive selection by searching for extreme distortions in the spatial distribution of the haplotype frequency spectrum along the genome relative to the genome-wide expectation taken as neutrality. Furthermore, the method simultaneously infers two parameters of the sweep: the number of sweeping haplotypes and the “width” of the sweep, which is related to the strength and timing of selection. We demonstrate that this method outperforms the leading haplotype-based selection statistics, though strong signals in low-recombination regions merit extra scrutiny. As a positive control, we apply it to two well-studied human populations from the 1000 Genomes Project and examine haplotype frequency spectrum patterns at the LCT and MHC loci. We also apply it to a data set of brown rats sampled in NYC and identify genes related to olfactory perception. To facilitate use of this method, we have implemented it in user-friendly open source software. Identifying regions of the genome that contain adaptive variation is of fundamental interest in evolutionary biology, providing insight into an organism’s history and biology. When positive selection is recent or ongoing, we expect to find genomic patterns such as high frequency haplotypes and low genetic diversity in the vicinity of the adaptive locus. Here we develop a statistic to identify these regions based on distortions of the haplotype frequency spectrum from a background distribution. We evaluate the performance of this statistic under numerous realistic settings of interest to empiricists and demonstrate its superior performance relative to other haplotype-based selection statistics. We also apply this statistic to real population-genetic data. As a positive control, we explore two well-studied loci, LCT and MHC, in a European and an African human population that show strong evidence for selection. We also apply this statistic to the genomes of an urban brown rat population, where we uncover evidence for adaptation in olfactory perception genes. We release user-friendly software implementing this statistic.
Collapse
Affiliation(s)
- Michael DeGiorgio
- Department of Electrical Engineering and Computer Science, Florida Atlantic University, Boca Raton, Florida, United States of America
- * E-mail: (MD); (ZAS)
| | - Zachary A. Szpiech
- Department of Biology, Pennsylvania State University, University Park, Pennsylvania, United States of America
- Institute for Computational and Data Sciences, Pennsylvania State University, University Park, Pennsylvania, United States of America
- * E-mail: (MD); (ZAS)
| |
Collapse
|
13
|
Ordoñez D, Bohórquez MD, Avendaño C, Patarroyo MA. Comparing Class II MHC DRB3 Diversity in Colombian Simmental and Simbrah Cattle Across Worldwide Bovine Populations. Front Genet 2022; 13:772885. [PMID: 35186024 PMCID: PMC8854852 DOI: 10.3389/fgene.2022.772885] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/08/2021] [Accepted: 01/17/2022] [Indexed: 11/22/2022] Open
Abstract
The major histocompatibility complex (MHC) exerts great influence on responses to infectious diseases and vaccination due to its fundamental role in the adaptive immune system. Knowledge about MHC polymorphism distribution among breeds can provide insights into cattle evolution and diversification as well as population-based immune response variability, thus guiding further studies. Colombian Simmental and Simbrah cattle’s BoLA-DRB3 genetic diversity was compared to that of taurine and zebuine breeds worldwide to estimate functional diversity. High allele richness was observed for Simmental and Simbrah cattle; nevertheless, high homozygosity was associated with individual low sequence variability in both the β1 domain and the peptide binding region (PBR), thereby implying reduced MHC-presented peptide repertoire size. There were strong signals of positive selection acting on BoLA-DRB3 in all populations, some of which were poorly structured and displayed common alleles accounting for their high genetic similarity. PBR sequence correlation analysis suggested that, except for a few populations exhibiting some divergence at PBR, global diversity regarding potential MHC-presented peptide repertoire could be similar for the cattle populations analyzed here, which points to the retention of functional diversity in spite of the selective pressures imposed by breeding.
Collapse
Affiliation(s)
- Diego Ordoñez
- Animal Science Faculty, Universidad de Ciencias Aplicadas y Ambientales (U.D.C.A), Bogotá, Colombia
- PhD Program in Tropical Health and Development, Universidad de Salamanca, Salamanca, Spain
| | - Michel David Bohórquez
- Molecular Biology and Immunology Department, Fundación Instituto de Inmunología de Colombia (FIDIC), Bogotá, Colombia
- MSc Program in Microbiology, Universidad Nacional de Colombia, Bogotá, Colombia
| | - Catalina Avendaño
- Animal Science Faculty, Universidad de Ciencias Aplicadas y Ambientales (U.D.C.A), Bogotá, Colombia
| | - Manuel Alfonso Patarroyo
- Molecular Biology and Immunology Department, Fundación Instituto de Inmunología de Colombia (FIDIC), Bogotá, Colombia
- Health Sciences Division, Main Campus, Universidad Santo Tomás, Bogotá, Colombia
- Microbiology Department, Faculty of Medicine, Universidad Nacional de Colombia, Bogotá, Colombia
- *Correspondence: Manuel Alfonso Patarroyo,
| |
Collapse
|
14
|
Termignoni-Garcia F, Kirchman JJ, Clark J, Edwards SV. Comparative Population Genomics of Cryptic Speciation and Adaptive Divergence in Bicknell's and Gray-Cheeked Thrushes (Aves: Catharus bicknelli and Catharus minimus). Genome Biol Evol 2022; 14:evab255. [PMID: 34999784 PMCID: PMC8743040 DOI: 10.1093/gbe/evab255] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 11/10/2021] [Indexed: 02/07/2023] Open
Abstract
Cryptic speciation may occur when reproductive isolation is recent or the accumulation of morphological differences between sister lineages is slowed by stabilizing selection preventing phenotypic differentiation. In North America, Bicknell's Thrush (Catharus bicknelli) and its sister species, the Gray-cheeked Thrush (Catharus minimus), are parapatrically breeding migratory songbirds, distinguishable in nature only by subtle differences in song and coloration, and were recognized as distinct species only in the 1990s. Previous molecular studies have estimated that the species diverged approximately 120,000-420,000 YBP and found very low levels of introgression despite their similarity and sympatry in the spring (prebreeding) migration. To further clarify the history, genetic divergence, genomic structure, and adaptive processes in C. bicknelli and C. minimus, we sequenced and assembled high-coverage reference genomes of both species and resequenced genomes from population samples of C. bicknelli, C. minimus, and two individuals of the Swainson's Thrush (Catharus ustulatus). The genome of C. bicknelli exhibits markedly higher abundances of transposable elements compared with other Catharus and chicken. Demographic and admixture analyses confirm moderate genome-wide differentiation (Fst ≈ 0.10) and limited gene flow between C. bicknelli and C. minimus, but suggest a more recent divergence than estimates based on mtDNA. We find evidence of rapid evolution of the Z-chromosome and elevated divergence consistent with natural selection on genomic regions near genes involved with neuronal processes in C. bicknelli. These genomes are a useful resource for future investigations of speciation, migration, and adaptation in Catharus thrushes.
Collapse
Affiliation(s)
- Flavia Termignoni-Garcia
- Department of Organismic and Evolutionary Biology and Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts, USA
| | | | - Johnathan Clark
- Department of Organismic and Evolutionary Biology and Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts, USA
| | - Scott V Edwards
- Department of Organismic and Evolutionary Biology and Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts, USA
| |
Collapse
|
15
|
Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids. Proc Natl Acad Sci U S A 2021; 118:2102467118. [PMID: 34429360 PMCID: PMC8536349 DOI: 10.1073/pnas.2102467118] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
Numerous animal lineages have maternally inherited symbionts that are required for host reproduction and growth. Endosymbionts also pose a risk to their hosts because of the mutational decay of their genomes through genetic drift or to selfish mutations that favor symbiont fitness over host fitness. One model for heritable endosymbiosis is the association of aphids with their obligate bacterial symbiont, Buchnera We experimentally established heteroplasmic pea aphid matrilines containing pairs of closely related Buchnera haplotypes and used deep sequencing of diagnostic markers to measure haplotype frequencies in successive host generations. These frequencies were used to estimate the effective population size of Buchnera within hosts (i.e., the transmission bottleneck size) and the extent of within-host selection. The within-host effective population size was in the range of 10 to 20, indicating a strong potential for genetic drift and fixation of deleterious mutations. Remarkably, closely related haplotypes were subject to strong within-host selection, with selection coefficients as high as 0.5 per aphid generation. In one case, the direction of selection depended on the thermal environment and went in the same direction as between-host selection. In another, a new mutant haplotype had a strong within-host advantage under both environments but had no discernible effect on host-level fitness under laboratory conditions. Thus, within-host selection can be strong, resulting in a rapid fixation of mutations with little impact on host-level fitness. Together, these results show that within-host selection can drive evolution of an obligate symbiont, accelerating sequence evolution.
Collapse
|
16
|
The importance of species interactions in eco-evolutionary community dynamics under climate change. Nat Commun 2021; 12:4759. [PMID: 34362916 PMCID: PMC8346578 DOI: 10.1038/s41467-021-24977-x] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/29/2020] [Accepted: 07/14/2021] [Indexed: 02/07/2023] Open
Abstract
Eco-evolutionary dynamics are essential in shaping the biological response of communities to ongoing climate change. Here we develop a spatially explicit eco-evolutionary framework which features more detailed species interactions, integrating evolution and dispersal. We include species interactions within and between trophic levels, and additionally, we incorporate the feature that species' interspecific competition might change due to increasing temperatures and affect the impact of climate change on ecological communities. Our modeling framework captures previously reported ecological responses to climate change, and also reveals two key results. First, interactions between trophic levels as well as temperature-dependent competition within a trophic level mitigate the negative impact of climate change on biodiversity, emphasizing the importance of understanding biotic interactions in shaping climate change impact. Second, our trait-based perspective reveals a strong positive relationship between the within-community variation in preferred temperatures and the capacity to respond to climate change. Temperature-dependent competition consistently results both in higher trait variation and more responsive communities to altered climatic conditions. Our study demonstrates the importance of species interactions in an eco-evolutionary setting, further expanding our knowledge of the interplay between ecological and evolutionary processes.
Collapse
|
17
|
Scobeyeva VA, Artyushin IV, Krinitsina AA, Nikitin PA, Antipin MI, Kuptsov SV, Belenikin MS, Omelchenko DO, Logacheva MD, Konorov EA, Samoilov AE, Speranskaya AS. Gene Loss, Pseudogenization in Plastomes of Genus Allium ( Amaryllidaceae), and Putative Selection for Adaptation to Environmental Conditions. Front Genet 2021; 12:674783. [PMID: 34306019 PMCID: PMC8296844 DOI: 10.3389/fgene.2021.674783] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2021] [Accepted: 06/15/2021] [Indexed: 01/07/2023] Open
Abstract
Amaryllidaceae is a large family with more than 1,600 species, belonging to 75 genera. The largest genus—Allium—is vast, comprising about a thousand species. Allium species (as well as other members of the Amaryllidaceae) are widespread and diversified, they are adapted to a wide range of habitats from shady forests to open habitats like meadows, steppes, and deserts. The genes present in chloroplast genomes (plastomes) play fundamental roles for the photosynthetic plants. Plastome traits could thus be associated with geophysical abiotic characteristics of habitats. Most chloroplast genes are highly conserved and are used as phylogenetic markers for many families of vascular plants. Nevertheless, some studies revealed signatures of positive selection in chloroplast genes of many plant families including Amaryllidaceae. We have sequenced plastomes of the following nine Allium (tribe Allieae of Allioideae) species: A. zebdanense, A. moly, A. victorialis, A. macleanii, A. nutans, A. obliquum, A. schoenoprasum, A. pskemense, A. platyspathum, A. fistulosum, A. semenovii, and Nothoscordum bivalve (tribe Leucocoryneae of Allioideae). We compared our data with previously published plastomes and provided our interpretation of Allium plastome genes’ annotations because we found some noteworthy inconsistencies with annotations previously reported. For Allium species we estimated the integral evolutionary rate, counted SNPs and indels per nucleotide position as well as compared pseudogenization events in species of three main phylogenetic lines of genus Allium to estimate whether they are potentially important for plant physiology or just follow the phylogenetic pattern. During examination of the 38 species of Allium and the 11 of other Amaryllidaceae species we found that rps16, rps2, infA, ccsA genes have lost their functionality multiple times in different species (regularly evolutionary events), while the pseudogenization of other genes was stochastic events. We found that the “normal” or “pseudo” state of rps16, rps2, infA, ccsA genes correlates well with the evolutionary line of genus the species belongs to. The positive selection in various NADH dehydrogenase (ndh) genes as well as in matK, accD, and some others were found. Taking into account known mechanisms of coping with excessive light by cyclic electron transport, we can hypothesize that adaptive evolution in genes, coding subunits of NADH-plastoquinone oxidoreductase could be driven by abiotic factors of alpine habitats, especially by intensive light and UV radiation.
Collapse
Affiliation(s)
- Victoria A Scobeyeva
- Department of Evolution, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia.,Department of Molecular and Biological Physics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Ilya V Artyushin
- Department of Vertebrate Zoology, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Anastasiya A Krinitsina
- Department of Higher Plants, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Pavel A Nikitin
- Faculty of Bioengineering and Bioinformatics, Lomonosov Moscow State University, Moscow, Russia
| | - Maxim I Antipin
- Botanical Garden, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Sergei V Kuptsov
- Botanical Garden, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia
| | - Maxim S Belenikin
- Department of Molecular and Biological Physics, Moscow Institute of Physics and Technology, Dolgoprudny, Russia
| | - Denis O Omelchenko
- Laboratory of Plant Genomics, Institute for Information Transmission Problems, Moscow, Russia
| | - Maria D Logacheva
- Center of Life Sciences, Skolkovo Institute of Science and Technology, Moscow, Russia
| | - Evgenii A Konorov
- Laboratory of Animal Genetics, Vavilov Institute of General Genetics, Russian Academy of Science (RAS), Moscow, Russia
| | - Andrey E Samoilov
- Group of Genomics and Postgenomic Technologies, Central Research Institute of Epidemiology, Moscow, Russia
| | - Anna S Speranskaya
- Department of Higher Plants, Faculty of Biology, Lomonosov Moscow State University, Moscow, Russia.,Group of Genomics and Postgenomic Technologies, Central Research Institute of Epidemiology, Moscow, Russia
| |
Collapse
|
18
|
Wang XJ, Barrett SCH, Zhong L, Wu ZK, Li DZ, Wang H, Zhou W. The Genomic Selfing Syndrome Accompanies the Evolutionary Breakdown of Heterostyly. Mol Biol Evol 2021; 38:168-180. [PMID: 32761213 PMCID: PMC7782863 DOI: 10.1093/molbev/msaa199] [Citation(s) in RCA: 17] [Impact Index Per Article: 5.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/28/2022] Open
Abstract
The evolutionary transition from outcrossing to selfing can have important genomic consequences. Decreased effective population size and the reduced efficacy of selection are predicted to play an important role in the molecular evolution of the genomes of selfing species. We investigated evidence for molecular signatures of the genomic selfing syndrome using 66 species of Primula including distylous (outcrossing) and derived homostylous (selfing) taxa. We complemented our comparative analysis with a microevolutionary study of P. chungensis, which is polymorphic for mating system and consists of both distylous and homostylous populations. We generated chloroplast and nuclear genomic data sets for distylous, homostylous, and distylous–homostylous species and identified patterns of nonsynonymous to synonymous divergence (dN/dS) and polymorphism (πN/πS) in species or lineages with contrasting mating systems. Our analysis of coding sequence divergence and polymorphism detected strongly reduced genetic diversity and heterozygosity, decreased efficacy of purifying selection, purging of large-effect deleterious mutations, and lower rates of adaptive evolution in samples from homostylous compared with distylous populations, consistent with theoretical expectations of the genomic selfing syndrome. Our results demonstrate that self-fertilization is a major driver of molecular evolutionary processes with genomic signatures of selfing evident in both old and relatively young homostylous populations.
Collapse
Affiliation(s)
- Xin-Jia Wang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Spencer C H Barrett
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Canada
| | - Li Zhong
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,University of Chinese Academy of Sciences, Beijing, China
| | - Zhi-Kun Wu
- Department of Pharmacy, Guizhou University of Traditional Chinese Medicine, Guiyang, Guizhou, China
| | - De-Zhu Li
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Hong Wang
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| | - Wei Zhou
- CAS Key Laboratory for Plant Diversity and Biogeography of East Asia, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China.,Plant Germplasm and Genomics Center, Germplasm Bank of Wild Species, Kunming Institute of Botany, Chinese Academy of Sciences, Kunming, Yunnan, China
| |
Collapse
|
19
|
Wang P, Burley JT, Liu Y, Chang J, Chen D, Lu Q, Li SH, Zhou X, Edwards S, Zhang Z. Genomic Consequences of Long-Term Population Decline in Brown Eared Pheasant. Mol Biol Evol 2021; 38:263-273. [PMID: 32853368 PMCID: PMC7783171 DOI: 10.1093/molbev/msaa213] [Citation(s) in RCA: 25] [Impact Index Per Article: 8.3] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/20/2022] Open
Abstract
Population genetic theory and empirical evidence indicate that deleterious alleles can be purged in small populations. However, this viewpoint remains controversial. It is unclear whether natural selection is powerful enough to purge deleterious mutations when wild populations continue to decline. Pheasants are terrestrial birds facing a long-term risk of extinction as a result of anthropogenic perturbations and exploitation. Nevertheless, there are scant genomics resources available for conservation management and planning. Here, we analyzed comparative population genomic data for the three extant isolated populations of Brown eared pheasant (Crossoptilon mantchuricum) in China. We showed that C. mantchuricum has low genome-wide diversity and a contracting effective population size because of persistent declines over the past 100,000 years. We compared genome-wide variation in C. mantchuricum with that of its closely related sister species, the Blue eared pheasant (C. auritum) for which the conservation concern is low. There were detrimental genetic consequences across all C. mantchuricum genomes including extended runs of homozygous sequences, slow rates of linkage disequilibrium decay, excessive loss-of-function mutations, and loss of adaptive genetic diversity at the major histocompatibility complex region. To the best of our knowledge, this study is the first to perform a comprehensive conservation genomic analysis on this threatened pheasant species. Moreover, we demonstrated that natural selection may not suffice to purge deleterious mutations in wild populations undergoing long-term decline. The findings of this study could facilitate conservation planning for threatened species and help recover their population size.
Collapse
Affiliation(s)
- Pengcheng Wang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing, China.,Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China.,Department of Organismic and Evolutionary Biology/Museum of Comparative Zoology, Harvard University, Cambridge, MA
| | - John T Burley
- Department of Ecology and Evolutionary Biology/Institute at Brown for Environment and Society, Brown University, Providence, RI
| | - Yang Liu
- State Key Laboratory of Biocontrol, School of Ecology, Sun Yat-Sen University, Guangzhou, China
| | - Jiang Chang
- State Key Laboratory of Environmental Criteria and Risk Assessment, Chinese Research Academy of Environmental Sciences, Beijing, China
| | - De Chen
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing, China
| | - Qi Lu
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing, China
| | - Shou-Hsien Li
- School of Life Science, National Taiwan Normal University, Taipei, Taiwan, China
| | - Xuming Zhou
- Key Laboratory of Animal Ecology and Conservation Biology, Institute of Zoology, Chinese Academy of Sciences, Beijing, China
| | - Scott Edwards
- Department of Organismic and Evolutionary Biology/Museum of Comparative Zoology, Harvard University, Cambridge, MA
| | - Zhengwang Zhang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing, China
| |
Collapse
|
20
|
A comparative analysis of SLA-DRB1 genetic diversity in Colombian (creoles and commercial line) and worldwide swine populations. Sci Rep 2021; 11:4340. [PMID: 33619347 PMCID: PMC7900169 DOI: 10.1038/s41598-021-83637-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/01/2020] [Accepted: 01/18/2021] [Indexed: 12/30/2022] Open
Abstract
Analysing pig class II mayor histocompatibility complex (MHC) molecules is mainly related to antigen presentation. Identifying frequently-occurring alleles in pig populations is an important aspect to be considered when developing peptide-based vaccines. Colombian creole pig populations have had to adapt to local conditions since entering Colombia; a recent census has shown low amounts of pigs which is why they are considered protected by the Colombian government. Commercial hybrids are more attractive regarding production. This research has been aimed at describing the allele distribution of Colombian pigs from diverse genetic backgrounds and comparing Colombian SLA-DRB1 locus diversity to that of internationally reported populations. Twenty SLA-DRB1 alleles were identified in the six populations analysed here using sequence-based typing. The amount of alleles ranged from six (Manta and Casco Mula) to nine (San Pedreño). Only one allele (01:02) having > 5% frequency was shared by all three commercial line populations. Allele 02:01:01 was shared by five populations (around > 5% frequency). Global FST indicated that pig populations were clearly structured, as 20.6% of total allele frequency variation was explained by differences between populations (FST = 0.206). This study’s results confirmed that the greatest diversity occurred in wild boars, thereby contrasting with low diversity in domestic pig populations.
Collapse
|
21
|
Gladieux P, De Bellis F, Hann-Soden C, Svedberg J, Johannesson H, Taylor JW. Neurospora from Natural Populations: Population Genomics Insights into the Life History of a Model Microbial Eukaryote. Methods Mol Biol 2021; 2090:313-336. [PMID: 31975173 DOI: 10.1007/978-1-0716-0199-0_13] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 02/09/2023]
Abstract
The ascomycete filamentous fungus Neurospora crassa played a historic role in experimental biology and became a model system for genetic research. Stimulated by a systematic effort to collect wild strains initiated by Stanford geneticist David Perkins, the genus Neurospora has also become a basic model for the study of evolutionary processes, speciation, and population biology. In this chapter, we will first trace the history that brought Neurospora into the era of population genomics. We will then cover the major contributions of population genomic investigations using Neurospora to our understanding of microbial biogeography and speciation, and review recent work using population genomics and genome-wide association mapping that illustrates the unique potential of Neurospora as a model for identifying the genetic basis of (potentially adaptive) phenotypes in filamentous fungi. The advent of population genomics has contributed to firmly establish Neurospora as a complete model system and we hope our review will entice biologists to include Neurospora in their research.
Collapse
Affiliation(s)
- Pierre Gladieux
- UMR BGPI, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, Montpellier, France.
| | - Fabien De Bellis
- UMR AGAP, Univ Montpellier, CIRAD, INRA, Montpellier SupAgro, Montpellier, France
| | - Christopher Hann-Soden
- Department of Plant and Microbial Biology, University of California, Berkeley, Berkeley, CA, USA
| | - Jesper Svedberg
- Department of Organismal Biology, Uppsala University, Uppsala, Sweden
| | - Hanna Johannesson
- Department of Organismal Biology, Uppsala University, Uppsala, Sweden
| | - John W Taylor
- Department of Plant and Microbial Biology, University of California, Berkeley, Berkeley, CA, USA
| |
Collapse
|
22
|
Leroy T, Rougemont Q. Introduction to Population Genomics Methods. Methods Mol Biol 2021; 2222:287-324. [PMID: 33301100 DOI: 10.1007/978-1-0716-0997-2_16] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/12/2022]
Abstract
High-throughput sequencing technologies have provided an unprecedented opportunity to study the different evolutionary forces that have shaped present-day patterns of genetic diversity, with important implications for many directions in plant biology research. To manage such massive quantities of sequencing data, biologists, however, need new additional skills in informatics and statistics. In this chapter, our objective is to introduce population genomics methods to beginners following a learning-by-doing strategy in order to help the reader to analyze the sequencing data by themselves. Conducted analyses cover several main areas of evolutionary biology, such as an initial description of the evolutionary history of a given species or the identification of genes targeted by natural or artificial selection. In addition to the practical advices, we performed re-analyses of two cases studies with different kind of data: a domesticated cereal (African rice) and a non-domesticated tree species (sessile oak). All the code needed to replicate this work is publicly available on github ( https://github.com/ThibaultLeroyFr/Intro2PopGenomics/ ).
Collapse
Affiliation(s)
- Thibault Leroy
- Montpellier Institute of Evolutionary Sciences (ISEM), Université de Montpellier, Montpellier, France. .,Department of Botany and Biodiversity Research, University of Vienna, Vienna, Austria.
| | - Quentin Rougemont
- Département de Biologie, Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Quebec, QC, Canada
| |
Collapse
|
23
|
Kutschera VE, Poelstra JW, Botero-Castro F, Dussex N, Gemmell NJ, Hunt GR, Ritchie MG, Rutz C, Wiberg RAW, Wolf JBW. Purifying Selection in Corvids Is Less Efficient on Islands. Mol Biol Evol 2020; 37:469-474. [PMID: 31633794 PMCID: PMC6993847 DOI: 10.1093/molbev/msz233] [Citation(s) in RCA: 14] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
Abstract
Theory predicts that deleterious mutations accumulate more readily in small populations. As a consequence, mutation load is expected to be elevated in species where life-history strategies and geographic or historical contingencies reduce the number of reproducing individuals. Yet, few studies have empirically tested this prediction using genome-wide data in a comparative framework. We collected whole-genome sequencing data for 147 individuals across seven crow species (Corvus spp.). For each species, we estimated the distribution of fitness effects of deleterious mutations and compared it with proxies of the effective population size Ne. Island species with comparatively smaller geographic range sizes had a significantly increased mutation load. These results support the view that small populations have an elevated risk of mutational meltdown, which may contribute to the higher extinction rates observed in island species.
Collapse
Affiliation(s)
- Verena E Kutschera
- Department of Evolutionary Biology, Uppsala University, Uppsala, Sweden.,Science for Life Laboratory, Department of Biochemistry and Biophysics, Stockholm University, Stockholm, Sweden
| | | | - Fidel Botero-Castro
- Division of Evolutionary Biology, Faculty of Biology, LMU Munich, Planegg-Martinsried, Germany
| | - Nicolas Dussex
- Department of Anatomy, University of Otago, Dunedin, New Zealand.,Department of Bioinformatics and Genetics, Swedish Museum of Natural History, Stockholm, Sweden
| | - Neil J Gemmell
- Department of Anatomy, University of Otago, Dunedin, New Zealand
| | | | - Michael G Ritchie
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, United Kingdom
| | - Christian Rutz
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, United Kingdom
| | - R Axel W Wiberg
- Centre for Biological Diversity, School of Biology, University of St Andrews, St Andrews, United Kingdom.,Department of Environmental Sciences, Evolutionary Biology, University of Basel, Basel, Switzerland
| | - Jochen B W Wolf
- Department of Evolutionary Biology, Uppsala University, Uppsala, Sweden.,Division of Evolutionary Biology, Faculty of Biology, LMU Munich, Planegg-Martinsried, Germany
| |
Collapse
|
24
|
Kurahashi R, Tanaka SI, Takano K. Highly active enzymes produced by directed evolution with stability-based selection. Enzyme Microb Technol 2020; 140:109626. [DOI: 10.1016/j.enzmictec.2020.109626] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.3] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/09/2020] [Revised: 06/12/2020] [Accepted: 06/12/2020] [Indexed: 12/22/2022]
|
25
|
Dayan A, Yeheskel A, Lamed R, Fleminger G, Ashur-Fabian O. Dihydrolipoamide dehydrogenase moonlighting activity as a DNA chelating agent. Proteins 2020; 89:21-28. [PMID: 32761961 DOI: 10.1002/prot.25991] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2020] [Revised: 06/14/2020] [Accepted: 07/26/2020] [Indexed: 12/14/2022]
Abstract
Dihydrolipoamide dehydrogenase (DLDH) is a mitochondrial enzyme that comprises an essential component of the pyruvate dehydrogenase complex. Lines of evidence have shown that many dehydrogenases possess unrelated actions known as moonlightings in addition to their oxidoreductase activity. As part of these activities, we have demonstrated that DLDH binds TiO2 as well as produces reactive oxygen species (ROS). This ROS production capability was harnessed for cancer therapy via integrin-mediated drug-delivery of RGD-modified DLDH (DLDHRGD ), leading to apoptotic cell death. In these experiments, DLDHRGD not only accumulated in the cytosol but also migrated to the cell nuclei, suggesting a potential DNA-binding capability of this enzyme. To explore this interaction under cell-free conditions, we have analyzed DLDH binding to phage lambda (λ) DNA by gel-shift assays and analytic ultracentrifugation, showing complex formation between the two, which led to full coverage of the DNA molecule with DLDH molecules. DNA binding did not affect DLDH enzymatic activity, indicating that there are neither conformational changes nor active site hindering in DLDH upon DNA-binding. A Docking algorithm for prediction of protein-DNA complexes, Paradoc, identified a putative DNA binding site at the C-terminus of DLDH. Our finding that TiO2 -bound DLDH failed to form a complex with DNA suggests partial overlapping between the two sites. To conclude, DLDH binding to DNA presents a novel moonlight activity which may be used for DNA alkylating in cancer treatment.
Collapse
Affiliation(s)
- Avraham Dayan
- The Shmunis School of Biomedicine and Cancer Research, Tel Aviv University, Ramat Aviv, Israel
| | - Adva Yeheskel
- Bioinformatics Unit, George S. Wise Faculty of Life Sciences, Tel Aviv University, Ramat Aviv, Israel
| | - Raphael Lamed
- The Shmunis School of Biomedicine and Cancer Research, Tel Aviv University, Ramat Aviv, Israel
| | - Gideon Fleminger
- The Shmunis School of Biomedicine and Cancer Research, Tel Aviv University, Ramat Aviv, Israel
| | - Osnat Ashur-Fabian
- The Department of Human Molecular Genetics and Biochemistry, Sackler School of Medicine, Tel Aviv University, Ramat Aviv, Israel
- Translational Oncology Laboratory, Meir Medical Center, Kfar-Saba, Israel
| |
Collapse
|
26
|
Chen Z, Ai F, Zhang J, Ma X, Yang W, Wang W, Su Y, Wang M, Yang Y, Mao K, Wang Q, Lascoux M, Liu J, Ma T. Survival in the Tropics despite isolation, inbreeding and asexual reproduction: insights from the genome of the world's southernmost poplar (Populus ilicifolia). THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2020; 103:430-442. [PMID: 32168389 DOI: 10.1111/tpj.14744] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 07/15/2019] [Revised: 02/04/2020] [Accepted: 03/02/2020] [Indexed: 05/16/2023]
Abstract
Species are becoming extinct at unprecedented rates as a consequence of human activity. Hence it is important to understand the evolutionary dynamics of species with already small population sizes. Populus ilicifolia is a vulnerable poplar species that is isolated from other poplar species and is uniquely adapted to the Tropics. It has a very limited size, reproduces partly clonally and is therefore an excellent case study for conservation genomics. We present here the first annotated draft genome of P. ilicifolia, characterize genome-wide patterns of polymorphisms and compare those to other poplar species with larger natural ranges. P. ilicifolia experienced a more prolonged and severe decline of effective population size (Ne ) and signs of genetic erosion than any other poplar species with which it was compared. At present, the species has the lowest genome-wide genetic diversity, the highest abundance of long runs of homozygosity, high inbreeding levels as well as a high overall accumulation of deleterious variants. However, more effective purging of severely deleterious variants and adaptation to the Tropics may have contributed to its survival. Hence, in spite of its limited genetic variation, it is certainly worth pursuing the conservation efforts of this unique species.
Collapse
Affiliation(s)
- Zeyuan Chen
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Fandi Ai
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Junlin Zhang
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Xinzhi Ma
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Wenlu Yang
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Weiwei Wang
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Yutao Su
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Mingcheng Wang
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Yongzhi Yang
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - Kangshan Mao
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| | - Qingfeng Wang
- Key Laboratory of Aquatic Botany and Watershed Ecology, The Chinese Academy of Sciences, Wuhan, 430074, Hubei, China
| | - Martin Lascoux
- Department of Ecology and Genetics, Evolutionary Biology Centre, Uppsala University, Norbyvägen, 18D 75326, Uppsala, Sweden
| | - Jianquan Liu
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
- State Key Laboratory of Grassland Agro-Ecosystem, Institute of Innovation Ecology & College of Life Sciences, Lanzhou University, Lanzhou, 730000, China
| | - Tao Ma
- Key Laboratory of Bio-Resource and Eco-Environment of Ministry of Education, College of Life Sciences, State Key Laboratory of Hydraulics and Mountain River Engineering, Sichuan University, Chengdu, 610065, China
| |
Collapse
|
27
|
Abstract
Biological systems are dynamic and display heterogeneity at all levels. Ubiquitous heterogeneity, here called for poikilosis, is an integral and important property of organisms and in molecules, systems and processes within them. Traditionally, heterogeneity in biology and experiments has been considered as unwanted noise, here poikilosis is shown to be the normal state. Acceptable variation ranges are called as lagom. Non-lagom, variations that are too extensive, have negative effects, which influence interconnected levels and once the variation is large enough cause a disease and can lead even to death. Poikilosis has numerous applications and consequences e.g. for how to design, analyze and report experiments, how to develop and apply prediction and modelling methods, and in diagnosis and treatment of diseases. Poikilosis-aware new and practical definitions are provided for life, death, senescence, disease, and lagom. Poikilosis is the first new unifying theory in biology since evolution and should be considered in every scientific study.
Collapse
Affiliation(s)
- Mauno Vihinen
- Department of Experimental Medical Science, Lund University, Lund, 22184, Sweden
| |
Collapse
|
28
|
Abstract
Biological systems are dynamic and display heterogeneity at all levels. Ubiquitous heterogeneity, here called for poikilosis, is an integral and important property of organisms and in molecules, systems and processes within them. Traditionally, heterogeneity in biology and experiments has been considered as unwanted noise, here poikilosis is shown to be the normal state. Acceptable variation ranges are called as lagom. Non-lagom, variations that are too extensive, have negative effects, which influence interconnected levels and once the variation is large enough cause a disease and can lead even to death. Poikilosis has numerous applications and consequences e.g. for how to design, analyze and report experiments, how to develop and apply prediction and modelling methods, and in diagnosis and treatment of diseases. Poikilosis-aware new and practical definitions are provided for life, death, senescence, disease, and lagom. Poikilosis is the first new unifying theory in biology since evolution and should be considered in every scientific study.
Collapse
Affiliation(s)
- Mauno Vihinen
- Department of Experimental Medical Science, Lund University, Lund, 22184, Sweden
| |
Collapse
|
29
|
Abstract
Neutral models of evolution assume the absence of natural selection. Formerly confined to ecology and evolutionary biology, neutral models are spreading. In recent years they've been applied to explaining the diversity of baby names, scientific citations, cryptocurrencies, pot decorations, literary lexica, tumour variants and much more besides. Here, we survey important neutral models and highlight their similarities. We investigate the most widely used tests of neutrality, show that they are weak and suggest more powerful methods. We conclude by discussing the role of neutral models in the explanation of diversity. We suggest that the ability of neutral models to fit low-information distributions should not be taken as evidence for the absence of selection. Nevertheless, many studies, in increasingly diverse fields, make just such claims. We call this tendency 'neutral syndrome'.
Collapse
|
30
|
Popitsch N, Huber CD, Buchumenski I, Eisenberg E, Jantsch M, von Haeseler A, Gallach M. A-to-I RNA Editing Uncovers Hidden Signals of Adaptive Genome Evolution in Animals. Genome Biol Evol 2020; 12:345-357. [PMID: 32145015 PMCID: PMC7186786 DOI: 10.1093/gbe/evaa046] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/03/2020] [Indexed: 02/06/2023] Open
Abstract
In animals, the most common type of RNA editing is the deamination of adenosines (A) into inosines (I). Because inosines basepair with cytosines (C), they are interpreted as guanosines (G) by the cellular machinery and genomically encoded G alleles at edited sites mimic the function of edited RNAs. The contribution of this hardwiring effect on genome evolution remains obscure. We looked for population genomics signatures of adaptive evolution associated with A-to-I RNA edited sites in humans and Drosophila melanogaster. We found that single nucleotide polymorphisms at edited sites occur 3 (humans) to 15 times (Drosophila) more often than at unedited sites, the nucleotide G is virtually the unique alternative allele at edited sites and G alleles segregate at higher frequency at edited sites than at unedited sites. Our study reveals that a significant fraction of coding synonymous and nonsynonymous as well as silent and intergenic A-to-I RNA editing sites are likely adaptive in the distantly related human and Drosophila lineages.
Collapse
Affiliation(s)
- Niko Popitsch
- Oxford NIHR Biomedical Research Center, Wellcome Trust Center for Human Genetics, University of Oxford, Oxford, United Kingdom
- Institute of Molecular Biotechnology (IMBA), Vienna BioCenter (VBC), Vienna, Austria
| | - Christian D Huber
- Australian Centre for Ancient DNA, The University of Adelaide, Adelaide, South Australia, Australia
| | - Ilana Buchumenski
- The Mina and Everard Goodman Faculty of Life Sciences, Bar-Ilan University, Ramat Gan, Israel
| | - Eli Eisenberg
- Raymond and Beverly Sackler School of Physics and Astronomy and Sagol School of Neuroscience, Tel Aviv University, Tel Aviv, Israel
| | - Michael Jantsch
- Department for Cell and Developmental Biology, Center for Anatomy and Cell Biology, Medical University of Vienna, Vienna, Austria
- Department for Medical Biochemistry, Max F. Perutz Laboratories, Medical University of Vienna, Vienna, Austria
| | - Arndt von Haeseler
- Bioinformatics and Computational Biology, Faculty of Computer Science, University of Vienna, Vienna, Austria
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna and Medical University of Vienna, Vienna, Austria
| | - Miguel Gallach
- Center for Integrative Bioinformatics Vienna, Max F. Perutz Laboratories, University of Vienna and Medical University of Vienna, Vienna, Austria
- iLabSystems, C/Alicante, Castellón, Spain
| |
Collapse
|
31
|
Allen B, Sample C, Jencks R, Withers J, Steinhagen P, Brizuela L, Kolodny J, Parke D, Lippner G, Dementieva YA. Transient amplifiers of selection and reducers of fixation for death-Birth updating on graphs. PLoS Comput Biol 2020; 16:e1007529. [PMID: 31951612 PMCID: PMC6968840 DOI: 10.1371/journal.pcbi.1007529] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2019] [Accepted: 10/30/2019] [Indexed: 11/30/2022] Open
Abstract
The spatial structure of an evolving population affects the balance of natural selection versus genetic drift. Some structures amplify selection, increasing the role that fitness differences play in determining which mutations become fixed. Other structures suppress selection, reducing the effect of fitness differences and increasing the role of random chance. This phenomenon can be modeled by representing spatial structure as a graph, with individuals occupying vertices. Births and deaths occur stochastically, according to a specified update rule. We study death-Birth updating: An individual is chosen to die and then its neighbors compete to reproduce into the vacant spot. Previous numerical experiments suggested that amplifiers of selection for this process are either rare or nonexistent. We introduce a perturbative method for this problem for weak selection regime, meaning that mutations have small fitness effects. We show that fixation probability under weak selection can be calculated in terms of the coalescence times of random walks. This result leads naturally to a new definition of effective population size. Using this and other methods, we uncover the first known examples of transient amplifiers of selection (graphs that amplify selection for a particular range of fitness values) for the death-Birth process. We also exhibit new families of “reducers of fixation”, which decrease the fixation probability of all mutations, whether beneficial or deleterious. Natural selection is often thought of as “survival of the fittest”, but random chance plays a significant role in which mutations persist and which are eliminated. The balance of selection versus randomness is affected by spatial structure—how individuals are arranged within their habitat. Some structures amplify the effects of selection, so that only the fittest mutations are likely to persist. Others suppress the effects of selection, making the survival of genes primarily a matter of random chance. We study this question using a mathematical model called the “death-Birth process”. Previous studies have found that spatial structure rarely, if ever, amplifies selection for this process. Here we report that spatial structure can indeed amplify selection, at least for mutations with small fitness effects. We also identify structures that reduce the spread of any new mutation, whether beneficial or deleterious. Our work introduces new mathematical techniques for assessing how population structure affects natural selection.
Collapse
Affiliation(s)
- Benjamin Allen
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
- * E-mail:
| | - Christine Sample
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - Robert Jencks
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - James Withers
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - Patricia Steinhagen
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - Lori Brizuela
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - Joshua Kolodny
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - Darren Parke
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| | - Gabor Lippner
- Department of Mathematics, Northeastern University, Boston, Massachusetts, United States of America
| | - Yulia A. Dementieva
- Department of Mathematics, Emmanuel College, Boston, Massachusetts, United States of America
| |
Collapse
|
32
|
Analysis of complete chloroplast genomes of Curcuma and the contribution to phylogeny and adaptive evolution. Gene 2020; 732:144355. [PMID: 31935501 DOI: 10.1016/j.gene.2020.144355] [Citation(s) in RCA: 23] [Impact Index Per Article: 5.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/11/2019] [Revised: 01/07/2020] [Accepted: 01/08/2020] [Indexed: 01/13/2023]
Abstract
Curcuma is an important member of Zingiberaceae. Many species of this genus are widely used in traditional medicine and have important cultural value in East Asia. Among them, C. longa is considered to be the main source of curcumin and has a very wide range of uses. The rapid development of molecular phylogeny has deepened our understanding of taxonomy and evolution of Curcuma. However, little is known about the chloroplast genome phylogeny and the genetic bases of adaptative evolution. In this work, we sequenced the complete chloroplast genome of 4 Curcuma species. Curcuma chloroplast genomes showed highly conserved structures and the length ranged from 159,423 bp to 152,723 bp. A total of 133 genes were observed. Multiple repeats and simple sequence repeats (SSRs) were detected. By comparing with related species, 7 highly variable regions were identified as potential specific DNA barcodes for species identification. Phylogenetic analysis of complete plastome sequences and specific data sets revealed discordance with expected genus boundary. Chloroplast phylogenetic relationships were better predicted by geography than by morphological and nuclear DNA, indicating a substantial existence of introgression. 9 genes were proved to have high posteriori probability in positive selection analysis, and 4 of them (psbA, psbD, PetA and rbcL) closely related to photosynthesis, implying that chloroplast genes may had undergone positive selection pressure in evolution. These results are of great significance for us to understand the genetic basis, phylogeny and adaptive evolution of Curcuma chloroplast.
Collapse
|
33
|
Beichman AC, Koepfli KP, Li G, Murphy W, Dobrynin P, Kliver S, Tinker MT, Murray MJ, Johnson J, Lindblad-Toh K, Karlsson EK, Lohmueller KE, Wayne RK. Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter. Mol Biol Evol 2019; 36:2631-2655. [PMID: 31212313 PMCID: PMC7967881 DOI: 10.1093/molbev/msz101] [Citation(s) in RCA: 37] [Impact Index Per Article: 7.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.
Collapse
Affiliation(s)
- Annabel C Beichman
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA
| | - Klaus-Peter Koepfli
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, DC
- Institute of Molecular and Cellular Biology, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Gang Li
- College of Life Science, Shaanxi Normal University, Xi’an, Shaanxi, China
| | - William Murphy
- Department of Veterinary Integrative Biosciences, Texas A&M University, College Station, TX
| | - Pasha Dobrynin
- Center for Species Survival, Smithsonian Conservation Biology Institute, National Zoological Park, Washington, DC
- Institute of Molecular and Cellular Biology, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Sergei Kliver
- Institute of Molecular and Cellular Biology, Siberian Branch of the Russian Academy of Sciences, Novosibirsk, Russian Federation
| | - Martin T Tinker
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, CA
| | | | - Jeremy Johnson
- Vertebrate Genome Biology, Broad Institute of MIT and Harvard, Cambridge, MA
| | - Kerstin Lindblad-Toh
- Vertebrate Genome Biology, Broad Institute of MIT and Harvard, Cambridge, MA
- Science for Life Laboratory, Department of Medical Biochemistry and Microbiology, Uppsala University, Uppsala, Sweden
| | - Elinor K Karlsson
- Vertebrate Genome Biology, Broad Institute of MIT and Harvard, Cambridge, MA
- Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA
| | - Kirk E Lohmueller
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA
- Interdepartmental Program in Bioinformatics, University of California, Los Angeles, CA
- Department of Human Genetics, David Geffen School of Medicine, University of California, Los Angeles, CA
| | - Robert K Wayne
- Department of Ecology and Evolutionary Biology, University of California, Los Angeles, CA
| |
Collapse
|
34
|
Braun EL. An evolutionary model motivated by physicochemical properties of amino acids reveals variation among proteins. Bioinformatics 2019; 34:i350-i356. [PMID: 29950007 PMCID: PMC6022633 DOI: 10.1093/bioinformatics/bty261] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
Motivation The relative rates of amino acid interchanges over evolutionary time are likely to vary among proteins. Variation in those rates has the potential to reveal information about constraints on proteins. However, the most straightforward model that could be used to estimate relative rates of amino acid substitution is parameter-rich and it is therefore impractical to use for this purpose. Results A six-parameter model of amino acid substitution that incorporates information about the physicochemical properties of amino acids was developed. It showed that amino acid side chain volume, polarity and aromaticity have major impacts on protein evolution. It also revealed variation among proteins in the relative importance of those properties. The same general approach can be used to improve the fit of empirical models such as the commonly used PAM and LG models. Availability and implementation Perl code and test data are available from https://github.com/ebraun68/sixparam. Supplementary information Supplementary data are available at Bioinformatics online.
Collapse
Affiliation(s)
- Edward L Braun
- Department of Biology and Genetics Institute, University of Florida, Gainesville, FL, USA
| |
Collapse
|
35
|
Abstract
For nearly a century adaptive landscapes have provided overviews of the evolutionary process and yet they remain metaphors. We redefine adaptive landscapes in terms of biological processes rather than descriptive phenomenology. We focus on the underlying mechanisms that generate emergent properties such as epistasis, dominance, trade-offs and adaptive peaks. We illustrate the utility of landscapes in predicting the course of adaptation and the distribution of fitness effects. We abandon aged arguments concerning landscape ruggedness in favor of empirically determining landscape architecture. In so doing, we transform the landscape metaphor into a scientific framework within which causal hypotheses can be tested.
Collapse
Affiliation(s)
- Xiao Yi
- BioTechnology Institute, University of Minnesota, St. Paul, MN
| | - Antony M Dean
- BioTechnology Institute, University of Minnesota, St. Paul, MN.,Department of Ecology, Evolution, and Behavior, University of Minnesota, St. Paul, MN
| |
Collapse
|
36
|
Bolívar P, Guéguen L, Duret L, Ellegren H, Mugal CF. GC-biased gene conversion conceals the prediction of the nearly neutral theory in avian genomes. Genome Biol 2019; 20:5. [PMID: 30616647 PMCID: PMC6322265 DOI: 10.1186/s13059-018-1613-z] [Citation(s) in RCA: 15] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/12/2018] [Accepted: 12/17/2018] [Indexed: 11/10/2022] Open
Abstract
BACKGROUND The nearly neutral theory of molecular evolution predicts that the efficacy of natural selection increases with the effective population size. This prediction has been verified by independent observations in diverse taxa, which show that life-history traits are strongly correlated with measures of the efficacy of selection, such as the dN/dS ratio. Surprisingly, avian taxa are an exception to this theory because correlations between life-history traits and dN/dS are apparently absent. Here we explore the role of GC-biased gene conversion on estimates of substitution rates as a potential driver of these unexpected observations. RESULTS We analyze the relationship between dN/dS estimated from alignments of 47 avian genomes and several proxies for effective population size. To distinguish the impact of GC-biased gene conversion from selection, we use an approach that accounts for non-stationary base composition and estimate dN/dS separately for changes affected or unaffected by GC-biased gene conversion. This analysis shows that the impact of GC-biased gene conversion on substitution rates can explain the lack of correlations between life-history traits and dN/dS. Strong correlations between life-history traits and dN/dS are recovered after accounting for GC-biased gene conversion. The correlations are robust to variation in base composition and genomic location. CONCLUSIONS Our study shows that gene sequence evolution across a wide range of avian lineages meets the prediction of the nearly neutral theory, the efficacy of selection increases with effective population size. Moreover, our study illustrates that accounting for GC-biased gene conversion is important to correctly estimate the strength of selection.
Collapse
Affiliation(s)
- Paulina Bolívar
- Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Laurent Guéguen
- Laboratoire de Biologie et Biométrie Évolutive CNRS UMR 5558, Université Claude Bernard Lyon 1, Lyon, France
| | - Laurent Duret
- Laboratoire de Biologie et Biométrie Évolutive CNRS UMR 5558, Université Claude Bernard Lyon 1, Lyon, France
| | - Hans Ellegren
- Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Carina F. Mugal
- Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| |
Collapse
|
37
|
Herman A, Brandvain Y, Weagley J, Jeffery WR, Keene AC, Kono TJY, Bilandžija H, Borowsky R, Espinasa L, O'Quin K, Ornelas-García CP, Yoshizawa M, Carlson B, Maldonado E, Gross JB, Cartwright RA, Rohner N, Warren WC, McGaugh SE. The role of gene flow in rapid and repeated evolution of cave-related traits in Mexican tetra, Astyanax mexicanus. Mol Ecol 2018; 27:4397-4416. [PMID: 30252986 PMCID: PMC6261294 DOI: 10.1111/mec.14877] [Citation(s) in RCA: 101] [Impact Index Per Article: 16.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/22/2018] [Revised: 08/08/2018] [Accepted: 08/19/2018] [Indexed: 12/13/2022]
Abstract
Understanding the molecular basis of repeatedly evolved phenotypes can yield key insights into the evolutionary process. Quantifying gene flow between populations is especially important in interpreting mechanisms of repeated phenotypic evolution, and genomic analyses have revealed that admixture occurs more frequently between diverging lineages than previously thought. In this study, we resequenced 47 whole genomes of the Mexican tetra from three cave populations, two surface populations and outgroup samples. We confirmed that cave populations are polyphyletic and two Astyanax mexicanus lineages are present in our data set. The two lineages likely diverged much more recently than previous mitochondrial estimates of 5-7 mya. Divergence of cave populations from their phylogenetically closest surface population likely occurred between ~161 and 191 k generations ago. The favoured demographic model for most population pairs accounts for divergence with secondary contact and heterogeneous gene flow across the genome, and we rigorously identified gene flow among all lineages sampled. Therefore, the evolution of cave-related traits occurred more rapidly than previously thought, and trogolomorphic traits are maintained despite gene flow with surface populations. The recency of these estimated divergence events suggests that selection may drive the evolution of cave-derived traits, as opposed to disuse and drift. Finally, we show that a key trogolomorphic phenotype QTL is enriched for genomic regions with low divergence between caves, suggesting that regions important for cave phenotypes may be transferred between caves via gene flow. Our study shows that gene flow must be considered in studies of independent, repeated trait evolution.
Collapse
Affiliation(s)
- Adam Herman
- Plant and Microbial Biology, Gortner Lab, University of Minnesota, Saint Paul, Minnesota
- Department of Molecular Biology, Rudjer Boskovic Institute, Zagreb, Croatia
| | - Yaniv Brandvain
- Plant and Microbial Biology, Gortner Lab, University of Minnesota, Saint Paul, Minnesota
| | - James Weagley
- Ecology, Evolution, and Behavior, Gortner Lab, University of Minnesota, Saint Paul, Minnesota
| | - William R Jeffery
- Department of Biology, University of Maryland, College Park, Maryland
| | - Alex C Keene
- Department of Biological Sciences, Florida Atlantic University, Jupiter, Florida
| | - Thomas J Y Kono
- Minnesota Supercomputing Institute, University of Minnesota, Minneapolis, Minnesota
| | - Helena Bilandžija
- Department of Molecular Biology, Rudjer Boskovic Institute, Zagreb, Croatia
- Department of Biology, University of Maryland, College Park, Maryland
| | | | - Luis Espinasa
- School of Science, Marist College, Poughkeepsie, New York
| | - Kelly O'Quin
- Department of Biology, Centre College, Danville, Kentucky
| | - Claudia P Ornelas-García
- Departamento de Zoología, Instituto de Biología, Universidad Nacional Autónoma de México, Coyoacán, Mexico
| | - Masato Yoshizawa
- Department of Biology, University of Hawai'i at Mānoa, Honolulu, Hawaii
| | - Brian Carlson
- Department of Biology, College of Wooster, Wooster, Ohio
| | - Ernesto Maldonado
- Unidad Académica de Sistemas Arrecifales, Instituto de Ciencias del Mar y Limnología, Universidad Nacional Autónoma de México, Puerto Morelos, Mexico
| | - Joshua B Gross
- Department of Biological Sciences, University of Cincinnati, Cincinnati, Ohio
| | - Reed A Cartwright
- The Biodesign Institute, Arizona State University, Tempe, Arizona
- School of Life Sciences, Arizona State University, Tempe, Arizona
| | - Nicolas Rohner
- Stowers Institute for Medical Research, Kansas City, Missouri
- Department of Molecular and Integrative Physiology, The University of Kansas Medical Center, Kansas City, Kansas
| | - Wesley C Warren
- McDonnell Genome Institute, Washington University, St Louis, Missouri
| | - Suzanne E McGaugh
- Department of Molecular Biology, Rudjer Boskovic Institute, Zagreb, Croatia
| |
Collapse
|
38
|
Ezaz T, Edwards SV. Editorial: Evolutionary Feedbacks Between Population Biology and Genome Architecture. Front Genet 2018; 9:329. [PMID: 30186309 PMCID: PMC6112050 DOI: 10.3389/fgene.2018.00329] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Key Words] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/16/2018] [Accepted: 07/31/2018] [Indexed: 01/11/2023] Open
Affiliation(s)
- Tariq Ezaz
- Institute for Applied Ecology, University of Canberra Canberra, ACT, Australia
| | - Scott V Edwards
- Department of Organismic and Evolutionary Biology and Museum of Comparative Zoology, Harvard University Cambridge, MA, United States
| |
Collapse
|
39
|
The direction of protein evolution is destined by the stability. Biochimie 2018; 150:100-109. [DOI: 10.1016/j.biochi.2018.05.006] [Citation(s) in RCA: 13] [Impact Index Per Article: 2.2] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/22/2018] [Accepted: 05/15/2018] [Indexed: 01/29/2023]
|
40
|
Mitaku S, Sawada R. Biological meaning of "habitable zone" in nucleotide composition space. Biophys Physicobiol 2018; 15:75-85. [PMID: 29892513 PMCID: PMC5992858 DOI: 10.2142/biophysico.15.0_75] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/20/2017] [Accepted: 03/17/2018] [Indexed: 12/01/2022] Open
Abstract
Organisms generally display two contrasting properties: large biodiversity and a uniform state of "life". In this study, we focused on the question of how genome sequences describe "life" where a large number of biomolecules are harmonized. We analyzed the whole genome sequence of 2664 organisms, paying attention to the nucleotide composition which is an intensive parameter from the genome sequence. The results showed that all organisms were plotted in narrow regions of the nucleotide composition space of the first and second letters of the codon. Since all genome sequences overlap irrespective of the living environment, it can be called a "habitable zone". The habitable zone deviates by 500 times the standard deviation from the nucleotide composition expected from the random sequence, indicating that unexpectedly rare sequences are realized. Furthermore, we found that the habitable zones at the first and second letters of the codon serve as the background mechanisms for the functional network of biological systems. The habitable zone at the second letter of the codon controls the formation of transmembrane regions and the habitable zone at the first letter controls the formation of molecular recognition unit. These analyses showed that the habitable zone of the nucleotide composition space and the exquisite arrangement of amino acids in the codon table are conjugated to form biological systems. Finally, we discussed the evolution of the higher order of genome sequences.
Collapse
Affiliation(s)
- Shigeki Mitaku
- Emeritus Professor of Nagoya University, Kokubunji, Tokyo 185-0021, Japan
| | - Ryusuke Sawada
- Division of System Cohort, Medical Institute of Bioregulation, Kyushu University, Fukuoka 812-8582, Japan
| |
Collapse
|
41
|
Protein Evolution is Potentially Governed by Protein Stability: Directed Evolution of an Esterase from the Hyperthermophilic Archaeon Sulfolobus tokodaii. J Mol Evol 2018; 86:283-292. [DOI: 10.1007/s00239-018-9843-y] [Citation(s) in RCA: 11] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2018] [Accepted: 04/18/2018] [Indexed: 11/27/2022]
|
42
|
Abstract
The rice blast fungus Magnaporthe oryzae (syn., Pyricularia oryzae) is both a threat to global food security and a model for plant pathology. Molecular pathologists need an accurate understanding of the origins and line of descent of M. oryzae populations in order to identify the genetic and functional bases of pathogen adaptation and to guide the development of more effective control strategies. We used a whole-genome sequence analysis of samples from different times and places to infer details about the genetic makeup of M. oryzae from a global collection of isolates. Analyses of population structure identified six lineages within M. oryzae, including two pandemic on japonica and indica rice, respectively, and four lineages with more restricted distributions. Tip-dating calibration indicated that M. oryzae lineages separated about a millennium ago, long after the initial domestication of rice. The major lineage endemic to continental Southeast Asia displayed signatures of sexual recombination and evidence of DNA acquisition from multiple lineages. Tests for weak natural selection revealed that the pandemic spread of clonal lineages entailed an evolutionary “cost,” in terms of the accumulation of deleterious mutations. Our findings reveal the coexistence of multiple endemic and pandemic lineages with contrasting population and genetic characteristics within a widely distributed pathogen. The rice blast fungus Magnaporthe oryzae (syn., Pyricularia oryzae) is a textbook example of a rapidly adapting pathogen, and it is responsible for one of the most damaging diseases of rice. Improvements in our understanding of Magnaporthe oryzae’s diversity and evolution are required to guide the development of more effective control strategies. We used genome sequencing data for samples from around the world to infer the evolutionary history of M. oryzae. We found that M. oryzae diversified about 1,000 years ago, separating into six main lineages: two pandemic on japonica and indica rice, respectively, and four with more restricted distributions. We also found that a lineage endemic to continental Southeast Asia displayed signatures of sexual recombination and the acquisition of genetic material from multiple lineages. This work provides a population-level genomic framework for defining molecular markers for the control of rice blast and investigations of the molecular basis of differences in pathogenicity between M. oryzae lineages.
Collapse
|
43
|
Platt A, Weber CC, Liberles DA. Protein evolution depends on multiple distinct population size parameters. BMC Evol Biol 2018; 18:17. [PMID: 29422024 PMCID: PMC5806465 DOI: 10.1186/s12862-017-1085-x] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/25/2017] [Accepted: 11/20/2017] [Indexed: 01/08/2023] Open
Abstract
That population size affects the fate of new mutations arising in genomes, modulating both how frequently they arise and how efficiently natural selection is able to filter them, is well established. It is therefore clear that these distinct roles for population size that characterize different processes should affect the evolution of proteins and need to be carefully defined. Empirical evidence is consistent with a role for demography in influencing protein evolution, supporting the idea that functional constraints alone do not determine the composition of coding sequences. Given that the relationship between population size, mutant fitness and fixation probability has been well characterized, estimating fitness from observed substitutions is well within reach with well-formulated models. Molecular evolution research has, therefore, increasingly begun to leverage concepts from population genetics to quantify the selective effects associated with different classes of mutation. However, in order for this type of analysis to provide meaningful information about the intra- and inter-specific evolution of coding sequences, a clear definition of concepts of population size, what they influence, and how they are best parameterized is essential. Here, we present an overview of the many distinct concepts that “population size” and “effective population size” may refer to, what they represent for studying proteins, and how this knowledge can be harnessed to produce better specified models of protein evolution.
Collapse
Affiliation(s)
- Alexander Platt
- Department of Biology and Center for Computational Genetics and Genomics, Temple University, Philadelphia, 19121, USA
| | - Claudia C Weber
- Department of Biology and Center for Computational Genetics and Genomics, Temple University, Philadelphia, 19121, USA
| | - David A Liberles
- Department of Biology and Center for Computational Genetics and Genomics, Temple University, Philadelphia, 19121, USA.
| |
Collapse
|
44
|
Nunez JCB, Elyanow RG, Ferranti DA, Rand DM. Population Genomics and Biogeography of the Northern Acorn Barnacle (Semibalanus balanoides) Using Pooled Sequencing Approaches. POPULATION GENOMICS 2018. [DOI: 10.1007/13836_2018_58] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/08/2023]
|
45
|
Vahdati AR, Wagner A. Population Size Affects Adaptation in Complex Ways: Simulations on Empirical Adaptive Landscapes. Evol Biol 2017. [DOI: 10.1007/s11692-017-9440-9] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/28/2022]
|
46
|
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
|
47
|
Ostrander EA, Wayne RK, Freedman AH, Davis BW. Demographic history, selection and functional diversity of the canine genome. Nat Rev Genet 2017; 18:705-720. [DOI: 10.1038/nrg.2017.67] [Citation(s) in RCA: 101] [Impact Index Per Article: 14.4] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
|
48
|
Evolutionary coincidence of adaptive changes in exuperantia and the emergence of bicoid in Cyclorrhapha (Diptera). Dev Genes Evol 2017; 227:355-365. [PMID: 28894941 PMCID: PMC5597691 DOI: 10.1007/s00427-017-0594-3] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2017] [Accepted: 08/20/2017] [Indexed: 11/15/2022]
Abstract
The great radiation in the infraorder Cyclorrhapha involved several morphological and molecular changes, including important changes in anterior egg development. During Drosophila oogenesis, exuperantia (exu) is critical for localizing bicoid (bcd) messenger RNA (mRNA) to the anterior region of the oocyte. Because it is phylogenetically older than bcd, which is exclusive to Cyclorrhapha, we hypothesize that exu has undergone adaptive changes to enable this new function. Although exu has been well studied in Drosophila, there is no functional or transcriptional information about it in any other Diptera. Here, we investigate exu in the South American fruit fly Anastrepha fraterculus, a Cyclorrhapha of great agricultural importance that have lost bcd, aiming to understand the evolution of exu in this infraorder. We assessed its pattern of gene expression in A. fraterculus by analyzing transcriptomes from cephalic and reproductive tissues. A combination of next-generation data with classical sequencing procedures enabled identification of the structure of exu and its alternative transcripts in this species. In addition to the sex-specific isoforms described for Drosophila, we found that not only exu is expressed in heads, but this is mediated by two transcripts with a specific 5′UTR exon—likely a result from usage of a third promoter. Furthermore, we tested the hypothesis that exu is evolving under positive selection in Cyclorrhapha after divergence from lower Diptera. We found evidence of positive selection at two important exu domains, EXO-like and SAM-like, both involved with mRNA binding during bcd mRNA localization in Drosophila, which could reflect its cooptation for the new function of bcd mRNA localization in Cyclorrhapha.
Collapse
|
49
|
Vahdati AR, Sprouffske K, Wagner A. Effect of Population Size and Mutation Rate on the Evolution of RNA Sequences on an Adaptive Landscape Determined by RNA Folding. Int J Biol Sci 2017; 13:1138-1151. [PMID: 29104505 PMCID: PMC5666329 DOI: 10.7150/ijbs.19436] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/01/2017] [Accepted: 07/05/2017] [Indexed: 02/04/2023] Open
Abstract
The dynamics of populations evolving on an adaptive landscape depends on multiple factors, including the structure of the landscape, the rate of mutations, and effective population size. Existing theoretical work often makes ad hoc and simplifying assumptions about landscape structure, whereas experimental work can vary important parameters only to a limited extent. We here overcome some of these limitations by simulating the adaptive evolution of RNA molecules, whose fitness is determined by the thermodynamics of RNA secondary structure folding. We study the influence of mutation rates and population sizes on final mean population fitness, on the substitution rates of mutations, and on population diversity. We show that evolutionary dynamics cannot be understood as a function of mutation rate µ, population size N, or population mutation rate Nµ alone. For example, at a given mutation rate, clonal interference prevents the fixation of beneficial mutations as population size increases, but larger populations still arrive at a higher mean fitness. In addition, at the highest population mutation rates we study, mean final fitness increases with population size, because small populations are driven to low fitness by the relatively higher incidence of mutations they experience. Our observations show that mutation rate and population size can interact in complex ways to influence the adaptive dynamics of a population on a biophysically motivated fitness landscape.
Collapse
Affiliation(s)
- Ali R Vahdati
- Institute of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland.,The Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Kathleen Sprouffske
- Institute of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland.,The Swiss Institute of Bioinformatics, Lausanne, Switzerland
| | - Andreas Wagner
- Institute of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland.,The Swiss Institute of Bioinformatics, Lausanne, Switzerland.,The Santa Fe Institute, Santa Fe, USA
| |
Collapse
|
50
|
Szövényi P, Ullrich KK, Rensing SA, Lang D, van Gessel N, Stenøien HK, Conti E, Reski R. Selfing in Haploid Plants and Efficacy of Selection: Codon Usage Bias in the Model Moss Physcomitrella patens. Genome Biol Evol 2017; 9:1528-1546. [PMID: 28549175 PMCID: PMC5507605 DOI: 10.1093/gbe/evx098] [Citation(s) in RCA: 14] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/25/2017] [Indexed: 12/15/2022] Open
Abstract
A long-term reduction in effective population size will lead to major shift in genome evolution. In particular, when effective population size is small, genetic drift becomes dominant over natural selection. The onset of self-fertilization is one evolutionary event considerably reducing effective size of populations. Theory predicts that this reduction should be more dramatic in organisms capable for haploid than for diploid selfing. Although theoretically well-grounded, this assertion received mixed experimental support. Here, we test this hypothesis by analyzing synonymous codon usage bias of genes in the model moss Physcomitrella patens frequently undergoing haploid selfing. In line with population genetic theory, we found that the effect of natural selection on synonymous codon usage bias is very weak. Our conclusion is supported by four independent lines of evidence: 1) Very weak or nonsignificant correlation between gene expression and codon usage bias, 2) no increased codon usage bias in more broadly expressed genes, 3) no evidence that codon usage bias would constrain synonymous and nonsynonymous divergence, and 4) predominant role of genetic drift on synonymous codon usage predicted by a model-based analysis. These findings show striking similarity to those observed in AT-rich genomes with weak selection for optimal codon usage and GC content overall. Our finding is in contrast to a previous study reporting adaptive codon usage bias in the moss P. patens.
Collapse
Affiliation(s)
- Péter Szövényi
- Department of Systematic and Evolutionary Botany, University of Zurich, Switzerland
| | - Kristian K. Ullrich
- Plant Cell Biology, Faculty of Biology, University of Marburg, Germany
- Present address: Max-Planck-Insitut für Evolutionsbiologie, Plön, Germany
| | - Stefan A. Rensing
- Plant Cell Biology, Faculty of Biology, University of Marburg, Germany
- BIOSS—Centre for Biological Signalling Studies, University of Freiburg, Germany
| | - Daniel Lang
- Plant Genome and Systems Biology, Helmholtz Zentrum München, German Research Center for Environmental Health, Neuherberg, Germany
| | - Nico van Gessel
- Plant Biotechnology, Faculty of Biology, University of Freiburg, Germany
| | | | - Elena Conti
- Department of Systematic and Evolutionary Botany, University of Zurich, Switzerland
| | - Ralf Reski
- BIOSS—Centre for Biological Signalling Studies, University of Freiburg, Germany
- Plant Biotechnology, Faculty of Biology, University of Freiburg, Germany
| |
Collapse
|