1
|
Zuo B, Nneji LM, Sun YB. Comparative genomics reveals insights into anuran genome size evolution. BMC Genomics 2023; 24:379. [PMID: 37415107 DOI: 10.1186/s12864-023-09499-8] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2023] [Accepted: 06/30/2023] [Indexed: 07/08/2023] Open
Abstract
BACKGROUND Amphibians, particularly anurans, display an enormous variation in genome size. Due to the unavailability of whole genome datasets in the past, the genomic elements and evolutionary causes of anuran genome size variation are poorly understood. To address this, we analyzed whole-genome sequences of 14 anuran species ranging in size from 1.1 to 6.8 Gb. By annotating multiple genomic elements, we investigated the genomic correlates of anuran genome size variation and further examined whether the genome size relates to habitat types. RESULTS Our results showed that intron expansions or contraction and Transposable Elements (TEs) diversity do not contribute significantly to genome size variation. However, the recent accumulation of transposable elements (TEs) and the lack of deletion of ancient TEs primarily accounted for the evolution of anuran genome sizes. Our study showed that the abundance and density of simple repeat sequences positively correlate with genome size. Ancestral state reconstruction revealed that genome size exhibits a taxon-specific pattern of evolution, with families Bufonidae and Pipidae experiencing extreme genome expansion and contraction events, respectively. Our result showed no relationship between genome size and habitat types, although large genome-sized species are predominantly found in humid habitats. CONCLUSIONS Overall, our study identified the genomic element and their evolutionary dynamics accounting for anuran genome size variation, thus paving a path to a greater understanding of the size evolution of the genome in amphibians.
Collapse
Affiliation(s)
- Bin Zuo
- Ministry of Education Key Laboratory for Transboundary Ecosecurity of Southwest China, Yunnan Key Laboratory of Plant Reproductive Adaptation and Evolutionary Ecology, Institute of Biodiversity, School of Ecology and Environmental Science, Yunnan University, Kunming, Yunnan, 650504, China
| | - Lotanna Micah Nneji
- Department of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ, 08544, USA
| | - Yan-Bo Sun
- Ministry of Education Key Laboratory for Transboundary Ecosecurity of Southwest China, Yunnan Key Laboratory of Plant Reproductive Adaptation and Evolutionary Ecology, Institute of Biodiversity, School of Ecology and Environmental Science, Yunnan University, Kunming, Yunnan, 650504, China.
- Laboratory for Conservation and Utilization of Bio-resources, Yunnan University, Kunming, 650091, China.
| |
Collapse
|
2
|
Kushwaha B, Nagpure NS, Srivastava S, Pandey M, Kumar R, Raizada S, Agarwal S, Singh M, Basheer VS, Kumar RG, Das P, Das SP, Patnaik S, Bit A, Srivastava SK, Vishwakarma AL, Joshi CG, Kumar D, Jena JK. Genome size estimation and its associations with body length, chromosome number and evolution in teleost fishes. Gene 2023; 864:147294. [PMID: 36858189 DOI: 10.1016/j.gene.2023.147294] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2022] [Revised: 02/07/2023] [Accepted: 02/15/2023] [Indexed: 03/02/2023]
Abstract
Precise estimation of genome size (GS) is vital for various genomic studies, such as deciding genome sequencing depth, genome assembly, biodiversity documentation, evolution, genetic disorders studies, duplication events etc. Animal Genome Size Database provides GS of over 2050 fish species, which ranges from 0.35 pg in pufferfish (Tetraodon nigroviridis) to 132.83 pg in marbled lungfish (Protopterus aethiopicus). The GS of majority of the fishes inhabiting waters of Indian subcontinent are still missing. In present study, we estimated GS of 51 freshwater teleost (31 commercially important, 7 vulnerable and 13 ornamental species) that ranged from 0.58 pg in banded gourami (Trichogaster fasciata) to 1.92 pg in scribbled goby (Awaous grammepomus). Substantial variation in GS was observed within the same fish orders (0.64-1.45 pg in cypriniformes, 0.70-1.41 pg in siluriformes and 0.58-1.92 pg in perciformes). We examined the relationship between the GS, chromosome number and body length across all the fishes. Body length was found to be associated with GS, whereas no relationship was noticed between the GS and the chromosome number. The analysis using ancestral information revealed haploid chromosome number 25, 27 and 24 for the most recent common ancestor of cypriniformes, siluriformes and perciformes, respectively. The study led to generation of new records on GS of 43 fish species and revalidated records for 8 species. The finding is valuable resource for further research in the areas of fish genomics, molecular ecology and evolutionary conservation genetics.
Collapse
Affiliation(s)
- Basdeo Kushwaha
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India.
| | - Naresh S Nagpure
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| | - Shreya Srivastava
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| | - Manmohan Pandey
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| | - Ravindra Kumar
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India.
| | - Sudhir Raizada
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| | - Suyash Agarwal
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| | - Mahender Singh
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| | - Valaparamail S Basheer
- PMFGR Division, ICAR-National Bureau of Fish Genetic Resources, CMFRI Campus, Ernakulam North, P.O. Kochi, 682 018 Kerala, India
| | - Rahul G Kumar
- PMFGR Division, ICAR-National Bureau of Fish Genetic Resources, CMFRI Campus, Ernakulam North, P.O. Kochi, 682 018 Kerala, India
| | - Paramananda Das
- ICAR-Central Institute of Freshwater Aquaculture, Kausalyanga, Bhubaneswar, 751 002 Odisha, India
| | - Sofia P Das
- ICAR-Central Institute of Freshwater Aquaculture, Kausalyanga, Bhubaneswar, 751 002 Odisha, India
| | - Siddhi Patnaik
- ICAR-Central Institute of Freshwater Aquaculture, Kausalyanga, Bhubaneswar, 751 002 Odisha, India
| | - Amrita Bit
- ICAR-Central Institute of Freshwater Aquaculture, Kausalyanga, Bhubaneswar, 751 002 Odisha, India
| | - Satish Kumar Srivastava
- Experimental Field Centre, ICAR-Directorate of Coldwater Fisheries Research, Champawat, 262 523 Uttarakhand, India
| | - Achchhe L Vishwakarma
- Flow Cytometry Lab, SAIF Division, CSIR-Central Drug Research Institute, Lucknow, 226 031 Uttar Pradesh, India
| | - Chaitanya G Joshi
- Department of Animal Biotechnology, College of Veterinary Science and Animal Husbandry, Anand Agricultural University, Anand, Gujarat 388 001, India
| | - Dinesh Kumar
- Centre for Agricultural Bio-informatics, ICAR-Indian Agricultural Statistics Research Institute, Library Avenue, New Delhi 110 012, India
| | - Joy K Jena
- ICAR- National Bureau of Fish Genetic Resources, Canal Ring Road, P.O. Dilkusha, Lucknow, 226 002 Uttar Pradesh, India
| |
Collapse
|
3
|
On the Base Composition of Transposable Elements. Int J Mol Sci 2022; 23:ijms23094755. [PMID: 35563146 PMCID: PMC9099904 DOI: 10.3390/ijms23094755] [Citation(s) in RCA: 5] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/24/2022] [Revised: 04/22/2022] [Accepted: 04/23/2022] [Indexed: 01/27/2023] Open
Abstract
Transposable elements exhibit a base composition that is often different from the genomic average and from hosts’ genes. The most common compositional bias is towards Adenosine and Thymine, although this bias is not universal, and elements with drastically different base composition can coexist within the same genome. The AT-richness of transposable elements is apparently maladaptive because it results in poor transcription and sub-optimal translation of proteins encoded by the elements. The cause(s) of this unusual base composition remain unclear and have yet to be investigated. Here, I review what is known about the nucleotide content of transposable elements and how this content can affect the genome of their host as well as their own replication. The compositional bias of transposable elements could result from several non-exclusive processes including horizontal transfer, mutational bias, and selection. It appears that mutation alone cannot explain the high AT-content of transposons and that selection plays a major role in the evolution of the compositional bias. The reason why selection would favor a maladaptive nucleotide content remains however unexplained and is an area of investigation that clearly deserves attention.
Collapse
|
4
|
Lehmann R, Kovařík A, Ocalewicz K, Kirtiklis L, Zuccolo A, Tegner JN, Wanzenböck J, Bernatchez L, Lamatsch DK, Symonová R. DNA Transposon Expansion is Associated with Genome Size Increase in Mudminnows. Genome Biol Evol 2021; 13:6380143. [PMID: 34599322 PMCID: PMC8557787 DOI: 10.1093/gbe/evab228] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 09/27/2021] [Indexed: 12/20/2022] Open
Abstract
Genome sizes of eukaryotic organisms vary substantially, with whole-genome duplications (WGD) and transposable element expansion acting as main drivers for rapid genome size increase. The two North American mudminnows, Umbra limi and Umbra pygmaea, feature genomes about twice the size of their sister lineage Esocidae (e.g., pikes and pickerels). However, it is unknown whether all Umbra species share this genome expansion and which causal mechanisms drive this expansion. Using flow cytometry, we find that the genome of the European mudminnow is expanded similarly to both North American species, ranging between 4.5 and 5.4 pg per diploid nucleus. Observed blocks of interstitially located telomeric repeats in U. limi suggest frequent Robertsonian rearrangements in its history. Comparative analyses of transcriptome and genome assemblies show that the genome expansion in Umbra is driven by the expansion of DNA transposon and unclassified repeat sequences without WGD. Furthermore, we find a substantial ongoing expansion of repeat sequences in the Alaska blackfish Dallia pectoralis, the closest relative to the family Umbridae, which might mark the beginning of a similar genome expansion. Our study suggests that the genome expansion in mudminnows, driven mainly by transposon expansion, but not WGD, occurred before the separation into the American and European lineage.
Collapse
Affiliation(s)
- Robert Lehmann
- Division of Biological and Environmental Sciences & Engineering, Computer, Electrical and Mathematical Sciences and Engineering Division, King Abdullah University of Science and Technology, Thuwal, Kingdom of Saudi Arabia
| | - Aleš Kovařík
- Laboratory of Molecular Epigenetics, Institute of Biophysics, Czech Academy of Science, Brno, Czech Republic
| | - Konrad Ocalewicz
- Department of Marine Biology and Ecology, Institute of Oceanography, Faculty of Oceanography and Geography, University of Gdansk, Gdansk, Poland
| | - Lech Kirtiklis
- Department of Zoology, Faculty of Biology and Biotechnology, University of Warmia and Mazury, Olsztyn, Poland
| | - Andrea Zuccolo
- Center for Desert Agriculture, Biological and Environmental Sciences & Engineering Division (BESE), King Abdullah University of Science and Technology, Thuwal, Kingdom of Saudi Arabia.,Institute of Life Sciences, Scuola Superiore Sant'Anna, Pisa, Italy
| | - Jesper N Tegner
- Division of Biological and Environmental Sciences & Engineering, Computer, Electrical and Mathematical Sciences and Engineering Division, King Abdullah University of Science and Technology, Thuwal, Kingdom of Saudi Arabia
| | - Josef Wanzenböck
- Research Department for Limnology Mondsee, University of Innsbruck, Mondsee, Austria
| | - Louis Bernatchez
- Department of Biology, IBIS (Institut de Biologie Intégrative et des Systèmes), Université Laval, Québec, QC, Canada
| | - Dunja K Lamatsch
- Research Department for Limnology Mondsee, University of Innsbruck, Mondsee, Austria
| | - Radka Symonová
- Department of Bioinformatics, Wissenschaftzentrum Weihenstephan, Technische Universität München, Freising, Germany.,Department of Biology, Faculty of Biology, University of Hradec Kralove, Czech Republic
| |
Collapse
|
5
|
Loewenthal G, Rapoport D, Avram O, Moshe A, Wygoda E, Itzkovitch A, Israeli O, Azouri D, Cartwright RA, Mayrose I, Pupko T. A probabilistic model for indel evolution: differentiating insertions from deletions. Mol Biol Evol 2021; 38:5769-5781. [PMID: 34469521 PMCID: PMC8662616 DOI: 10.1093/molbev/msab266] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/27/2022] Open
Abstract
Insertions and deletions (indels) are common molecular evolutionary events. However, probabilistic models for indel evolution are under-developed due to their computational complexity. Here, we introduce several improvements to indel modeling: 1) While previous models for indel evolution assumed that the rates and length distributions of insertions and deletions are equal, here we propose a richer model that explicitly distinguishes between the two; 2) we introduce numerous summary statistics that allow approximate Bayesian computation-based parameter estimation; 3) we develop a method to correct for biases introduced by alignment programs, when inferring indel parameters from empirical data sets; and 4) using a model-selection scheme, we test whether the richer model better fits biological data compared with the simpler model. Our analyses suggest that both our inference scheme and the model-selection procedure achieve high accuracy on simulated data. We further demonstrate that our proposed richer model better fits a large number of empirical data sets and that, for the majority of these data sets, the deletion rate is higher than the insertion rate.
Collapse
Affiliation(s)
- Gil Loewenthal
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Dana Rapoport
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Oren Avram
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Asher Moshe
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Elya Wygoda
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Alon Itzkovitch
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Omer Israeli
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Dana Azouri
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel.,School of Plant Sciences and Food Security, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Reed A Cartwright
- The Biodesign Institute, Arizona State University, Tempe, Arizona, USA.,School of Life Sciences, Arizona State University, Tempe, Arizona, USA
| | - Itay Mayrose
- School of Plant Sciences and Food Security, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| | - Tal Pupko
- The Shmunis School of Biomedicine and Cancer Research, George S. Wise Faculty of Life Sciences, Tel Aviv University, Tel Aviv 69978, Israel
| |
Collapse
|
6
|
Wang X, Morton JA, Pellicer J, Leitch IJ, Leitch AR. Genome downsizing after polyploidy: mechanisms, rates and selection pressures. THE PLANT JOURNAL : FOR CELL AND MOLECULAR BIOLOGY 2021; 107:1003-1015. [PMID: 34077584 DOI: 10.1111/tpj.15363] [Citation(s) in RCA: 30] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/11/2021] [Revised: 05/07/2021] [Accepted: 05/13/2021] [Indexed: 05/20/2023]
Abstract
An analysis of over 10 000 plant genome sizes (GSs) indicates that most species have smaller genomes than expected given the incidence of polyploidy in their ancestries, suggesting selection for genome downsizing. However, comparing ancestral GS with the incidence of ancestral polyploidy suggests that the rate of DNA loss following polyploidy is likely to have been very low (4-70 Mb/million years, 4-482 bp/generation). This poses a problem. How might such small DNA losses be visible to selection, overcome the power of genetic drift and drive genome downsizing? Here we explore that problem, focussing on the role that double-strand break (DSB) repair pathways (non-homologous end joining and homologous recombination) may have played. We also explore two hypotheses that could explain how selection might favour genome downsizing following polyploidy: to reduce (i) nitrogen (N) and phosphate (P) costs associated with nucleic acid synthesis in the nucleus and the transcriptome and (ii) the impact of scaling effects of GS on cell size, which influences CO2 uptake and water loss. We explore the hypothesis that losses of DNA must be fastest in early polyploid generations. Alternatively, if DNA loss is a more continuous process over evolutionary time, then we propose it is a byproduct of selection elsewhere, such as limiting the damaging activity of repetitive DNA. If so, then the impact of GS on photosynthesis, water use efficiency and/or nutrient costs at the nucleus level may be emergent properties, which have advantages, but not ones that could have been selected for over generational timescales.
Collapse
Affiliation(s)
- Xiaotong Wang
- Royal Botanic Gardens, Kew, Surrey, TW9 3AB, UK
- Queen Mary University of London, Mile End Road, London, E1 4NS, UK
| | - Joseph A Morton
- Royal Botanic Gardens, Kew, Surrey, TW9 3AB, UK
- Queen Mary University of London, Mile End Road, London, E1 4NS, UK
| | - Jaume Pellicer
- Royal Botanic Gardens, Kew, Surrey, TW9 3AB, UK
- Institut Botànic de Barcelona (IBB, CSIC-Ajuntament de Barcelona), Passeig del Migdia sn, Barcelona, 08038, Spain
| | | | - Andrew R Leitch
- Queen Mary University of London, Mile End Road, London, E1 4NS, UK
| |
Collapse
|
7
|
Sellis D, Guérin F, Arnaiz O, Pett W, Lerat E, Boggetto N, Krenek S, Berendonk T, Couloux A, Aury JM, Labadie K, Malinsky S, Bhullar S, Meyer E, Sperling L, Duret L, Duharcourt S. Massive colonization of protein-coding exons by selfish genetic elements in Paramecium germline genomes. PLoS Biol 2021; 19:e3001309. [PMID: 34324490 PMCID: PMC8354472 DOI: 10.1371/journal.pbio.3001309] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/11/2021] [Revised: 08/10/2021] [Accepted: 06/04/2021] [Indexed: 11/18/2022] Open
Abstract
Ciliates are unicellular eukaryotes with both a germline genome and a somatic genome in the same cytoplasm. The somatic macronucleus (MAC), responsible for gene expression, is not sexually transmitted but develops from a copy of the germline micronucleus (MIC) at each sexual generation. In the MIC genome of Paramecium tetraurelia, genes are interrupted by tens of thousands of unique intervening sequences called internal eliminated sequences (IESs), which have to be precisely excised during the development of the new MAC to restore functional genes. To understand the evolutionary origin of this peculiar genomic architecture, we sequenced the MIC genomes of 9 Paramecium species (from approximately 100 Mb in Paramecium aurelia species to >1.5 Gb in Paramecium caudatum). We detected several waves of IES gains, both in ancestral and in more recent lineages. While the vast majority of IESs are single copy in present-day genomes, we identified several families of mobile IESs, including nonautonomous elements acquired via horizontal transfer, which generated tens to thousands of new copies. These observations provide the first direct evidence that transposable elements can account for the massive proliferation of IESs in Paramecium. The comparison of IESs of different evolutionary ages indicates that, over time, IESs shorten and diverge rapidly in sequence while they acquire features that allow them to be more efficiently excised. We nevertheless identified rare cases of IESs that are under strong purifying selection across the aurelia clade. The cases examined contain or overlap cellular genes that are inactivated by excision during development, suggesting conserved regulatory mechanisms. Similar to the evolution of introns in eukaryotes, the evolution of Paramecium IESs highlights the major role played by selfish genetic elements in shaping the complexity of genome architecture and gene expression. A comparative genomics study of nine Paramecium species reveals successful invasion of genes by transposable elements in their germline genomes, showing that the internal eliminated sequences (IESs) followed an evolutionary trajectory remarkably similar to that of spliceosomal introns.
Collapse
Affiliation(s)
- Diamantis Sellis
- Université de Lyon, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR 5558, Villeurbanne, France
| | - Frédéric Guérin
- Université de Paris, CNRS, Institut Jacques Monod, Paris, France
| | - Olivier Arnaiz
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), Gif-sur-Yvette, France
| | - Walker Pett
- Université de Lyon, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR 5558, Villeurbanne, France
| | - Emmanuelle Lerat
- Université de Lyon, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR 5558, Villeurbanne, France
| | - Nicole Boggetto
- Université de Paris, CNRS, Institut Jacques Monod, Paris, France
| | - Sascha Krenek
- TU Dresden, Institute of Hydrobiology, Dresden, Germany
| | | | - Arnaud Couloux
- Génomique Métabolique, Genoscope, Institut de biologie François Jacob, CEA, CNRS, Université d’Évry, Université Paris-Saclay, Evry, France
| | - Jean-Marc Aury
- Génomique Métabolique, Genoscope, Institut de biologie François Jacob, CEA, CNRS, Université d’Évry, Université Paris-Saclay, Evry, France
| | - Karine Labadie
- Genoscope, Institut de biologie François-Jacob, Commissariat à l’Energie Atomique (CEA), Université Paris-Saclay, Evry, France
| | - Sophie Malinsky
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), Ecole Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
- Université de Paris, Paris, France
| | - Simran Bhullar
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), Ecole Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
| | - Eric Meyer
- Institut de Biologie de l’Ecole Normale Supérieure (IBENS), Ecole Normale Supérieure, CNRS, INSERM, Université PSL, Paris, France
| | - Linda Sperling
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), Gif-sur-Yvette, France
| | - Laurent Duret
- Université de Lyon, CNRS, Laboratoire de Biométrie et Biologie Evolutive UMR 5558, Villeurbanne, France
- * E-mail: (LD); (SD)
| | - Sandra Duharcourt
- Université de Paris, CNRS, Institut Jacques Monod, Paris, France
- * E-mail: (LD); (SD)
| |
Collapse
|
8
|
Liu Q, Deng Y, Song A, Xiang Y, Chen D, Wei L. Comparative analysis of mite genomes reveals positive selection for diet adaptation. Commun Biol 2021; 4:668. [PMID: 34083730 PMCID: PMC8175442 DOI: 10.1038/s42003-021-02173-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/14/2020] [Accepted: 04/30/2021] [Indexed: 02/06/2023] Open
Abstract
Diet is a powerful evolutionary force for species adaptation and diversification. Acari is one of the most abundant clades of Arachnida, exhibiting diverse dietary types, while the underlying genetic adaptive mechanisms are not fully understood. Based on comparative analyses of 15 Acari genomes, we found genetic bases for three specialized diets. Herbivores experienced stronger selection pressure than other groups; the olfactory genes and gene families involving metabolizing toxins showed strong adaptive signals. Genes and gene families related to anticoagulation, detoxification, and haemoglobin digestion were found to be under strong selection pressure or significantly expanded in the blood-feeding species. Lipid metabolism genes have a faster evolutionary rate and been subjected to greater selection pressures in fat-feeding species; one positively selected site in the fatty-acid amide hydrolases 2 gene was identified. Our research provides a new perspective for the evolution of Acari and offers potential target loci for novel pesticide development.
Collapse
Affiliation(s)
- Qiong Liu
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China
| | - Yuhua Deng
- Clinical Research Institute, The First People's Hospital of Foshan, Foshan, China
| | - An Song
- ShaanXi JunDa Forensic Medicine Expertise Station, The Fourth Military Medical University, Xi'an, China
| | - Yifan Xiang
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China
| | - De Chen
- MOE Key Laboratory for Biodiversity Science and Ecological Engineering, College of Life Sciences, Beijing Normal University, Beijing, China.
| | - Lai Wei
- State Key Laboratory of Ophthalmology, Zhongshan Ophthalmic Center, Sun Yat-sen University, Guangzhou, China.
| |
Collapse
|
9
|
Stelzer CP, Pichler M, Hatheuer A. Linking genome size variation to population phenotypic variation within the rotifer, Brachionus asplanchnoidis. Commun Biol 2021; 4:596. [PMID: 34011946 PMCID: PMC8134563 DOI: 10.1038/s42003-021-02131-z] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Accepted: 04/19/2021] [Indexed: 11/17/2022] Open
Abstract
Eukaryotic organisms usually contain much more genomic DNA than expected from their biological complexity. In explaining this pattern, selection-based hypotheses suggest that genome size evolves through selection acting on correlated life history traits, implicitly assuming the existence of phenotypic effects of (extra) genomic DNA that are independent of its information content. Here, we present conclusive evidence of such phenotypic effects within a well-mixed natural population that shows heritable variation in genome size. We found that genome size is positively correlated with body size, egg size, and embryonic development time in a population of the monogonont rotifer Brachionus asplanchnoidis. The effect on embryonic development time was mediated partly by an indirect effect (via egg size), and a direct effect, the latter indicating an increased replication cost of the larger amounts of DNA during mitosis. Our results suggest that selection-based change of genome size can operate in this population, provided it is strong enough to overcome drift or mutational change of genome size.
Collapse
Affiliation(s)
| | - Maria Pichler
- University of Innsbruck, Mondseestr. 9, 5310, Mondsee, Austria
| | - Anita Hatheuer
- University of Innsbruck, Mondseestr. 9, 5310, Mondsee, Austria
| |
Collapse
|
10
|
Lauria Sneideman MP, Meller VH. Drosophila Satellite Repeats at the Intersection of Chromatin, Gene Regulation and Evolution. PROGRESS IN MOLECULAR AND SUBCELLULAR BIOLOGY 2021; 60:1-26. [PMID: 34386870 DOI: 10.1007/978-3-030-74889-0_1] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 01/24/2023]
Abstract
Satellite repeats make up a large fraction of the genomes of many higher eukaryotes. Until recently these sequences were viewed as molecular parasites with few functions. Drosophila melanogaster and related species have a wealth of diverse satellite repeats. Comparative studies of Drosophilids have been instrumental in understanding how these rapidly evolving sequences change and move. Remarkably, satellite repeats have been found to modulate gene expression and mediate genetic conflicts between chromosomes and between closely related fly species. This suggests that satellites play a key role in speciation. We have taken advantage of the depth of research on satellite repeats in flies to review the known functions of these sequences and consider their central role in evolution and gene expression.
Collapse
Affiliation(s)
| | - Victoria H Meller
- Department of Biological Sciences, Wayne State University, Detroit, MI, USA.
| |
Collapse
|
11
|
Zhang SJ, Liu L, Yang R, Wang X. Genome Size Evolution Mediated by Gypsy Retrotransposons in Brassicaceae. GENOMICS PROTEOMICS & BIOINFORMATICS 2020; 18:321-332. [PMID: 33137519 PMCID: PMC7801240 DOI: 10.1016/j.gpb.2018.07.009] [Citation(s) in RCA: 13] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/09/2018] [Revised: 05/31/2018] [Accepted: 07/16/2018] [Indexed: 12/29/2022]
Abstract
The dynamic activity of transposable elements (TEs) contributes to the vast diversity of genome size and architecture among plants. Here, we examined the genomic distribution and transposition activity of long terminal repeat retrotransposons (LTR-RTs) in Arabidopsis thaliana (Ath) and three of its relatives, Arabidopsis lyrata (Aly), Eutrema salsugineum (Esa), and Schrenkiella parvula (Spa), in Brassicaceae. Our analyses revealed the distinct evolutionary dynamics of Gypsyretrotransposons, which reflects the different patterns of genome size changes of the four species over the past million years. The rate of Gypsy transposition in Aly is approximately five times more rapid than that of Ath and Esa, suggesting an expanding Aly genome. Gypsy insertions in Esa are strictly confined to pericentromeric heterochromatin and associated with dramatic centromere expansion. In contrast, Gypsy insertions in Spa have been largely suppressed over the last million years, likely as a result of a combination of an inherent molecular mechanism of preferential DNA removal and purifying selection at Gypsy elements. Additionally, species-specific clades of Gypsy elements shaped the distinct genome architectures of Aly and Esa.
Collapse
Affiliation(s)
- Shi-Jian Zhang
- Department of Crop Genomics and Bioinformatics, College of Agronomy and Biotechnology, National Maize Improvement Center of China, China Agricultural University, Beijing 100094, China
| | - Lei Liu
- Beijing Key Laboratory of Plant Resources Research and Development, School of Sciences, Beijing Technology and Business University, Beijing 100048, China
| | - Ruolin Yang
- College of Life Sciences, Northwest A&F University, Yangling 712100, China.
| | - Xiangfeng Wang
- Department of Crop Genomics and Bioinformatics, College of Agronomy and Biotechnology, National Maize Improvement Center of China, China Agricultural University, Beijing 100094, China.
| |
Collapse
|
12
|
Tracking the Distribution and Burst of Nuclear Mitochondrial DNA Sequences (NUMTs) in Fig Wasp Genomes. INSECTS 2020; 11:insects11100680. [PMID: 33036463 PMCID: PMC7600805 DOI: 10.3390/insects11100680] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 09/04/2020] [Revised: 09/29/2020] [Accepted: 10/06/2020] [Indexed: 11/17/2022]
Abstract
Simple Summary Nuclear mitochondrial DNA sequences (NUMTs), which result from the insertion of exogenous mtDNA into the nuclear genome, are widely distributed in eukaryotes. However, how NUMTs are inserted into the nuclear genome and their post-insertion fates remain a mystery. Previous studies have suggested that Hymenoptera may be a group rich in NUMTs, which will be helpful to study the biological issues of NUMTs. We here select 11 species of fig wasps (Chalcidoidea, Hymenoptera) to analyze the distribution and evolution of NUMTs at the genomic level. The results show that the distributions of NUMTs are species- or lineage-specific. Furthermore, genomic environmental factors such as genome size, the damage-prone regions, and the mode of TE dynamics can determine the insertion and post-insertion fate of NUMTs. Especially because of TEs, the fragmentation and duplication of NUMTs, and thus their burst, are common. This is a relatively comprehensive investigation of the specific distribution of NUMTs and its influencing factors. Our study will help people to understand the evolution of exogenous fragments in the nuclear genome. Abstract Mitochondrial DNA sequences can be transferred into the nuclear genome, giving rise to nuclear mitochondrial DNA sequences (NUMTs). NUMTs have been described in numerous eukaryotes. However, the studies on the distribution of NUMTs and its influencing factors are still inadequate and even controversial. Previous studies have suggested that Hymenoptera may be a group rich in NUMTs, in which we selected 11 species of fig wasps (Chalcidoidea, Hymenoptera) to analyze the distribution and evolution of NUMTs at the genomic level. The results showed that the contents of NUMTs varied greatly in these species, and bursts of NUMTs existed in some species or lineages. Further detailed analyses showed that the large number of NUMTs might be related to the large genomes; NUMTs tended to be inserted into unstable regions of the genomes; and the inserted NUMTs might also be affected by transposable elements (TEs) in the neighbors, leading to fragmentations and duplications, followed by bursts of NUMTs. In summary, our results suggest that a variety of genomic environmental factors can determine the insertion and post-insertion fate of NUMTs, resulting in their species- or lineage-specific distribution patterns, and that studying the evolution of NUMTs can provide good evidence and theoretical basis for exploring the dynamics of exogenous DNA entering into the nuclear genome.
Collapse
|
13
|
Abbate L, Mercati F, Di Noto G, Heuertz M, Carimi F, Fatta del Bosco S, Schicchi R. Genetic Distinctiveness Highlights the Conservation Value of a Sicilian Manna Ash Germplasm Collection Assigned to Fraxinus angustifolia (Oleaceae). PLANTS 2020; 9:plants9081035. [PMID: 32824084 PMCID: PMC7463994 DOI: 10.3390/plants9081035] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/24/2020] [Revised: 08/10/2020] [Accepted: 08/12/2020] [Indexed: 11/16/2022]
Abstract
The cosmopolitan genus Fraxinus comprises about 40 species occupying several habitats in the Northern Hemisphere. With some species hybridizing and sharing genetic variants, questions remain on the species assignment of germplasm within the genus Fraxinus despite numerous species-specific assessments. A multidisciplinary approach was employed to provide a definitive insight into the genetics of an endangered Fraxinus “manna ash” collection, located in a rich plant biodiversity hotspot of the Madonie Mountains (Sicily). Although the collection size was small, genetic diversity, assessed by chloroplast (cpSSR) and nuclear (nSSR) microsatellites (SSR—Simple Sequence Repeats), allowed identifying three different chloroplast haplotypes, with one (H5) dominant, and several polymorphic loci, able to discriminate most of the local accessions studied. Molecular data were linked to cytofluorimetric and phenotypic evaluations and, contrary to popular belief that manna ash is Fraxinus ornus L., the germplasm currently used for manna production belongs to Fraxinus angustifolia Vahl. Interestingly, joint analysis of our genetic panel with a large European dataset of Fraxinus spp. suggested the presence of a possible glacial refuge in Sicily, confirming its importance as biodiversity source. Our results will be helpful for the design of long-term conservation programs for genetic resources, such as in situ and ex situ conservation, seed collection and tree reintroduction.
Collapse
Affiliation(s)
- Loredana Abbate
- Institute of Biosciences and Bioresources (IBBR), National Research Council, Corso Calatafimi 414, 90129 Palermo, Italy; (L.A.); (F.C.); (S.F.d.B.)
| | - Francesco Mercati
- Institute of Biosciences and Bioresources (IBBR), National Research Council, Corso Calatafimi 414, 90129 Palermo, Italy; (L.A.); (F.C.); (S.F.d.B.)
- Correspondence: ; Tel.: +39-091-6574578
| | - Giuseppe Di Noto
- Department of Agricultural, Food and Forestry Sciences (SAAF), University of Palermo, Via Archirafi 38, 90123 Palermo, Italy; (G.D.N.); (R.S.)
| | - Myriam Heuertz
- Institut National de Recherche Pour l’agriculture, l’alimentation et l’environnement (INRAE), Univ. Bordeaux, BIOGECO, 69 route d’Arcachon, F-33610 Cestas, France;
| | - Francesco Carimi
- Institute of Biosciences and Bioresources (IBBR), National Research Council, Corso Calatafimi 414, 90129 Palermo, Italy; (L.A.); (F.C.); (S.F.d.B.)
| | - Sergio Fatta del Bosco
- Institute of Biosciences and Bioresources (IBBR), National Research Council, Corso Calatafimi 414, 90129 Palermo, Italy; (L.A.); (F.C.); (S.F.d.B.)
| | - Rosario Schicchi
- Department of Agricultural, Food and Forestry Sciences (SAAF), University of Palermo, Via Archirafi 38, 90123 Palermo, Italy; (G.D.N.); (R.S.)
| |
Collapse
|
14
|
Hjelmen CE, Parrott JJ, Srivastav SP, McGuane AS, Ellis LL, Stewart AD, Johnston JS, Tarone AM. Effect of Phenotype Selection on Genome Size Variation in Two Species of Diptera. Genes (Basel) 2020; 11:genes11020218. [PMID: 32093067 PMCID: PMC7074110 DOI: 10.3390/genes11020218] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/30/2020] [Revised: 02/11/2020] [Accepted: 02/15/2020] [Indexed: 11/16/2022] Open
Abstract
Genome size varies widely across organisms yet has not been found to be related to organismal complexity in eukaryotes. While there is no evidence for a relationship with complexity, there is evidence to suggest that other phenotypic characteristics, such as nucleus size and cell-cycle time, are associated with genome size, body size, and development rate. However, what is unknown is how the selection for divergent phenotypic traits may indirectly affect genome size. Drosophila melanogaster were selected for small and large body size for up to 220 generations, while Cochliomyia macellaria were selected for 32 generations for fast and slow development. Size in D. melanogaster significantly changed in terms of both cell-count and genome size in isolines, but only the cell-count changed in lines which were maintained at larger effective population sizes. Larger genome sizes only occurred in a subset of D. melanogaster isolines originated from flies selected for their large body size. Selection for development time did not change average genome size yet decreased the within-population variation in genome size with increasing generations of selection. This decrease in variation and convergence on a similar mean genome size was not in correspondence with phenotypic variation and suggests stabilizing selection on genome size in laboratory conditions.
Collapse
Affiliation(s)
- Carl E. Hjelmen
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
- Department of Biology, Texas A&M University, College Station, TX 77843, USA
- Correspondence: or
| | - Jonathan J. Parrott
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
- School of Mathematical and Natural Sciences, Arizona State University, Glendale, AZ 85306, USA
| | - Satyam P. Srivastav
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY 14853, USA
| | - Alexander S. McGuane
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
- Harris County Institute of Forensic Sciences, 1861 Old Spanish Trail, Houston, TX 77054, USA
| | - Lisa L. Ellis
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
- Department of Biology, Houston Baptist University, Houston, TX 77074, USA
| | | | - J. Spencer Johnston
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
| | - Aaron M. Tarone
- Department of Entomology, Texas A&M University, College Station, TX 77843, USA; (J.J.P.); (S.P.S.); (A.S.M.); (L.L.E.); (J.S.J.); (A.M.T.)
| |
Collapse
|
15
|
Carrubba A, Abbate L, Sarno M, Sunseri F, Mauceri A, Lupini A, Mercati F. Characterization of Sicilian rosemary (Rosmarinus officinalis L.) germplasm through a multidisciplinary approach. PLANTA 2020; 251:37. [PMID: 31907671 DOI: 10.1007/s00425-019-03327-8] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/23/2019] [Accepted: 12/09/2019] [Indexed: 06/10/2023]
Abstract
In Sicily, small differences exist between wild and cultivated rosemary biotypes; VOCs and genetic profiles may be a useful tool to distinguish them. A germplasm collection of Rosmarinus officinalis L. was harvested from 15 locations in Sicily. Eleven wild and four cultivated populations were collected and, due to the surveyed area covered, they can be considered as a representative panel of Sicilian genetic background of the species. Ex situ plant collection was transferred to the field cultivation in homogeneous conditions for characterizing through a multidisciplinary approach. The study included morphological traits observations (growth habitus, flower color, number and size of leaves, length and number of internodes), VOC profiles using HS-SPME, genome size by flow cytometry analysis, and genetic characterization by means of DNA and nuclear microsatellite (nSSR) investigation. To detect any pattern within- and among-populations variability, all morphological and chemical data were submitted to ANOVA, while clustering and structure population analysis were carried out using genetic profiles. The present work allowed us to distinguish rather well between wild and cultivated genotypes and to underline the biodiversity richness among rosemary Sicilian germplasm, never highlighted, useful for future breeding programs addressed to exploit this important resource.
Collapse
Affiliation(s)
- Alessandra Carrubba
- Department of Agriculture, Food and Forest Sciences, University of Palermo, Palermo, Italy.
| | - Loredana Abbate
- Institute of Biosciences and Bioresources (IBBR), National Research Council, Palermo, Italy
| | - Mauro Sarno
- Department of Agriculture, Food and Forest Sciences, University of Palermo, Palermo, Italy
| | - Francesco Sunseri
- Dipartimento AGRARIA, Località Feo di Vito snc, 89121, Reggio Calabria, Italy
| | - Antonio Mauceri
- Dipartimento AGRARIA, Località Feo di Vito snc, 89121, Reggio Calabria, Italy
| | - Antonio Lupini
- Dipartimento AGRARIA, Località Feo di Vito snc, 89121, Reggio Calabria, Italy
| | - Francesco Mercati
- Institute of Biosciences and Bioresources (IBBR), National Research Council, Palermo, Italy
| |
Collapse
|
16
|
Legrand S, Caron T, Maumus F, Schvartzman S, Quadrana L, Durand E, Gallina S, Pauwels M, Mazoyer C, Huyghe L, Colot V, Hanikenne M, Castric V. Differential retention of transposable element-derived sequences in outcrossing Arabidopsis genomes. Mob DNA 2019; 10:30. [PMID: 31346350 PMCID: PMC6636163 DOI: 10.1186/s13100-019-0171-6] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/02/2019] [Accepted: 06/28/2019] [Indexed: 12/20/2022] Open
Abstract
Background Transposable elements (TEs) are genomic parasites with major impacts on host genome architecture and host adaptation. A proper evaluation of their evolutionary significance has been hampered by the paucity of short scale phylogenetic comparisons between closely related species. Here, we characterized the dynamics of TE accumulation at the micro-evolutionary scale by comparing two closely related plant species, Arabidopsis lyrata and A. halleri. Results Joint genome annotation in these two outcrossing species confirmed that both contain two distinct populations of TEs with either 'recent' or 'old' insertion histories. Identification of rare segregating insertions suggests that diverse TE families contribute to the ongoing dynamics of TE accumulation in the two species. Orthologous TE fragments (i.e. those that have been maintained in both species), tend to be located closer to genes than those that are retained in one species only. Compared to non-orthologous TE insertions, those that are orthologous tend to produce fewer short interfering RNAs, are less heavily methylated when found within or adjacent to genes and these tend to have lower expression levels. These findings suggest that long-term retention of TE insertions reflects their frequent acquisition of adaptive roles and/or the deleterious effects of removing nearly neutral TE insertions when they are close to genes. Conclusion Our results indicate a rapid evolutionary dynamics of the TE landscape in these two outcrossing species, with an important input of a diverse set of new insertions with variable propensity to resist deletion.
Collapse
Affiliation(s)
- Sylvain Legrand
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Thibault Caron
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Florian Maumus
- 2URGI, INRA, Université Paris-Saclay, 78026 Versailles, France
| | - Sol Schvartzman
- 3InBioS - PhytoSystems, Functional Genomics and Plant Molecular Imaging, University of Liège, 4000 Liège, Belgium
| | - Leandro Quadrana
- 4IBENS, Département de Biologie, Ecole Normale Supérieure, CNRS, Inserm, PSL Research University, F-75005 Paris, France
| | - Eléonore Durand
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Sophie Gallina
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Maxime Pauwels
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Clément Mazoyer
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Lucie Huyghe
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| | - Vincent Colot
- 4IBENS, Département de Biologie, Ecole Normale Supérieure, CNRS, Inserm, PSL Research University, F-75005 Paris, France
| | - Marc Hanikenne
- 3InBioS - PhytoSystems, Functional Genomics and Plant Molecular Imaging, University of Liège, 4000 Liège, Belgium
| | - Vincent Castric
- 1Univ. Lille, CNRS, UMR 8198 - Evo-Eco-Paleo, F-59000 Lille, France
| |
Collapse
|
17
|
Rubenstein DR, Ågren JA, Carbone L, Elde NC, Hoekstra HE, Kapheim KM, Keller L, Moreau CS, Toth AL, Yeaman S, Hofmann HA. Coevolution of Genome Architecture and Social Behavior. Trends Ecol Evol 2019; 34:844-855. [PMID: 31130318 DOI: 10.1016/j.tree.2019.04.011] [Citation(s) in RCA: 22] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/05/2019] [Revised: 04/03/2019] [Accepted: 04/17/2019] [Indexed: 01/02/2023]
Abstract
Although social behavior can have a strong genetic component, it can also result in selection on genome structure and function, thereby influencing the evolution of the genome itself. Here we explore the bidirectional links between social behavior and genome architecture by considering variation in social and/or mating behavior among populations (social polymorphisms) and across closely related species. We propose that social behavior can influence genome architecture via associated demographic changes due to social living. We establish guidelines to exploit emerging whole-genome sequences using analytical approaches that examine genome structure and function at different levels (regulatory vs structural variation) from the perspective of both molecular biology and population genetics in an ecological context.
Collapse
Affiliation(s)
- Dustin R Rubenstein
- Columbia University, Department of Ecology, Evolution, and Environmental Biology and Center for Integrative Animal Behavior, New York, NY 10027, USA.
| | - J Arvid Ågren
- Harvard University, Department of Organismic and Evolutionary Biology, Cambridge, MA 02138, USA
| | - Lucia Carbone
- Oregon Health & Science University, Department of Medicine, KCVI, Portland, OR 97239, USA; Oregon National Primate Research Center, Division of Genetics, Beaverton, OR 97006, USA
| | - Nels C Elde
- University of Utah School of Medicine, Department of Human Genetics, Salt Lake City, UT 84112, USA
| | - Hopi E Hoekstra
- Harvard University, Department of Organismic and Evolutionary Biology, Cambridge, MA 02138, USA; Harvard University, Howard Hughes Medical Institute, Department of Molecular and Cellular Biology, Cambridge, MA 02138, USA
| | - Karen M Kapheim
- Utah State University, Department of Biology, Logan, UT 84322, USA
| | - Laurent Keller
- University of Lausanne, Department of Ecology and Evolution, Biophore, UNIL, 1015 Lausanne, Switzerland
| | - Corrie S Moreau
- Cornell University, Departments of Entomology and Ecology and Evolutionary Biology, Ithaca, NY 14850, USA
| | - Amy L Toth
- Iowa State University, Department of Ecology, Evolution, and Organismal Biology and Department of Entomology, Ames, IA 50011, USA
| | - Sam Yeaman
- University of Calgary, Department of Biological Sciences, Calgary, AB T2N 1N4, Canada
| | - Hans A Hofmann
- The University of Texas at Austin, Department of Integrative Biology and Institute for Cellular and Molecular Biology, 2415 Speedway C-0990, Austin, TX 78712, USA.
| |
Collapse
|
18
|
|
19
|
Hjelmen CE, Garrett MA, Holmes VR, Mynes M, Piron E, Johnston JS. Genome Size Evolution within and between the Sexes. J Hered 2019; 110:219-228. [PMID: 30476187 DOI: 10.1093/jhered/esy063] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/26/2018] [Accepted: 11/21/2018] [Indexed: 02/02/2023] Open
Abstract
Genome sizes are known to vary between closely related species, but the patterns behind this variation have yet to be fully understood. Although this variation has been evaluated between species and within sexes, unknown is the extent to which this variation is driven by differentiation in sex chromosomes. To address this longstanding question, we examine the mode and tempo of genome size evolution for a total of 87 species of Drosophilidae, estimating and updating male genome size values for 44 of these species. We compare the evolution of genome size within each sex to the evolution of the differences between the sexes. Utilizing comparative phylogenetic methods, we find that male and female genome size evolution is largely a neutral process, reflective of phylogenetic relatedness between species, which supports the newly proposed accordion model for genome size change. When similarly analyzed, the difference between the sexes due to heteromorphic sex chromosomes is a dynamic process; the male-female genome size difference increases with time with or without known neo-Y events or complete loss of the Y. Observed instances of rapid change match theoretical expectations and known neo-Y and Y loss events in individual species.
Collapse
Affiliation(s)
- Carl E Hjelmen
- Department of Entomology, Texas A&M University, College Station, TX
| | - Margaret A Garrett
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, TX
| | - V Renee Holmes
- Department of Entomology, Texas A&M University, College Station, TX
| | - Melissa Mynes
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, TX
| | - Elizabeth Piron
- Department of Biochemistry and Biophysics, Texas A&M University, College Station, TX
| | | |
Collapse
|
20
|
Blankers T, Oh KP, Bombarely A, Shaw KL. The Genomic Architecture of a Rapid Island Radiation: Recombination Rate Variation, Chromosome Structure, and Genome Assembly of the Hawaiian Cricket Laupala. Genetics 2018; 209:1329-1344. [PMID: 29875253 PMCID: PMC6063224 DOI: 10.1534/genetics.118.300894] [Citation(s) in RCA: 27] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/07/2018] [Accepted: 06/03/2018] [Indexed: 12/30/2022] Open
Abstract
Phenotypic evolution and speciation depend on recombination in many ways. Within populations, recombination can promote adaptation by bringing together favorable mutations and decoupling beneficial and deleterious alleles. As populations diverge, crossing over can give rise to maladapted recombinants and impede or reverse diversification. Suppressed recombination due to genomic rearrangements, modifier alleles, and intrinsic chromosomal properties may offer a shield against maladaptive gene flow eroding coadapted gene complexes. Both theoretical and empirical results support this relationship. However, little is known about this relationship in the context of behavioral isolation, where coevolving signals and preferences are the major hybridization barrier. Here we examine the genomic architecture of recently diverged, sexually isolated Hawaiian swordtail crickets (Laupala). We assemble a de novo genome and generate three dense linkage maps from interspecies crosses. In line with expectations based on the species' recent divergence and successful interbreeding in the laboratory, the linkage maps are highly collinear and show no evidence for large-scale chromosomal rearrangements. Next, the maps were used to anchor the assembly to pseudomolecules and estimate recombination rates across the genome to test the hypothesis that loci involved in behavioral isolation (song and preference divergence) are in regions of low interspecific recombination. Contrary to our expectations, the genomic region where a male song and female preference QTL colocalize is not associated with particularly low recombination rates. This study provides important novel genomic resources for an emerging evolutionary genetics model system and suggests that trait-preference coevolution is not necessarily facilitated by locally suppressed recombination.
Collapse
Affiliation(s)
- Thomas Blankers
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853
| | - Kevin P Oh
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853
| | - Aureliano Bombarely
- Department of Horticulture, Virginia Polytechnic Institute and State University, Blacksburg, Virginia 24061
| | - Kerry L Shaw
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853
| |
Collapse
|
21
|
Robertson HM, Baits RL, Walden KK, Wada‐Katsumata A, Schal C. Enormous expansion of the chemosensory gene repertoire in the omnivorous German cockroach Blattella germanica. JOURNAL OF EXPERIMENTAL ZOOLOGY. PART B, MOLECULAR AND DEVELOPMENTAL EVOLUTION 2018; 330:265-278. [PMID: 29566459 PMCID: PMC6175461 DOI: 10.1002/jez.b.22797] [Citation(s) in RCA: 52] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 01/23/2018] [Accepted: 02/27/2018] [Indexed: 01/26/2023]
Abstract
The acquisition of genome sequences from a wide range of insects and other arthropods has revealed a broad positive correlation between the complexity of their chemical ecology and the size of their chemosensory gene repertoire. The German cockroach Blattella germanica is an extreme omnivore and has the largest chemosensory gene repertoire known for an arthropod, exceeding even the highly polyphagous spider mite Tetranychus urticae. While the Odorant Receptor family is not particularly large, with 123 genes potentially encoding 134 receptors (105 intact), the Gustatory Receptor family is greatly expanded to 431 genes potentially encoding 545 receptors (483 intact), the largest known for insects and second only to the spider mite. The Ionotropic Receptor family of olfactory and gustatory receptors is vastly expanded to at least 897 genes (604 intact), the largest size known in arthropods, far surpassing the 150 known from the dampwood termite Zootermopsis nevadensis. Commensurately, the Odorant Binding Protein family is expanded to the largest known for insects at 109 genes (all intact). Comparison with the far more specialized, but phylogenetically related termite, within the Dictyoptera, reveals considerable gene losses from the termite, and massive species-specific gene expansions in the cockroach. The cockroach has lost function of 11%-41% of these three chemoreceptor gene families to pseudogenization, and most of these are young events, implying rapid turnover of genes along with these major expansions, presumably in response to changes in its chemical ecology.
Collapse
Affiliation(s)
- Hugh M. Robertson
- Department of EntomologyUniversity of Illinois at Urbana‐ChampaignUrbanaIllinoisUSA
| | - Rachel L. Baits
- Department of EntomologyUniversity of Illinois at Urbana‐ChampaignUrbanaIllinoisUSA
| | - Kimberly K.O. Walden
- Department of EntomologyUniversity of Illinois at Urbana‐ChampaignUrbanaIllinoisUSA
| | - Ayako Wada‐Katsumata
- Department of Entomology and Plant PathologyNorth Carolina State UniversityRaleighNorth CarolinaUSA
| | - Coby Schal
- Department of Entomology and Plant PathologyNorth Carolina State UniversityRaleighNorth CarolinaUSA
| |
Collapse
|
22
|
Bilinski P, Albert PS, Berg JJ, Birchler JA, Grote MN, Lorant A, Quezada J, Swarts K, Yang J, Ross-Ibarra J. Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays. PLoS Genet 2018; 14:e1007162. [PMID: 29746459 PMCID: PMC5944917 DOI: 10.1371/journal.pgen.1007162] [Citation(s) in RCA: 61] [Impact Index Per Article: 10.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2017] [Accepted: 12/20/2017] [Indexed: 12/03/2022] Open
Abstract
While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes. Genome size in plants can vary by orders of magnitude, but this variation has long been considered to be of little functional consequence. Studying three independent adaptations to high altitude in Zea mays, we find that genome size experiences parallel pressures from natural selection, causing a reduction in genome size with increasing altitude. Though reductions in overall repetitive content are responsible for the genome size change, we find that only those individual loci contributing most to the variation in genome size are individually targeted by selection. To identify the phenotype influenced by genome size, we study how variation in genome size within a single wild population impacts leaf growth and cell division. We find that genome size variation correlates negatively with the rate of cell division, suggesting that individuals with larger genomes require longer to complete a mitotic cycle. Finally, we reanalyze data from maize inbreds to show that faster cell division is correlated with earlier flowering, connecting observed variation in genome size to an important adaptive phenotype.
Collapse
Affiliation(s)
- Paul Bilinski
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
- Research Group for Ancient Genomics and Evolution, Department of Molecular Biology, Max Planck Institute for Developmental Biology, Tuebingen, Germany
- * E-mail: (PB); (JRI)
| | - Patrice S. Albert
- Division of Biological Sciences, University of Missouri, Columbia, Missouri, United States of America
| | - Jeremy J. Berg
- Center for Population Biology, University of California, Davis, Davis, California, United States of America
- Department of Evolution and Ecology, University of California, Davis, Davis, California, United States of America
| | - James A. Birchler
- Division of Biological Sciences, University of Missouri, Columbia, Missouri, United States of America
| | - Mark N. Grote
- Department of Anthropology, University of California, Davis, Davis, California, United States of America
| | - Anne Lorant
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
| | - Juvenal Quezada
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
| | - Kelly Swarts
- Research Group for Ancient Genomics and Evolution, Department of Molecular Biology, Max Planck Institute for Developmental Biology, Tuebingen, Germany
| | - Jinliang Yang
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
- Department of Agronomy and Horticulture, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
| | - Jeffrey Ross-Ibarra
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
- Center for Population Biology, University of California, Davis, Davis, California, United States of America
- Genome Center, University of California, Davis, Davis, California, United States of America
- * E-mail: (PB); (JRI)
| |
Collapse
|
23
|
Bilinski P, Albert PS, Berg JJ, Birchler JA, Grote MN, Lorant A, Quezada J, Swarts K, Yang J, Ross-Ibarra J. Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays. PLoS Genet 2018; 14:e1007162. [PMID: 29746459 DOI: 10.1371/journal.pgen.1007162.g001] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/10/2017] [Accepted: 12/20/2017] [Indexed: 05/23/2023] Open
Abstract
While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes.
Collapse
Affiliation(s)
- Paul Bilinski
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
- Research Group for Ancient Genomics and Evolution, Department of Molecular Biology, Max Planck Institute for Developmental Biology, Tuebingen, Germany
| | - Patrice S Albert
- Division of Biological Sciences, University of Missouri, Columbia, Missouri, United States of America
| | - Jeremy J Berg
- Center for Population Biology, University of California, Davis, Davis, California, United States of America
- Department of Evolution and Ecology, University of California, Davis, Davis, California, United States of America
| | - James A Birchler
- Division of Biological Sciences, University of Missouri, Columbia, Missouri, United States of America
| | - Mark N Grote
- Department of Anthropology, University of California, Davis, Davis, California, United States of America
| | - Anne Lorant
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
| | - Juvenal Quezada
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
| | - Kelly Swarts
- Research Group for Ancient Genomics and Evolution, Department of Molecular Biology, Max Planck Institute for Developmental Biology, Tuebingen, Germany
| | - Jinliang Yang
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
- Department of Agronomy and Horticulture, University of Nebraska-Lincoln, Lincoln, Nebraska, United States of America
| | - Jeffrey Ross-Ibarra
- Department of Plant Sciences, University of California, Davis, Davis, California, United States of America
- Center for Population Biology, University of California, Davis, Davis, California, United States of America
- Genome Center, University of California, Davis, Davis, California, United States of America
| |
Collapse
|
24
|
Estimation of universal and taxon-specific parameters of prokaryotic genome evolution. PLoS One 2018; 13:e0195571. [PMID: 29652909 PMCID: PMC5898727 DOI: 10.1371/journal.pone.0195571] [Citation(s) in RCA: 5] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/16/2017] [Accepted: 03/23/2018] [Indexed: 01/13/2023] Open
Abstract
The results of our recent study on mathematical modeling of microbial genome evolution indicate that, on average, genomes of bacteria and archaea evolve in the regime of mutation-selection balance defined by positive selection coefficients associated with gene acquisition that is counter-acted by the intrinsic deletion bias. This analysis was based on the strong assumption that parameters of genome evolution are universal across the diversity of bacteria and archaea, and yielded extremely low values of the selection coefficient. Here we further refine the modeling approach by taking into account evolutionary factors specific for individual groups of microbes using two independent fitting strategies, an ad hoc hard fitting scheme and a mixture model. The resulting estimate of the mean selection coefficient of s∼10−10 associated with the gain of one gene implies that, on average, acquisition of a gene is beneficial, and that microbial genomes typically evolve under a weak selection regime that might transition to strong selection in highly abundant organisms with large effective population sizes. The apparent selective pressure towards larger genomes is balanced by the deletion bias, which is estimated to be consistently greater than unity for all analyzed groups of microbes. The estimated values of s are more realistic than the lower values obtained previously, indicating that global and group-specific evolutionary factors synergistically affect microbial genome evolution that seems to be driven primarily by adaptation to existence in diverse niches.
Collapse
|
25
|
Turnell BR, Shaw KL, Reeve HK. Modeling strategic sperm allocation: Tailoring the predictions to the species. Evolution 2018; 72:414-425. [PMID: 29331038 DOI: 10.1111/evo.13423] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/23/2016] [Revised: 12/06/2017] [Accepted: 12/14/2017] [Indexed: 12/01/2022]
Abstract
Two major challenges exist when empirically testing the predictions of sperm allocation theory. First, the study species must adhere to the assumptions of the model being tested. Unfortunately, the common assumption of sperm allocation models that females mate a maximum of once or twice does not hold for many, if not most, multiply and sequentially mating animals. Second, a model's parameters, which dictate its predictions, must be measured in the study species. Common examples of such parameters, female mating frequency and sperm precedence patterns, are unknown for many species used in empirical tests. Here, we present a broadly applicable model, appropriate for multiply, sequentially mating animals, and test it in three species for which data on all the relevant parameter values are available. The model predicts that relative allocation to virgin females, compared to nonvirgins, depends on the interaction between female mating rate and the sperm precedence pattern: relative allocation to virgins increases with female mating rate under first-male precedence, while the opposite is true under later-male precedence. Our model is moderately successful in predicting actual allocation patterns in the three species, including a cricket in which we measured the parameter values and performed an empirical test of allocation.
Collapse
Affiliation(s)
- Biz R Turnell
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853.,Current Address: Department of Biology, Technische Universität Dresden, 01062 Dresden, Germany
| | - Kerry L Shaw
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853
| | - H Kern Reeve
- Department of Neurobiology and Behavior, Cornell University, Ithaca, New York 14853
| |
Collapse
|
26
|
Chong RA, Mueller RL. Polymorphic duplicate genes and persistent non-coding sequences reveal heterogeneous patterns of mitochondrial DNA loss in salamanders. BMC Genomics 2017; 18:992. [PMID: 29281973 PMCID: PMC5745709 DOI: 10.1186/s12864-017-4358-2] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/24/2016] [Accepted: 11/29/2017] [Indexed: 12/03/2022] Open
Abstract
Background Mitochondria are the site of the citric acid cycle and oxidative phosphorylation (OXPHOS). In metazoans, the mitochondrial genome is a small, circular molecule averaging 16.5 kb in length. Despite evolutionarily conserved gene content, metazoan mitochondrial genomes show a diversity of gene orders most commonly explained by the duplication-random loss (DRL) model. In the DRL model, (1) a sequence of genes is duplicated in tandem, (2) one paralog sustains a loss-of-function mutation, resulting in selection to retain the other copy, and (3) the non-functional paralog is eventually deleted from the genome. Despite its apparent role in generating mitochondrial gene order diversity, little is known about the tempo and mode of random gene loss after duplication events. Here, we determine mitochondrial gene order across the salamander genus Aneides, which was previously shown to include at least two DRL-mediated rearrangement events. We then analyze these gene orders in a phylogenetic context to reveal patterns of DNA loss after mitochondrial gene duplication. Results We identified two separate duplication events that resulted in mitochondrial gene rearrangements in Aneides; one occurred at the base of the clade tens of millions of years ago, while the other occurred much more recently (i.e. within a single species), resulting in gene order polymorphism and paralogs that are readily identifiable. We demonstrate that near-complete removal of duplicate rRNA genes has occurred since the recent duplication event, whereas duplicate protein-coding genes persist as pseudogenes and duplicate tRNAs persist as functionally intact paralogs. In addition, we show that non-coding DNA duplicated at the base of the clade has persisted across species for tens of millions of years. Conclusions The evolutionary history of the mitochondrial genome, from its inception as a bacterial endosymbiont, includes massive genomic reduction. Consistent with this overall trend, selection for efficiency of mitochondrial replication and transcription has been hypothesized to favor elimination of extra sequence. Our results, however, suggest that there may be no strong disadvantage to extraneous sequences in salamander mitochondrial genomes, although duplicate rRNA genes may be deleterious. Electronic supplementary material The online version of this article (10.1186/s12864-017-4358-2) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Rebecca A Chong
- Department of Biology, Colorado State University, Fort Collins, CO, 80523-1878, USA. .,Present address: Department of Biology, University of Hawai'i at Mānoa, Honolulu, HI, 96822, USA.
| | | |
Collapse
|
27
|
Pai AA, Henriques T, McCue K, Burkholder A, Adelman K, Burge CB. The kinetics of pre-mRNA splicing in the Drosophila genome and the influence of gene architecture. eLife 2017; 6:32537. [PMID: 29280736 PMCID: PMC5762160 DOI: 10.7554/elife.32537] [Citation(s) in RCA: 44] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/05/2017] [Accepted: 12/22/2017] [Indexed: 12/28/2022] Open
Abstract
Production of most eukaryotic mRNAs requires splicing of introns from pre-mRNA. The splicing reaction requires definition of splice sites, which are initially recognized in either intron-spanning (‘intron definition’) or exon-spanning (‘exon definition’) pairs. To understand how exon and intron length and splice site recognition mode impact splicing, we measured splicing rates genome-wide in Drosophila, using metabolic labeling/RNA sequencing and new mathematical models to estimate rates. We found that the modal intron length range of 60–70 nt represents a local maximum of splicing rates, but that much longer exon-defined introns are spliced even faster and more accurately. We observed unexpectedly low variation in splicing rates across introns in the same gene, suggesting the presence of gene-level influences, and we identified multiple gene level variables associated with splicing rate. Together our data suggest that developmental and stress response genes may have preferentially evolved exon definition in order to enhance the rate or accuracy of splicing.
Collapse
Affiliation(s)
- Athma A Pai
- Departments of Biology and Biological Engineering, Massachusetts Institute of Technology, Cambridge, United States
| | - Telmo Henriques
- Epigenetics and Stem Cell Biology Laboratory, National Institute of Environmental Health Sciences, Research Triangle, United States.,Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
| | - Kayla McCue
- Program in Computational and Systems Biology, Massachusetts Institute of Technology, Cambridge, United States
| | - Adam Burkholder
- Center for Integrative Bioinformatics, National Institute of Environmental Health Sciences, Research Triangle, United States
| | - Karen Adelman
- Epigenetics and Stem Cell Biology Laboratory, National Institute of Environmental Health Sciences, Research Triangle, United States.,Department of Biological Chemistry and Molecular Pharmacology, Harvard Medical School, Boston, United States
| | - Christopher B Burge
- Departments of Biology and Biological Engineering, Massachusetts Institute of Technology, Cambridge, United States.,Program in Computational and Systems Biology, Massachusetts Institute of Technology, Cambridge, United States
| |
Collapse
|
28
|
An Important Role for Purifying Selection in Archaeal Genome Evolution. mSystems 2017; 2:mSystems00112-17. [PMID: 29085915 PMCID: PMC5655593 DOI: 10.1128/msystems.00112-17] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/22/2017] [Accepted: 10/06/2017] [Indexed: 12/20/2022] Open
Abstract
As the null hypothesis of genome evolution, population genetic theory suggests that selection strength controls genome size. Through the process of genetic drift, this theory predicts that compact genomes are maintained by strong purifying selection while complex genomes are enabled by weak purifying selection. It offers a unifying framework that explains why prokaryotic genomes are much smaller than their eukaryotic counterparts. However, recent findings suggest that bigger prokaryotic genomes appear to experience stronger purifying selection, indicating that purifying selection may not dominate prokaryotic genome evolution. Since archaeal genomes were underrepresented in those studies, generalization of the conclusions to both archaeal and bacterial genomes may not be warranted. In this study, we revisited this matter by focusing on archaeal and bacterial genomes separately. We found that bigger bacterial genomes indeed experienced stronger purifying selection, but the opposite was observed in archaeal genomes. This new finding would predict an enrichment of noncoding sequences in large archaeal genomes, which was confirmed by an analysis of coding density. In contrast, coding density remained stable regardless of bacterial genome size. In conclusion, this study suggests that purifying selection may play a more important role in archaeal genome evolution than previously hypothesized, indicating that there could be a major difference between the evolutionary regimes of Archaea and Bacteria. IMPORTANCE The evolution of genome complexity is a fundamental question in biology. A hallmark of eukaryotic genome complexity is that larger genomes tend to have more noncoding sequences, which are believed to be minimal in archaeal and bacterial genomes. However, we found that archaeal genomes also possessed this eukaryotic feature while bacterial genomes did not. This could be predicted from our analysis on genetic drift, which showed a relaxation of purifying selection in larger archaeal genomes, also a eukaryotic feature. In contrast, the opposite was evident in bacterial genomes.
Collapse
|
29
|
Rodrigo G. Evolutionary impact of copy number variation rates. BMC Res Notes 2017; 10:393. [PMID: 28797272 PMCID: PMC5553737 DOI: 10.1186/s13104-017-2741-3] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2017] [Accepted: 08/08/2017] [Indexed: 02/17/2023] Open
Abstract
Objective Copy number variation is now recognized as one of the major sources of genetic variation among individuals in natural populations of any species. However, the relevance of these unexpected observations goes beyond diagnosing high diversity. Results Here, it is argued that the molecular rates of copy number variation, mainly the deletion rate upon variation, determine the evolutionary road of the genome regarding size. Genetic drift will govern this process only if the effective population size is lower than the inverse of the deletion rate. Otherwise, natural selection will do.
Collapse
Affiliation(s)
- Guillermo Rodrigo
- Institute for Integrative Systems Biology, CSIC-UV, 46980, Paterna, Spain.
| |
Collapse
|
30
|
Rapid and Recent Evolution of LTR Retrotransposons Drives Rice Genome Evolution During the Speciation of AA-Genome Oryza Species. G3-GENES GENOMES GENETICS 2017; 7:1875-1885. [PMID: 28413161 PMCID: PMC5473765 DOI: 10.1534/g3.116.037572] [Citation(s) in RCA: 38] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/20/2022]
Abstract
The dynamics of long terminal repeat (LTR) retrotransposons and their contribution to genome evolution during plant speciation have remained largely unanswered. Here, we perform a genome-wide comparison of all eight Oryza AA-genome species, and identify 3911 intact LTR retrotransposons classified into 790 families. The top 44 most abundant LTR retrotransposon families show patterns of rapid and distinct diversification since the species split over the last ∼4.8 MY (million years). Phylogenetic and read depth analyses of 11 representative retrotransposon families further provide a comprehensive evolutionary landscape of these changes. Compared with Ty1-copia, independent bursts of Ty3-gypsy retrotransposon expansions have occurred with the three largest showing signatures of lineage-specific evolution. The estimated insertion times of 2213 complete retrotransposons from the top 23 most abundant families reveal divergent life histories marked by speedy accumulation, decline, and extinction that differed radically between species. We hypothesize that this rapid evolution of LTR retrotransposons not only divergently shaped the architecture of rice genomes but also contributed to the process of speciation and diversification of rice.
Collapse
|
31
|
Vu GTH, Cao HX, Reiss B, Schubert I. Deletion-bias in DNA double-strand break repair differentially contributes to plant genome shrinkage. THE NEW PHYTOLOGIST 2017; 214:1712-1721. [PMID: 28245065 DOI: 10.1111/nph.14490] [Citation(s) in RCA: 23] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/23/2016] [Accepted: 01/22/2017] [Indexed: 06/06/2023]
Abstract
In order to prevent genome instability, cells need to be protected by a number of repair mechanisms, including DNA double-strand break (DSB) repair. The extent to which DSB repair, biased towards deletions or insertions, contributes to evolutionary diversification of genome size is still under debate. We analyzed mutation spectra in Arabidopsis thaliana and in barley (Hordeum vulgare) by PacBio sequencing of three DSB-targeted loci each, uncovering repair via gene conversion, single strand annealing (SSA) or nonhomologous end-joining (NHEJ). Furthermore, phylogenomic comparisons between A. thaliana and two related species were used to detect naturally occurring deletions during Arabidopsis evolution. Arabidopsis thaliana revealed significantly more and larger deletions after DSB repair than barley, and barley displayed more and larger insertions. Arabidopsis displayed a clear net loss of DNA after DSB repair, mainly via SSA and NHEJ. Barley revealed a very weak net loss of DNA, apparently due to less active break-end resection and easier copying of template sequences into breaks. Comparative phylogenomics revealed several footprints of SSA in the A. thaliana genome. Quantitative assessment of DNA gain and loss through DSB repair processes suggests deletion-biased DSB repair causing ongoing genome shrinking in A. thaliana, whereas genome size in barley remains nearly constant.
Collapse
Affiliation(s)
- Giang T H Vu
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), D-06466, Gatersleben, Germany
| | - Hieu X Cao
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), D-06466, Gatersleben, Germany
| | - Bernd Reiss
- Max Planck Institute for Plant Breeding Research, 50829, Köln, Germany
| | - Ingo Schubert
- Leibniz Institute of Plant Genetics and Crop Plant Research (IPK), D-06466, Gatersleben, Germany
| |
Collapse
|
32
|
Xiao-Ming Z, Junrui W, Li F, Sha L, Hongbo P, Lan Q, Jing L, Yan S, Weihua Q, Lifang Z, Yunlian C, Qingwen Y. Inferring the evolutionary mechanism of the chloroplast genome size by comparing whole-chloroplast genome sequences in seed plants. Sci Rep 2017; 7:1555. [PMID: 28484234 PMCID: PMC5431534 DOI: 10.1038/s41598-017-01518-5] [Citation(s) in RCA: 62] [Impact Index Per Article: 8.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/17/2016] [Accepted: 03/29/2017] [Indexed: 01/08/2023] Open
Abstract
The chloroplast genome originated from photosynthetic organisms and has retained the core genes that mainly encode components of photosynthesis. However, the causes of variations in chloroplast genome size in seed plants have only been thoroughly analyzed within small subsets of spermatophytes. In this study, we conducted the first comparative analysis on a large scale to examine the relationship between sequence characteristics and genome size in 272 seed plants based on cross-species and phylogenetic signal analysis. Our results showed that inverted repeat regions, large or small single copies, intergenic regions, and gene number can be attributed to the variations in chloroplast genome size among closely related species. However, chloroplast gene length underwent evolution affecting chloroplast genome size in seed plants irrespective of whether phylogenetic information was incorporated. Among chloroplast genes, atpA, accD and ycf1 account for 13% of the variation in genome size, and the average Ka/Ks values of homologous pairs of the three genes are larger than 1. The relationship between chloroplast genome size and gene length might be affected by selection during the evolution of spermatophytes. The variation in chloroplast genome size may influence energy generation and ecological strategy in seed plants.
Collapse
Affiliation(s)
- Zheng Xiao-Ming
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Wang Junrui
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Feng Li
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Liu Sha
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Pang Hongbo
- College of Chemistry and Life Science, Shenyang Normal University, Shenyang, 110034, China
| | - Qi Lan
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Li Jing
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Sun Yan
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Qiao Weihua
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Zhang Lifang
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Cheng Yunlian
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China
| | - Yang Qingwen
- National Key Facility for Crop Gene Resources and Genetic Improvement, Institute of Crop Sciences, Chinese Academy of Agricultural Sciences, Beijing, 100081, China.
| |
Collapse
|
33
|
Hjelmen CE, Johnston JS. The mode and tempo of genome size evolution in the subgenus Sophophora. PLoS One 2017; 12:e0173505. [PMID: 28267812 PMCID: PMC5340367 DOI: 10.1371/journal.pone.0173505] [Citation(s) in RCA: 10] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/09/2016] [Accepted: 02/21/2017] [Indexed: 01/05/2023] Open
Abstract
Genome size varies widely across organisms, with no apparent tie to organismal complexity. While genome size is inherited, there is no established evolutionary model for this trait. Hypotheses have been postulated for the observed variation in genome sizes across species, most notably the effective population size hypothesis, the mutational equilibrium hypothesis, and the adaptive hypothesis. While much data has been collected on genome size, the above hypotheses have largely ignored impacts from phylogenetic relationships. In order to test these competing hypotheses, genome sizes of 87 Sophophora species were analyzed in a comparative phylogenetic approach using Pagel’s parameters of evolution, Blomberg’s K, Abouheif’s Cmean and Moran’s I. In addition to testing the mode and rate of genome size evolution in Sophophora species, the effect of number of taxa on detection of phylogenetic signal was analyzed for each of these comparative phylogenetic methods. Sophophora genome size was found to be dependent on the phylogeny, indicating that evolutionary time was important for predicting the variation among species. Genome size was found to evolve gradually on branches of the tree, with a rapid burst of change early in the phylogeny. These results suggest that Sophophora genome size has experienced gradual changes, which support the largely theoretical mutational equilibrium hypothesis. While some methods (Abouheif’s Cmean and Moran’s I) were found to be affected by increasing taxa numbers, more commonly used methods (λ and Blomberg’s K) were found to have increasing reliability with increasing taxa number, with significantly more support with fifteen or more taxa. Our results suggest that these comparative phylogenetic methods, with adequate taxon sampling, can be a powerful way to uncover the enigma that is genome size variation through incorporation of phylogenetic relationships.
Collapse
Affiliation(s)
- Carl E. Hjelmen
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
- * E-mail:
| | - J. Spencer Johnston
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
| |
Collapse
|
34
|
Abstract
Genome size in mammals and birds shows remarkably little interspecific variation compared with other taxa. However, genome sequencing has revealed that many mammal and bird lineages have experienced differential rates of transposable element (TE) accumulation, which would be predicted to cause substantial variation in genome size between species. Thus, we hypothesize that there has been covariation between the amount of DNA gained by transposition and lost by deletion during mammal and avian evolution, resulting in genome size equilibrium. To test this model, we develop computational methods to quantify the amount of DNA gained by TE expansion and lost by deletion over the last 100 My in the lineages of 10 species of eutherian mammals and 24 species of birds. The results reveal extensive variation in the amount of DNA gained via lineage-specific transposition, but that DNA loss counteracted this expansion to various extents across lineages. Our analysis of the rate and size spectrum of deletion events implies that DNA removal in both mammals and birds has proceeded mostly through large segmental deletions (>10 kb). These findings support a unified "accordion" model of genome size evolution in eukaryotes whereby DNA loss counteracting TE expansion is a major determinant of genome size. Furthermore, we propose that extensive DNA loss, and not necessarily a dearth of TE activity, has been the primary force maintaining the greater genomic compaction of flying birds and bats relative to their flightless relatives.
Collapse
|
35
|
Abstract
Bacteria and archaea typically possess small genomes that are tightly packed with protein-coding genes. The compactness of prokaryotic genomes is commonly perceived as evidence of adaptive genome streamlining caused by strong purifying selection in large microbial populations. In such populations, even the small cost incurred by nonfunctional DNA because of extra energy and time expenditure is thought to be sufficient for this extra genetic material to be eliminated by selection. However, contrary to the predictions of this model, there exists a consistent, positive correlation between the strength of selection at the protein sequence level, measured as the ratio of nonsynonymous to synonymous substitution rates, and microbial genome size. Here, by fitting the genome size distributions in multiple groups of prokaryotes to predictions of mathematical models of population evolution, we show that only models in which acquisition of additional genes is, on average, slightly beneficial yield a good fit to genomic data. These results suggest that the number of genes in prokaryotic genomes reflects the equilibrium between the benefit of additional genes that diminishes as the genome grows and deletion bias (i.e., the rate of deletion of genetic material being slightly greater than the rate of acquisition). Thus, new genes acquired by microbial genomes, on average, appear to be adaptive. The tight spacing of protein-coding genes likely results from a combination of the deletion bias and purifying selection that efficiently eliminates nonfunctional, noncoding sequences.
Collapse
|
36
|
Zenil-Ferguson R, Ponciano JM, Burleigh JG. Evaluating the role of genome downsizing and size thresholds from genome size distributions in angiosperms. AMERICAN JOURNAL OF BOTANY 2016; 103:1175-1186. [PMID: 27206462 DOI: 10.3732/ajb.1500408] [Citation(s) in RCA: 23] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/17/2016] [Accepted: 01/07/2016] [Indexed: 06/05/2023]
Abstract
PREMISE OF THE STUDY Whole-genome duplications (WGDs) can rapidly increase genome size in angiosperms. Yet their mean genome size is not correlated with ploidy. We compared three hypotheses to explain the constancy of genome size means across ploidies. The genome downsizing hypothesis suggests that genome size will decrease by a given percentage after a WGD. The genome size threshold hypothesis assumes that taxa with large genomes or large monoploid numbers will fail to undergo or survive WGDs. Finally, the genome downsizing and threshold hypothesis suggests that both genome downsizing and thresholds affect the relationship between genome size means and ploidy. METHODS We performed nonparametric bootstrap simulations to compare observed angiosperm genome size means among species or genera against simulated genome sizes under the three different hypotheses. We evaluated the hypotheses using a decision theory approach and estimated the expected percentage of genome downsizing. KEY RESULTS The threshold hypothesis improves the approximations between mean genome size and simulated genome size. At the species level, the genome downsizing with thresholds hypothesis best explains the genome size means with a 15% genome downsizing percentage. In the genus level simulations, the monoploid number threshold hypothesis best explains the data. CONCLUSIONS Thresholds of genome size and monoploid number added to genome downsizing at species level simulations explain the observed means of angiosperm genome sizes, and monoploid number is important for determining the genome size mean at the genus level.
Collapse
Affiliation(s)
- Rosana Zenil-Ferguson
- Department of Biology, University of Florida, P. O. Box 118525, Gainesville, Florida 32611-8525 USA
| | - José M Ponciano
- Department of Biology, University of Florida, P. O. Box 118525, Gainesville, Florida 32611-8525 USA
| | - J Gordon Burleigh
- Department of Biology, University of Florida, P. O. Box 118525, Gainesville, Florida 32611-8525 USA
| |
Collapse
|
37
|
Gupta A, LaBar T, Miyagi M, Adami C. Evolution of Genome Size in Asexual Digital Organisms. Sci Rep 2016; 6:25786. [PMID: 27181837 PMCID: PMC4867773 DOI: 10.1038/srep25786] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Accepted: 04/22/2016] [Indexed: 12/21/2022] Open
Abstract
Genome sizes have evolved to vary widely, from 250 bases in viroids to 670 billion bases in some amoebas. This remarkable variation in genome size is the outcome of complex interactions between various evolutionary factors such as mutation rate and population size. While comparative genomics has uncovered how some of these evolutionary factors influence genome size, we still do not understand what drives genome size evolution. Specifically, it is not clear how the primordial mutational processes of base substitutions, insertions, and deletions influence genome size evolution in asexual organisms. Here, we use digital evolution to investigate genome size evolution by tracking genome edits and their fitness effects in real time. In agreement with empirical data, we find that mutation rate is inversely correlated with genome size in asexual populations. We show that at low point mutation rate, insertions are significantly more beneficial than deletions, driving genome expansion and the acquisition of phenotypic complexity. Conversely, the high mutational load experienced at high mutation rates inhibits genome growth, forcing the genomes to compress their genetic information. Our analyses suggest that the inverse relationship between mutation rate and genome size is a result of the tradeoff between evolving phenotypic innovation and limiting the mutational load.
Collapse
Affiliation(s)
- Aditi Gupta
- BEACON Center for the Study of Evolution in Action, Michigan State University, East Lansing, MI 48824, USA.,Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, MI 48824, USA
| | - Thomas LaBar
- BEACON Center for the Study of Evolution in Action, Michigan State University, East Lansing, MI 48824, USA.,Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, MI 48824, USA.,Program in Ecology, Evolutionary Biology, and Behavior, Michigan State University, East Lansing, MI 48824, USA
| | - Michael Miyagi
- Department of Integrative Biology, University of Texas at Austin, Austin, TX 78712, USA
| | - Christoph Adami
- BEACON Center for the Study of Evolution in Action, Michigan State University, East Lansing, MI 48824, USA.,Department of Microbiology and Molecular Genetics, Michigan State University, East Lansing, MI 48824, USA.,Program in Ecology, Evolutionary Biology, and Behavior, Michigan State University, East Lansing, MI 48824, USA.,Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA
| |
Collapse
|
38
|
Rius N, Guillén Y, Delprat A, Kapusta A, Feschotte C, Ruiz A. Exploration of the Drosophila buzzatii transposable element content suggests underestimation of repeats in Drosophila genomes. BMC Genomics 2016; 17:344. [PMID: 27164953 PMCID: PMC4862133 DOI: 10.1186/s12864-016-2648-8] [Citation(s) in RCA: 21] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/22/2015] [Accepted: 04/22/2016] [Indexed: 11/10/2022] Open
Abstract
Background Many new Drosophila genomes have been sequenced in recent years using new-generation sequencing platforms and assembly methods. Transposable elements (TEs), being repetitive sequences, are often misassembled, especially in the genomes sequenced with short reads. Consequently, the mobile fraction of many of the new genomes has not been analyzed in detail or compared with that of other genomes sequenced with different methods, which could shed light into the understanding of genome and TE evolution. Here we compare the TE content of three genomes: D. buzzatii st-1, j-19, and D. mojavensis. Results We have sequenced a new D. buzzatii genome (j-19) that complements the D. buzzatii reference genome (st-1) already published, and compared their TE contents with that of D. mojavensis. We found an underestimation of TE sequences in Drosophila genus NGS-genomes when compared to Sanger-genomes. To be able to compare genomes sequenced with different technologies, we developed a coverage-based method and applied it to the D. buzzatii st-1 and j-19 genome. Between 10.85 and 11.16 % of the D. buzzatii st-1 genome is made up of TEs, between 7 and 7,5 % of D. buzzatii j-19 genome, while TEs represent 15.35 % of the D. mojavensis genome. Helitrons are the most abundant order in the three genomes. Conclusions TEs in D. buzzatii are less abundant than in D. mojavensis, as expected according to the genome size and TE content positive correlation. However, TEs alone do not explain the genome size difference. TEs accumulate in the dot chromosomes and proximal regions of D. buzzatii and D. mojavensis chromosomes. We also report a significantly higher TE density in D. buzzatii and D. mojavensis X chromosomes, which is not expected under the current models. Our easy-to-use correction method allowed us to identify recently active families in D. buzzatii st-1 belonging to the LTR-retrotransposon superfamily Gypsy. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2648-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Nuria Rius
- Department de Genética i Microbiologia, Universitat Autònoma de Barcelona, Bellaterra (Barcelona), Spain.
| | - Yolanda Guillén
- Department de Genética i Microbiologia, Universitat Autònoma de Barcelona, Bellaterra (Barcelona), Spain
| | - Alejandra Delprat
- Department de Genética i Microbiologia, Universitat Autònoma de Barcelona, Bellaterra (Barcelona), Spain
| | - Aurélie Kapusta
- Department of Human Genetics, University of Utah School of Medicine, Salt Lake City, UT, USA
| | - Cédric Feschotte
- Department of Human Genetics, University of Utah School of Medicine, Salt Lake City, UT, USA
| | - Alfredo Ruiz
- Department de Genética i Microbiologia, Universitat Autònoma de Barcelona, Bellaterra (Barcelona), Spain
| |
Collapse
|
39
|
Canapa A, Barucca M, Biscotti MA, Forconi M, Olmo E. Transposons, Genome Size, and Evolutionary Insights in Animals. Cytogenet Genome Res 2016; 147:217-39. [PMID: 26967166 DOI: 10.1159/000444429] [Citation(s) in RCA: 83] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 12/03/2015] [Indexed: 11/19/2022] Open
Abstract
The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.
Collapse
Affiliation(s)
- Adriana Canapa
- Dipartimento di Scienze della Vita e dell'Ambiente, Universitx00E0; Politecnica delle Marche, Ancona, Italy
| | | | | | | | | |
Collapse
|
40
|
Craddock EM, Gall JG, Jonas M. Hawaiian Drosophila genomes: size variation and evolutionary expansions. Genetica 2016; 144:107-24. [PMID: 26790663 DOI: 10.1007/s10709-016-9882-5] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/01/2015] [Accepted: 01/09/2016] [Indexed: 01/24/2023]
Abstract
This paper reports genome sizes of one Hawaiian Scaptomyza and 16 endemic Hawaiian Drosophila species that include five members of the antopocerus species group, one member of the modified mouthpart group, and ten members of the picture wing clade. Genome size expansions have occurred independently multiple times among Hawaiian Drosophila lineages, and have resulted in an over 2.3-fold range of genome sizes among species, with the largest observed in Drosophila cyrtoloma (1C = 0.41 pg). We find evidence that these repeated genome size expansions were likely driven by the addition of significant amounts of heterochromatin and satellite DNA. For example, our data reveal that the addition of seven heterochromatic chromosome arms to the ancestral haploid karyotype, and a remarkable proportion of ~70 % satellite DNA, account for the greatly expanded size of the D. cyrtoloma genome. Moreover, the genomes of 13/17 Hawaiian picture wing species are composed of substantial proportions (22-70 %) of detectable satellites (all but one of which are AT-rich). Our results suggest that in this tightly knit group of recently evolved species, genomes have expanded, in large part, via evolutionary amplifications of satellite DNA sequences in centric and pericentric domains (especially of the X and dot chromosomes), which have resulted in longer acrocentric chromosomes or metacentrics with an added heterochromatic chromosome arm. We discuss possible evolutionary mechanisms that may have shaped these patterns, including rapid fixation of novel expanded genomes during founder-effect speciation.
Collapse
Affiliation(s)
- Elysse M Craddock
- Natural Sciences Building, Purchase College, State University of New York, 735 Anderson Hill Road, Purchase, NY, 10577, USA.
| | - Joseph G Gall
- Department of Embryology, Carnegie Institution for Science, Baltimore, MD, USA
| | - Mark Jonas
- Natural Sciences Building, Purchase College, State University of New York, 735 Anderson Hill Road, Purchase, NY, 10577, USA
| |
Collapse
|
41
|
Fierst JL, Willis JH, Thomas CG, Wang W, Reynolds RM, Ahearne TE, Cutter AD, Phillips PC. Reproductive Mode and the Evolution of Genome Size and Structure in Caenorhabditis Nematodes. PLoS Genet 2015; 11:e1005323. [PMID: 26114425 PMCID: PMC4482642 DOI: 10.1371/journal.pgen.1005323] [Citation(s) in RCA: 82] [Impact Index Per Article: 9.1] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/29/2014] [Accepted: 05/31/2015] [Indexed: 11/18/2022] Open
Abstract
The self-fertile nematode worms Caenorhabditis elegans, C. briggsae, and C. tropicalis evolved independently from outcrossing male-female ancestors and have genomes 20-40% smaller than closely related outcrossing relatives. This pattern of smaller genomes for selfing species and larger genomes for closely related outcrossing species is also seen in plants. We use comparative genomics, including the first high quality genome assembly for an outcrossing member of the genus (C. remanei) to test several hypotheses for the evolution of genome reduction under a change in mating system. Unlike plants, it does not appear that reductions in the number of repetitive elements, such as transposable elements, are an important contributor to the change in genome size. Instead, all functional genomic categories are lost in approximately equal proportions. Theory predicts that self-fertilization should equalize the effective population size, as well as the resulting effects of genetic drift, between the X chromosome and autosomes. Contrary to this, we find that the self-fertile C. briggsae and C. elegans have larger intergenic spaces and larger protein-coding genes on the X chromosome when compared to autosomes, while C. remanei actually has smaller introns on the X chromosome than either self-reproducing species. Rather than being driven by mutational biases and/or genetic drift caused by a reduction in effective population size under self reproduction, changes in genome size in this group of nematodes appear to be caused by genome-wide patterns of gene loss, most likely generated by genomic adaptation to self reproduction per se.
Collapse
Affiliation(s)
- Janna L. Fierst
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - John H. Willis
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - Cristel G. Thomas
- Department of Ecology and Evolutionary Biology and Centre for the Analysis of Genome Evolution and Function, University of Toronto, Ontario, Canada
| | - Wei Wang
- Department of Ecology and Evolutionary Biology and Centre for the Analysis of Genome Evolution and Function, University of Toronto, Ontario, Canada
| | - Rose M. Reynolds
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - Timothy E. Ahearne
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| | - Asher D. Cutter
- Department of Ecology and Evolutionary Biology and Centre for the Analysis of Genome Evolution and Function, University of Toronto, Ontario, Canada
| | - Patrick C. Phillips
- Institute of Ecology and Evolution, University of Oregon, Eugene, Oregon, United States of America
| |
Collapse
|
42
|
Rastogi S, Kalra A, Gupta V, Khan F, Lal RK, Tripathi AK, Parameswaran S, Gopalakrishnan C, Ramaswamy G, Shasany AK. Unravelling the genome of Holy basil: an "incomparable" "elixir of life" of traditional Indian medicine. BMC Genomics 2015; 16:413. [PMID: 26017011 PMCID: PMC4445982 DOI: 10.1186/s12864-015-1640-z] [Citation(s) in RCA: 35] [Impact Index Per Article: 3.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/01/2014] [Accepted: 05/18/2015] [Indexed: 01/22/2023] Open
Abstract
Background Ocimum sanctum L. (O. tenuiflorum) family-Lamiaceae is an important component of Indian tradition of medicine as well as culture around the world, and hence is known as “Holy basil” in India. This plant is mentioned in the ancient texts of Ayurveda as an “elixir of life” (life saving) herb and worshipped for over 3000 years due to its healing properties. Although used in various ailments, validation of molecules for differential activities is yet to be fully analyzed, as about 80 % of the patents on this plant are on extracts or the plant parts, and mainly focussed on essential oil components. With a view to understand the full metabolic potential of this plant whole nuclear and chloroplast genomes were sequenced for the first time combining the sequence data from 4 libraries and three NGS platforms. Results The saturated draft assembly of the genome was about 386 Mb, along with the plastid genome of 142,245 bp, turning out to be the smallest in Lamiaceae. In addition to SSR markers, 136 proteins were identified as homologous to five important plant genomes. Pathway analysis indicated an abundance of phenylpropanoids in O. sanctum. Phylogenetic analysis for chloroplast proteome placed Salvia miltiorrhiza as the nearest neighbor. Comparison of the chemical compounds and genes availability in O. sanctum and S. miltiorrhiza indicated the potential for the discovery of new active molecules. Conclusion The genome sequence and annotation of O. sanctum provides new insights into the function of genes and the medicinal nature of the metabolites synthesized in this plant. This information is highly beneficial for mining biosynthetic pathways for important metabolites in related species. Electronic supplementary material The online version of this article (doi:10.1186/s12864-015-1640-z) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Shubhra Rastogi
- Biotechnology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, U.P., India.
| | - Alok Kalra
- Microbial Technology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, Uttar Pradesh, India.
| | - Vikrant Gupta
- Biotechnology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, U.P., India.
| | - Feroz Khan
- Metabolic and Structural Biology Department, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, U.P., India.
| | - Raj Kishori Lal
- Genetics and Plant Breeding Division, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, U.P., India.
| | - Anil Kumar Tripathi
- Biotechnology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, U.P., India.
| | - Sriram Parameswaran
- Research and Development Unit, Genotypic Technology Private Limited, Bangalore, Karnataka, 560094, India.
| | - Chellappa Gopalakrishnan
- Research and Development Unit, Genotypic Technology Private Limited, Bangalore, Karnataka, 560094, India.
| | - Gopalakrishna Ramaswamy
- Research and Development Unit, Genotypic Technology Private Limited, Bangalore, Karnataka, 560094, India.
| | - Ajit Kumar Shasany
- Biotechnology Division, CSIR-Central Institute of Medicinal and Aromatic Plants, P.O. CIMAP, Lucknow, 226015, U.P., India.
| |
Collapse
|
43
|
Kern AD, Barbash DA, Chang Mell J, Hupalo D, Jensen A. Highly constrained intergenic Drosophila ultraconserved elements are candidate ncRNAs. Genome Biol Evol 2015; 7:689-98. [PMID: 25618141 PMCID: PMC5322558 DOI: 10.1093/gbe/evv011] [Citation(s) in RCA: 14] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
Eukaryotes contain short (∼80-200 bp) regions that have few or no substitutions among species that represent hundreds of millions of years of evolutionary divergence. These ultraconserved elements (UCEs) are candidates for containing essential functions, but their biological roles remain largely unknown. Here, we report the discovery and characterization of UCEs from 12 sequenced Drosophila species. We identified 98 elements ≥80 bp long with very high conservation across the Drosophila phylogeny. Population genetic analyses reveal that these UCEs are not present in mutational cold spots. Instead we infer that they experience a level of selective constraint almost 10-fold higher compared with missense mutations in protein-coding sequences, which is substantially higher than that observed previously for human UCEs. About one-half of these Drosophila UCEs overlap the transcribed portion of genes, with many of those that are within coding sequences likely to correspond to sites of ADAR-dependent RNA editing. For the remaining UCEs that are in nongenic regions, we find that many are potentially capable of forming RNA secondary structures. Among ten chosen for further analysis, we discovered that the majority are transcribed in multiple tissues of Drosophila melanogaster. We conclude that Drosophila species are rich with UCEs and that many of them may correspond to novel noncoding RNAs.
Collapse
Affiliation(s)
| | - Daniel A Barbash
- Department of Molecular Biology and Genetics, Cornell University
| | - Joshua Chang Mell
- Department of Microbiology and Immunology, Drexel University College of Medicine
| | - Daniel Hupalo
- Department of Biology, Dartmouth College, Hanover, New Hampshire
| | - Amanda Jensen
- Department of Biology, Dartmouth College, Hanover, New Hampshire
| |
Collapse
|
44
|
Frahry MB, Sun C, Chong RA, Mueller RL. Low levels of LTR retrotransposon deletion by ectopic recombination in the gigantic genomes of salamanders. J Mol Evol 2015; 80:120-9. [PMID: 25608479 DOI: 10.1007/s00239-014-9663-7] [Citation(s) in RCA: 26] [Impact Index Per Article: 2.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2014] [Accepted: 12/30/2014] [Indexed: 11/25/2022]
Abstract
Across the tree of life, species vary dramatically in nuclear genome size. Mutations that add or remove sequences from genomes-insertions or deletions, or indels-are the ultimate source of this variation. Differences in the tempo and mode of insertion and deletion across taxa have been proposed to contribute to evolutionary diversity in genome size. Among vertebrates, most of the largest genomes are found within the salamanders, an amphibian clade with genome sizes ranging from ~14 to ~120 Gb. Salamander genomes have been shown to experience slower rates of DNA loss through small (i.e., <30 bp) deletions than do other vertebrate genomes. However, no studies have addressed DNA loss from salamander genomes resulting from larger deletions. Here, we focus on one type of large deletion-ectopic-recombination-mediated removal of LTR retrotransposon sequences. In ectopic recombination, double-strand breaks are repaired using a "wrong" (i.e., ectopic, or non-allelic) template sequence-typically another locus of similar sequence. When breaks occur within the LTR portions of LTR retrotransposons, ectopic-recombination-mediated repair can produce deletions that remove the internal transposon sequence and the equivalent of one of the two LTR sequences. These deletions leave a signature in the genome-a solo LTR sequence. We compared levels of solo LTRs in the genomes of four salamander species with levels present in five vertebrates with smaller genomes. Our results demonstrate that salamanders have low levels of solo LTRs, suggesting that ectopic-recombination-mediated deletion of LTR retrotransposons occurs more slowly than in other vertebrates with smaller genomes.
Collapse
Affiliation(s)
- Matthew Blake Frahry
- Department of Biology, Colorado State University, Fort Collins, CO, 80523-1878, USA
| | | | | | | |
Collapse
|
45
|
Sloan DB, Wu Z. History of plastid DNA insertions reveals weak deletion and at mutation biases in angiosperm mitochondrial genomes. Genome Biol Evol 2014; 6:3210-21. [PMID: 25416619 PMCID: PMC4986453 DOI: 10.1093/gbe/evu253] [Citation(s) in RCA: 44] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/31/2023] Open
Abstract
Angiosperm mitochondrial genomes exhibit many unusual properties, including heterogeneous nucleotide composition and exceptionally large and variable genome sizes. Determining the role of nonadaptive mechanisms such as mutation bias in shaping the molecular evolution of these unique genomes has proven challenging because their dynamic structures generally prevent identification of homologous intergenic sequences for comparative analyses. Here, we report an analysis of angiosperm mitochondrial DNA sequences that are derived from inserted plastid DNA (mtpts). The availability of numerous completely sequenced plastid genomes allows us to infer the evolutionary history of these insertions, including the specific nucleotide substitutions and indels that have occurred because their incorporation into the mitochondrial genome. Our analysis confirmed that many mtpts have a complex history, including frequent gene conversion and multiple examples of horizontal transfer between divergent angiosperm lineages. Nevertheless, it is clear that the majority of extant mtpt sequence in angiosperms is the product of recent transfer (or gene conversion) and is subject to rapid loss/deterioration, suggesting that most mtpts are evolving relatively free from functional constraint. The evolution of mtpt sequences reveals a pattern of biased mutational input in angiosperm mitochondrial genomes, including an excess of small deletions over insertions and a skew toward nucleotide substitutions that increase AT content. However, these mutation biases are far weaker than have been observed in many other cellular genomes, providing insight into some of the notable features of angiosperm mitochondrial architecture, including the retention of large intergenic regions and the relatively neutral GC content found in these regions.
Collapse
Affiliation(s)
- Daniel B Sloan
- Department of Biology, Colorado State University, Fort Collins
| | - Zhiqiang Wu
- Department of Biology, Colorado State University, Fort Collins
| |
Collapse
|
46
|
Metcalfe CJ, Casane D. Accommodating the load: The transposable element content of very large genomes. Mob Genet Elements 2014; 3:e24775. [PMID: 24616835 PMCID: PMC3943481 DOI: 10.4161/mge.24775] [Citation(s) in RCA: 28] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/22/2013] [Revised: 04/20/2013] [Accepted: 04/22/2013] [Indexed: 01/31/2023] Open
Abstract
Very large genomes, that is, those above 20 Gb, are rare but widely distributed throughout the eukaryotes. They are found within the diatoms, dinoflagellates, metazoans and green plants, but so far have not been found in the excavates. There is a known positive correlation between genome size and the proportion of the genome composed of transposable elements (TEs). Very large genomes may therefore be expected to be almost entirely composed of TEs. Of the large genomes examined, in the angiosperms, gymnosperms and the dinoflagellates only a small portion of the genome was identified as TEs, most of these genomes were unidentified and may be novel or diverse TEs. In the salamanders and lungfish, 25 to 47% of the genome were identifiable retrotransposons, that is, TEs that copy themselves before insertion. However, the predominant class of TEs found in the lungfish was not the same as that found in the salamanders. The little data we have at the moment suggests therefore that the diversity and abundance of TEs is variable between taxa with large genomes, similar to patterns found in taxa with smaller genomes. Based on results from the human genome, we suggest that the ‘missing’ portion of the lungfish and salamander genomes are old, highly divergent, and therefore inactive copies of TEs. The data available indicate that, unlike plants with large genomes, neither the lungfish nor the salamanders show an increased risk of extinction. Based on a slow rate of DNA loss in salamanders it has been suggested that the large salamander genome is the result of run-away genome expansion involving genome size increases via TE proliferation associated with reduced recombination rate. We know of no studies on DNA loss or recombination rates in lungfish genomes, however a similar scenario could describe the process of genome expansion in the lungfish. A series of waves of TE transposition and sequence decay would describe the pattern of TE content seen in both the lungfish and the salamanders. The lungfish and salamanders, therefore, may accommodate their large load of TEs because these TEs have accumulated gradually over a long period of time and have been subject to inactivation and decay.
Collapse
Affiliation(s)
- Cushla J Metcalfe
- Instituto de Biociências; Universidade de São Paulo; Cidade Universitária; São Paulo, Brazil
| | - Didier Casane
- Laboratoire Evolution Génomes et Spéciation; UPR9034 CNRS; Gif-sur-Yvette, France ; Université Paris Diderot; Sorbonne Paris Cité, France
| |
Collapse
|
47
|
Song H, Moulton MJ, Whiting MF. Rampant nuclear insertion of mtDNA across diverse lineages within Orthoptera (Insecta). PLoS One 2014; 9:e110508. [PMID: 25333882 PMCID: PMC4204883 DOI: 10.1371/journal.pone.0110508] [Citation(s) in RCA: 36] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/18/2014] [Accepted: 09/23/2014] [Indexed: 02/06/2023] Open
Abstract
Nuclear mitochondrial pseudogenes (numts) are non-functional fragments of mtDNA inserted into the nuclear genome. Numts are prevalent across eukaryotes and a positive correlation is known to exist between the number of numts and the genome size. Most numt surveys have relied on model organisms with fully sequenced nuclear genomes, but such analyses have limited utilities for making a generalization about the patterns of numt accumulation for any given clade. Among insects, the order Orthoptera is known to have the largest nuclear genome and it is also reported to include several species with a large number of numts. In this study, we use Orthoptera as a case study to document the diversity and abundance of numts by generating numts of three mitochondrial loci across 28 orthopteran families, representing the phylogenetic diversity of the order. We discover that numts are rampant in all lineages, but there is no discernable and consistent pattern of numt accumulation among different lineages. Likewise, we do not find any evidence that a certain mitochondrial gene is more prone to nuclear insertion than others. We also find that numt insertion must have occurred continuously and frequently throughout the diversification of Orthoptera. Although most numts are the result of recent nuclear insertion, we find evidence of very ancient numt insertion shared by highly divergent families dating back to the Jurassic period. Finally, we discuss several factors contributing to the extreme prevalence of numts in Orthoptera and highlight the importance of exploring the utility of numts in evolutionary studies.
Collapse
Affiliation(s)
- Hojun Song
- Department of Biology, University of Central Florida, Orlando, Florida, United States of America
| | - Matthew J. Moulton
- Department of Human Genetics, University of Utah, Salt Lake City, Utah, United States of America
- Department of Biology and M. L. Bean Museum, Brigham Young University, Provo, Utah, United States of America
| | - Michael F. Whiting
- Department of Biology and M. L. Bean Museum, Brigham Young University, Provo, Utah, United States of America
| |
Collapse
|
48
|
Fischer S, Bernard S, Beslon G, Knibbe C. A model for genome size evolution. Bull Math Biol 2014; 76:2249-91. [PMID: 25142746 PMCID: PMC4153982 DOI: 10.1007/s11538-014-9997-8] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/15/2013] [Accepted: 07/15/2014] [Indexed: 01/27/2023]
Abstract
We present a model for genome size evolution that takes into account both local mutations such as small insertions and small deletions, and large chromosomal rearrangements such as duplications and large deletions. We introduce the possibility of undergoing several mutations within one generation. The model, albeit minimalist, reveals a non-trivial spontaneous dynamics of genome size: in the absence of selection, an arbitrary large part of genomes remains beneath a finite size, even for a duplication rate 2.6-fold higher than the rate of large deletions, and even if there is also a systematic bias toward small insertions compared to small deletions. Specifically, we show that the condition of existence of an asymptotic stationary distribution for genome size non-trivially depends on the rates and mean sizes of the different mutation types. We also give upper bounds for the median and other quantiles of the genome size distribution, and argue that these bounds cannot be overcome by selection. Taken together, our results show that the spontaneous dynamics of genome size naturally prevents it from growing infinitely, even in cases where intuition would suggest an infinite growth. Using quantitative numerical examples, we show that, in practice, a shrinkage bias appears very quickly in genomes undergoing mutation accumulation, even though DNA gains and losses appear to be perfectly symmetrical at first sight. We discuss this spontaneous dynamics in the light of the other evolutionary forces proposed in the literature and argue that it provides them a stability-related size limit below which they can act.
Collapse
Affiliation(s)
- Stephan Fischer
- INSA-Lyon, Inria, CNRS, LIRIS, UMR5205, 69621, Villeurbanne, France,
| | | | | | | |
Collapse
|
49
|
Ellis LL, Huang W, Quinn AM, Ahuja A, Alfrejd B, Gomez FE, Hjelmen CE, Moore KL, Mackay TFC, Johnston JS, Tarone AM. Intrapopulation genome size variation in D. melanogaster reflects life history variation and plasticity. PLoS Genet 2014; 10:e1004522. [PMID: 25057905 PMCID: PMC4109859 DOI: 10.1371/journal.pgen.1004522] [Citation(s) in RCA: 49] [Impact Index Per Article: 4.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2013] [Accepted: 04/24/2014] [Indexed: 12/01/2022] Open
Abstract
We determined female genome sizes using flow cytometry for 211 Drosophila melanogaster sequenced inbred strains from the Drosophila Genetic Reference Panel, and found significant conspecific and intrapopulation variation in genome size. We also compared several life history traits for 25 lines with large and 25 lines with small genomes in three thermal environments, and found that genome size as well as genome size by temperature interactions significantly correlated with survival to pupation and adulthood, time to pupation, female pupal mass, and female eclosion rates. Genome size accounted for up to 23% of the variation in developmental phenotypes, but the contribution of genome size to variation in life history traits was plastic and varied according to the thermal environment. Expression data implicate differences in metabolism that correspond to genome size variation. These results indicate that significant genome size variation exists within D. melanogaster and this variation may impact the evolutionary ecology of the species. Genome size variation accounts for a significant portion of life history variation in an environmentally dependent manner, suggesting that potential fitness effects associated with genome size variation also depend on environmental conditions. Genome size evolution is ubiquitous, and–even after decades of research–mysterious. There are two major classes of hypotheses regarding genome size evolution, those that attribute its causes to evolutionarily neutral processes and those that suggest a role for selection. Numerous correlations between genome size and fitness-related phenotypes have been documented, suggesting selection could play a role in genome size evolution. Unfortunately, many of the effects in those studies are confounded with factors that could provide alternative explanations. Here, we show that 211 inbred strains of Drosophila melanogaster exhibit abundant variation in genome size, which correlates with life history traits in a temperature-dependent manner. Gene expression analyses suggest a role for differences in metabolism between strains with large and small genomes. Thus, there is genetic variation in genome size within D. melanogaster, and this variation is connected to variation in environmentally dependent life history traits. These observations indicate that selection is indeed a potential mechanism by which genome size can evolve. Our results also suggest that higher levels of genetic architecture may explain some of the genetic contribution to biologically important complex traits and raise the possibility that nucleotide quantity can contribute to phenotype in addition to quality.
Collapse
Affiliation(s)
- Lisa L. Ellis
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
| | - Wen Huang
- Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina, United States of America
| | - Andrew M. Quinn
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
- University of Texas Health Science Center, Houston, Texas, United States of America
| | - Astha Ahuja
- Baylor Human Genome Sequencing Center, Houston, Texas, United States of America
| | - Ben Alfrejd
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
| | - Francisco E. Gomez
- Department of Soil and Crop Science, Texas A&M University, College Station, Texas, United States of America
| | - Carl E. Hjelmen
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
| | - Kristi L. Moore
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
| | - Trudy F. C. Mackay
- Department of Biological Sciences, North Carolina State University, Raleigh, North Carolina, United States of America
| | - J. Spencer Johnston
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
| | - Aaron M. Tarone
- Department of Entomology, Texas A&M University, College Station, Texas, United States of America
- * E-mail:
| |
Collapse
|
50
|
Abstract
Mutation is the ultimate source of genetic variation. The most direct and unbiased method of studying spontaneous mutations is via mutation accumulation (MA) lines. Until recently, MA experiments were limited by the cost of sequencing and thus provided us with small numbers of mutational events and therefore imprecise estimates of rates and patterns of mutation. We used whole-genome sequencing to identify nearly 1,000 spontaneous mutation events accumulated over ∼311,000 generations in 145 diploid MA lines of the budding yeast Saccharomyces cerevisiae. MA experiments are usually assumed to have negligible levels of selection, but even mild selection will remove strongly deleterious events. We take advantage of such patterns of selection and show that mutation classes such as indels and aneuploidies (especially monosomies) are proportionately much more likely to contribute mutations of large effect. We also provide conservative estimates of indel, aneuploidy, environment-dependent dominant lethal, and recessive lethal mutation rates. To our knowledge, for the first time in yeast MA data, we identified a sufficiently large number of single-nucleotide mutations to measure context-dependent mutation rates and were able to (i) confirm strong AT bias of mutation in yeast driven by high rate of mutations from C/G to T/A and (ii) detect a higher rate of mutation at C/G nucleotides in two specific contexts consistent with cytosine methylation in S. cerevisiae.
Collapse
|