1
|
Rosser N, Seixas F, Queste LM, Cama B, Mori-Pezo R, Kryvokhyzha D, Nelson M, Waite-Hudson R, Goringe M, Costa M, Elias M, Mendes Eleres de Figueiredo C, Freitas AVL, Joron M, Kozak K, Lamas G, Martins ARP, McMillan WO, Ready J, Rueda-Muñoz N, Salazar C, Salazar P, Schulz S, Shirai LT, Silva-Brandão KL, Mallet J, Dasmahapatra KK. Hybrid speciation driven by multilocus introgression of ecological traits. Nature 2024; 628:811-817. [PMID: 38632397 PMCID: PMC11041799 DOI: 10.1038/s41586-024-07263-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2023] [Accepted: 03/01/2024] [Indexed: 04/19/2024]
Abstract
Hybridization allows adaptations to be shared among lineages and may trigger the evolution of new species1,2. However, convincing examples of homoploid hybrid speciation remain rare because it is challenging to demonstrate that hybridization was crucial in generating reproductive isolation3. Here we combine population genomic analysis with quantitative trait locus mapping of species-specific traits to examine a case of hybrid speciation in Heliconius butterflies. We show that Heliconius elevatus is a hybrid species that is sympatric with both parents and has persisted as an independently evolving lineage for at least 180,000 years. This is despite pervasive and ongoing gene flow with one parent, Heliconius pardalinus, which homogenizes 99% of their genomes. The remaining 1% introgressed from the other parent, Heliconius melpomene, and is scattered widely across the H. elevatus genome in islands of divergence from H. pardalinus. These islands contain multiple traits that are under disruptive selection, including colour pattern, wing shape, host plant preference, sex pheromones and mate choice. Collectively, these traits place H. elevatus on its own adaptive peak and permit coexistence with both parents. Our results show that speciation was driven by introgression of ecological traits, and that speciation with gene flow is possible with a multilocus genetic architecture.
Collapse
Affiliation(s)
- Neil Rosser
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA.
- Department of Biology, University of York, York, UK.
| | - Fernando Seixas
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
| | | | - Bruna Cama
- Department of Biology, University of York, York, UK
| | - Ronald Mori-Pezo
- URKU Estudios Amazónicos, Tarapoto, Perú
- Universidad Nacional Autónoma de Alto Amazona, Yurimaguas, Perú
| | - Dmytro Kryvokhyzha
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA
- Department of Clinical Sciences, Lund University Diabetes Centre, Malmö, Sweden
| | | | | | - Matt Goringe
- Department of Biology, University of York, York, UK
| | | | - Marianne Elias
- Institut Systématique, Evolution, Biodiversité, UMR 7205 MNHN-CNRS-EPHE-UPMC Sorbonne Universités, Muséum National d'Histoire Naturelle, Paris, France
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Clarisse Mendes Eleres de Figueiredo
- Institute for Biological Sciences, Federal University of Pará (UFPA), Belém, Brazil
- Centre for Advanced Studies of Biodiversity (CEABIO), Belém, Brazil
| | - André Victor Lucci Freitas
- Departamento de Biologia Animal and Museu de Diversidade Biológica, Instituto de Biologia, Universidade Estadual de Campinas, São Paulo, Brazil
| | - Mathieu Joron
- Centre d'Ecologie Fonctionnelle et Evolutive, UMR 5175 CNRS, Université de Montpellier-Université Paul Valéry Montpellier-EPHE, Montpellier, France
| | - Krzysztof Kozak
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Gerardo Lamas
- Museo de Historia Natural, Universidad Nacional Mayor de San Marcos, Lima, Peru
| | | | - W Owen McMillan
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Jonathan Ready
- Institute for Biological Sciences, Federal University of Pará (UFPA), Belém, Brazil
- Centre for Advanced Studies of Biodiversity (CEABIO), Belém, Brazil
| | - Nicol Rueda-Muñoz
- Biology Program, Faculty of Natural Sciences, Universidad del Rosario, Bogotá, Colombia
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences, Universidad del Rosario, Bogotá, Colombia
| | - Patricio Salazar
- Ecology and Evolutionary Biology, School of Biosciences, University of Sheffield, Sheffield, UK
| | - Stefan Schulz
- Institut für Organische Chemie, Technische Universität Braunschweig, Braunschweig, Germany
| | - Leila T Shirai
- Departamento de Biologia Animal and Museu de Diversidade Biológica, Instituto de Biologia, Universidade Estadual de Campinas, São Paulo, Brazil
| | - Karina L Silva-Brandão
- Leibniz Institute for the Analysis of Biodiversity Change, Museum de Natur Hamburg Zoology, Hamburg, Germany
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, USA.
| | - Kanchon K Dasmahapatra
- Department of Biology, University of York, York, UK
- Leverhulme Centre for Anthropocene Biodiversity, Department of Biology, University of York, York, UK
| |
Collapse
|
2
|
Wright CJ, Stevens L, Mackintosh A, Lawniczak M, Blaxter M. Comparative genomics reveals the dynamics of chromosome evolution in Lepidoptera. Nat Ecol Evol 2024; 8:777-790. [PMID: 38383850 PMCID: PMC11009112 DOI: 10.1038/s41559-024-02329-4] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/09/2023] [Accepted: 01/12/2024] [Indexed: 02/23/2024]
Abstract
Chromosomes are a central unit of genome organization. One-tenth of all described species on Earth are butterflies and moths, the Lepidoptera, which generally possess 31 chromosomes. However, some species display dramatic variation in chromosome number. Here we analyse 210 chromosomally complete lepidopteran genomes and show that the chromosomes of extant lepidopterans are derived from 32 ancestral linkage groups, which we term Merian elements. Merian elements have remained largely intact through 250 million years of evolution and diversification. Against this stable background, eight lineages have undergone extensive reorganization either through numerous fissions or a combination of fusion and fission events. Outside these lineages, fusions are rare and fissions are rarer still. Fusions often involve small, repeat-rich Merian elements and the sex-linked element. Our results reveal the constraints on genome architecture in Lepidoptera and provide a deeper understanding of chromosomal rearrangements in eukaryotic genome evolution.
Collapse
Affiliation(s)
| | - Lewis Stevens
- Tree of Life, Wellcome Sanger Institute, Cambridge, UK
| | | | | | - Mark Blaxter
- Tree of Life, Wellcome Sanger Institute, Cambridge, UK.
| |
Collapse
|
3
|
Briscoe AD. A genetic cause of male mate preference. Science 2024; 383:1290-1291. [PMID: 38513043 DOI: 10.1126/science.ado4079] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 03/23/2024]
Abstract
A gene for mate preference has been shared between hybridizing butterfly species.
Collapse
Affiliation(s)
- Adriana D Briscoe
- Department of Ecology and Evolutionary Biology, University of California, Irvine, CA, USA
| |
Collapse
|
4
|
Orteu A, Kucka M, Gordon IJ, Ng’iru I, van der Heijden ESM, Talavera G, Warren IA, Collins S, ffrench-Constant RH, Martins DJ, Chan YF, Jiggins CD, Martin SH. Transposable Element Insertions Are Associated with Batesian Mimicry in the Pantropical Butterfly Hypolimnas misippus. Mol Biol Evol 2024; 41:msae041. [PMID: 38401262 PMCID: PMC10924252 DOI: 10.1093/molbev/msae041] [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: 08/04/2023] [Revised: 02/14/2024] [Accepted: 02/16/2024] [Indexed: 02/26/2024] Open
Abstract
Hypolimnas misippus is a Batesian mimic of the toxic African Queen butterfly (Danaus chrysippus). Female H. misippus butterflies use two major wing patterning loci (M and A) to imitate three color morphs of D. chrysippus found in different regions of Africa. In this study, we examine the evolution of the M locus and identify it as an example of adaptive atavism. This phenomenon involves a morphological reversion to an ancestral character that results in an adaptive phenotype. We show that H. misippus has re-evolved an ancestral wing pattern present in other Hypolimnas species, repurposing it for Batesian mimicry of a D. chrysippus morph. Using haplotagging, a linked-read sequencing technology, and our new analytical tool, Wrath, we discover two large transposable element insertions located at the M locus and establish that these insertions are present in the dominant allele responsible for producing mimetic phenotype. By conducting a comparative analysis involving additional Hypolimnas species, we demonstrate that the dominant allele is derived. This suggests that, in the derived allele, the transposable elements disrupt a cis-regulatory element, leading to the reversion to an ancestral phenotype that is then utilized for Batesian mimicry of a distinct model, a different morph of D. chrysippus. Our findings present a compelling instance of convergent evolution and adaptive atavism, in which the same pattern element has independently evolved multiple times in Hypolimnas butterflies, repeatedly playing a role in Batesian mimicry of diverse model species.
Collapse
Affiliation(s)
- Anna Orteu
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, UK
- Tree of Life Programme, Wellcome Sanger Institute, Hinxton, UK
| | - Marek Kucka
- Friedrich Miescher Laboratory of the Max Planck Society, Tübingen, Germany
| | - Ian J Gordon
- Centre of Excellence in Biodiversity, University of Rwanda, Huye, Rwanda
| | - Ivy Ng’iru
- Mpala Research Centre, Nanyuki 10400, Laikipia, Kenya
- School of Biosciences, Cardiff University, Cardiff CF 10 3AX, UK
- UK Centre for Ecology and Hydrology, Wallingford OX10 8BB, UK
| | - Eva S M van der Heijden
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, UK
- Tree of Life Programme, Wellcome Sanger Institute, Hinxton, UK
| | - Gerard Talavera
- Institut Botànic de Barcelona (IBB), CSIC-CMCNB, Barcelona, Catalonia, Spain
| | - Ian A Warren
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, UK
| | - Steve Collins
- African Butterfly Research Institute, Nairobi, Kenya
| | | | - Dino J Martins
- Turkana Basin Institute, Stony Brook University, Stony Brook, NY 11794, USA
| | | | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, UK
| | - Simon H Martin
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, UK
| |
Collapse
|
5
|
Laetsch DR, Bisschop G, Martin SH, Aeschbacher S, Setter D, Lohse K. Demographically explicit scans for barriers to gene flow using gIMble. PLoS Genet 2023; 19:e1010999. [PMID: 37816069 PMCID: PMC10610087 DOI: 10.1371/journal.pgen.1010999] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/06/2023] [Revised: 10/27/2023] [Accepted: 09/25/2023] [Indexed: 10/12/2023] Open
Abstract
Identifying regions of the genome that act as barriers to gene flow between recently diverged taxa has remained challenging given the many evolutionary forces that generate variation in genetic diversity and divergence along the genome, and the stochastic nature of this variation. Progress has been impeded by a conceptual and methodological divide between analyses that infer the demographic history of speciation and genome scans aimed at identifying locally maladaptive alleles i.e. genomic barriers to gene flow. Here we implement genomewide IM blockwise likelihood estimation (gIMble), a composite likelihood approach for the quantification of barriers, that bridges this divide. This analytic framework captures background selection and selection against barriers in a model of isolation with migration (IM) as heterogeneity in effective population size (Ne) and effective migration rate (me), respectively. Variation in both effective demographic parameters is estimated in sliding windows via pre-computed likelihood grids. gIMble includes modules for pre-processing/filtering of genomic data and performing parametric bootstraps using coalescent simulations. To demonstrate the new approach, we analyse data from a well-studied pair of sister species of tropical butterflies with a known history of post-divergence gene flow: Heliconius melpomene and H. cydno. Our analyses uncover both large-effect barrier loci (including well-known wing-pattern genes) and a genome-wide signal of a polygenic barrier architecture.
Collapse
Affiliation(s)
- Dominik R. Laetsch
- Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, United Kingdom
| | - Gertjan Bisschop
- Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, United Kingdom
| | - Simon H. Martin
- Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, United Kingdom
| | - Simon Aeschbacher
- Department of Evolutionary Biology and Environmental Studies, University of Zurich, Zurich, Switzerland
| | - Derek Setter
- Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, United Kingdom
| | - Konrad Lohse
- Institute of Ecology and Evolution, University of Edinburgh, Edinburgh, United Kingdom
| |
Collapse
|
6
|
Näsvall K, Boman J, Höök L, Vila R, Wiklund C, Backström N. Nascent evolution of recombination rate differences as a consequence of chromosomal rearrangements. PLoS Genet 2023; 19:e1010717. [PMID: 37549188 PMCID: PMC10434929 DOI: 10.1371/journal.pgen.1010717] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/25/2023] [Revised: 08/17/2023] [Accepted: 07/17/2023] [Indexed: 08/09/2023] Open
Abstract
Reshuffling of genetic variation occurs both by independent assortment of chromosomes and by homologous recombination. Such reshuffling can generate novel allele combinations and break linkage between advantageous and deleterious variants which increases both the potential and the efficacy of natural selection. Here we used high-density linkage maps to characterize global and regional recombination rate variation in two populations of the wood white butterfly (Leptidea sinapis) that differ considerably in their karyotype as a consequence of at least 27 chromosome fissions and fusions. The recombination data were compared to estimates of genetic diversity and measures of selection to assess the relationship between chromosomal rearrangements, crossing over, maintenance of genetic diversity and adaptation. Our data show that the recombination rate is influenced by both chromosome size and number, but that the difference in the number of crossovers between karyotypes is reduced as a consequence of a higher frequency of double crossovers in larger chromosomes. As expected from effects of selection on linked sites, we observed an overall positive association between recombination rate and genetic diversity in both populations. Our results also revealed a significant effect of chromosomal rearrangements on the rate of intergenic diversity change between populations, but limited effects on polymorphisms in coding sequence. We conclude that chromosomal rearrangements can have considerable effects on the recombination landscape and consequently influence both maintenance of genetic diversity and efficiency of selection in natural populations.
Collapse
Affiliation(s)
- Karin Näsvall
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, Uppsala, Sweden
| | - Jesper Boman
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, Uppsala, Sweden
| | - Lars Höök
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, Uppsala, Sweden
| | - Roger Vila
- Butterfly Diversity and Evolution Lab, Institut de Biologia Evolutiva (CSIC-Univ. Pompeu Fabra), Barcelona, Spain
| | - Christer Wiklund
- Department of Zoology: Division of Ecology, Stockholm University, Stockholm, Sweden
| | - Niclas Backström
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, Uppsala, Sweden
| |
Collapse
|
7
|
Zhang L, Chaturvedi S, Nice CC, Lucas LK, Gompert Z. Population genomic evidence of selection on structural variants in a natural hybrid zone. Mol Ecol 2023; 32:1497-1514. [PMID: 35398939 DOI: 10.1111/mec.16469] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/17/2022] [Revised: 03/18/2022] [Accepted: 03/28/2022] [Indexed: 12/01/2022]
Abstract
Structural variants (SVs) can promote speciation by directly causing reproductive isolation or by suppressing recombination across large genomic regions. Whereas examples of each mechanism have been documented, systematic tests of the role of SVs in speciation are lacking. Here, we take advantage of long-read (Oxford nanopore) whole-genome sequencing and a hybrid zone between two Lycaeides butterfly taxa (L. melissa and Jackson Hole Lycaeides) to comprehensively evaluate genome-wide patterns of introgression for SVs and relate these patterns to hypotheses about speciation. We found >100,000 SVs segregating within or between the two hybridizing species. SVs and SNPs exhibited similar levels of genetic differentiation between species, with the exception of inversions, which were more differentiated. We detected credible variation in patterns of introgression among SV loci in the hybrid zone, with 562 of 1419 ancestry-informative SVs exhibiting genomic clines that deviated from null expectations based on genome-average ancestry. Overall, hybrids exhibited a directional shift towards Jackson Hole Lycaeides ancestry at SV loci, consistent with the hypothesis that these loci experienced more selection on average than SNP loci. Surprisingly, we found that deletions, rather than inversions, showed the highest skew towards excess ancestry from Jackson Hole Lycaeides. Excess Jackson Hole Lycaeides ancestry in hybrids was also especially pronounced for Z-linked SVs and inversions containing many genes. In conclusion, our results show that SVs are ubiquitous and suggest that SVs in general, but especially deletions, might disproportionately affect hybrid fitness and thus contribute to reproductive isolation.
Collapse
Affiliation(s)
- Linyi Zhang
- Department of Ecology and Evolutionary Biology, University of Toronto, Toronto, Ontario, Canada.,Department of Biology, Utah State University, Logan, Utah, USA
| | - Samridhi Chaturvedi
- Department of Integrative Biology, University of California, Berkeley, California, USA
| | - Chris C Nice
- Department of Biology, Texas State University, San Marcos, Texas, USA
| | - Lauren K Lucas
- Department of Biology, Utah State University, Logan, Utah, USA
| | | |
Collapse
|
8
|
Predicting recombination frequency from map distance. Heredity (Edinb) 2023; 130:114-121. [PMID: 36566319 PMCID: PMC9981558 DOI: 10.1038/s41437-022-00585-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2022] [Revised: 12/08/2022] [Accepted: 12/08/2022] [Indexed: 12/25/2022] Open
Abstract
Map distance is one of the key measures in genetics and indicates the expected number of crossovers between two loci. Map distance is estimated from the observed recombination frequency using mapping functions, the most widely used of those, Haldane and Kosambi, being developed at the time when the number of markers was low and unobserved crossovers had a substantial effect on the recombination fractions. In contemporary high-density marker data, the probability of multiple crossovers between adjacent loci is negligible and different mapping functions yield the same result, that is, the recombination frequency between adjacent loci is equal to the map distance in Morgans. However, high-density linkage maps contain an interpretation problem: the map distance over a long interval is additive and its association with recombination frequency is not defined. Here, we demonstrate with high-density linkage maps from humans and stickleback fishes that the inverses of Haldane's and Kosambi's mapping functions systematically underpredict recombination frequencies from map distance. To remedy this, we formulate a piecewise function that yields more accurate predictions of recombination frequency from map distance. Our results demonstrate that the association between map distance and recombination frequency is context-dependent and without a universal solution.
Collapse
|
9
|
Bastide H, López-Villavicencio M, Ogereau D, Lledo J, Dutrillaux AM, Debat V, Llaurens V. Genome assembly of 3 Amazonian Morpho butterfly species reveals Z-chromosome rearrangements between closely related species living in sympatry. Gigascience 2022; 12:giad033. [PMID: 37216769 PMCID: PMC10202424 DOI: 10.1093/gigascience/giad033] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/25/2022] [Revised: 02/13/2023] [Accepted: 05/12/2023] [Indexed: 05/24/2023] Open
Abstract
The genomic processes enabling speciation and species coexistence in sympatry are still largely unknown. Here we describe the whole-genome sequencing and assembly of 3 closely related species from the butterfly genus Morpho: Morpho achilles (Linnaeus, 1758), Morpho helenor (Cramer, 1776), and Morpho deidamia (Höbner, 1819). These large blue butterflies are emblematic species of the Amazonian rainforest. They live in sympatry in a wide range of their geographical distribution and display parallel diversification of dorsal wing color pattern, suggesting local mimicry. By sequencing, assembling, and annotating their genomes, we aim at uncovering prezygotic barriers preventing gene flow between these sympatric species. We found a genome size of 480 Mb for the 3 species and a chromosomal number ranging from 2n = 54 for M. deidamia to 2n = 56 for M. achilles and M. helenor. We also detected inversions on the sex chromosome Z that were differentially fixed between species, suggesting that chromosomal rearrangements may contribute to their reproductive isolation. The annotation of their genomes allowed us to recover in each species at least 12,000 protein-coding genes and to discover duplications of genes potentially involved in prezygotic isolation like genes controlling color discrimination (L-opsin). Altogether, the assembly and the annotation of these 3 new reference genomes open new research avenues into the genomic architecture of speciation and reinforcement in sympatry, establishing Morpho butterflies as a new eco-evolutionary model.
Collapse
Affiliation(s)
| | - Manuela López-Villavicencio
- Institut de Systématique, Evolution et Biodiversité (UMR 7205 CNRS/MNHN/SU/EPHE/UA), Muséum National d’Histoire Naturelle–CP50, 75005 Paris, France
| | | | - Joanna Lledo
- GeT-PlaGe, Bât G2, INRAe, 31326 Castanet-Tolosan Cedex, France
| | - Anne-Marie Dutrillaux
- Institut de Systématique, Evolution et Biodiversité (UMR 7205 CNRS/MNHN/SU/EPHE/UA), Muséum National d’Histoire Naturelle–CP50, 75005 Paris, France
| | - Vincent Debat
- Institut de Systématique, Evolution et Biodiversité (UMR 7205 CNRS/MNHN/SU/EPHE/UA), Muséum National d’Histoire Naturelle–CP50, 75005 Paris, France
| | - Violaine Llaurens
- Institut de Systématique, Evolution et Biodiversité (UMR 7205 CNRS/MNHN/SU/EPHE/UA), Muséum National d’Histoire Naturelle–CP50, 75005 Paris, France
| |
Collapse
|
10
|
Schuldiner‐Harpaz T, Merrill RM, Jiggins CD. Evolution of physical linkage between loci controlling ecological traits and mating preferences. J Evol Biol 2022; 35:1537-1547. [PMID: 36196988 PMCID: PMC9827829 DOI: 10.1111/jeb.14105] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/06/2022] [Revised: 07/26/2022] [Accepted: 08/08/2022] [Indexed: 01/12/2023]
Abstract
Coupling of multiple barriers to gene-flow, such as divergent local adaptation and reproductive isolation, facilitates speciation. However, alleles at loci that contribute to barrier effects can be dissociated by recombination. Models of linkage between diverging alleles often consider elements that reduce recombination, such as chromosomal inversions and alleles that modify recombination rate between existing loci. In contrast, here, we consider the evolution of linkage due to the close proximity of loci on the same chromosome. Examples of such physical linkage exist in several species, but in other cases, strong associations are maintained without physical linkage. We use an individual-based model to study the conditions under which the physical linkage between loci controlling ecological traits and mating preferences might be expected to evolve. We modelled a single locus controlling an ecological trait that acts also as a mating cue. Mating preferences are controlled by multiple loci, formed by mutations that are randomly placed in the "genome", within varying distances from the ecological trait locus, allowing us to examine which genomic architectures spread across the population. Our model reveals that stronger physical linkage is favoured when mating preferences and selection are weaker. Under such conditions mating among divergent phenotypes is more frequent, and matching ecological trait and mating preference alleles are more likely to become dissociated by recombination, favouring the evolution of genetic linkage. While most theoretical studies on clustering of divergent loci focus on how physical linkage influences speciation, we show how physical linkage itself can arise, establishing conditions that can favour speciation.
Collapse
|
11
|
Shipilina D, Näsvall K, Höök L, Vila R, Talavera G, Backström N. Linkage mapping and genome annotation give novel insights into gene family expansions and regional recombination rate variation in the painted lady (Vanessa cardui) butterfly. Genomics 2022; 114:110481. [PMID: 36115505 DOI: 10.1016/j.ygeno.2022.110481] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/04/2022] [Revised: 09/01/2022] [Accepted: 09/10/2022] [Indexed: 01/14/2023]
Abstract
Characterization of gene family expansions and crossing over is crucial for understanding how organisms adapt to the environment. Here, we develop a high-density linkage map and detailed genome annotation of the painted lady butterfly (Vanessa cardui) - a non-diapausing, highly polyphagous species famous for its long-distance migratory behavior and almost cosmopolitan distribution. Our results reveal a complex interplay between regional recombination rate variation, gene duplications and transposable element activity shaping the genome structure of the painted lady. We identify several lineage specific gene family expansions. Their functions are mainly associated with protein and fat metabolism, detoxification, and defense against infection - critical processes for the painted lady's unique life-history. Furthermore, the detailed recombination maps allow us to characterize the regional recombination landscape, data that reveal a strong effect of chromosome size on the recombination rate, a limited impact of GC-biased gene conversion and a positive association between recombination and short interspersed elements.
Collapse
Affiliation(s)
- Daria Shipilina
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden; Swedish Collegium for Advanced Study, Thunbergsvägen 2, 75236 Uppsala, Sweden.
| | - Karin Näsvall
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Lars Höök
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| | - Roger Vila
- The Butterfly Diversity and Evolution Lab, Institut de Biologia Evolutiva, Passeig Martim de la Barceloneta 37-49, 08003 Barcelona, Spain
| | - Gerard Talavera
- Institut Botànic de Barcelona (IBB), CSIC-Ajuntament de Barcelona, Passeig del Migdia s/n, 08038 Barcelona, Spain
| | - Niclas Backström
- Evolutionary Biology Program, Department of Ecology and Genetics, Uppsala University, Norbyvägen 18D, 75236 Uppsala, Sweden
| |
Collapse
|
12
|
Ruggieri AA, Livraghi L, Lewis JJ, Evans E, Cicconardi F, Hebberecht L, Ortiz-Ruiz Y, Montgomery SH, Ghezzi A, Rodriguez-Martinez JA, Jiggins CD, McMillan WO, Counterman BA, Papa R, Van Belleghem SM. A butterfly pan-genome reveals that a large amount of structural variation underlies the evolution of chromatin accessibility. Genome Res 2022; 32:1862-1875. [PMID: 36109150 PMCID: PMC9712634 DOI: 10.1101/gr.276839.122] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/16/2022] [Accepted: 09/13/2022] [Indexed: 01/16/2023]
Abstract
Despite insertions and deletions being the most common structural variants (SVs) found across genomes, not much is known about how much these SVs vary within populations and between closely related species, nor their significance in evolution. To address these questions, we characterized the evolution of indel SVs using genome assemblies of three closely related Heliconius butterfly species. Over the relatively short evolutionary timescales investigated, up to 18.0% of the genome was composed of indels between two haplotypes of an individual Heliconius charithonia butterfly and up to 62.7% included lineage-specific SVs between the genomes of the most distant species (11 Mya). Lineage-specific sequences were mostly characterized as transposable elements (TEs) inserted at random throughout the genome and their overall distribution was similarly affected by linked selection as single nucleotide substitutions. Using chromatin accessibility profiles (i.e., ATAC-seq) of head tissue in caterpillars to identify sequences with potential cis-regulatory function, we found that out of the 31,066 identified differences in chromatin accessibility between species, 30.4% were within lineage-specific SVs and 9.4% were characterized as TE insertions. These TE insertions were localized closer to gene transcription start sites than expected at random and were enriched for sites with significant resemblance to several transcription factor binding sites with known function in neuron development in Drosophila We also identified 24 TE insertions with head-specific chromatin accessibility. Our results show high rates of structural genome evolution that were previously overlooked in comparative genomic studies and suggest a high potential for structural variation to serve as raw material for adaptive evolution.
Collapse
Affiliation(s)
- Angelo A Ruggieri
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan PR 00931, Puerto Rico
| | - Luca Livraghi
- Department of Biological Sciences, The George Washington University, Washington, DC 20052, USA
- Smithsonian Tropical Research Institute, Apartado 0843-03092 Panamá, Panama
| | - James J Lewis
- Department of Zoology, University of British Columbia, Vancouver, British Columbia V6T 1Z4, Canada
| | - Elizabeth Evans
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan PR 00931, Puerto Rico
| | - Francesco Cicconardi
- School of Biological Sciences, Bristol University, Bristol BS8 1QU, United Kingdom
| | - Laura Hebberecht
- School of Biological Sciences, Bristol University, Bristol BS8 1QU, United Kingdom
| | - Yadira Ortiz-Ruiz
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan PR 00931, Puerto Rico
- Molecular Sciences and Research Center, University of Puerto Rico, San Juan 00926, Puerto Rico
| | - Stephen H Montgomery
- School of Biological Sciences, Bristol University, Bristol BS8 1QU, United Kingdom
| | - Alfredo Ghezzi
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan PR 00931, Puerto Rico
| | | | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, United Kingdom
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Apartado 0843-03092 Panamá, Panama
| | - Brian A Counterman
- Department of Biological Sciences, Auburn University, Auburn, Alabama 36849, USA
| | - Riccardo Papa
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan PR 00931, Puerto Rico
- Molecular Sciences and Research Center, University of Puerto Rico, San Juan 00926, Puerto Rico
| | - Steven M Van Belleghem
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan PR 00931, Puerto Rico
- Ecology, Evolution and Conservation Biology, Biology Department, KU Leuven, 3000 Leuven, Belgium
| |
Collapse
|
13
|
Guo R, Papanicolaou A, Fritz ML. Validation of reference-assisted assembly using existing and novel Heliothine genomes. Genomics 2022; 114:110441. [PMID: 35931274 DOI: 10.1016/j.ygeno.2022.110441] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/29/2021] [Revised: 07/19/2022] [Accepted: 07/29/2022] [Indexed: 11/16/2022]
Abstract
Chloridea subflexa and Chloridea virescens are a pair of closely related noctuid species exhibiting pheromone-based sexual isolation and divergent host plant preferences. We produced a novel Illumina short read C. subflexa genome assembly and an improved C. virescens genome assembly, which offer opportunities to study the genomic basis for evolutionarily important traits in this lepidopteran family with few genomic resources. We then examined the feasibility of reference-assisted assembly, an approach that leverages existing high quality genomic resources for genome improvement in closely related taxa and applied it to our Heliothine genomes. Our work demonstrates that reference-assisted assembly has the potential to enhance contiguity and completeness of existing insect genomic resources with minimal additional laboratory costs. We conclude by discussing both the potential and pitfalls of reference-assisted assembly according to the intended downstream assembly application.
Collapse
Affiliation(s)
- Rong Guo
- Department of Entomology, University of Maryland, College Park, MD 20742, USA; Computational Biology, Bioinformatics and Genomics Program, Department of Biological Sciences, University of Maryland, College Park, MD 20742, USA
| | - Alexie Papanicolaou
- Hawkesbury Institute for the Environment, Western Sydney University, Richmond, NSW 2753, Australia.
| | - Megan L Fritz
- Department of Entomology, University of Maryland, College Park, MD 20742, USA; Computational Biology, Bioinformatics and Genomics Program, Department of Biological Sciences, University of Maryland, College Park, MD 20742, USA.
| |
Collapse
|
14
|
Jay P, Leroy M, Le Poul Y, Whibley A, Arias M, Chouteau M, Joron M. Association mapping of colour variation in a butterfly provides evidence that a supergene locks together a cluster of adaptive loci. Philos Trans R Soc Lond B Biol Sci 2022; 377:20210193. [PMID: 35694756 PMCID: PMC9189503 DOI: 10.1098/rstb.2021.0193] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022] Open
Abstract
Supergenes are genetic architectures associated with discrete and concerted variation in multiple traits. It has long been suggested that supergenes control these complex polymorphisms by suppressing recombination between sets of coadapted genes. However, because recombination suppression hinders the dissociation of the individual effects of genes within supergenes, there is still little evidence that supergenes evolve by tightening linkage between coadapted genes. Here, combining a landmark-free phenotyping algorithm with multivariate genome-wide association studies, we dissected the genetic basis of wing pattern variation in the butterfly Heliconius numata. We show that the supergene controlling the striking wing pattern polymorphism displayed by this species contains several independent loci associated with different features of wing patterns. The three chromosomal inversions of this supergene suppress recombination between these loci, supporting the hypothesis that they may have evolved because they captured beneficial combinations of alleles. Some of these loci are, however, associated with colour variations only in a subset of morphs where the phenotype is controlled by derived inversion forms, indicating that they were recruited after the formation of the inversions. Our study shows that supergenes and clusters of adaptive loci in general may form via the evolution of chromosomal rearrangements suppressing recombination between co-adapted loci but also via the subsequent recruitment of linked adaptive mutations. This article is part of the theme issue 'Genomic architecture of supergenes: causes and evolutionary consequences'.
Collapse
Affiliation(s)
- Paul Jay
- CEFE, Université de Montpellier, CNRS, EPHE, IRD, 34293 Montpellier cedex 5, France
| | - Manon Leroy
- CEFE, Université de Montpellier, CNRS, EPHE, IRD, 34293 Montpellier cedex 5, France
| | - Yann Le Poul
- CEFE, Université de Montpellier, CNRS, EPHE, IRD, 34293 Montpellier cedex 5, France
| | - Annabel Whibley
- School of Biological Sciences, University of Auckland, Auckland 1010, New Zealand
| | - Mónica Arias
- CIRAD, UMR PHIM, F-34398 Montpellier, France.,PHIM, Univ Montpellier, CIRAD, INRAE, Institut Agro, IRD, CEDEX 5, 34398 Montpellier, France
| | - Mathieu Chouteau
- CEFE, Université de Montpellier, CNRS, EPHE, IRD, 34293 Montpellier cedex 5, France.,LEEISA, USR 63456, Université de Guyane, CNRS, IFREMER, 275 route de Montabo, 797334 Cayenne, French Guiana
| | - Mathieu Joron
- CEFE, Université de Montpellier, CNRS, EPHE, IRD, 34293 Montpellier cedex 5, France
| |
Collapse
|
15
|
Smolander OP, Blande D, Ahola V, Rastas P, Tanskanen J, Kammonen JI, Oostra V, Pellegrini L, Ikonen S, Dallas T, DiLeo MF, Duplouy A, Duru IC, Halimaa P, Kahilainen A, Kuwar SS, Kärenlampi SO, Lafuente E, Luo S, Makkonen J, Nair A, de la Paz Celorio-Mancera M, Pennanen V, Ruokolainen A, Sundell T, Tervahauta AI, Twort V, van Bergen E, Österman-Udd J, Paulin L, Frilander MJ, Auvinen P, Saastamoinen M. Improved chromosome-level genome assembly of the Glanville fritillary butterfly (Melitaea cinxia) integrating Pacific Biosciences long reads and a high-density linkage map. Gigascience 2022; 11:6505122. [PMID: 35022701 PMCID: PMC8756199 DOI: 10.1093/gigascience/giab097] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2020] [Revised: 05/03/2021] [Accepted: 12/14/2021] [Indexed: 11/12/2022] Open
Abstract
BACKGROUND The Glanville fritillary (Melitaea cinxia) butterfly is a model system for metapopulation dynamics research in fragmented landscapes. Here, we provide a chromosome-level assembly of the butterfly's genome produced from Pacific Biosciences sequencing of a pool of males, combined with a linkage map from population crosses. RESULTS The final assembly size of 484 Mb is an increase of 94 Mb on the previously published genome. Estimation of the completeness of the genome with BUSCO indicates that the genome contains 92-94% of the BUSCO genes in complete and single copies. We predicted 14,810 genes using the MAKER pipeline and manually curated 1,232 of these gene models. CONCLUSIONS The genome and its annotated gene models are a valuable resource for future comparative genomics, molecular biology, transcriptome, and genetics studies on this species.
Collapse
Affiliation(s)
- Olli-Pekka Smolander
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland.,Department of Chemistry and Biotechnology, Tallinn University of Technology, 12618 Tallinn, Estonia
| | - Daniel Blande
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Virpi Ahola
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland.,Ming Wai Lau Centre for Reparative Medicine, Karolinska Institutet, 171 77 Stockholm, Hong Kong
| | - Pasi Rastas
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | | | - Juhana I Kammonen
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Vicencio Oostra
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland.,Department of Evolution, Ecology and Behaviour, University of Liverpool, Liverpool CH64 7TE, UK
| | - Lorenzo Pellegrini
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Suvi Ikonen
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Tad Dallas
- Department of Biological Sciences, Louisiana State University, Baton Rouge, LA 70803, USA
| | - Michelle F DiLeo
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Anne Duplouy
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland.,Department of Biology, Lund University, 223 62 Lund, Sweden
| | - Ilhan Cem Duru
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Pauliina Halimaa
- Department of Environmental and Biological Sciences, University of Eastern Finland, 70211 KUOPIO, Finland
| | - Aapo Kahilainen
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Suyog S Kuwar
- Department of Entomology and Nematology, University of Florida, Gainesville, FL 32611-0620, USA.,Department of Zoology, Loknete Vyankatrao Hiray Arts, Science & Commerce College, 422003, Maharashtra, India
| | - Sirpa O Kärenlampi
- Department of Environmental and Biological Sciences, University of Eastern Finland, 70211 KUOPIO, Finland
| | - Elvira Lafuente
- Department of Aquatic Ecology, Swiss Federal Institute of Aquatic Science and Technology, CH-8600 Dübendorf, Switzerland
| | - Shiqi Luo
- College of Plant Protection, China Agricultural University, Beijing 100083, China
| | - Jenny Makkonen
- Department of Environmental and Biological Sciences, University of Eastern Finland, 70211 KUOPIO, Finland
| | - Abhilash Nair
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | | | - Ville Pennanen
- Viikki Plant Science Centre, Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Annukka Ruokolainen
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Tarja Sundell
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Arja I Tervahauta
- Department of Environmental and Biological Sciences, University of Eastern Finland, 70211 KUOPIO, Finland
| | - Victoria Twort
- Department of Biology, Lund University, 223 62 Lund, Sweden
| | - Erik van Bergen
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Janina Österman-Udd
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland
| | - Lars Paulin
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Mikko J Frilander
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Petri Auvinen
- Institute of Biotechnology, University of Helsinki, 00790 Helsinki, Finland
| | - Marjo Saastamoinen
- Organismal and Evolutionary Biology Research Programme, University of Helsinki, 00014 Helsinki, Finland.,Helsinki Institute of Life Science (HiLIFE), University of Helsinki, 00014 Helsinki, Finland
| |
Collapse
|
16
|
Thawornwattana Y, Seixas FA, Yang Z, Mallet J. OUP accepted manuscript. Syst Biol 2022; 71:1159-1177. [PMID: 35169847 PMCID: PMC9366460 DOI: 10.1093/sysbio/syac009] [Citation(s) in RCA: 10] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Revised: 02/01/2022] [Accepted: 02/08/2022] [Indexed: 11/21/2022] Open
Abstract
Introgressive hybridization plays a key role in adaptive evolution and species diversification in many groups of species. However, frequent hybridization and gene flow between species make estimation of the species phylogeny and key population parameters challenging. Here, we show that by accounting for phasing and using full-likelihood methods, introgression histories and population parameters can be estimated reliably from whole-genome sequence data. We employ the multispecies coalescent (MSC) model with and without gene flow to infer the species phylogeny and cross-species introgression events using genomic data from six members of the erato-sara clade of Heliconius butterflies. The methods naturally accommodate random fluctuations in genealogical history across the genome due to deep coalescence. To avoid heterozygote phasing errors in haploid sequences commonly produced by genome assembly methods, we process and compile unphased diploid sequence alignments and use analytical methods to average over uncertainties in heterozygote phase resolution. There is robust evidence for introgression across the genome, both among distantly related species deep in the phylogeny and between sister species in shallow parts of the tree. We obtain chromosome-specific estimates of key population parameters such as introgression directions, times and probabilities, as well as species divergence times and population sizes for modern and ancestral species. We confirm ancestral gene flow between the sara clade and an ancestral population of Heliconius telesiphe, a likely hybrid speciation origin for Heliconius hecalesia, and gene flow between the sister species Heliconius erato and Heliconius himera. Inferred introgression among ancestral species also explains the history of two chromosomal inversions deep in the phylogeny of the group. This study illustrates how a full-likelihood approach based on the MSC makes it possible to extract rich historical information of species divergence and gene flow from genomic data. [3s; bpp; gene flow; Heliconius; hybrid speciation; introgression; inversion; multispecies coalescent]
Collapse
Affiliation(s)
- Yuttapong Thawornwattana
- Correspondence to be sent to: Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA; E-mail: ; (Y.T. and J.M.); Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK; E-mail: (Z.Y.)
| | - Fernando A Seixas
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | - Ziheng Yang
- Correspondence to be sent to: Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA; E-mail: ; (Y.T. and J.M.); Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK; E-mail: (Z.Y.)
| | - James Mallet
- Correspondence to be sent to: Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA; E-mail: ; (Y.T. and J.M.); Department of Genetics, Evolution and Environment, University College London, London WC1E 6BT, UK; E-mail: (Z.Y.)
| |
Collapse
|
17
|
Ogilvie JG, Van Belleghem S, Range R, Papa R, McMillan OW, Chouteau M, Counterman BA. Balanced polymorphisms and their divergence in a Heliconius butterfly. Ecol Evol 2021; 11:18319-18330. [PMID: 35003675 PMCID: PMC8717333 DOI: 10.1002/ece3.8423] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2021] [Revised: 11/15/2021] [Accepted: 11/19/2021] [Indexed: 11/10/2022] Open
Abstract
The evolution of mimicry in similarly defended prey is well described by the Müllerian mimicry theory, which predicts the convergence of warning patterns in order to gain the most protection from predators. However, despite this prediction, we can find great diversity of color patterns among Müllerian mimics such as Heliconius butterflies in the neotropics. Furthermore, some species have evolved the ability to maintain multiple distinct warning patterns in single populations, a phenomenon known as polymorphic mimicry. The adaptive benefit of these polymorphisms is questionable since variation from the most common warning patterns is expected to be disadvantageous as novel signals are punished by predators naive to them. In this study, we use artificial butterfly models throughout Central and South America to characterize the selective pressures maintaining polymorphic mimicry in Heliconius doris. Our results highlight the complexity of positive frequency-dependent selection, the principal selective pressure driving convergence among Müllerian mimics, and its impacts on interspecific variation of mimetic warning coloration. We further show how this selection regime can both limit and facilitate the diversification of mimetic traits.
Collapse
Affiliation(s)
- James G. Ogilvie
- Department of Biological SciencesAuburn UniversityAuburnAlabamaUSA
- Smithsonian Tropical Research InstitutePanama CityPanama
| | | | - Ryan Range
- Department of Biological SciencesAuburn UniversityAuburnAlabamaUSA
| | - Riccardo Papa
- Department of BiologyUniversity of Puerto RicoRio PiedrasPuerto Rico
| | | | - Mathieu Chouteau
- Laboratoire Écologie, Évolution, Interactions des Systèmes Amazoniens (LEEISA)Université de GuyaneCNRSIFREMERCayenneFrench Guiana
| | | |
Collapse
|
18
|
Ferrari G, Atmore LM, Jentoft S, Jakobsen KS, Makowiecki D, Barrett JH, Star B. An accurate assignment test for extremely low-coverage whole-genome sequence data. Mol Ecol Resour 2021; 22:1330-1344. [PMID: 34779123 DOI: 10.1111/1755-0998.13551] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/23/2021] [Revised: 10/28/2021] [Accepted: 11/02/2021] [Indexed: 11/28/2022]
Abstract
Genomic assignment tests can provide important diagnostic biological characteristics, such as population of origin or ecotype. Yet, assignment tests often rely on moderate- to high-coverage sequence data that can be difficult to obtain for fields such as molecular ecology and ancient DNA. We have developed a novel approach that efficiently assigns biologically relevant information (i.e., population identity or structural variants such as inversions) in extremely low-coverage sequence data. First, we generate databases from existing reference data using a subset of diagnostic single nucleotide polymorphisms (SNPs) associated with a biological characteristic. Low-coverage alignment files are subsequently compared to these databases to ascertain allelic state, yielding a joint probability for each association. To assess the efficacy of this approach, we assigned haplotypes and population identity in Heliconius butterflies, Atlantic herring, and Atlantic cod using chromosomal inversion sites and whole-genome data. We scored both modern and ancient specimens, including the first whole-genome sequence data recovered from ancient Atlantic herring bones. The method accurately assigns biological characteristics, including population membership, using extremely low-coverage data (as low as 0.0001x) based on genome-wide SNPs. This approach will therefore increase the number of samples in evolutionary, ecological and archaeological research for which relevant biological information can be obtained.
Collapse
Affiliation(s)
- Giada Ferrari
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Lane M Atmore
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Sissel Jentoft
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Kjetill S Jakobsen
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| | - Daniel Makowiecki
- Department of Environmental Archaeology and Human Paleoecology, Institute of Archaeology, Nicolaus Copernicus University, Torun, Poland
| | - James H Barrett
- McDonald Institute for Archaeological Research, Department of Archaeology, University of Cambridge, Cambridge, UK.,Department of Archaeology and Cultural History, NTNU University Museum, Trondheim, Norway
| | - Bastiaan Star
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo, Oslo, Norway
| |
Collapse
|
19
|
Rosser N, Edelman NB, Queste LM, Nelson M, Seixas F, Dasmahapatra KK, Mallet J. Complex basis of hybrid female sterility and Haldane's rule in Heliconius butterflies: Z-linkage and epistasis. Mol Ecol 2021; 31:959-977. [PMID: 34779079 DOI: 10.1111/mec.16272] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2021] [Revised: 09/14/2021] [Accepted: 09/27/2021] [Indexed: 12/19/2022]
Abstract
Hybrids between species are often sterile or inviable. Hybrid unfitness usually evolves first in the heterogametic sex-a pattern known as Haldane's rule. The genetics of Haldane's rule have been extensively studied in species where the male is the heterogametic (XX/XY) sex, but its basis in taxa where the female is heterogametic (ZW/ZZ), such as Lepidoptera and birds, is largely unknown. Here, we analyse a new case of female hybrid sterility between geographic subspecies of Heliconius pardalinus. The two subspecies mate freely in captivity, but female F1 hybrids in both directions of cross are sterile. Sterility is due to arrested development of oocytes after they become differentiated from nurse cells, but before yolk deposition. We backcrossed fertile male F1 hybrids to parental females and mapped quantitative trait loci (QTLs) for female sterility. We also identified genes differentially expressed in the ovary as a function of oocyte development. The Z chromosome has a major effect, similar to the 'large X effect' in Drosophila, with strong epistatic interactions between loci at either end of the Z chromosome, and between the Z chromosome and autosomal loci on chromosomes 8 and 20. By intersecting the list of genes within these QTLs with those differentially expressed in sterile and fertile hybrids, we identified three candidate genes with relevant phenotypes. This study is the first to characterize hybrid sterility using genome mapping in the Lepidoptera and shows that it is produced by multiple complex epistatic interactions often involving the sex chromosome, as predicted by the dominance theory of Haldane's rule.
Collapse
Affiliation(s)
- Neil Rosser
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA.,Department of Biology, University of York, York, UK
| | - Nathaniel B Edelman
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA.,Yale Institute for Biospheric Studies, Yale University, New Haven, Connecticut, USA.,Yale School for the Environment, Yale University, New Haven, Connecticut, USA
| | | | | | - Fernando Seixas
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
| | | | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
| |
Collapse
|
20
|
Blischak PD, Barker MS, Gutenkunst RN. Chromosome-scale inference of hybrid speciation and admixture with convolutional neural networks. Mol Ecol Resour 2021; 21:2676-2688. [PMID: 33682305 PMCID: PMC8675098 DOI: 10.1111/1755-0998.13355] [Citation(s) in RCA: 7] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/30/2020] [Revised: 01/26/2021] [Accepted: 02/05/2021] [Indexed: 11/30/2022]
Abstract
Inferring the frequency and mode of hybridization among closely related organisms is an important step for understanding the process of speciation and can help to uncover reticulated patterns of phylogeny more generally. Phylogenomic methods to test for the presence of hybridization come in many varieties and typically operate by leveraging expected patterns of genealogical discordance in the absence of hybridization. An important assumption made by these tests is that the data (genes or SNPs) are independent given the species tree. However, when the data are closely linked, it is especially important to consider their nonindependence. Recently, deep learning techniques such as convolutional neural networks (CNNs) have been used to perform population genetic inferences with linked SNPs coded as binary images. Here, we use CNNs for selecting among candidate hybridization scenarios using the tree topology (((P1 , P2 ), P3 ), Out) and a matrix of pairwise nucleotide divergence (dXY ) calculated in windows across the genome. Using coalescent simulations to train and independently test a neural network showed that our method, HyDe-CNN, was able to accurately perform model selection for hybridization scenarios across a wide breath of parameter space. We then used HyDe-CNN to test models of admixture in Heliconius butterflies, as well as comparing it to phylogeny-based introgression statistics. Given the flexibility of our approach, the dropping cost of long-read sequencing and the continued improvement of CNN architectures, we anticipate that inferences of hybridization using deep learning methods like ours will help researchers to better understand patterns of admixture in their study organisms.
Collapse
Affiliation(s)
- Paul D. Blischak
- Department of Ecology & Evolutionary Biology, University of Arizona, Tucson, AZ, 85721, USA
- Department of Molecular & Cellular Biology, University of Arizona, Tucson, AZ, 85721, USA
| | - Michael S. Barker
- Department of Ecology & Evolutionary Biology, University of Arizona, Tucson, AZ, 85721, USA
| | - Ryan N. Gutenkunst
- Department of Molecular & Cellular Biology, University of Arizona, Tucson, AZ, 85721, USA
| |
Collapse
|
21
|
Mongue AJ, Hansen ME, Walters JR. Support for faster and more adaptive Z chromosome evolution in two divergent lepidopteran lineages. Evolution 2021; 76:332-345. [PMID: 34463346 PMCID: PMC9291949 DOI: 10.1111/evo.14341] [Citation(s) in RCA: 12] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/14/2021] [Revised: 07/21/2021] [Accepted: 08/06/2021] [Indexed: 12/14/2022]
Abstract
The rateof divergence for Z or X chromosomes is usually observed to be greater than autosomes, but the proposed evolutionary causes for this pattern vary, as do empirical results from diverse taxa. Even among moths and butterflies (Lepidoptera), which generally share a single-origin Z chromosome, the handful of available studies give mixed support for faster or more adaptive evolution of the Z chromosome, depending on the species assayed. Here, we examine the molecular evolution of Z chromosomes in two additional lepidopteran species: the Carolina sphinx moth and the monarch butterfly, the latter of which possesses a recent chromosomal fusion yielding a segment of newly Z-linked DNA. We find evidence for both faster and more adaptive Z chromosome evolution in both species, although this effect is strongest in the neo-Z portion of the monarch sex chromosome. The neo-Z is less male-biased than expected of a Z chromosome, and unbiased and female-biased genes drive the signal for adaptive evolution here. Together these results suggest that male-biased gene accumulation and haploid selection have opposing effects on long-term rates of adaptation and may help explain the discrepancies in previous findings as well as the repeated evolution of neo-sex chromosomes in Lepidoptera.
Collapse
Affiliation(s)
- Andrew J Mongue
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas, 66045.,Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH93FL, United Kingdom
| | - Megan E Hansen
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas, 66045
| | - James R Walters
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas, 66045
| |
Collapse
|
22
|
Manthey JD, Klicka J, Spellman GM. The Genomic Signature of Allopatric Speciation in a Songbird Is Shaped by Genome Architecture (Aves: Certhia americana). Genome Biol Evol 2021; 13:evab120. [PMID: 34042960 PMCID: PMC8364988 DOI: 10.1093/gbe/evab120] [Citation(s) in RCA: 6] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/24/2021] [Indexed: 12/31/2022] Open
Abstract
The genomic signature of speciation with gene flow is often attributed to the strength of divergent selection and recombination rate in regions harboring targets for selection. In contrast, allopatric speciation provides a different geographic context and evolutionary scenario, whereby introgression is limited by isolation rather than selection against gene flow. Lacking shared divergent selection or selection against hybridization, we would predict the genomic signature of allopatric speciation would largely be shaped by genomic architecture-the nonrandom distribution of functional elements and chromosomal characteristics-through its role in affecting the processes of selection and drift. Here, we built and annotated a chromosome-scale genome assembly for a songbird (Passeriformes: Certhia americana). We show that the genomic signature of allopatric speciation between its two primary lineages is largely shaped by genomic architecture. Regionally, gene density and recombination rate variation explain a large proportion of variance in genomic diversity, differentiation, and divergence. We identified a heterogeneous landscape of selection and neutrality, with a large portion of the genome under the effects of indirect selection. We found higher proportions of small chromosomes under the effects of indirect selection, likely because they have relatively higher gene density. At the chromosome scale, differential genomic architecture of macro- and microchromosomes shapes the genomic signatures of speciation: chromosome size has: 1) a positive relationship with genetic differentiation, genetic divergence, rate of lineage sorting in the contact zone, and proportion neutral evolution and 2) a negative relationship with genetic diversity and recombination rate.
Collapse
Affiliation(s)
- Joseph D Manthey
- Department of Biological Sciences, Texas Tech University, Lubbock, Texas, USA
| | - John Klicka
- Burke Museum of Natural History, University of Washington, Seattle, Washington, USA
- Department of Biology, University of Washington, Seattle, Washington, USA
| | - Garth M Spellman
- Department of Zoology, Denver Museum of Nature & Science, Denver, Colorado, USA
| |
Collapse
|
23
|
How Important Are Structural Variants for Speciation? Genes (Basel) 2021; 12:genes12071084. [PMID: 34356100 PMCID: PMC8305853 DOI: 10.3390/genes12071084] [Citation(s) in RCA: 21] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Revised: 07/04/2021] [Accepted: 07/14/2021] [Indexed: 12/11/2022] Open
Abstract
Understanding the genetic basis of reproductive isolation is a central issue in the study of speciation. Structural variants (SVs); that is, structural changes in DNA, including inversions, translocations, insertions, deletions, and duplications, are common in a broad range of organisms and have been hypothesized to play a central role in speciation. Recent advances in molecular and statistical methods have identified structural variants, especially inversions, underlying ecologically important traits; thus, suggesting these mutations contribute to adaptation. However, the contribution of structural variants to reproductive isolation between species—and the underlying mechanism by which structural variants most often contribute to speciation—remain unclear. Here, we review (i) different mechanisms by which structural variants can generate or maintain reproductive isolation; (ii) patterns expected with these different mechanisms; and (iii) relevant empirical examples of each. We also summarize the available sequencing and bioinformatic methods to detect structural variants. Lastly, we suggest empirical approaches and new research directions to help obtain a more complete assessment of the role of structural variants in speciation.
Collapse
|
24
|
Haplotype tagging reveals parallel formation of hybrid races in two butterfly species. Proc Natl Acad Sci U S A 2021; 118:2015005118. [PMID: 34155138 PMCID: PMC8237668 DOI: 10.1073/pnas.2015005118] [Citation(s) in RCA: 33] [Impact Index Per Article: 11.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/18/2022] Open
Abstract
A defining goal in genetics is linking variation in DNA sequence to trait evolution between populations and, ultimately, species. Genome sequencing efficiently captures such variation but typically in millions of tiny fragments that omit haplotype or linkage information. We present “haplotagging,” a simple, rapid linked-read sequencing technique that allows high-throughput sequencing without sacrificing haplotype information. We validated this affordable approach for whole-genome haplotyping in large populations. We used haplotagging to investigate the rise of a novel hybrid morph in parallel hybrid zones of two comimetic Heliconius butterfly species in Ecuador. Our results reveal that strikingly parallel divergences in their genomes produced coordinated shifts in haplotype frequencies across the hybrid zone, giving rise to comimetic hybrid morphs in each species. Genetic variation segregates as linked sets of variants or haplotypes. Haplotypes and linkage are central to genetics and underpin virtually all genetic and selection analysis. Yet, genomic data often omit haplotype information due to constraints in sequencing technologies. Here, we present “haplotagging,” a simple, low-cost linked-read sequencing technique that allows sequencing of hundreds of individuals while retaining linkage information. We apply haplotagging to construct megabase-size haplotypes for over 600 individual butterflies (Heliconius erato and H. melpomene), which form overlapping hybrid zones across an elevational gradient in Ecuador. Haplotagging identifies loci controlling distinctive high- and lowland wing color patterns. Divergent haplotypes are found at the same major loci in both species, while chromosome rearrangements show no parallelism. Remarkably, in both species, the geographic clines for the major wing-pattern loci are displaced by 18 km, leading to the rise of a novel hybrid morph in the center of the hybrid zone. We propose that shared warning signaling (Müllerian mimicry) may couple the cline shifts seen in both species and facilitate the parallel coemergence of a novel hybrid morph in both comimetic species. Our results show the power of efficient haplotyping methods when combined with large-scale sequencing data from natural populations.
Collapse
|
25
|
Boman J, Mugal CF, Backström N. The Effects of GC-Biased Gene Conversion on Patterns of Genetic Diversity among and across Butterfly Genomes. Genome Biol Evol 2021; 13:evab064. [PMID: 33760095 PMCID: PMC8175052 DOI: 10.1093/gbe/evab064] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/22/2021] [Indexed: 12/28/2022] Open
Abstract
Recombination reshuffles the alleles of a population through crossover and gene conversion. These mechanisms have considerable consequences on the evolution and maintenance of genetic diversity. Crossover, for example, can increase genetic diversity by breaking the linkage between selected and nearby neutral variants. Bias in favor of G or C alleles during gene conversion may instead promote the fixation of one allele over the other, thus decreasing diversity. Mutation bias from G or C to A and T opposes GC-biased gene conversion (gBGC). Less recognized is that these two processes may-when balanced-promote genetic diversity. Here, we investigate how gBGC and mutation bias shape genetic diversity patterns in wood white butterflies (Leptidea sp.). This constitutes the first in-depth investigation of gBGC in butterflies. Using 60 resequenced genomes from six populations of three species, we find substantial variation in the strength of gBGC across lineages. When modeling the balance of gBGC and mutation bias and comparing analytical results with empirical data, we reject gBGC as the main determinant of genetic diversity in these butterfly species. As alternatives, we consider linked selection and GC content. We find evidence that high values of both reduce diversity. We also show that the joint effects of gBGC and mutation bias can give rise to a diversity pattern which resembles the signature of linked selection. Consequently, gBGC should be considered when interpreting the effects of linked selection on levels of genetic diversity.
Collapse
Affiliation(s)
- Jesper Boman
- Evolutionary Biology Program, Department of Ecology and Genetics (IEG), Uppsala University, Sweden
| | - Carina F Mugal
- Evolutionary Biology Program, Department of Ecology and Genetics (IEG), Uppsala University, Sweden
| | - Niclas Backström
- Evolutionary Biology Program, Department of Ecology and Genetics (IEG), Uppsala University, Sweden
| |
Collapse
|
26
|
Stankowski S, Ravinet M. Defining the speciation continuum. Evolution 2021; 75:1256-1273. [PMID: 33754340 DOI: 10.1111/evo.14215] [Citation(s) in RCA: 48] [Impact Index Per Article: 16.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/01/2020] [Revised: 03/01/2021] [Accepted: 03/13/2021] [Indexed: 02/06/2023]
Abstract
A primary roadblock to our understanding of speciation is that it usually occurs over a timeframe that is too long to study from start to finish. The idea of a speciation continuum provides something of a solution to this problem; rather than observing the entire process, we can simply reconstruct it from the multitude of speciation events that surround us. But what do we really mean when we talk about the speciation continuum, and can it really help us understand speciation? We explored these questions using a literature review and online survey of speciation researchers. Although most researchers were familiar with the concept and thought it was useful, our survey revealed extensive disagreement about what the speciation continuum actually tells us. This is due partly to the lack of a clear definition. Here, we provide an explicit definition that is compatible with the Biological Species Concept. That is, the speciation continuum is a continuum of reproductive isolation. After outlining the logic of the definition in light of alternatives, we explain why attempts to reconstruct the speciation process from present-day populations will ultimately fail. We then outline how we think the speciation continuum concept can continue to act as a foundation for understanding the continuum of reproductive isolation that surrounds us.
Collapse
Affiliation(s)
- Sean Stankowski
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, S10 2TN, United Kingdom.,Current Address: Sean Stankowski, IST Austria, Klosterneuburg, 3400, Austria
| | - Mark Ravinet
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Oslo, 0316, Norway.,School of Life Sciences, University of Nottingham, Nottingham, NG7 2RD, United Kingdom
| |
Collapse
|
27
|
Seixas FA, Edelman NB, Mallet J. Synteny-Based Genome Assembly for 16 Species of Heliconius Butterflies, and an Assessment of Structural Variation across the Genus. Genome Biol Evol 2021; 13:6207971. [PMID: 33792688 PMCID: PMC8290116 DOI: 10.1093/gbe/evab069] [Citation(s) in RCA: 9] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 03/29/2021] [Indexed: 12/11/2022] Open
Abstract
Heliconius butterflies (Lepidoptera: Nymphalidae) are a group of 48 neotropical species widely studied in evolutionary research. Despite the wealth of genomic data generated in past years, chromosomal level genome assemblies currently exist for only two species, Heliconius melpomene and Heliconius erato, each a representative of one of the two major clades of the genus. Here, we use these reference genomes to improve the contiguity of previously published draft genome assemblies of 16 Heliconius species. Using a reference-assisted scaffolding approach, we place and order the scaffolds of these genomes onto chromosomes, resulting in 95.7-99.9% of their genomes anchored to chromosomes. Genome sizes are somewhat variable among species (270-422 Mb) and in one small group of species (Heliconius hecale, Heliconius elevatus, and Heliconius pardalinus) expansions in genome size are driven mainly by repetitive sequences that map to four small regions in the H. melpomene reference genome. Genes from these repeat regions show an increase in exon copy number, an absence of internal stop codons, evidence of constraint on nonsynonymous changes, and increased expression, all of which suggest that at least some of the extra copies are functional. Finally, we conducted a systematic search for inversions and identified five moderately large inversions fixed between the two major Heliconius clades. We infer that one of these inversions was transferred by introgression between the lineages leading to the erato/sara and burneyi/doris clades. These reference-guided assemblies represent a major improvement in Heliconius genomic resources that enable further genetic and evolutionary discoveries in this genus.
Collapse
Affiliation(s)
- Fernando A Seixas
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
| | - Nathaniel B Edelman
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA.,Yale Institute for Biospheric Studies, Yale University, New Haven, Connecticut, USA
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
| |
Collapse
|
28
|
Van Belleghem SM, Lewis JJ, Rivera ES, Papa R. Heliconius butterflies: a window into the evolution and development of diversity. Curr Opin Genet Dev 2021; 69:72-81. [PMID: 33714874 DOI: 10.1016/j.gde.2021.01.010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/17/2020] [Revised: 01/14/2021] [Accepted: 01/25/2021] [Indexed: 10/21/2022]
Abstract
Butterflies have become prominent models for studying the evolution and development of phenotypic variation. In Heliconius, extraordinary within species divergence and between species convergence in wing color patterns has driven decades of comparative genetic studies. However, connecting genetic patterns of diversification to the molecular mechanisms of adaptation has remained elusive. Recent studies are bridging this gap between genome and function and have driven substantial advances in deciphering the genetic architecture of diversification in Heliconius. While only a handful of large-effect genes were initially identified in the diversification of Heliconius color patterns, recent experiments have begun to unravel the underlying gene regulatory networks and how these have evolved. These results reveal an evolutionary story of many interacting loci and partly independent genetic architectures that underlie convergent evolution.
Collapse
Affiliation(s)
| | - James J Lewis
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY, USA; Baker Institute for Animal Health, Cornell University, Ithaca, NY, USA
| | - Edgardo S Rivera
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan, Puerto Rico; Chairs of Biomaterials, University of Bayreuth, Bayreuth, Bayern, Germany
| | - Riccardo Papa
- Department of Biology, University of Puerto Rico-Rio Piedras, San Juan, Puerto Rico; Molecular Sciences and Research Center, University of Puerto Rico, San Juan, Puerto Rico.
| |
Collapse
|
29
|
Darragh K, Orteu A, Black D, Byers KJRP, Szczerbowski D, Warren IA, Rastas P, Pinharanda A, Davey JW, Fernanda Garza S, Abondano Almeida D, Merrill RM, McMillan WO, Schulz S, Jiggins CD. A novel terpene synthase controls differences in anti-aphrodisiac pheromone production between closely related Heliconius butterflies. PLoS Biol 2021; 19:e3001022. [PMID: 33465061 PMCID: PMC7815096 DOI: 10.1371/journal.pbio.3001022] [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: 06/12/2020] [Accepted: 11/30/2020] [Indexed: 02/07/2023] Open
Abstract
Plants and insects often use the same compounds for chemical communication, but not much is known about the genetics of convergent evolution of chemical signals. The terpene (E)-β-ocimene is a common component of floral scent and is also used by the butterfly Heliconius melpomene as an anti-aphrodisiac pheromone. While the biosynthesis of terpenes has been described in plants and microorganisms, few terpene synthases (TPSs) have been identified in insects. Here, we study the recent divergence of 2 species, H. melpomene and Heliconius cydno, which differ in the presence of (E)-β-ocimene; combining linkage mapping, gene expression, and functional analyses, we identify 2 novel TPSs. Furthermore, we demonstrate that one, HmelOS, is able to synthesise (E)-β-ocimene in vitro. We find no evidence for TPS activity in HcydOS (HmelOS ortholog of H. cydno), suggesting that the loss of (E)-β-ocimene in this species is the result of coding, not regulatory, differences. The TPS enzymes we discovered are unrelated to previously described plant and insect TPSs, demonstrating that chemical convergence has independent evolutionary origins. Plants and insects often use the same compounds for chemical communication, but little is known about the convergent evolution of such chemical signals. This study identifies a novel terpene synthase involved in production of an anti-aphrodisiac pheromone by the butterfly Heliconius melpomene. This enzyme is unrelated to other insect terpene synthases, providing evidence that the ability to synthesise terpenes has arisen multiple times independently within the insects.
Collapse
Affiliation(s)
- Kathy Darragh
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Smithsonian Tropical Research Institute, Panamá, Panamá
- * E-mail:
| | - Anna Orteu
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Daniella Black
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- School of Biology, University of St Andrews, St Andrews, United Kingdom
| | - Kelsey J. R. P. Byers
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Smithsonian Tropical Research Institute, Panamá, Panamá
| | - Daiane Szczerbowski
- Institute of Organic Chemistry, Department of Life Sciences, Technische Universität Braunschweig, Braunschweig, Germany
| | - Ian A. Warren
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Pasi Rastas
- Institute of Biotechnology, University of Helsinki, Helsinki, Finland
| | - Ana Pinharanda
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - John W. Davey
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | | | | | - Richard M. Merrill
- Smithsonian Tropical Research Institute, Panamá, Panamá
- Division of Evolutionary Biology, Ludwig-Maximilians-Universität München, Munich, Germany
| | | | - Stefan Schulz
- Institute of Organic Chemistry, Department of Life Sciences, Technische Universität Braunschweig, Braunschweig, Germany
| | - Chris D. Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Smithsonian Tropical Research Institute, Panamá, Panamá
| |
Collapse
|
30
|
Byers KJRP, Darragh K, Fernanda Garza S, Abondano Almeida D, Warren IA, Rastas PMA, Merrill RM, Schulz S, McMillan WO, Jiggins CD. Clustering of loci controlling species differences in male chemical bouquets of sympatric Heliconius butterflies. Ecol Evol 2021; 11:89-107. [PMID: 33437416 PMCID: PMC7790645 DOI: 10.1002/ece3.6947] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/23/2020] [Revised: 09/29/2020] [Accepted: 09/30/2020] [Indexed: 12/23/2022] Open
Abstract
The degree to which loci promoting reproductive isolation cluster in the genome-that is, the genetic architecture of reproductive isolation-can influence the tempo and mode of speciation. Tight linkage between these loci can facilitate speciation in the face of gene flow. Pheromones play a role in reproductive isolation in many Lepidoptera species, and the role of endogenously produced compounds as secondary metabolites decreases the likelihood of pleiotropy associated with many barrier loci. Heliconius butterflies use male sex pheromones to both court females (aphrodisiac wing pheromones) and ward off male courtship (male-transferred antiaphrodisiac genital pheromones), and it is likely that these compounds play a role in reproductive isolation between Heliconius species. Using a set of backcross hybrids between H. melpomene and H. cydno, we investigated the genetic architecture of putative male pheromone compound production. We found a set of 40 significant quantitative trait loci (QTL) representing 33 potential pheromone compounds. QTL clustered significantly on two chromosomes, chromosome 8 for genital compounds and chromosome 20 for wing compounds, and chromosome 20 was enriched for potential pheromone biosynthesis genes. There was minimal overlap between pheromone QTL and known QTL for mate choice and color pattern. Nonetheless, we did detect linkage between a QTL for wing androconial area and optix, a color pattern locus known to play a role in reproductive isolation in these species. This tight clustering of putative pheromone loci might contribute to coincident reproductive isolating barriers, facilitating speciation despite ongoing gene flow.
Collapse
Affiliation(s)
- Kelsey J. R. P. Byers
- Department of ZoologyUniversity of CambridgeCambridgeUK
- Smithsonian Tropical Research InstitutePanamaPanama
- Present address:
Department of Cell and Developmental BiologyJohn Innes CentreNorwichUK
| | - Kathy Darragh
- Department of ZoologyUniversity of CambridgeCambridgeUK
- Smithsonian Tropical Research InstitutePanamaPanama
- Present address:
Department of Evolution and EcologyUniversity of California DavisDavisCAUSA
| | - Sylvia Fernanda Garza
- Smithsonian Tropical Research InstitutePanamaPanama
- Present address:
Department of Collective BehaviourMax Planck Institute of Animal BehaviourKonstanzGermany
| | - Diana Abondano Almeida
- Smithsonian Tropical Research InstitutePanamaPanama
- Present address:
Institute for Ecology, Evolution and DiversityGoethe UniversitätFrankfurtGermany
| | - Ian A. Warren
- Department of ZoologyUniversity of CambridgeCambridgeUK
| | | | - Richard M. Merrill
- Smithsonian Tropical Research InstitutePanamaPanama
- Division of Evolutionary BiologyLudwig‐Maximilians‐Universität MünchenMunichGermany
| | - Stefan Schulz
- Institute of Organic ChemistryDepartment of Life SciencesTechnische Universität BraunschweigBraunschweigGermany
| | | | - Chris D. Jiggins
- Department of ZoologyUniversity of CambridgeCambridgeUK
- Smithsonian Tropical Research InstitutePanamaPanama
| |
Collapse
|
31
|
Bainbridge HE, Brien MN, Morochz C, Salazar PA, Rastas P, Nadeau NJ. Limited genetic parallels underlie convergent evolution of quantitative pattern variation in mimetic butterflies. J Evol Biol 2020; 33:1516-1529. [DOI: 10.1111/jeb.13704] [Citation(s) in RCA: 12] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/11/2020] [Revised: 08/05/2020] [Accepted: 09/04/2020] [Indexed: 01/28/2023]
Affiliation(s)
- Hannah E. Bainbridge
- Department of Animal and Plant Sciences The University of Sheffield Sheffield UK
| | - Melanie N. Brien
- Department of Animal and Plant Sciences The University of Sheffield Sheffield UK
| | - Carlos Morochz
- Biology & Research Department Mashpi Lodge Mashpi Ecuador
| | - Patricio A. Salazar
- Department of Animal and Plant Sciences The University of Sheffield Sheffield UK
| | - Pasi Rastas
- Institute of Biotechnology University of Helsinki Helsinki Finland
| | - Nicola J. Nadeau
- Department of Animal and Plant Sciences The University of Sheffield Sheffield UK
| |
Collapse
|
32
|
Rivera-Colón AG, Rochette NC, Catchen JM. Simulation with RADinitio improves RADseq experimental design and sheds light on sources of missing data. Mol Ecol Resour 2020; 21:363-378. [PMID: 32275349 DOI: 10.1111/1755-0998.13163] [Citation(s) in RCA: 20] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [Key Words] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/10/2019] [Accepted: 03/25/2020] [Indexed: 12/20/2022]
Abstract
Restriction-site associated DNA sequencing (RADseq) has become a powerful and versatile tool in modern population genomics, enabling large-scale evolutionary and genomic analyses in otherwise inaccessible biological systems. With its widespread use, different variants on the protocol have been developed to suit specific experimental needs. Researchers face the challenge of choosing the optimal molecular and sequencing protocols for their reduced representation experimental design, an often-complicated process. Strategic errors can lead to biased data generation that has reduced power to answer biological questions. Here, we present RADinitio, simulation software for the selection and optimization of RADseq experiments via the generation of sequencing data that behave similarly to empirical sources. RADinitio provides an evolutionary simulation of populations, implementation of various RADseq protocols with customizable parameters, and thorough assessment of missing data. We test the efficacy of the software using different RAD protocols across several organisms, highlighting the importance of protocol selection on the magnitude and quality of data acquired. Additionally, we test the effects of RAD library preparation and sequencing on allelic dropout, observing that library preparation and sequencing often contributes more to missing alleles than population-level variation.
Collapse
Affiliation(s)
- Angel G Rivera-Colón
- Department of Evolution, Ecology and Behavior, University of Illinois, Urbana, Illinois, USA
| | - Nicolas C Rochette
- Department of Evolution, Ecology and Behavior, University of Illinois, Urbana, Illinois, USA
| | - Julian M Catchen
- Department of Evolution, Ecology and Behavior, University of Illinois, Urbana, Illinois, USA
| |
Collapse
|
33
|
Darragh K, Montejo‐Kovacevich G, Kozak KM, Morrison CR, Figueiredo CME, Ready JS, Salazar C, Linares M, Byers KJRP, Merrill RM, McMillan WO, Schulz S, Jiggins CD. Species specificity and intraspecific variation in the chemical profiles of Heliconius butterflies across a large geographic range. Ecol Evol 2020; 10:3895-3918. [PMID: 32489619 PMCID: PMC7244815 DOI: 10.1002/ece3.6079] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2019] [Revised: 01/06/2020] [Accepted: 01/07/2020] [Indexed: 02/01/2023] Open
Abstract
In many animals, mate choice is important for the maintenance of reproductive isolation between species. Traits important for mate choice and behavioral isolation are predicted to be under strong stabilizing selection within species; however, such traits can also exhibit variation at the population level driven by neutral and adaptive evolutionary processes. Here, we describe patterns of divergence among androconial and genital chemical profiles at inter- and intraspecific levels in mimetic Heliconius butterflies. Most variation in chemical bouquets was found between species, but there were also quantitative differences at the population level. We found a strong correlation between interspecific chemical and genetic divergence, but this correlation varied in intraspecific comparisons. We identified "indicator" compounds characteristic of particular species that included compounds already known to elicit a behavioral response, suggesting an approach for identification of candidate compounds for future behavioral studies in novel systems. Overall, the strong signal of species identity suggests a role for these compounds in species recognition, but with additional potentially neutral variation at the population level.
Collapse
Affiliation(s)
- Kathy Darragh
- Department of ZoologyUniversity of CambridgeCambridgeUK
- Smithsonian Tropical Research InstitutePanama CityPanama
| | | | | | - Colin R. Morrison
- Smithsonian Tropical Research InstitutePanama CityPanama
- Department of Integrative BiologyThe University of Texas at AustinAustinTXUSA
| | | | - Jonathan S. Ready
- Institute for Biological SciencesUniversidade Federal do ParáBelémBrazil
| | - Camilo Salazar
- Biology ProgramFaculty of Natural Sciences and MathematicsUniversidad del RosarioBogotaColombia
| | - Mauricio Linares
- Biology ProgramFaculty of Natural Sciences and MathematicsUniversidad del RosarioBogotaColombia
| | - Kelsey J. R. P. Byers
- Department of ZoologyUniversity of CambridgeCambridgeUK
- Smithsonian Tropical Research InstitutePanama CityPanama
| | - Richard M. Merrill
- Smithsonian Tropical Research InstitutePanama CityPanama
- Division of Evolutionary BiologyFaculty of BiologyLudwig‐Maximilians‐Universität MünchenMunichGermany
| | | | - Stefan Schulz
- Institute of Organic ChemistryTechnische Universität BraunschweigBraunschweigGermany
| | - Chris D. Jiggins
- Department of ZoologyUniversity of CambridgeCambridgeUK
- Smithsonian Tropical Research InstitutePanama CityPanama
| |
Collapse
|
34
|
Martin SH, Singh KS, Gordon IJ, Omufwoko KS, Collins S, Warren IA, Munby H, Brattström O, Traut W, Martins DJ, Smith DAS, Jiggins CD, Bass C, ffrench-Constant RH. Whole-chromosome hitchhiking driven by a male-killing endosymbiont. PLoS Biol 2020; 18:e3000610. [PMID: 32108180 PMCID: PMC7046192 DOI: 10.1371/journal.pbio.3000610] [Citation(s) in RCA: 34] [Impact Index Per Article: 8.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/05/2019] [Accepted: 01/23/2020] [Indexed: 12/30/2022] Open
Abstract
Neo-sex chromosomes are found in many taxa, but the forces driving their emergence and spread are poorly understood. The female-specific neo-W chromosome of the African monarch (or queen) butterfly Danaus chrysippus presents an intriguing case study because it is restricted to a single 'contact zone' population, involves a putative colour patterning supergene, and co-occurs with infection by the male-killing endosymbiont Spiroplasma. We investigated the origin and evolution of this system using whole genome sequencing. We first identify the 'BC supergene', a broad region of suppressed recombination across nearly half a chromosome, which links two colour patterning loci. Association analysis suggests that the genes yellow and arrow in this region control the forewing colour pattern differences between D. chrysippus subspecies. We then show that the same chromosome has recently formed a neo-W that has spread through the contact zone within approximately 2,200 years. We also assembled the genome of the male-killing Spiroplasma, and find that it shows perfect genealogical congruence with the neo-W, suggesting that the neo-W has hitchhiked to high frequency as the male-killer has spread through the population. The complete absence of female crossing-over in the Lepidoptera causes whole-chromosome hitchhiking of a single neo-W haplotype, carrying a single allele of the BC supergene and dragging multiple non-synonymous mutations to high frequency. This has created a population of infected females that all carry the same recessive colour patterning allele, making the phenotypes of each successive generation highly dependent on uninfected male immigrants. Our findings show how hitchhiking can occur between the physically unlinked genomes of host and endosymbiont, with dramatic consequences.
Collapse
Affiliation(s)
- Simon H. Martin
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Kumar Saurabh Singh
- Centre for Ecology and Conservation, University of Exeter, Penryn Campus, Penryn, United Kingdom
| | - Ian J. Gordon
- Center of Excellence in Biodiversity and Natural Resource Management, University of Rwanda, Huye, Rwanda
| | - Kennedy Saitoti Omufwoko
- Department of Ecology and Evolutionary Biology, Princeton University, Princeton, United States of America
- Mpala Research Centre, Nanyuki, Kenya
| | - Steve Collins
- African Butterfly Research Institute, Nairobi, Kenya
| | - Ian A. Warren
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Hannah Munby
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Oskar Brattström
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Walther Traut
- Institut für Biologie, Universität Lübeck, Lübeck, Germany
| | - Dino J. Martins
- Department of Ecology and Evolutionary Biology, Princeton University, Princeton, United States of America
- Mpala Research Centre, Nanyuki, Kenya
| | | | - Chris D. Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Chris Bass
- Centre for Ecology and Conservation, University of Exeter, Penryn Campus, Penryn, United Kingdom
| | | |
Collapse
|
35
|
Moest M, Van Belleghem SM, James JE, Salazar C, Martin SH, Barker SL, Moreira GRP, Mérot C, Joron M, Nadeau NJ, Steiner FM, Jiggins CD. Selective sweeps on novel and introgressed variation shape mimicry loci in a butterfly adaptive radiation. PLoS Biol 2020; 18:e3000597. [PMID: 32027643 PMCID: PMC7029882 DOI: 10.1371/journal.pbio.3000597] [Citation(s) in RCA: 38] [Impact Index Per Article: 9.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/03/2019] [Revised: 02/19/2020] [Accepted: 01/15/2020] [Indexed: 11/21/2022] Open
Abstract
Natural selection leaves distinct signatures in the genome that can reveal the targets and history of adaptive evolution. By analysing high-coverage genome sequence data from 4 major colour pattern loci sampled from nearly 600 individuals in 53 populations, we show pervasive selection on wing patterns in the Heliconius adaptive radiation. The strongest signatures correspond to loci with the greatest phenotypic effects, consistent with visual selection by predators, and are found in colour patterns with geographically restricted distributions. These recent sweeps are similar between co-mimics and indicate colour pattern turn-over events despite strong stabilising selection. Using simulations, we compare sweep signatures expected under classic hard sweeps with those resulting from adaptive introgression, an important aspect of mimicry evolution in Heliconius butterflies. Simulated recipient populations show a distinct 'volcano' pattern with peaks of increased genetic diversity around the selected target, characteristic of sweeps of introgressed variation and consistent with diversity patterns found in some populations. Our genomic data reveal a surprisingly dynamic history of colour pattern selection and co-evolution in this adaptive radiation.
Collapse
Affiliation(s)
- Markus Moest
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Department of Ecology, University of Innsbruck, Innsbruck, Austria
| | - Steven M. Van Belleghem
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Department of Biology, University of Puerto Rico, Rio Piedras, Puerto Rico
| | - Jennifer E. James
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Department of Ecology and Evolutionary Biology, University of Arizona, Tucson, Arizona, United States of America
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences and Mathematics, Universidad del Rosario, Bogota D.C., Colombia
| | - Simon H. Martin
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, United Kingdom
| | - Sarah L. Barker
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Gilson R. P. Moreira
- Departamento de Zoologia, Universidade Federal do Rio Grande do Sul, Porto Alegre, Brazil
| | - Claire Mérot
- IBIS, Department of Biology, Université Laval, Québec, Canada
| | - Mathieu Joron
- Centre d'Ecologie Fonctionnelle et Evolutive, UMR 5175 CNRS—Université de Montpellier—Université Paul Valéry Montpellier—EPHE, Montpellier, France
| | - Nicola J. Nadeau
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, United Kingdom
| | | | - Chris D. Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
36
|
Byers KJRP, Darragh K, Musgrove J, Almeida DA, Garza SF, Warren IA, Rastas PM, Kučka M, Chan YF, Merrill RM, Schulz S, McMillan WO, Jiggins CD. A major locus controls a biologically active pheromone component in Heliconius melpomene. Evolution 2020; 74:349-364. [PMID: 31913497 PMCID: PMC7027519 DOI: 10.1111/evo.13922] [Citation(s) in RCA: 8] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/23/2019] [Revised: 12/20/2019] [Accepted: 01/03/2020] [Indexed: 12/23/2022]
Abstract
Understanding the production, response, and genetics of signals used in mate choice can inform our understanding of the evolution of both intraspecific mate choice and reproductive isolation. Sex pheromones are important for courtship and mate choice in many insects, but we know relatively little of their role in butterflies. The butterfly Heliconius melpomene uses a complex blend of wing androconial compounds during courtship. Electroantennography in H. melpomene and its close relative Heliconius cydno showed that responses to androconial extracts were not species specific. Females of both species responded equally strongly to extracts of both species, suggesting conservation of peripheral nervous system elements across the two species. Individual blend components provoked little to no response, with the exception of octadecanal, a major component of the H. melpomene blend. Supplementing octadecanal on the wings of octadecanal-rich H. melpomene males led to an increase in the time until mating, demonstrating the bioactivity of octadecanal in Heliconius. Using quantitative trait locus (QTL) mapping, we identified a single locus on chromosome 20 responsible for 41% of the parental species' difference in octadecanal production. This QTL does not overlap with any of the major wing color or mate choice loci, nor does it overlap with known regions of elevated or reduced FST . A set of 16 candidate fatty acid biosynthesis genes lies underneath the QTL. Pheromones in Heliconius carry information relevant for mate choice and are under simple genetic control, suggesting they could be important during speciation.
Collapse
Affiliation(s)
- Kelsey J. R. P. Byers
- Department of ZoologyUniversity of CambridgeCambridgeCB2 3EJUnited Kingdom
- Smithsonian Tropical Research InstitutePanamaPanama
| | - Kathy Darragh
- Department of ZoologyUniversity of CambridgeCambridgeCB2 3EJUnited Kingdom
- Smithsonian Tropical Research InstitutePanamaPanama
| | | | - Diana Abondano Almeida
- Smithsonian Tropical Research InstitutePanamaPanama
- Current address: Institute for Ecology, Evolution, and DiversityGoethe Universität60323FrankfurtGermany
| | - Sylvia Fernanda Garza
- Smithsonian Tropical Research InstitutePanamaPanama
- Current address: Department of Collective BehaviorMax Planck Institute of Animal Behavior78315KonstanzGermany
- Current address: Centre for the Advanced Study of Collective BehaviorUniversity of Konstanz78464KonstanzGermany
| | - Ian A. Warren
- Department of ZoologyUniversity of CambridgeCambridgeCB2 3EJUnited Kingdom
| | - Pasi M. Rastas
- Institute of BiotechnologyUniversity of Helsinki00014HelsinkiFinland
| | - Marek Kučka
- Friedrich Miescher LaboratoryMax Planck Society72076TübingenGermany
| | | | - Richard M. Merrill
- Division of Evolutionary BiologyLudwig‐Maximilians‐Universität München80539MunichGermany
| | - Stefan Schulz
- Department of Life SciencesInstitute of Organic Chemistry, Technische Universität Braunschweig38106BraunschweigGermany
| | | | - Chris D. Jiggins
- Department of ZoologyUniversity of CambridgeCambridgeCB2 3EJUnited Kingdom
- Smithsonian Tropical Research InstitutePanamaPanama
| |
Collapse
|
37
|
Li F, Zhao X, Li M, He K, Huang C, Zhou Y, Li Z, Walters JR. Insect genomes: progress and challenges. INSECT MOLECULAR BIOLOGY 2019; 28:739-758. [PMID: 31120160 DOI: 10.1111/imb.12599] [Citation(s) in RCA: 78] [Impact Index Per Article: 15.6] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/02/2017] [Revised: 03/22/2019] [Accepted: 05/14/2019] [Indexed: 05/24/2023]
Abstract
In the wake of constant improvements in sequencing technologies, numerous insect genomes have been sequenced. Currently, 1219 insect genome-sequencing projects have been registered with the National Center for Biotechnology Information, including 401 that have genome assemblies and 155 with an official gene set of annotated protein-coding genes. Comparative genomics analysis showed that the expansion or contraction of gene families was associated with well-studied physiological traits such as immune system, metabolic detoxification, parasitism and polyphagy in insects. Here, we summarize the progress of insect genome sequencing, with an emphasis on how this impacts research on pest control. We begin with a brief introduction to the basic concepts of genome assembly, annotation and metrics for evaluating the quality of draft assemblies. We then provide an overview of genome information for numerous insect species, highlighting examples from prominent model organisms, agricultural pests and disease vectors. We also introduce the major insect genome databases. The increasing availability of insect genomic resources is beneficial for developing alternative pest control methods. However, many opportunities remain for developing data-mining tools that make maximal use of the available insect genome resources. Although rapid progress has been achieved, many challenges remain in the field of insect genomics.
Collapse
Affiliation(s)
- F Li
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - X Zhao
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - M Li
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - K He
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - C Huang
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - Y Zhou
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - Z Li
- Ministry of Agriculture Key Lab of Molecular Biology of Crop Pathogens and Insects, Zhejiang University, Hangzhou, China
| | - J R Walters
- Ecology and Evolutionary Biology, University of Kansas, Lawrence, KS, USA
| |
Collapse
|
38
|
Catanach A, Crowhurst R, Deng C, David C, Bernatchez L, Wellenreuther M. The genomic pool of standing structural variation outnumbers single nucleotide polymorphism by threefold in the marine teleost Chrysophrys auratus. Mol Ecol 2019; 28:1210-1223. [PMID: 30770610 DOI: 10.1111/mec.15051] [Citation(s) in RCA: 35] [Impact Index Per Article: 7.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/07/2018] [Revised: 01/31/2019] [Accepted: 02/01/2019] [Indexed: 12/22/2022]
Abstract
Recent studies have highlighted an important role of structural variation (SV) in ecological and evolutionary processes, but few have studied nonmodel species in the wild. As part of our long-term research programme on the nonmodel teleost fish Australasian snapper (Chrysophrys auratus), we aim to build one of the first catalogues of genomic variants (SNPs and indels, and deletions, duplications and inversions) in fishes and evaluate overlap of genomic variants with regions under putative selection (Tajima's D and π), and coding sequences (genes). For this, we analysed six males and six females from three locations in New Zealand and generated a high-resolution genomic variation catalogue. We characterized 20,385 SVs and found they intersected with almost a third of all annotated genes. Together with small indels, SVs account for three times more variation in the genome in terms of bases affected compared to SNPs. We found that a sizeable portion of detected SVs was in the upper and lower genomic regions of Tajima's D and π, indicating that some of these have an effect on the phenotype. Together, these results shed light on the often neglected genomic variation that is produced by SVs and highlights the need to go beyond the mere measure of SNPs when investigating evolutionary processes, such as species diversification and adaptation.
Collapse
Affiliation(s)
- Andrew Catanach
- The New Zealand Institute for Plant & Food Research Ltd, Lincoln, New Zealand
| | - Ross Crowhurst
- The New Zealand Institute for Plant & Food Research Ltd, Auckland, New Zealand
| | - Cecilia Deng
- The New Zealand Institute for Plant & Food Research Ltd, Auckland, New Zealand
| | - Charles David
- The New Zealand Institute for Plant & Food Research Ltd, Lincoln, New Zealand
| | - Louis Bernatchez
- Institut de Biologie Intégrative et des Systèmes (IBIS), Université Laval, Québec City, Quebec, Canada
| | - Maren Wellenreuther
- The New Zealand Institute for Plant & Food Research Ltd, Nelson, New Zealand.,School of Biological Sciences, University of Auckland, Auckland, New Zealand
| |
Collapse
|
39
|
Edelman NB, Frandsen PB, Miyagi M, Clavijo B, Davey J, Dikow RB, García-Accinelli G, Van Belleghem SM, Patterson N, Neafsey DE, Challis R, Kumar S, Moreira GRP, Salazar C, Chouteau M, Counterman BA, Papa R, Blaxter M, Reed RD, Dasmahapatra KK, Kronforst M, Joron M, Jiggins CD, McMillan WO, Di Palma F, Blumberg AJ, Wakeley J, Jaffe D, Mallet J. Genomic architecture and introgression shape a butterfly radiation. Science 2019; 366:594-599. [PMID: 31672890 PMCID: PMC7197882 DOI: 10.1126/science.aaw2090] [Citation(s) in RCA: 247] [Impact Index Per Article: 49.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/11/2018] [Accepted: 09/16/2019] [Indexed: 12/26/2022]
Abstract
We used 20 de novo genome assemblies to probe the speciation history and architecture of gene flow in rapidly radiating Heliconius butterflies. Our tests to distinguish incomplete lineage sorting from introgression indicate that gene flow has obscured several ancient phylogenetic relationships in this group over large swathes of the genome. Introgressed loci are underrepresented in low-recombination and gene-rich regions, consistent with the purging of foreign alleles more tightly linked to incompatibility loci. Here, we identify a hitherto unknown inversion that traps a color pattern switch locus. We infer that this inversion was transferred between lineages by introgression and is convergent with a similar rearrangement in another part of the genus. These multiple de novo genome sequences enable improved understanding of the importance of introgression and selective processes in adaptive radiation.
Collapse
Affiliation(s)
- Nathaniel B Edelman
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA.
| | - Paul B Frandsen
- Department of Plant and Wildlife Sciences, Brigham Young University, Provo, UT 84602, USA
- Data Science Lab, Office of the Chief Information Officer, Smithsonian Institution, Washington, DC 20560, USA
| | - Miriam Miyagi
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | | | - John Davey
- Bioscience Technology Facility, Department of Biology, University of York, York YO10 5DD, UK
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, UK
| | - Rebecca B Dikow
- Data Science Lab, Office of the Chief Information Officer, Smithsonian Institution, Washington, DC 20560, USA
| | | | - Steven M Van Belleghem
- Department of Biology, University of Puerto Rico, Río Piedras Campus, San Juan, PR 00931-3360, Puerto Rico
| | - Nick Patterson
- Department of Human Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
- Broad Institute of MIT and Harvard, Cambridge, MA, 02142 USA
| | - Daniel E Neafsey
- Broad Institute of MIT and Harvard, Cambridge, MA, 02142 USA
- Harvard TH Chan School of Public Health, Boston, MA 02115, USA
| | - Richard Challis
- Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK
| | - Sujai Kumar
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh EH9 3JT, UK
| | - Gilson R P Moreira
- Departamento de Zoologia, Universidade Federal do Rio Grande do Sul, Porto Alegre, 91501-970 Brasil
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences and Mathematics, Universidad del Rosario, Carrera 24, No. 63C-69, Bogotá D.C. 111221, Colombia
| | - Mathieu Chouteau
- Laboratoire Ecologie, Evolution, Interactions des Systèmes Amazoniens (LEEISA), USR 3456, Université De Guyane, CNRS Guyane, 275 Route de Montabo, 97334 Cayenne, French Guiana
| | - Brian A Counterman
- Department of Biological Sciences, Mississippi State University, Starkville, MS 39762, USA
| | - Riccardo Papa
- Department of Biology, University of Puerto Rico, Río Piedras Campus, San Juan, PR 00931-3360, Puerto Rico
- Molecular Sciences and Research Center, University of Puerto Rico, San Juan, PR 00931-3360, Puerto Rico
| | - Mark Blaxter
- Wellcome Sanger Institute, Wellcome Genome Campus, Cambridge CB10 1SA, UK
| | - Robert D Reed
- Department of Ecology and Evolutionary Biology, Cornell University, Ithaca, NY 14853, USA
| | - Kanchon K Dasmahapatra
- Bioscience Technology Facility, Department of Biology, University of York, York YO10 5DD, UK
| | - Marcus Kronforst
- Department of Ecology and Evolution, University of Chicago, Chicago, IL 60637, USA
| | - Mathieu Joron
- CEFE, CNRS, Université de Montpellier, Université Paul Valéry Montpellier 3, EPHE, IRD, 34090 Montpellier, France
| | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge CB2 3EJ, UK
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Apartado 0843-03092 Panamá, Panama
| | | | - Andrew J Blumberg
- Department of Mathematics, University of Texas, Austin, TX 78712, USA
| | - John Wakeley
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA
| | - David Jaffe
- Broad Institute of MIT and Harvard, Cambridge, MA, 02142 USA
- 10x Genomics, Pleasanton, CA 94566, USA
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA 02138, USA.
| |
Collapse
|
40
|
Mackintosh A, Laetsch DR, Hayward A, Charlesworth B, Waterfall M, Vila R, Lohse K. The determinants of genetic diversity in butterflies. Nat Commun 2019; 10:3466. [PMID: 31371715 PMCID: PMC6672018 DOI: 10.1038/s41467-019-11308-4] [Citation(s) in RCA: 48] [Impact Index Per Article: 9.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/06/2019] [Accepted: 07/01/2019] [Indexed: 11/09/2022] Open
Abstract
Under the neutral theory, genetic diversity is expected to increase with population size. While comparative analyses have consistently failed to find strong relationships between census population size and genetic diversity, a recent study across animals identified a strong correlation between propagule size and genetic diversity, suggesting that r-strategists that produce many small offspring, have greater long-term population sizes. Here we compare genome-wide genetic diversity across 38 species of European butterflies (Papilionoidea), a group that shows little variation in reproductive strategy. We show that genetic diversity across butterflies varies over an order of magnitude and that this variation cannot be explained by differences in current abundance, propagule size, host or geographic range. Instead, neutral genetic diversity is negatively correlated with body size and positively with the length of the genetic map. This suggests that genetic diversity is determined both by differences in long-term population size and the effect of selection on linked sites.
Collapse
Affiliation(s)
- Alexander Mackintosh
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3FL, UK.
| | - Dominik R Laetsch
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Alexander Hayward
- Centre for Ecology and Conservation, University of Exeter, Penryn Campus, Cornwall, TR10 9FE, UK
| | - Brian Charlesworth
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Martin Waterfall
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3FL, UK
| | - Roger Vila
- Institut de Biologia Evolutiva (CSIC Universitat Pompeu Fabra), Passeig Marítim de la Barceloneta 37, ESP-08003, Barcelona, Spain
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, Edinburgh, EH9 3FL, UK.
| |
Collapse
|
41
|
Barth JMI, Villegas-Ríos D, Freitas C, Moland E, Star B, André C, Knutsen H, Bradbury I, Dierking J, Petereit C, Righton D, Metcalfe J, Jakobsen KS, Olsen EM, Jentoft S. Disentangling structural genomic and behavioural barriers in a sea of connectivity. Mol Ecol 2019; 28:1394-1411. [PMID: 30633410 PMCID: PMC6518941 DOI: 10.1111/mec.15010] [Citation(s) in RCA: 36] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/01/2018] [Revised: 12/21/2018] [Accepted: 12/27/2018] [Indexed: 12/17/2022]
Abstract
Genetic divergence among populations arises through natural selection or drift and is counteracted by connectivity and gene flow. In sympatric populations, isolating mechanisms are thus needed to limit the homogenizing effects of gene flow to allow for adaptation and speciation. Chromosomal inversions act as an important mechanism maintaining isolating barriers, yet their role in sympatric populations and divergence with gene flow is not entirely understood. Here, we revisit the question of whether inversions play a role in the divergence of connected populations of the marine fish Atlantic cod (Gadus morhua), by exploring a unique data set combining whole‐genome sequencing data and behavioural data obtained with acoustic telemetry. Within a confined fjord environment, we find three genetically differentiated Atlantic cod types belonging to the oceanic North Sea population, the western Baltic population and a local fjord‐type cod. Continuous behavioural tracking over 4 year revealed temporally stable sympatry of these types within the fjord. Despite overall weak genetic differentiation consistent with high levels of gene flow, we detected significant frequency shifts of three previously identified inversions, indicating an adaptive barrier to gene flow. In addition, behavioural data indicated that North Sea cod and individuals homozygous for the LG12 inversion had lower fitness in the fjord environment. However, North Sea and fjord‐type cod also occupy different depths, possibly contributing to prezygotic reproductive isolation and representing a behavioural barrier to gene flow. Our results provide the first insights into a complex interplay of genomic and behavioural isolating barriers in Atlantic cod and establish a new model system towards an understanding of the role of genomic structural variants in adaptation and diversification.
Collapse
Affiliation(s)
- Julia M I Barth
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway.,Zoological Institute, University of Basel, Basel, Switzerland
| | - David Villegas-Ríos
- Department of Ecology and Marine Resources, Mediterranean Institute for Advanced Studies, IMEDEA CSIC-UIB, Esporles, Spain.,Department of Ecology and Marine Resources, Institute of Marine Research, (IIM CSIC), Vigo, Spain
| | - Carla Freitas
- Institute for Marine Research, Flødevigen, Norway.,Centre for Coastal Research, University of Agder, Agder, Norway.,Oceanic Observatory of Madeira, Funchal, Portugal
| | - Even Moland
- Institute for Marine Research, Flødevigen, Norway.,Centre for Coastal Research, University of Agder, Agder, Norway
| | - Bastiaan Star
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Carl André
- Department of Marine Sciences - Tjärnö, University of Gothenburg, Gothenburg, Sweden
| | - Halvor Knutsen
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway.,Institute for Marine Research, Flødevigen, Norway.,Centre for Coastal Research, University of Agder, Agder, Norway
| | - Ian Bradbury
- Science Branch, Fisheries and Oceans Canada, St. John's, Newfoundland and Labrador, Canada
| | - Jan Dierking
- GEOMAR Helmholtz Centre for Ocean Research, Kiel, Germany
| | | | - David Righton
- Centre for Environment, Fisheries and Aquaculture Science (CEFAS), Lowestoft, UK
| | - Julian Metcalfe
- Centre for Environment, Fisheries and Aquaculture Science (CEFAS), Lowestoft, UK
| | - Kjetill S Jakobsen
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| | - Esben M Olsen
- Institute for Marine Research, Flødevigen, Norway.,Centre for Coastal Research, University of Agder, Agder, Norway
| | - Sissel Jentoft
- Centre for Ecological and Evolutionary Synthesis (CEES), Department of Biosciences, University of Oslo, Oslo, Norway
| |
Collapse
|
42
|
Faria R, Chaube P, Morales HE, Larsson T, Lemmon AR, Lemmon EM, Rafajlović M, Panova M, Ravinet M, Johannesson K, Westram AM, Butlin RK. Multiple chromosomal rearrangements in a hybrid zone between Littorina saxatilis ecotypes. Mol Ecol 2019; 28:1375-1393. [PMID: 30537056 PMCID: PMC6518922 DOI: 10.1111/mec.14972] [Citation(s) in RCA: 69] [Impact Index Per Article: 13.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/29/2018] [Revised: 10/07/2018] [Accepted: 10/08/2018] [Indexed: 12/22/2022]
Abstract
Both classical and recent studies suggest that chromosomal inversion polymorphisms are important in adaptation and speciation. However, biases in discovery and reporting of inversions make it difficult to assess their prevalence and biological importance. Here, we use an approach based on linkage disequilibrium among markers genotyped for samples collected across a transect between contrasting habitats to detect chromosomal rearrangements de novo. We report 17 polymorphic rearrangements in a single locality for the coastal marine snail, Littorina saxatilis. Patterns of diversity in the field and of recombination in controlled crosses provide strong evidence that at least the majority of these rearrangements are inversions. Most show clinal changes in frequency between habitats, suggestive of divergent selection, but only one appears to be fixed for different arrangements in the two habitats. Consistent with widespread evidence for balancing selection on inversion polymorphisms, we argue that a combination of heterosis and divergent selection can explain the observed patterns and should be considered in other systems spanning environmental gradients.
Collapse
Affiliation(s)
- Rui Faria
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | - Pragya Chaube
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK
| | - Hernán E Morales
- Department of Marine Sciences, Centre for Marine Evolutionary Biology, University of Gothenburg, Gothenburg, Sweden
| | - Tomas Larsson
- Department of Marine Sciences, Centre for Marine Evolutionary Biology, University of Gothenburg, Gothenburg, Sweden
| | - Alan R Lemmon
- Department of Scientific Computing, Florida State University, Tallahassee, Florida
| | - Emily M Lemmon
- Department of Biological Science, Florida State University, Tallahassee, Florida
| | - Marina Rafajlović
- Department of Marine Sciences, Centre for Marine Evolutionary Biology, University of Gothenburg, Gothenburg, Sweden
| | - Marina Panova
- Department of Marine Sciences at Tjärnö, Centre for Marine Evolutionary Biology, University of Gothenburg, Strömstad, Sweden
| | - Mark Ravinet
- Centre for Ecological and Evolutionary Synthesis, University of Oslo, Oslo, Norway
| | - Kerstin Johannesson
- Department of Marine Sciences at Tjärnö, Centre for Marine Evolutionary Biology, University of Gothenburg, Strömstad, Sweden
| | - Anja M Westram
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK.,IST Austria, Klosterneuburg, Austria
| | - Roger K Butlin
- Department of Animal and Plant Sciences, University of Sheffield, Sheffield, UK.,Department of Marine Sciences at Tjärnö, Centre for Marine Evolutionary Biology, University of Gothenburg, Strömstad, Sweden
| |
Collapse
|
43
|
Merrill RM, Rastas P, Martin SH, Melo MC, Barker S, Davey J, McMillan WO, Jiggins CD. Genetic dissection of assortative mating behavior. PLoS Biol 2019; 17:e2005902. [PMID: 30730873 PMCID: PMC6366751 DOI: 10.1371/journal.pbio.2005902] [Citation(s) in RCA: 52] [Impact Index Per Article: 10.4] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/03/2018] [Accepted: 11/06/2018] [Indexed: 12/26/2022] Open
Abstract
The evolution of new species is made easier when traits under divergent ecological selection are also mating cues. Such ecological mating cues are now considered more common than previously thought, but we still know little about the genetic changes underlying their evolution or more generally about the genetic basis for assortative mating behaviors. Both tight physical linkage and the existence of large-effect preference loci will strengthen genetic associations between behavioral and ecological barriers, promoting the evolution of assortative mating. The warning patterns of Heliconius melpomene and H. cydno are under disruptive selection due to increased predation of nonmimetic hybrids and are used during mate recognition. We carried out a genome-wide quantitative trait locus (QTL) analysis of preference behaviors between these species and showed that divergent male preference has a simple genetic basis. We identify three QTLs that together explain a large proportion (approximately 60%) of the difference in preference behavior observed between the parental species. One of these QTLs is just 1.2 (0-4.8) centiMorgans (cM) from the major color pattern gene optix, and, individually, all three have a large effect on the preference phenotype. Genomic divergence between H. cydno and H. melpomene is high but broadly heterogenous, and admixture is reduced at the preference-optix color pattern locus but not the other preference QTLs. The simple genetic architecture we reveal will facilitate the evolution and maintenance of new species despite ongoing gene flow by coupling behavioral and ecological aspects of reproductive isolation.
Collapse
Affiliation(s)
- Richard M. Merrill
- Division of Evolutionary Biology, Ludwig-Maximilians-Universität, München, Germany
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Smithsonian Tropical Research Institute, Panama City, Panama
| | - Pasi Rastas
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Simon H. Martin
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - Maria C. Melo
- Smithsonian Tropical Research Institute, Panama City, Panama
- IST Austria, Klosterburg, Austria
| | - Sarah Barker
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - John Davey
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
- Department of Biology, University of York, York, United Kingdom
| | | | - Chris D. Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
44
|
Martin SH, Davey JW, Salazar C, Jiggins CD. Recombination rate variation shapes barriers to introgression across butterfly genomes. PLoS Biol 2019; 17:e2006288. [PMID: 30730876 PMCID: PMC6366726 DOI: 10.1371/journal.pbio.2006288] [Citation(s) in RCA: 163] [Impact Index Per Article: 32.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/08/2018] [Accepted: 01/07/2019] [Indexed: 12/30/2022] Open
Abstract
Hybridisation and introgression can dramatically alter the relationships among groups of species, leading to phylogenetic discordance across the genome and between populations. Introgression can also erode species differences over time, but selection against introgression at certain loci acts to maintain postmating species barriers. Theory predicts that species barriers made up of many loci throughout the genome should lead to a broad correlation between introgression and recombination rate, which determines the extent to which selection on deleterious foreign alleles will affect neutral alleles at physically linked loci. Here, we describe the variation in genealogical relationships across the genome among three species of Heliconius butterflies: H. melpomene (mel), H. cydno (cyd), and H. timareta (tim), using whole genomes of 92 individuals, and ask whether this variation can be explained by heterogeneous barriers to introgression. We find that species relationships vary predictably at the chromosomal scale. By quantifying recombination rate and admixture proportions, we then show that rates of introgression are predicted by variation in recombination rate. This implies that species barriers are highly polygenic, with selection acting against introgressed alleles across most of the genome. In addition, long chromosomes, which have lower recombination rates, produce stronger barriers on average than short chromosomes. Finally, we find a consistent difference between two species pairs on either side of the Andes, which suggests differences in the architecture of the species barriers. Our findings illustrate how the combined effects of hybridisation, recombination, and natural selection, acting at multitudes of loci over long periods, can dramatically sculpt the phylogenetic relationships among species.
Collapse
Affiliation(s)
- Simon H. Martin
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| | - John W. Davey
- Department of Biology, University of York, York, United Kingdom
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences and Mathematics, Universidad del Rosario, Bogota, Colombia
| | - Chris D. Jiggins
- Department of Zoology, University of Cambridge, Cambridge, United Kingdom
| |
Collapse
|
45
|
Doellman MM, Feder JL. Genomic transitions during host race and species formation. CURRENT OPINION IN INSECT SCIENCE 2019; 31:84-92. [PMID: 31109679 DOI: 10.1016/j.cois.2018.11.006] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 11/03/2018] [Revised: 11/19/2018] [Accepted: 11/20/2018] [Indexed: 06/09/2023]
Abstract
Darwin recognized species as discontinuous, yet considered them to be formed by an incremental process of natural selection. Recent theoretical work on 'genome-wide congealing' is bridging this gap between the gradualism of divergent selection and rapid genome-wide divergence, particularly during ecological speciation-with-gene-flow. Host races and species of phytophagous insects, displaying a spectrum of divergence and gene flow among member taxa, provide model systems for testing predicted non-linear transitions from 'genic' divergence at a few uncoupled loci to 'genomic' divergence with genome-wide coupling of selected loci and strong reproductive isolation. Integrating across natural history, genomics, and evolutionary theory, emerging research suggests a tipping point from 'genic' to 'genomic' divergence between host races and species, during both sympatric speciation and secondary contact.
Collapse
Affiliation(s)
- Meredith M Doellman
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN 46556, USA.
| | - Jeffrey L Feder
- Department of Biological Sciences, University of Notre Dame, Notre Dame, IN 46556, USA
| |
Collapse
|
46
|
Kolora SRR, Weigert A, Saffari A, Kehr S, Walter Costa MB, Spröer C, Indrischek H, Chintalapati M, Lohse K, Doose G, Overmann J, Bunk B, Bleidorn C, Grimm-Seyfarth A, Henle K, Nowick K, Faria R, Stadler PF, Schlegel M. Divergent evolution in the genomes of closely related lacertids, Lacerta viridis and L. bilineata, and implications for speciation. Gigascience 2019; 8:giy160. [PMID: 30535196 PMCID: PMC6381762 DOI: 10.1093/gigascience/giy160] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 09/19/2018] [Accepted: 11/29/2018] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND Lacerta viridis and Lacerta bilineata are sister species of European green lizards (eastern and western clades, respectively) that, until recently, were grouped together as the L. viridis complex. Genetic incompatibilities were observed between lacertid populations through crossing experiments, which led to the delineation of two separate species within the L. viridis complex. The population history of these sister species and processes driving divergence are unknown. We constructed the first high-quality de novo genome assemblies for both L. viridis and L. bilineata through Illumina and PacBio sequencing, with annotation support provided from transcriptome sequencing of several tissues. To estimate gene flow between the two species and identify factors involved in reproductive isolation, we studied their evolutionary history, identified genomic rearrangements, detected signatures of selection on non-coding RNA, and on protein-coding genes. FINDINGS Here we show that gene flow was primarily unidirectional from L. bilineata to L. viridis after their split at least 1.15 million years ago. We detected positive selection of the non-coding repertoire; mutations in transcription factors; accumulation of divergence through inversions; selection on genes involved in neural development, reproduction, and behavior, as well as in ultraviolet-response, possibly driven by sexual selection, whose contribution to reproductive isolation between these lacertid species needs to be further evaluated. CONCLUSION The combination of short and long sequence reads resulted in one of the most complete lizard genome assemblies. The characterization of a diverse array of genomic features provided valuable insights into the demographic history of divergence among European green lizards, as well as key species differences, some of which are candidates that could have played a role in speciation. In addition, our study generated valuable genomic resources that can be used to address conservation-related issues in lacertids.
Collapse
Affiliation(s)
- Sree Rohit Raj Kolora
- German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Deutscher Platz 5e, Leipzig, 04103, Germany
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstrasse 16-18, Leipzig, 04107, Germany
- Molecular Evolution and Systematics of Animals, Institute of Biology, University of Leipzig, Talstrasse 33, Leipzig, 04103, Germany
| | - Anne Weigert
- Molecular Evolution and Systematics of Animals, Institute of Biology, University of Leipzig, Talstrasse 33, Leipzig, 04103, Germany
- Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, Leipzig, 04103, Germany
| | - Amin Saffari
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstrasse 16-18, Leipzig, 04107, Germany
- Human Biology Group, Institute for Zoology, Department of Biology, Chemistry and Pharmacy, Freie Universität Berlin, Königin-Luise-Straße 1–3, Berlin, D-14195, Germany
| | - Stephanie Kehr
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstrasse 16-18, Leipzig, 04107, Germany
| | - Maria Beatriz Walter Costa
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstrasse 16-18, Leipzig, 04107, Germany
- Embrapa Agroenergia, Parque Estacaeo Biologica (PqEB), Asa Norte, Brasilia/DF, 70770-901, Brazil
| | - Cathrin Spröer
- Department of Microbial Ecology and Diversity Research, Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Inhoffenstrasse 7B, Braunschweig, 38124, Germany
| | - Henrike Indrischek
- Max Planck Institute of Molecular Cell Biology and Genetics, Pfotenhauerstrasse 108, Dresden, 01307, Germany
- Max Planck Institute for Physics of Complex Systems, Noethnitzerstrasse 38, 01187 Dresden, Germany
- Center for Systems Biology Dresden, Pfotenhauerstrasse 108, 01397 Dresden, Germany
| | - Manjusha Chintalapati
- Max Planck Institute for Evolutionary Anthropology, Deutscher Platz 6, Leipzig, 04103, Germany
| | - Konrad Lohse
- Institute of Evolutionary Biology, University of Edinburgh, King's Buildings, Charlotte Auerbach Road, Edinburgh, EH9 3FL, United Kingdom
| | - Gero Doose
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstrasse 16-18, Leipzig, 04107, Germany
| | - Jörg Overmann
- German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Deutscher Platz 5e, Leipzig, 04103, Germany
- Department of Microbial Ecology and Diversity Research, Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Inhoffenstrasse 7B, Braunschweig, 38124, Germany
| | - Boyke Bunk
- Department of Microbial Ecology and Diversity Research, Leibniz Institute DSMZ-German Collection of Microorganisms and Cell Cultures, Inhoffenstrasse 7B, Braunschweig, 38124, Germany
| | - Christoph Bleidorn
- German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Deutscher Platz 5e, Leipzig, 04103, Germany
- Department of Animal Evolution and Biodiversity, University of Göttingen, Untere Karspüle 2, Göttingen, 37073, Germany
- Museo Nacional de Ciencias Naturales, Spanish National Research Council (CSIC), Madrid, 28006, Spain
| | - Annegret Grimm-Seyfarth
- Department of Conservation Biology, UFZ - Helmholtz Center for Environmental Research, Permoserstrasse 15, Leipzig, 04318, Germany
- Plant Ecology and Nature Conservation, University of Potsdam, Am Mühlenberg 3, Potsdam, 14476, Germany
| | - Klaus Henle
- Department of Conservation Biology, UFZ - Helmholtz Center for Environmental Research, Permoserstrasse 15, Leipzig, 04318, Germany
| | - Katja Nowick
- Human Biology Group, Institute for Zoology, Department of Biology, Chemistry and Pharmacy, Freie Universität Berlin, Königin-Luise-Straße 1–3, Berlin, D-14195, Germany
| | - Rui Faria
- Department of Animal and Plant Sciences, Alfred Building, University of Sheffield, Western Bank, Sheffield, S10 2TN, United Kingdom
| | - Peter F Stadler
- German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Deutscher Platz 5e, Leipzig, 04103, Germany
- Bioinformatics Group, Department of Computer Science, and Interdisciplinary Center for Bioinformatics, Universität Leipzig, Härtelstrasse 16-18, Leipzig, 04107, Germany
- Competence Center for Scalable Data Services and Solutions Dresden/Leipzig, Universität Leipzig, Augustusplatz 12, Leipzig, 04107, Germany
- Max-Planck-Institute for Mathematics in the Sciences, Inselstrasse 22, Leipzig, 04103, Germany
- Fraunhofer Institut Für Zelltherapie Und Immunologie, Perlickstrasse 1, Leipzig, 04103, Germany
- Department of Theoretical Chemistry, University of Vienna, Währinger strasse 17, Wien, 1090, Austria
- Center for non-Coding RNA in Technology and Health, University of Copenhagen, Gronnegardsvej 3, Frederiksberg C, 1870, Denmark
- Santa Fe Institute, 1399 Hyde Park Road, Santa Fe, New Mexico, 87501, USA
| | - Martin Schlegel
- German Centre for Integrative Biodiversity Research (iDiv) Halle-Jena-Leipzig, Deutscher Platz 5e, Leipzig, 04103, Germany
- Molecular Evolution and Systematics of Animals, Institute of Biology, University of Leipzig, Talstrasse 33, Leipzig, 04103, Germany
| |
Collapse
|
47
|
Cocker JM, Wright J, Li J, Swarbreck D, Dyer S, Caccamo M, Gilmartin PM. Primula vulgaris (primrose) genome assembly, annotation and gene expression, with comparative genomics on the heterostyly supergene. Sci Rep 2018; 8:17942. [PMID: 30560928 PMCID: PMC6299000 DOI: 10.1038/s41598-018-36304-4] [Citation(s) in RCA: 21] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/15/2018] [Accepted: 11/14/2018] [Indexed: 11/24/2022] Open
Abstract
Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species.
Collapse
Affiliation(s)
- Jonathan M Cocker
- School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ, United Kingdom.,Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, United Kingdom
| | - Jonathan Wright
- Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, United Kingdom
| | - Jinhong Li
- School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ, United Kingdom.,Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, United Kingdom
| | - David Swarbreck
- Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, United Kingdom
| | - Sarah Dyer
- National Institute for Agricultural Botany, Huntingdon Road, Cambridge, CB3 0LE, United Kingdom
| | - Mario Caccamo
- National Institute for Agricultural Botany, Huntingdon Road, Cambridge, CB3 0LE, United Kingdom
| | - Philip M Gilmartin
- School of Biological Sciences, University of East Anglia, Norwich Research Park, Norwich, NR4 7TJ, United Kingdom. .,Earlham Institute, Norwich Research Park, Norwich, NR4 7UZ, United Kingdom.
| |
Collapse
|
48
|
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
|
49
|
Fuller ZL, Leonard CJ, Young RE, Schaeffer SW, Phadnis N. Ancestral polymorphisms explain the role of chromosomal inversions in speciation. PLoS Genet 2018; 14:e1007526. [PMID: 30059505 PMCID: PMC6085072 DOI: 10.1371/journal.pgen.1007526] [Citation(s) in RCA: 48] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/17/2018] [Revised: 08/09/2018] [Accepted: 06/29/2018] [Indexed: 01/28/2023] Open
Abstract
Understanding the role of chromosomal inversions in speciation is a fundamental problem in evolutionary genetics. Here, we perform a comprehensive reconstruction of the evolutionary histories of the chromosomal inversions in Drosophila persimilis and D. pseudoobscura. We provide a solution to the puzzling origins of the selfish Sex-Ratio arrangement in D. persimilis and uncover surprising patterns of phylogenetic discordance on this chromosome. These patterns show that, contrary to widely held views, all fixed chromosomal inversions between D. persimilis and D. pseudoobscura were already present in their ancestral population long before the species split. Our results suggest that patterns of higher genomic divergence and an association of reproductive isolation genes with chromosomal inversions may be a direct consequence of incomplete lineage sorting of ancestral polymorphisms. These findings force a reconsideration of the role of chromosomal inversions in speciation, not as protectors of existing hybrid incompatibilities, but as fertile grounds for their formation.
Collapse
Affiliation(s)
- Zachary L. Fuller
- Department of Biology, Erwin W. Mueller Laboratories, The Pennsylvania State University, University Park, PA, United States of America
| | | | - Randee E. Young
- Department of Biology, University of Utah, Salt Lake City, UT, United States of America
| | - Stephen W. Schaeffer
- Department of Biology, Erwin W. Mueller Laboratories, The Pennsylvania State University, University Park, PA, United States of America
| | - Nitin Phadnis
- Department of Biology, University of Utah, Salt Lake City, UT, United States of America
| |
Collapse
|
50
|
Van Belleghem SM, Baquero M, Papa R, Salazar C, McMillan WO, Counterman BA, Jiggins CD, Martin SH. Patterns of Z chromosome divergence among Heliconius species highlight the importance of historical demography. Mol Ecol 2018; 27:3852-3872. [PMID: 29569384 PMCID: PMC6151167 DOI: 10.1111/mec.14560] [Citation(s) in RCA: 43] [Impact Index Per Article: 7.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/21/2017] [Revised: 02/14/2018] [Accepted: 02/20/2018] [Indexed: 12/31/2022]
Abstract
Sex chromosomes are disproportionately involved in reproductive isolation and adaptation. In support of such a “large‐X” effect, genome scans between recently diverged populations and species pairs often identify distinct patterns of divergence on the sex chromosome compared to autosomes. When measures of divergence between populations are higher on the sex chromosome compared to autosomes, such patterns could be interpreted as evidence for faster divergence on the sex chromosome, that is “faster‐X”, barriers to gene flow on the sex chromosome. However, demographic changes can strongly skew divergence estimates and are not always taken into consideration. We used 224 whole‐genome sequences representing 36 populations from two Heliconius butterfly clades (H. erato and H. melpomene) to explore patterns of Z chromosome divergence. We show that increased divergence compared to equilibrium expectations can in many cases be explained by demographic change. Among Heliconius erato populations, for instance, population size increase in the ancestral population can explain increased absolute divergence measures on the Z chromosome compared to the autosomes, as a result of increased ancestral Z chromosome genetic diversity. Nonetheless, we do identify increased divergence on the Z chromosome relative to the autosomes in parapatric or sympatric species comparisons that imply postzygotic reproductive barriers. Using simulations, we show that this is consistent with reduced gene flow on the Z chromosome, perhaps due to greater accumulation of incompatibilities. Our work demonstrates the importance of taking demography into account to interpret patterns of divergence on the Z chromosome, but nonetheless provides evidence to support the Z chromosome as a strong barrier to gene flow in incipient Heliconius butterfly species.
Collapse
Affiliation(s)
- Steven M Van Belleghem
- Department of Zoology, University of Cambridge, Cambridge, UK.,Department of Biological Sciences, Mississippi State University, Mississippi State, MS, USA.,Department of Biology, Center for Applied Tropical Ecology and Conservation, University of Puerto Rico, Rio Piedras, Puerto Rico.,Smithsonian Tropical Research Institute, Apartado, Panamá, Panama
| | - Margarita Baquero
- Department of Biological Sciences, Mississippi State University, Mississippi State, MS, USA
| | - Riccardo Papa
- Department of Biology, Center for Applied Tropical Ecology and Conservation, University of Puerto Rico, Rio Piedras, Puerto Rico
| | - Camilo Salazar
- Biology Program, Faculty of Natural Sciences and Mathematics, Universidad del Rosario, Carrera, Bogota, Colombia
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Apartado, Panamá, Panama
| | - Brian A Counterman
- Department of Biological Sciences, Mississippi State University, Mississippi State, MS, USA
| | - Chris D Jiggins
- Department of Zoology, University of Cambridge, Cambridge, UK
| | - Simon H Martin
- Department of Zoology, University of Cambridge, Cambridge, UK
| |
Collapse
|