1
|
Hopkins BR, Angus-Henry A, Kim BY, Carlisle JA, Thompson A, Kopp A. Decoupled evolution of the Sex Peptide gene family and Sex Peptide Receptor in Drosophilidae. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.06.29.547128. [PMID: 37425821 PMCID: PMC10327216 DOI: 10.1101/2023.06.29.547128] [Citation(s) in RCA: 2] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 07/11/2023]
Abstract
Across internally fertilising species, males transfer ejaculate proteins that trigger wide-ranging changes in female behaviour and physiology. Much theory has been developed to explore the drivers of ejaculate protein evolution. The accelerating availability of high-quality genomes now allows us to test how these proteins are evolving at fine taxonomic scales. Here, we use genomes from 264 species to chart the evolutionary history of Sex Peptide (SP), a potent regulator of female post-mating responses in Drosophila melanogaster. We infer that SP first evolved in the Drosophilinae subfamily and has followed markedly different evolutionary trajectories in different lineages. Outside of the Sophophora-Lordiphosa, SP exists largely as a single-copy gene with independent losses in several lineages. Within the Sophophora-Lordiphosa, the SP gene family has repeatedly and independently expanded. Up to seven copies, collectively displaying extensive sequence variation, are present in some species. Despite these changes, SP expression remains restricted to the male reproductive tract. Alongside, we document considerable interspecific variation in the presence and morphology of seminal microcarriers that, despite the critical role SP plays in microcarrier assembly in D. melanogaster, appear to be independent of changes in the presence/absence or sequence of SP. We end by providing evidence that SP's evolution is decoupled from that of its receptor, SPR, in which we detect no evidence of correlated diversifying selection. Collectively, our work describes the divergent evolutionary trajectories that a novel gene has taken following its origin and finds a surprisingly weak coevolutionary signal between a supposedly sexually antagonistic protein and its receptor.
Collapse
Affiliation(s)
- Ben R. Hopkins
- Department of Evolution and Ecology, University of California – Davis, CA, USA
| | - Aidan Angus-Henry
- Department of Evolution and Ecology, University of California – Davis, CA, USA
| | | | - Jolie A. Carlisle
- Department of Molecular Biology and Genetics, Cornell University, Ithaca, NY, USA
| | - Ammon Thompson
- Department of Evolution and Ecology, University of California – Davis, CA, USA
| | - Artyom Kopp
- Department of Evolution and Ecology, University of California – Davis, CA, USA
| |
Collapse
|
2
|
Karageorgiou C, Tarrío R, Rodríguez-Trelles F. The Cyclically Seasonal Drosophila subobscura Inversion O 7 Originated From Fragile Genomic Sites and Relocated Immunity and Metabolic Genes. Front Genet 2020; 11:565836. [PMID: 33193649 PMCID: PMC7584159 DOI: 10.3389/fgene.2020.565836] [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] [Received: 05/26/2020] [Accepted: 09/09/2020] [Indexed: 11/28/2022] Open
Abstract
Chromosome inversions are important contributors to standing genetic variation in Drosophila subobscura. Presently, the species is experiencing a rapid replacement of high-latitude by low-latitude inversions associated with global warming. Yet not all low-latitude inversions are correlated with the ongoing warming trend. This is particularly unexpected in the case of O7 because it shows a regular seasonal cycle that peaks in summer and rose with a heatwave. The inconsistent behavior of O7 across components of the ambient temperature suggests that is causally more complex than simply due to temperature alone. In order to understand the dynamics of O7, high-quality genomic data are needed to determine both the breakpoints and the genetic content. To fill this gap, here we generated a PacBio long read-based chromosome-scale genome assembly, from a highly homozygous line made isogenic for an O3 + 4 + 7 chromosome. Then we isolated the complete continuous sequence of O7 by conserved synteny analysis with the available reference genome. Main findings include the following: (i) the assembled O7 inversion stretches 9.936 Mb, containing > 1,000 annotated genes; (ii) O7 had a complex origin, involving multiple breaks associated with non-B DNA-forming motifs, formation of a microinversion, and ectopic repair in trans with the two homologous chromosomes; (iii) the O7 breakpoints carry a pre-inversion record of fragility, including a sequence insertion, and transposition with later inverted duplication of an Attacin immunity gene; and (iv) the O7 inversion relocated the major insulin signaling forkhead box subgroup O (foxo) gene in tight linkage with its antagonistic regulatory partner serine/threonine-protein kinase B (Akt1) and disrupted concerted evolution of the two inverted Attacin duplicates, reattaching them to dFOXO metabolic enhancers. Our findings suggest that O7 exerts antagonistic pleiotropic effects on reproduction and immunity, setting a framework to understand its relationship with climate change. Furthermore, they are relevant for fragility in genome rearrangement evolution and for current views on the contribution of breakage versus repair in shaping inversion-breakpoint junctions.
Collapse
Affiliation(s)
- Charikleia Karageorgiou
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Barcelona, Spain
| | - Rosa Tarrío
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Barcelona, Spain
| | - Francisco Rodríguez-Trelles
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Barcelona, Spain
| |
Collapse
|
3
|
Chromosome-Level Assembly of Drosophila bifasciata Reveals Important Karyotypic Transition of the X Chromosome. G3-GENES GENOMES GENETICS 2020; 10:891-897. [PMID: 31969429 PMCID: PMC7056972 DOI: 10.1534/g3.119.400922] [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] [Indexed: 12/22/2022]
Abstract
The Drosophila obscura species group is one of the most studied clades of Drosophila and harbors multiple distinct karyotypes. Here we present a de novo genome assembly and annotation of D. bifasciata, a species which represents an important subgroup for which no high-quality chromosome-level genome assembly currently exists. We combined long-read sequencing (Nanopore) and Hi-C scaffolding to achieve a highly contiguous genome assembly approximately 193 Mb in size, with repetitive elements constituting 30.1% of the total length. Drosophila bifasciata harbors four large metacentric chromosomes and the small dot, and our assembly contains each chromosome in a single scaffold, including the highly repetitive pericentromeres, which were largely composed of Jockey and Gypsy transposable elements. We annotated a total of 12,821 protein-coding genes and comparisons of synteny with D. athabasca orthologs show that the large metacentric pericentromeric regions of multiple chromosomes are conserved between these species. Importantly, Muller A (X chromosome) was found to be metacentric in D. bifasciata and the pericentromeric region appears homologous to the pericentromeric region of the fused Muller A-AD (XL and XR) of pseudoobscura/affinis subgroup species. Our finding suggests a metacentric ancestral X fused to a telocentric Muller D and created the large neo-X (Muller A-AD) chromosome ∼15 MYA. We also confirm the fusion of Muller C and D in D. bifasciata and show that it likely involved a centromere-centromere fusion.
Collapse
|
4
|
Orengo DJ, Puerma E, Cereijo U, Aguadé M. The molecular genealogy of sequential overlapping inversions implies both homologous chromosomes of a heterokaryotype in an inversion origin. Sci Rep 2019; 9:17009. [PMID: 31740730 PMCID: PMC6861252 DOI: 10.1038/s41598-019-53582-8] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/12/2019] [Accepted: 10/14/2019] [Indexed: 11/25/2022] Open
Abstract
Cytological and molecular studies have revealed that inversion chromosomal polymorphism is widespread across taxa and that inversions are among the most common structural changes fixed between species. Two major mechanisms have been proposed for the origin of inversions considering that breaks occur at either repetitive or non-homologous sequences. While inversions originating through the first mechanism might have a multiple origin, those originating through the latter mechanism would have a unique origin. Variation at regions flanking inversion breakpoints can be informative on the origin and history of inversions given the reduced recombination in heterokaryotypes. Here, we have analyzed nucleotide variation at a fragment flanking the most centromere-proximal shared breakpoint of several sequential overlapping inversions of the E chromosome of Drosophila subobscura —inversions E1, E2, E9 and E3. The molecular genealogy inferred from variation at this shared fragment does not exhibit the branching pattern expected according to the sequential origin of inversions. The detected discordance between the molecular and cytological genealogies has led us to consider a novel possibility for the origin of an inversion, and more specifically that one of these inversions originated on a heterokaryotype for chromosomal arrangements. Based on this premise, we propose three new models for inversions origin.
Collapse
Affiliation(s)
- Dorcas J Orengo
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, i Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - Eva Puerma
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, i Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain
| | - Unai Cereijo
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, i Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain.,Centre for Research in Agricultural Genomics, CSIC-IRTA-UAB-UB, Campus UAB, Bellaterra (Cerdanyola del Vallès), 08193, Barcelona, Spain
| | - Montserrat Aguadé
- Departament de Genètica, Microbiologia i Estadística, Facultat de Biologia, i Institut de Recerca de la Biodiversitat (IRBio), Universitat de Barcelona, Barcelona, Spain.
| |
Collapse
|
5
|
Bracewell R, Chatla K, Nalley MJ, Bachtrog D. Dynamic turnover of centromeres drives karyotype evolution in Drosophila. eLife 2019; 8:e49002. [PMID: 31524597 PMCID: PMC6795482 DOI: 10.7554/elife.49002] [Citation(s) in RCA: 42] [Impact Index Per Article: 8.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/03/2019] [Accepted: 09/12/2019] [Indexed: 12/21/2022] Open
Abstract
Centromeres are the basic unit for chromosome inheritance, but their evolutionary dynamics is poorly understood. We generate high-quality reference genomes for multiple Drosophila obscura group species to reconstruct karyotype evolution. All chromosomes in this lineage were ancestrally telocentric and the creation of metacentric chromosomes in some species was driven by de novo seeding of new centromeres at ancestrally gene-rich regions, independently of chromosomal rearrangements. The emergence of centromeres resulted in a drastic size increase due to repeat accumulation, and dozens of genes previously located in euchromatin are now embedded in pericentromeric heterochromatin. Metacentric chromosomes secondarily became telocentric in the pseudoobscura subgroup through centromere repositioning and a pericentric inversion. The former (peri)centric sequences left behind shrunk dramatically in size after their inactivation, yet contain remnants of their evolutionary past, including increased repeat-content and heterochromatic environment. Centromere movements are accompanied by rapid turnover of the major satellite DNA detected in (peri)centromeric regions.
Collapse
Affiliation(s)
- Ryan Bracewell
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyUnited States
| | - Kamalakar Chatla
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyUnited States
| | - Matthew J Nalley
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyUnited States
| | - Doris Bachtrog
- Department of Integrative BiologyUniversity of California, BerkeleyBerkeleyUnited States
| |
Collapse
|
6
|
Taverner AM, Yang L, Barile ZJ, Lin B, Peng J, Pinharanda AP, Rao AS, Roland BP, Talsma AD, Wei D, Petschenka G, Palladino MJ, Andolfatto P. Adaptive substitutions underlying cardiac glycoside insensitivity in insects exhibit epistasis in vivo. eLife 2019; 8:48224. [PMID: 31453806 PMCID: PMC6733596 DOI: 10.7554/elife.48224] [Citation(s) in RCA: 18] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/06/2019] [Accepted: 08/24/2019] [Indexed: 01/20/2023] Open
Abstract
Predicting how species will respond to selection pressures requires understanding the factors that constrain their evolution. We use genome engineering of Drosophila to investigate constraints on the repeated evolution of unrelated herbivorous insects to toxic cardiac glycosides, which primarily occurs via a small subset of possible functionally-relevant substitutions to Na+,K+-ATPase. Surprisingly, we find that frequently observed adaptive substitutions at two sites, 111 and 122, are lethal when homozygous and adult heterozygotes exhibit dominant neural dysfunction. We identify a phylogenetically correlated substitution, A119S, that partially ameliorates the deleterious effects of substitutions at 111 and 122. Despite contributing little to cardiac glycoside-insensitivity in vitro, A119S, like substitutions at 111 and 122, substantially increases adult survivorship upon cardiac glycoside exposure. Our results demonstrate the importance of epistasis in constraining adaptive paths. Moreover, by revealing distinct effects of substitutions in vitro and in vivo, our results underscore the importance of evaluating the fitness of adaptive substitutions and their interactions in whole organisms.
Collapse
Affiliation(s)
- Andrew M Taverner
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, United States
| | - Lu Yang
- Department of Ecology and Evolutionary Biology, Princeton University, Princeton, United States
| | - Zachary J Barile
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, United States.,Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, United States
| | - Becky Lin
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, United States.,Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, United States
| | - Julie Peng
- Lewis-Sigler Institute for Integrative Genomics, Princeton University, Princeton, United States
| | - Ana P Pinharanda
- Department of Biological Sciences, Columbia University, New York, United States
| | - Arya S Rao
- Department of Biological Sciences, Columbia University, New York, United States
| | - Bartholomew P Roland
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, United States.,Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, United States
| | - Aaron D Talsma
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, United States.,Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, United States
| | - Daniel Wei
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, United States.,Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, United States
| | - Georg Petschenka
- Institute for Insect Biotechnology, Justus-Liebig-Universität Gießen, Hesse, Germany
| | - Michael J Palladino
- Department of Pharmacology and Chemical Biology, University of Pittsburgh, Pittsburgh, United States.,Pittsburgh Institute for Neurodegenerative Diseases (PIND), University of Pittsburgh School of Medicine, Pittsburgh, United States
| | - Peter Andolfatto
- Department of Biological Sciences, Columbia University, New York, United States
| |
Collapse
|
7
|
Karageorgiou C, Gámez-Visairas V, Tarrío R, Rodríguez-Trelles F. Long-read based assembly and synteny analysis of a reference Drosophila subobscura genome reveals signatures of structural evolution driven by inversions recombination-suppression effects. BMC Genomics 2019; 20:223. [PMID: 30885123 PMCID: PMC6423853 DOI: 10.1186/s12864-019-5590-8] [Citation(s) in RCA: 9] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/18/2018] [Accepted: 03/06/2019] [Indexed: 01/06/2023] Open
Abstract
BACKGROUND Drosophila subobscura has long been a central model in evolutionary genetics. Presently, its use is hindered by the lack of a reference genome. To bridge this gap, here we used PacBio long-read technology, together with the available wealth of genetic marker information, to assemble and annotate a high-quality nuclear and complete mitochondrial genome for the species. With the obtained assembly, we performed the first synteny analysis of genome structure evolution in the subobscura subgroup. RESULTS We generated a highly-contiguous ~ 129 Mb-long nuclear genome, consisting of six pseudochromosomes corresponding to the six chromosomes of a female haploid set, and a complete 15,764 bp-long mitogenome, and provide an account of their numbers and distributions of codifying and repetitive content. All 12 identified paracentric inversion differences in the subobscura subgroup would have originated by chromosomal breakage and repair, with some associated duplications, but no evidence of direct gene disruptions by the breakpoints. Between lineages, inversion fixation rates were 10 times higher in continental D. subobscura than in the two small oceanic-island endemics D. guanche and D. madeirensis. Within D. subobscura, we found contrasting ratios of chromosomal divergence to polymorphism between the A sex chromosome and the autosomes. CONCLUSIONS We present the first high-quality, long-read sequencing of a D. subobscura genome. Our findings generally support genome structure evolution in this species being driven indirectly, through the inversions' recombination-suppression effects in maintaining sets of adaptive alleles together in the face of gene flow. The resources developed will serve to further establish the subobscura subgroup as model for comparative genomics and evolutionary indicator of global change.
Collapse
Affiliation(s)
- Charikleia Karageorgiou
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Bellaterra, Barcelona, Spain
| | - Víctor Gámez-Visairas
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Bellaterra, Barcelona, Spain
| | - Rosa Tarrío
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Bellaterra, Barcelona, Spain
| | - Francisco Rodríguez-Trelles
- Grup de Genòmica, Bioinformàtica i Biologia Evolutiva (GGBE), Departament de Genètica i de Microbiologia, Universitat Autonòma de Barcelona, Bellaterra, Barcelona, Spain
| |
Collapse
|
8
|
The molecular characterization of fixed inversions breakpoints unveils the ancestral character of the Drosophila guanche chromosomal arrangements. Sci Rep 2019; 9:1706. [PMID: 30737415 PMCID: PMC6368638 DOI: 10.1038/s41598-018-37121-5] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/06/2018] [Accepted: 12/04/2018] [Indexed: 12/20/2022] Open
Abstract
Cytological studies revealed that the number of chromosomes and their organization varies across species. The increasing availability of whole genome sequences of multiple species across specific phylogenies has confirmed and greatly extended these cytological observations. In the Drosophila genus, the ancestral karyotype consists of five rod-like acrocentric chromosomes (Muller elements A to E) and one dot-like chromosome (element F), each exhibiting a generally conserved gene content. Chromosomal fusions and paracentric inversions are thus the major contributors, respectively, to chromosome number variation among species and to gene order variation within chromosomal element. The subobscura cluster of Drosophila consists in three species that retain the genus ancestral karyotype and differ by a reduced number of fixed inversions. Here, we have used cytological information and the D. guanche genome sequence to identify and molecularly characterize the breakpoints of inversions that became fixed since the D. guanche-D. subobscura split. Our results have led us to propose a modified version of the D. guanche cytological map of its X chromosome, and to establish that (i) most inversions became fixed in the D. subobscura lineage and (ii) the order in which the four X chromosome overlapping inversions occurred and became fixed.
Collapse
|