1
|
Bock DG, Baeckens S, Kolbe JJ, Losos JB. When adaptation is slowed down: Genomic analysis of evolutionary stasis in thermal tolerance during biological invasion in a novel climate. Mol Ecol 2024; 33:e17075. [PMID: 37489260 DOI: 10.1111/mec.17075] [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: 01/30/2023] [Revised: 06/25/2023] [Accepted: 07/04/2023] [Indexed: 07/26/2023]
Abstract
Research conducted during the past two decades has demonstrated that biological invasions are excellent models of rapid evolution. Even so, characteristics of invasive populations such as a short time for recombination to assemble optimal combinations of alleles may occasionally limit adaptation to new environments. Here, we investigated such genetic constraints to adaptation in the invasive brown anole (Anolis sagrei)-a tropical ectotherm that was introduced to the southeastern United States, a region with a much colder climate than in its native Caribbean range. We examined thermal physiology for 30 invasive populations and tested for a climatic cline in cold tolerance. Also, we used genomics to identify mechanisms that may limit adaptation. We found no support for a climatic cline, indicating that thermal tolerance did not shift adaptively. Concomitantly, population genomic results were consistent with the occurrence of recombination cold spots that comprise more than half of the genome and maintain long-range associations among alleles in invasive populations. These genomic regions overlap with both candidate thermal tolerance loci that we identified using a standard genome-wide association test. Moreover, we found that recombination cold spots do not have a large contribution to population differentiation in the invasive range, contrary to observations in the native range. We suggest that limited recombination is constraining the contribution of large swaths of the genome to adaptation in invasive brown anoles. Our study provides an example of evolutionary stasis during invasion and highlights the possibility that reduced recombination occasionally slows down adaptation in invasive populations.
Collapse
Affiliation(s)
- Dan G Bock
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
- Department of Biology, Washington University, St. Louis, Missouri, USA
| | - Simon Baeckens
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
- Evolution and Optics of Nanostructures Lab, Department of Biology, Ghent University, Ghent, Belgium
- Functional Morphology Lab, Department of Biology, University of Antwerp, Antwerp, Belgium
| | - Jason J Kolbe
- Department of Biological Sciences, University of Rhode Island, Kingston, Rhode Island, USA
| | - Jonathan B Losos
- Department of Biology, Washington University, St. Louis, Missouri, USA
| |
Collapse
|
2
|
Mirchandani CD, Shultz AJ, Thomas GWC, Smith SJ, Baylis M, Arnold B, Corbett-Detig R, Enbody E, Sackton TB. A Fast, Reproducible, High-throughput Variant Calling Workflow for Population Genomics. Mol Biol Evol 2024; 41:msad270. [PMID: 38069903 PMCID: PMC10764099 DOI: 10.1093/molbev/msad270] [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: 06/22/2023] [Revised: 10/27/2023] [Accepted: 11/22/2023] [Indexed: 01/05/2024] Open
Abstract
The increasing availability of genomic resequencing data sets and high-quality reference genomes across the tree of life present exciting opportunities for comparative population genomic studies. However, substantial challenges prevent the simple reuse of data across different studies and species, arising from variability in variant calling pipelines, data quality, and the need for computationally intensive reanalysis. Here, we present snpArcher, a flexible and highly efficient workflow designed for the analysis of genomic resequencing data in nonmodel organisms. snpArcher provides a standardized variant calling pipeline and includes modules for variant quality control, data visualization, variant filtering, and other downstream analyses. Implemented in Snakemake, snpArcher is user-friendly, reproducible, and designed to be compatible with high-performance computing clusters and cloud environments. To demonstrate the flexibility of this pipeline, we applied snpArcher to 26 public resequencing data sets from nonmammalian vertebrates. These variant data sets are hosted publicly to enable future comparative population genomic analyses. With its extensibility and the availability of public data sets, snpArcher will contribute to a broader understanding of genetic variation across species by facilitating the rapid use and reuse of large genomic data sets.
Collapse
Affiliation(s)
- Cade D Mirchandani
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Allison J Shultz
- Ornithology Department, Natural History Museum of Los Angeles County, Los Angeles, CA 90007, USA
| | | | - Sara J Smith
- Informatics Group, Harvard University, Cambridge, MA, USA
- Biology, Mount Royal University, Calgary, AB T3E 6K6, Canada
| | - Mara Baylis
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Brian Arnold
- Department of Ecology and Evolutionary Biology, Princeton University, Princeton, NJ, USA
- Center for Statistics and Machine Learning, Princeton University, Princeton, NJ, USA
| | - Russ Corbett-Detig
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
- Genomics Institute, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | - Erik Enbody
- Department of Biomolecular Engineering, University of California Santa Cruz, Santa Cruz, CA 95064, USA
| | | |
Collapse
|
3
|
Farleigh K, Ascanio A, Farleigh ME, Schield DR, Card DC, Leal M, Castoe TA, Jezkova T, Rodríguez-Robles JA. Signals of differential introgression in the genome of natural hybrids of Caribbean anoles. Mol Ecol 2023; 32:6000-6017. [PMID: 37861454 DOI: 10.1111/mec.17170] [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: 12/28/2021] [Revised: 08/30/2023] [Accepted: 10/03/2023] [Indexed: 10/21/2023]
Abstract
Hybridization facilitates recombination between divergent genetic lineages and can be shaped by both neutral and selective processes. Upon hybridization, loci with no net fitness effects introgress randomly from parental species into the genomes of hybrid individuals. Conversely, alleles from one parental species at some loci may provide a selective advantage to hybrids, resulting in patterns of introgression that do not conform to random expectations. We investigated genomic patterns of differential introgression in natural hybrids of two species of Caribbean anoles, Anolis pulchellus and A. krugi in Puerto Rico. Hybrids exhibit A. pulchellus phenotypes but possess A. krugi mitochondrial DNA, originated from multiple, independent hybridization events, and appear to have replaced pure A. pulchellus across a large area in western Puerto Rico. Combining genome-wide SNP datasets with bioinformatic methods to identify signals of differential introgression in hybrids, we demonstrate that the genomes of hybrids are dominated by pulchellus-derived alleles and show only 10%-20% A. krugi ancestry. The majority of A. krugi loci in hybrids exhibit a signal of non-random differential introgression and include loci linked to genes involved in development and immune function. Three of these genes (delta like canonical notch ligand 1, jagged1 and notch receptor 1) affect cell differentiation and growth and interact with mitochondrial function. Our results suggest that differential non-random introgression for a subset of loci may be driven by selection favouring the inheritance of compatible mitochondrial and nuclear-encoded genes in hybrids.
Collapse
Affiliation(s)
- Keaka Farleigh
- Department of Biology, Miami University, Oxford, Ohio, USA
| | | | | | - Drew R Schield
- Department of Ecology and Evolutionary Biology, University of Colorado, Boulder, Colorado, USA
| | - Daren C Card
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, Massachusetts, USA
- Museum of Comparative Zoology, Harvard University, Cambridge, Massachusetts, USA
| | - Manuel Leal
- Division of Biological Sciences, University of Missouri, Columbia, Missouri, USA
| | - Todd A Castoe
- Department of Biology, University of Texas, Arlington, Arlington, Texas, USA
| | - Tereza Jezkova
- Department of Biology, Miami University, Oxford, Ohio, USA
| | | |
Collapse
|
4
|
Card DC, Van Camp AG, Santonastaso T, Jensen-Seaman MI, Anthony NM, Edwards SV. Structure and evolution of the squamate major histocompatibility complex as revealed by two Anolis lizard genomes. Front Genet 2022; 13:979746. [PMID: 36425073 PMCID: PMC9679377 DOI: 10.3389/fgene.2022.979746] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Accepted: 10/20/2022] [Indexed: 11/10/2022] Open
Abstract
The major histocompatibility complex (MHC) is an important genomic region for adaptive immunity and has long been studied in ecological and evolutionary contexts, such as disease resistance and mate and kin selection. The MHC has been investigated extensively in mammals and birds but far less so in squamate reptiles, the third major radiation of amniotes. We localized the core MHC genomic region in two squamate species, the green anole (Anolis carolinensis) and brown anole (A. sagrei), and provide the first detailed characterization of the squamate MHC, including the presence and ordering of known MHC genes in these species and comparative assessments of genomic structure and composition in MHC regions. We find that the Anolis MHC, located on chromosome 2 in both species, contains homologs of many previously-identified mammalian MHC genes in a single core MHC region. The repetitive element composition in anole MHC regions was similar to those observed in mammals but had important distinctions, such as higher proportions of DNA transposons. Moreover, longer introns and intergenic regions result in a much larger squamate MHC region (11.7 Mb and 24.6 Mb in the green and brown anole, respectively). Evolutionary analyses of MHC homologs of anoles and other representative amniotes uncovered generally monophyletic relationships between species-specific homologs and a loss of the peptide-binding domain exon 2 in one of two mhc2β gene homologs of each anole species. Signals of diversifying selection in each anole species was evident across codons of mhc1, many of which appear functionally relevant given known structures of this protein from the green anole, chicken, and human. Altogether, our investigation fills a major gap in understanding of amniote MHC diversity and evolution and provides an important foundation for future squamate-specific or vertebrate-wide investigations of the MHC.
Collapse
Affiliation(s)
- Daren C. Card
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, United States
- Museum of Comparative Zoology, Harvard University, Cambridge, MA, United States
- *Correspondence: Daren C. Card,
| | - Andrew G. Van Camp
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, United States
- Museum of Comparative Zoology, Harvard University, Cambridge, MA, United States
| | - Trenten Santonastaso
- Department of Biological Sciences, University of New Orleans, New Orleans, LA, United States
| | | | - Nicola M. Anthony
- Department of Biological Sciences, University of New Orleans, New Orleans, LA, United States
| | - Scott V. Edwards
- Department of Organismic and Evolutionary Biology, Harvard University, Cambridge, MA, United States
- Museum of Comparative Zoology, Harvard University, Cambridge, MA, United States
| |
Collapse
|
5
|
Chromosome-scale genome assembly of the brown anole (Anolis sagrei), an emerging model species. Commun Biol 2022; 5:1126. [PMID: 36284162 PMCID: PMC9596491 DOI: 10.1038/s42003-022-04074-5] [Citation(s) in RCA: 9] [Impact Index Per Article: 4.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/27/2021] [Accepted: 10/06/2022] [Indexed: 12/12/2022] Open
Abstract
Rapid technological improvements are democratizing access to high quality, chromosome-scale genome assemblies. No longer the domain of only the most highly studied model organisms, now non-traditional and emerging model species can be genome-enabled using a combination of sequencing technologies and assembly software. Consequently, old ideas built on sparse sampling across the tree of life have recently been amended in the face of genomic data drawn from a growing number of high-quality reference genomes. Arguably the most valuable are those long-studied species for which much is already known about their biology; what many term emerging model species. Here, we report a highly complete chromosome-scale genome assembly for the brown anole, Anolis sagrei – a lizard species widely studied across a variety of disciplines and for which a high-quality reference genome was long overdue. This assembly exceeds the vast majority of existing reptile and snake genomes in contiguity (N50 = 253.6 Mb) and annotation completeness. Through the analysis of this genome and population resequence data, we examine the history of repetitive element accumulation, identify the X chromosome, and propose a hypothesis for the evolutionary history of fusions between autosomes and the X that led to the sex chromosomes of A. sagrei. A highly-complete chromosome-scale genome assembly of the brown anole, Anolis sagrei, provides insight into the evolution of sex chromosomes and is a crucial resource for this model lizard species.
Collapse
|
6
|
Provost K, Shue SY, Forcellati M, Smith BT. The Genomic Landscapes of Desert Birds Form over Multiple Time Scales. Mol Biol Evol 2022; 39:6711078. [PMID: 36134537 PMCID: PMC9577548 DOI: 10.1093/molbev/msac200] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/15/2022] Open
Abstract
Spatial models show that genetic differentiation between populations can be explained by factors ranging from geographic distance to environmental resistance across the landscape. However, genomes exhibit a landscape of differentiation, indicating that multiple processes may mediate divergence in different portions of the genome. We tested this idea by comparing alternative geographic predctors of differentiation in ten bird species that co-occur in Sonoran and Chihuahuan Deserts of North America. Using population-level genomic data, we described the genomic landscapes across species and modeled conditions that represented historical and contemporary mechanisms. The characteristics of genomic landscapes differed across species, influenced by varying levels of population structuring and admixture between deserts, and the best-fit models contrasted between the whole genome and partitions along the genome. Both historical and contemporary mechanisms were important in explaining genetic distance, but particularly past and current environments, suggesting that genomic evolution was modulated by climate and habitat There were also different best-ftit models across genomic partitions of the data, indicating that these regions capture different evolutionary histories. These results show that the genomic landscape of differentiation can be associated with alternative geographic factors operating on different portions of the genome, which reflect how heterogeneous patterns of genetic differentiation can evolve across species and genomes.
Collapse
Affiliation(s)
| | - Stephanie Yun Shue
- Bergen County Academies, Hackensack, NJ, USA,Biological Sciences, University of California Berkeley, Berkeley, CA, USA
| | - Meghan Forcellati
- Bergen County Academies, Hackensack, NJ, USA,Ecology, Evolution, and Environmental Biology, Columbia University, New York, NY, USA
| | - Brian Tilston Smith
- Department of Ornithology, American Museum of Natural History, New York, NY, USA
| |
Collapse
|
7
|
Lyman RA, Edwards CE. Revisiting the comparative phylogeography of unglaciated eastern North America: 15 years of patterns and progress. Ecol Evol 2022; 12:e8827. [PMID: 35475178 PMCID: PMC9019306 DOI: 10.1002/ece3.8827] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2021] [Revised: 03/17/2022] [Accepted: 03/18/2022] [Indexed: 11/28/2022] Open
Abstract
In a landmark comparative phylogeographic study, “Comparative phylogeography of unglaciated eastern North America,” Soltis et al. (Molecular Ecology, 2006, 15, 4261) identified geographic discontinuities in genetic variation shared across taxa occupying unglaciated eastern North America and proposed several common biogeographical discontinuities related to past climate fluctuations and geographic barriers. Since 2006, researchers have published many phylogeographical studies and achieved many advances in genotyping and analytical techniques; however, it is unknown how this work has changed our understanding of the factors shaping the phylogeography of eastern North American taxa. We analyzed 184 phylogeographical studies of eastern North American taxa published between 2007 and 2019 to evaluate: (1) the taxonomic focus of studies and whether a previously detected taxonomic bias towards studies focused on vertebrates has changed over time, (2) the extent to which studies have adopted genotyping technologies that improve the resolution of genetic groups (i.e., NGS DNA sequencing) and analytical approaches that facilitate hypothesis‐testing (i.e., divergence time estimation and niche modeling), and (3) whether new studies support the hypothesized biogeographic discontinuities proposed by Soltis et al. (Molecular Ecology, 2006, 15, 4261) or instead support new, previously undetected discontinuities. We observed little change in taxonomic focus over time, with studies still biased toward vertebrates. Although many technological and analytical advances became available during the period, uptake was slow and they were employed in only a small proportion of studies. We found variable support for previously identified discontinuities and identified one new recurrent discontinuity. However, the limited resolution and taxonomic breadth of many studies hindered our ability to clarify the most important climatological or geographical factors affecting taxa in the region. Broadening the taxonomic focus to include more non‐vertebrate taxa, employing technologies that improve genetic resolution, and using analytical approaches that improve hypothesis testing are necessary to strengthen our inference of the forces shaping the phylogeography of eastern North America.
Collapse
Affiliation(s)
- Rachel Ann Lyman
- Ecology, Evolution, and Population Biology Program Washington University in St. Louis St. Louis Missouri USA
- Center for Conservation and Sustainable Development Missouri Botanical Garden St. Louis Missouri USA
| | - Christine E. Edwards
- Center for Conservation and Sustainable Development Missouri Botanical Garden St. Louis Missouri USA
| |
Collapse
|
8
|
Versoza CJ, Rivera JA, Rosenblum EB, Vital-García C, Hews DK, Pfeifer SP. The recombination landscapes of spiny lizards (genus Sceloporus). G3 GENES|GENOMES|GENETICS 2022; 12:6433156. [PMID: 34878100 PMCID: PMC9210290 DOI: 10.1093/g3journal/jkab402] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 07/07/2021] [Accepted: 11/14/2021] [Indexed: 11/16/2022]
Abstract
Despite playing a critical role in evolutionary processes and outcomes, relatively little is known about rates of recombination in the vast majority of species, including squamate reptiles—the second largest order of extant vertebrates, many species of which serve as important model organisms in evolutionary and ecological studies. This paucity of data has resulted in limited resolution on questions related to the causes and consequences of rate variation between species and populations, the determinants of within-genome rate variation, as well as the general tempo of recombination rate evolution on this branch of the tree of life. In order to address these questions, it is thus necessary to begin broadening our phylogenetic sampling. We here provide the first fine-scale recombination maps for two species of spiny lizards, Sceloporus jarrovii and Sceloporus megalepidurus, which diverged at least 12 Mya. As might be expected from similarities in karyotype, population-scaled recombination landscapes are largely conserved on the broad-scale. At the same time, considerable variation exists at the fine-scale, highlighting the importance of incorporating species-specific recombination maps in future population genomic studies.
Collapse
Affiliation(s)
- Cyril J Versoza
- School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ 85281, USA
| | - Julio A Rivera
- School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
| | - Erica Bree Rosenblum
- Department of Environmental Science, Policy and Management, University of California, Berkeley, Berkeley, CA 94720, USA
| | - Cuauhcihuatl Vital-García
- Departamento de Ciencias Veterinarias, Programa de Maestría en Ciencia Animal, Universidad Autónoma de Ciudad Juárez México, Chihuahua 32315, Mexico
| | - Diana K Hews
- Department of Biology, Indiana State University, Terre Haute, IN 47809, USA
| | - Susanne P Pfeifer
- School of Life Sciences, Arizona State University, Tempe, AZ 85281, USA
- Center for Evolution and Medicine, Arizona State University, Tempe, AZ 85281, USA
- Center for Mechanisms of Evolution, Arizona State University, Tempe, AZ 85281, USA
| |
Collapse
|
9
|
Hench K, Helmkampf M, McMillan WO, Puebla O. Rapid radiation in a highly diverse marine environment. Proc Natl Acad Sci U S A 2022; 119:e2020457119. [PMID: 35042790 PMCID: PMC8794831 DOI: 10.1073/pnas.2020457119] [Citation(s) in RCA: 4] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/30/2020] [Accepted: 11/18/2021] [Indexed: 12/16/2022] Open
Abstract
Rapid diversification is often observed when founding species invade isolated or newly formed habitats that provide ecological opportunity for adaptive radiation. However, most of the Earth's diversity arose in diverse environments where ecological opportunities appear to be more constrained. Here, we present a striking example of a rapid radiation in a highly diverse marine habitat. The hamlets, a group of reef fishes from the wider Caribbean, have radiated into a stunning diversity of color patterns but show low divergence across other ecological axes. Although the hamlet lineage is ∼26 My old, the radiation appears to have occurred within the last 10,000 generations in a burst of diversification that ranks among the fastest in fishes. As such, the hamlets provide a compelling backdrop to uncover the genomic elements associated with phenotypic diversification and an excellent opportunity to build a broader comparative framework for understanding the drivers of adaptive radiation. The analysis of 170 genomes suggests that color pattern diversity is generated by different combinations of alleles at a few large-effect loci. Such a modular genomic architecture of diversification has been documented before in Heliconius butterflies, capuchino finches, and munia finches, three other tropical radiations that took place in highly diverse and complex environments. The hamlet radiation also occurred in a context of high effective population size, which is typical of marine populations. This allows for the accumulation of new variants through mutation and the retention of ancestral genetic variation, both of which appear to be important in this radiation.
Collapse
Affiliation(s)
- Kosmas Hench
- Ecology Department, Leibniz Centre for Tropical Marine Research, 28359 Bremen, Germany;
| | - Martin Helmkampf
- Ecology Department, Leibniz Centre for Tropical Marine Research, 28359 Bremen, Germany
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Apartado Postal 0843-03092, Republic of Panama
| | - Oscar Puebla
- Ecology Department, Leibniz Centre for Tropical Marine Research, 28359 Bremen, Germany;
- Smithsonian Tropical Research Institute, Apartado Postal 0843-03092, Republic of Panama
- Institute for Chemistry and Biology of the Marine Environment, 26111 Oldenburg, Germany
- Marine Evolutionary Ecology, GEOMAR Helmholtz Centre for Ocean Research Kiel, 24105 Kiel, Germany
| |
Collapse
|
10
|
Wollenberg Valero KC, Garcia-Porta J, Irisarri I, Feugere L, Bates A, Kirchhof S, Jovanović Glavaš O, Pafilis P, Samuel SF, Müller J, Vences M, Turner AP, Beltran-Alvarez P, Storey KB. Functional genomics of abiotic environmental adaptation in lacertid lizards and other vertebrates. J Anim Ecol 2021; 91:1163-1179. [PMID: 34695234 DOI: 10.1111/1365-2656.13617] [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: 03/05/2021] [Accepted: 09/27/2021] [Indexed: 11/27/2022]
Abstract
Understanding the genomic basis of adaptation to different abiotic environments is important in the context of climate change and resulting short-term environmental fluctuations. Using functional and comparative genomics approaches, we here investigated whether signatures of genomic adaptation to a set of environmental parameters are concentrated in specific subsets of genes and functions in lacertid lizards and other vertebrates. We first identify 200 genes with signatures of positive diversifying selection from transcriptomes of 24 species of lacertid lizards and demonstrate their involvement in physiological and morphological adaptations to climate. To understand how functionally similar these genes are to previously predicted candidate functions for climate adaptation and to compare them with other vertebrate species, we then performed a meta-analysis of 1,100 genes under selection obtained from -omics studies in vertebrate species adapted to different abiotic factors. We found that the vertebrate gene set formed a tightly connected interactome, which was to 23% enriched in previously predicted functions of adaptation to climate, and to a large part (18%) involved in organismal stress response. We found a much higher degree of identical genes being repeatedly selected among different animal groups (43.6%), and of functional similarity and post-translational modifications than expected by chance, and no clear functional division between genes used for ectotherm and endotherm physiological strategies. In total, 171 out of 200 genes of Lacertidae were part of this network. These results highlight an important role of a comparatively small set of genes and their functions in environmental adaptation and narrow the set of candidate pathways and markers to be used in future research on adaptation and stress response related to climate change.
Collapse
Affiliation(s)
| | - Joan Garcia-Porta
- Department of Biology, Washington University in St. Louis, St. Louis, MO, USA
| | - Iker Irisarri
- Department of Applied Bioinformatics, Institute for Microbiology and Genetics, University of Göttingen, Göttingen, Germany.,Campus Institut Data Science (CIDAS), Göttingen, Germany
| | - Lauric Feugere
- Department of Biological and Marine Sciences, University of Hull, Kingston-Upon-Hull, UK
| | - Adam Bates
- Department of Biological and Marine Sciences, University of Hull, Kingston-Upon-Hull, UK
| | - Sebastian Kirchhof
- Museum für Naturkunde, Leibniz Institute for Evolution and Biodiversity Science, Berlin, Germany.,New York University Abu Dhabi, Abu Dhabi, United Arab Emirates
| | | | - Panayiotis Pafilis
- Section of Zoology and Marine Biology, Department of Biology, National and Kapodistrian University of Athens, Athens, Greece
| | - Sabrina F Samuel
- Department of Biomedical Sciences, University of Hull, Kingston-Upon-Hull, UK
| | - Johannes Müller
- Museum für Naturkunde, Leibniz Institute for Evolution and Biodiversity Science, Berlin, Germany
| | - Miguel Vences
- Zoological Institute, Braunschweig University of Technology, Braunschweig, Germany
| | - Alexander P Turner
- Department of Computer Science, University of Nottingham, Nottingham, UK
| | | | | |
Collapse
|
11
|
Alam O, Gutaker RM, Wu CC, Hicks KA, Bocinsky K, Castillo CC, Acabado S, Fuller D, d'Alpoim Guedes JA, Hsing YI, Purugganan MD. Genome analysis traces regional dispersal of rice in Taiwan and Southeast Asia. Mol Biol Evol 2021; 38:4832-4846. [PMID: 34240169 PMCID: PMC8557449 DOI: 10.1093/molbev/msab209] [Citation(s) in RCA: 10] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/13/2022] Open
Abstract
The dispersal of rice (Oryza sativa) following domestication influenced massive social and cultural changes across South, East, and Southeast Asia. The history of dispersal across islands of Southeast Asia, and the role of Taiwan and the Austronesian expansion in this process remain largely unresolved. Here, we reconstructed the routes of dispersal of O. sativa ssp. japonica rice through Taiwan and the northern Philippines using whole-genome re-sequencing of indigenous rice landraces coupled with archaeological and paleoclimate data. Our results indicate that japonica rice found in the northern Philippines diverged from Indonesian landraces as early as 3500 BP. In contrast, rice cultivated by the indigenous peoples of the Taiwanese mountains has complex origins. It comprises two distinct populations, each best explained as a result of admixture between temperate japonica that presumably came from northeast Asia, and tropical japonica from the northern Philippines and mainland Southeast Asia respectively. We find that the temperate japonica component of these indigenous Taiwan populations diverged from northeast Asia subpopulations at about 2600 BP, while gene flow from the northern Philippines occurred before ∼1300 years BP. This coincides with a period of intensified trade established across the South China Sea. Finally, we find evidence for positive selection acting on distinct genomic regions in different rice subpopulations, indicating local adaptation associated with the spread of japonica rice.
Collapse
Affiliation(s)
- Ornob Alam
- Center for Genomics and Systems Biology, New York University, New York, NY 10003 USA
| | - Rafal M Gutaker
- Center for Genomics and Systems Biology, New York University, New York, NY 10003 USA.,Royal Botanic Garden, Kew, Richmond, London, TW9 3AE UK
| | - Cheng-Chieh Wu
- Institute of Plant and Microbial Biology, Academia Sinica, Nankang, Taipei 115, Taiwan.,Institute of Plant Biology, National Taiwan University, Taipei 106, Taiwan
| | - Karen A Hicks
- Department of Biology, Kenyon College, Gambier, Ohio 43022 USA
| | | | | | - Stephen Acabado
- Department of Anthropology, University of California, Los Angeles, CA USA
| | - Dorian Fuller
- Institute of Archaeology, University College London, London, United Kingdom.,School of Cultural Heritage, North-West University, Xi'an, China
| | - Jade A d'Alpoim Guedes
- Department of Anthropology and Scripps Institution of Oceanography, University of California, San Diego, CA, USA
| | - Yue-Ie Hsing
- Institute of Plant and Microbial Biology, Academia Sinica, Nankang, Taipei 115, Taiwan
| | - Michael D Purugganan
- Center for Genomics and Systems Biology, New York University, New York, NY 10003 USA.,Institute for the Study of the Ancient World, New York University, New York, NY 10028 USA
| |
Collapse
|
12
|
Araya-Donoso R, San Juan E, Tamburrino Í, Lamborot M, Veloso C, Véliz D. Integrating genetics, physiology and morphology to study desert adaptation in a lizard species. J Anim Ecol 2021; 91:1148-1162. [PMID: 34048024 DOI: 10.1111/1365-2656.13546] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/28/2020] [Accepted: 05/24/2021] [Indexed: 11/28/2022]
Abstract
Integration of multiple approaches is key to understand the evolutionary processes of local adaptation and speciation. Reptiles have successfully colonized desert environments, that is, extreme and arid conditions that constitute a strong selective pressure on organisms. Here, we studied genomic, physiological and morphological variations of the lizard Liolaemus fuscus to detect adaptations to the Atacama Desert. By comparing populations of L. fuscus inhabiting the Atacama Desert with populations from the Mediterranean forests from central Chile, we aimed at characterizing features related to desert adaptation. We combined ddRAD sequencing with physiological (evaporative water loss, metabolic rate and selected temperature) and morphological (linear and geometric morphometrics) measurements. We integrated the genomic and phenotypic data using redundancy analyses. Results showed strong genetic divergence, along with a high number of fixed loci between desert and forest populations. Analyses detected 110 fixed and 30 outlier loci located within genes, from which 43 were in coding regions, and 12 presented non-synonymous mutations. The candidate genes were associated with cellular membrane and development. Desert lizards presented lower evaporative water loss than those from the forest. Morphological data showed that desert lizards had smaller body size, different allometry, larger eyeballs and more dorsoventrally compressed heads. Our results suggest incipient speciation between desert and forest populations. The adaptive signal must be cautiously interpreted since genetic drift could also contribute to the divergence pattern. Nonetheless, we propose water and resource availability, and changes in habitat structure, as the most relevant challenges for desert reptiles. This study provides insights of the mechanisms that allow speciation as well as desert adaptation in reptiles at multiple levels, and highlights the benefit of integrating independent evidence.
Collapse
Affiliation(s)
- Raúl Araya-Donoso
- Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile.,Núcleo Milenio de Ecología y Manejo Sustentable de Islas Oceánicas (ESMOI), Departamento de Biología Marina, Universidad Católica del Norte, Coquimbo, Chile.,School of Life Sciences, Arizona State University, Tempe, AZ, USA
| | - Esteban San Juan
- Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
| | - Ítalo Tamburrino
- Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
| | - Madeleine Lamborot
- Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
| | - Claudio Veloso
- Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile
| | - David Véliz
- Departamento de Ciencias Ecológicas, Facultad de Ciencias, Universidad de Chile, Santiago, Chile.,Núcleo Milenio de Ecología y Manejo Sustentable de Islas Oceánicas (ESMOI), Departamento de Biología Marina, Universidad Católica del Norte, Coquimbo, Chile
| |
Collapse
|
13
|
Emami-Khoyi A, Knapp IS, Monsanto DM, van Vuuren BJ, Toonen RJ, Teske PR. Genomic divergence and differential gene expression between crustacean ecotypes across a marine thermal gradient. Mar Genomics 2021; 58:100847. [PMID: 33637426 DOI: 10.1016/j.margen.2021.100847] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/10/2020] [Revised: 12/15/2020] [Accepted: 01/29/2021] [Indexed: 11/30/2022]
Abstract
Environmental gradients between marine biogeographical provinces separate distinct faunal communities. However, the absence of absolute dispersal barriers allows numerous species to occur on both sides of such boundaries. While the regional populations of such widespread species are often morphologically indistinguishable from each other, genetic evidence suggests that they represent unique ecotypes, and likely even cryptic species, that may be uniquely adapted to their local environment. Here, we explored genomic divergence in four sympatric southern African decapod crustaceans whose ranges span the boundary between the cool-temperate west coast (south-eastern Atlantic) and the warm-temperate south coast (south-western Indian Ocean) near the southern tip of the African continent. Using genome-wide data, we found that all four species comprise distinct west coast and south coast ecotypes, with molecular dating suggesting divergence during the Pleistocene. Transcriptomic data from the hepatopancreas of twelve specimens of one of these species, the mudprawn Upogebia africana, which were exposed to either 10 °C or 20 °C, showed a clear difference in gene expression profiles between the west- and south coast ecotypes. This difference was particularly clear at 10 °C, where individuals from the south coast experienced a 'transcriptomic shock'. This low temperature is more typical of the west coast during upwelling events, and the physiological stress experienced by the south coast ecotype under such conditions may explain its absence from that coastline. Our results shed new light on the processes involved in driving genomic divergence and incipient speciation along coastlines with porous dispersal barriers.
Collapse
Affiliation(s)
- Arsalan Emami-Khoyi
- Centre for Ecological Genomics and Wildlife Conservation, Department of Zoology, University of Johannesburg, Auckland Park 2006, South Africa.
| | - Ingrid S Knapp
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, Kāne'ohe, Hawai'i, Honolulu, HI, USA
| | - Daniela M Monsanto
- Centre for Ecological Genomics and Wildlife Conservation, Department of Zoology, University of Johannesburg, Auckland Park 2006, South Africa
| | - Bettine Jansen van Vuuren
- Centre for Ecological Genomics and Wildlife Conservation, Department of Zoology, University of Johannesburg, Auckland Park 2006, South Africa
| | - Robert J Toonen
- Hawai'i Institute of Marine Biology, University of Hawai'i at Mānoa, Kāne'ohe, Hawai'i, Honolulu, HI, USA
| | - Peter R Teske
- Centre for Ecological Genomics and Wildlife Conservation, Department of Zoology, University of Johannesburg, Auckland Park 2006, South Africa.
| |
Collapse
|
14
|
Momigliano P, Florin AB, Merilä J. Biases in Demographic Modeling Affect Our Understanding of Recent Divergence. Mol Biol Evol 2021; 38:2967-2985. [PMID: 33624816 PMCID: PMC8233503 DOI: 10.1093/molbev/msab047] [Citation(s) in RCA: 14] [Impact Index Per Article: 4.7] [Reference Citation Analysis] [Abstract] [Key Words] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Testing among competing demographic models of divergence has become an important component of evolutionary research in model and non-model organisms. However, the effect of unaccounted demographic events on model choice and parameter estimation remains largely unexplored. Using extensive simulations, we demonstrate that under realistic divergence scenarios, failure to account for population size (Ne) changes in daughter and ancestral populations leads to strong biases in divergence time estimates as well as model choice. We illustrate these issues reconstructing the recent demographic history of North Sea and Baltic Sea turbots (Scophthalmus maximus) by testing 16 isolation with migration (IM) and 16 secondary contact (SC) scenarios, modeling changes in Ne as well as the effects of linked selection and barrier loci. Failure to account for changes in Ne resulted in selecting SC models with long periods of strict isolation and divergence times preceding the formation of the Baltic Sea. In contrast, models accounting for Ne changes suggest recent (<6 kya) divergence with constant gene flow. We further show how interpreting genomic landscapes of differentiation can help discerning among competing models. For example, in the turbot data, islands of differentiation show signatures of recent selective sweeps, rather than old divergence resisting secondary introgression. The results have broad implications for the study of population divergence by highlighting the potential effects of unmodeled changes in Ne on demographic inference. Tested models should aim at representing realistic divergence scenarios for the target taxa, and extreme caution should always be exercised when interpreting results of demographic modeling.
Collapse
Affiliation(s)
- Paolo Momigliano
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, University of Helsinki, Helsinki, Finland
| | - Ann-Britt Florin
- Department of Aquatic Resources, Institute of Coastal Research, Swedish University of Agricultural Sciences, Öregrund, Sweden
| | - Juha Merilä
- Ecological Genetics Research Unit, Organismal and Evolutionary Biology Research Programme, University of Helsinki, Helsinki, Finland.,Division of Ecology and Biodiversity, Faculty of Science, The University of Hong Kong, Hong Kong SAR
| |
Collapse
|
15
|
Svardal H, Salzburger W, Malinsky M. Genetic Variation and Hybridization in Evolutionary Radiations of Cichlid Fishes. Annu Rev Anim Biosci 2020; 9:55-79. [PMID: 33197206 DOI: 10.1146/annurev-animal-061220-023129] [Citation(s) in RCA: 9] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/09/2022]
Abstract
Evolutionary radiations are responsible for much of the variation in biodiversity across taxa. Cichlid fishes are well known for spectacular evolutionary radiations, as they have repeatedly evolved into large and phenotypically diverse arrays of species. Cichlid genomes carry signatures of past events and, at the same time, are the substrate for ongoing evolution. We survey genome-wide data and the available literature covering 438 cichlid populations (412 species) across multiple radiations to synthesize information about patterns and sharing of genetic variation. Nucleotide diversity within species is low in cichlids, with 92% of surveyed populations having less diversity than the median value found in other vertebrates. Divergence within radiations is also low, and a large proportion of variation is shared among species due to incomplete lineage sorting and widespread hybridization. Population genetics therefore provides a suitable conceptual framework for evolutionary genomic studies of cichlid radiations. We focus in detail on the roles of hybridization in shaping the patterns of genetic variation and in promoting cichlid diversification.
Collapse
Affiliation(s)
- Hannes Svardal
- Department of Biology, University of Antwerp, 2020 Antwerp, Belgium; .,Naturalis Biodiversity Center, 2333 Leiden, The Netherlands
| | - Walter Salzburger
- Zoological Institute, University of Basel, 4051 Basel, Switzerland; ,
| | - Milan Malinsky
- Zoological Institute, University of Basel, 4051 Basel, Switzerland; ,
| |
Collapse
|
16
|
Bourgeois Y, Ruggiero RP, Hariyani I, Boissinot S. Disentangling the determinants of transposable elements dynamics in vertebrate genomes using empirical evidences and simulations. PLoS Genet 2020; 16:e1009082. [PMID: 33017388 PMCID: PMC7561263 DOI: 10.1371/journal.pgen.1009082] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/23/2020] [Revised: 10/15/2020] [Accepted: 08/25/2020] [Indexed: 12/14/2022] Open
Abstract
The interactions between transposable elements (TEs) and their hosts constitute one of the most profound co-evolutionary processes found in nature. The population dynamics of TEs depends on factors specific to each TE families, such as the rate of transposition and insertional preference, the demographic history of the host and the genomic landscape. How these factors interact has yet to be investigated holistically. Here we are addressing this question in the green anole (Anolis carolinensis) whose genome contains an extraordinary diversity of TEs (including non-LTR retrotransposons, SINEs, LTR-retrotransposons and DNA transposons). We observed a positive correlation between recombination rate and frequency of TEs and densities for LINEs, SINEs and DNA transposons. For these elements, there was a clear impact of demography on TE frequency and abundance, with a loss of polymorphic elements and skewed frequency spectra in recently expanded populations. On the other hand, some LTR-retrotransposons displayed patterns consistent with a very recent phase of intense amplification. To determine how demography, genomic features and intrinsic properties of TEs interact we ran simulations using SLiM3. We determined that i) short TE insertions are not strongly counter-selected, but long ones are, ii) neutral demographic processes, linked selection and preferential insertion may explain positive correlations between average TE frequency and recombination, iii) TE insertions are unlikely to have been massively recruited in recent adaptation. We demonstrate that deterministic and stochastic processes have different effects on categories of TEs and that a combination of empirical analyses and simulations can disentangle these mechanisms.
Collapse
Affiliation(s)
- Yann Bourgeois
- School of Biological Sciences, University of Portsmouth, Portsmouth, United Kingdom
- New York University Abu Dhabi, Saadiyat Island Campus, Abu Dhabi, United Arab Emirates
- * E-mail: (YB); (SB)
| | - Robert P. Ruggiero
- New York University Abu Dhabi, Saadiyat Island Campus, Abu Dhabi, United Arab Emirates
- Department of Biology, Southeast Missouri State University, Cape Girardeau, MO, United States of America
| | - Imtiyaz Hariyani
- New York University Abu Dhabi, Saadiyat Island Campus, Abu Dhabi, United Arab Emirates
| | - Stéphane Boissinot
- New York University Abu Dhabi, Saadiyat Island Campus, Abu Dhabi, United Arab Emirates
- * E-mail: (YB); (SB)
| |
Collapse
|
17
|
Zhang J, Nie C, Li X, Ning Z, Chen Y, Jia Y, Han J, Wang L, Lv X, Yang W, Qu L. Genome-Wide Population Genetic Analysis of Commercial, Indigenous, Game, and Wild Chickens Using 600K SNP Microarray Data. Front Genet 2020; 11:543294. [PMID: 33101376 PMCID: PMC7545075 DOI: 10.3389/fgene.2020.543294] [Citation(s) in RCA: 26] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/16/2020] [Accepted: 09/02/2020] [Indexed: 12/13/2022] Open
Abstract
Following chicken domestication, diversified chicken breeds were developed by both natural and artificial selection, which led to the accumulation of abundant genetic and phenotypic variations, making chickens an ideal genetic research model. To better understand the genetic structure of chicken breeds under different selection pressures, we genotyped various chicken populations with specific selection targets, including indigenous, commercial, gamecock, and wild ancestral chickens, using the 600K SNP array. We analyzed the population structure, genetic relationships, run of homozygosity (ROH), effective population number (Ne), and other genetic parameters. The wild ancestral population, red junglefowl (RJF), possessed the highest diversity, in comparison with all other domesticated populations, which was supported by linkage disequilibrium decay (LD), effective population number, and ROH analyses. The gamecock breeds, which were subjected to stronger male-biased selection for fighting-related traits, also presented higher variation than the commercial and indigenous breeds. Admixture analysis also indicated that game breed is a relatively independent branch of Chinese local breeds. Following intense selection for reproductive and productive traits, the commercial lines showed the least diversity. We also observed that the European local chickens had lower genetic variation than the Chinese local breeds, which could be attributed to the shorter history of the European breed. ROH were present in a breed specific manner and 191 ROH island were detected on four groups (commercial, local, game and wild chickens). These ROH islands were involved in egg production, growth and silky feathers and other traits. Moreover, we estimated the effective sex ratio of these breeds to demonstrate the change in the ratio of the two sexes. We found that commercial chickens had a greater sex imbalance between females and males. The commercial lines showed the highest female-to-male ratios. Interestingly, RJF comprised a greater proportion of males than females. Our results show the population genetics of chickens under selection pressures, and can aid in the development of better conservation strategies for different chicken breeds.
Collapse
Affiliation(s)
- Jinxin Zhang
- Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Changsheng Nie
- Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Xinghua Li
- Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Zhonghua Ning
- Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| | - Yu Chen
- Beijing Municipal General Station of Animal Science, Beijing, China
| | - Yaxiong Jia
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Jianlin Han
- Institute of Animal Sciences, Chinese Academy of Agricultural Sciences, Beijing, China
| | - Liang Wang
- Beijing Municipal General Station of Animal Science, Beijing, China
| | - Xueze Lv
- Beijing Municipal General Station of Animal Science, Beijing, China
| | - Weifang Yang
- Beijing Municipal General Station of Animal Science, Beijing, China
| | - Lujiang Qu
- Department of Animal Genetics and Breeding, National Engineering Laboratory for Animal Breeding, College of Animal Science and Technology, China Agricultural University, Beijing, China
| |
Collapse
|
18
|
Bourgeois Y, Boissinot S. Selection at behavioural, developmental and metabolic genes is associated with the northward expansion of a successful tropical colonizer. Mol Ecol 2019; 28:3523-3543. [PMID: 31233650 DOI: 10.1111/mec.15162] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/03/2018] [Accepted: 05/28/2019] [Indexed: 02/06/2023]
Abstract
What makes a species able to colonize novel environments? This question is key to understand the dynamics of adaptive radiations and ecological niche shifts, but the mechanisms that underlie expansion into novel habitats remain poorly understood at a genomic scale. Lizards from the genus Anolis are typically tropical, and the green anole (Anolis carolinensis) constitutes an exception since it expanded into temperate North America from subtropical Florida. Thus, we used the green anole as a model to investigate signatures of selection associated with colonization of a new environment, namely temperate North America. To this end, we analysed 29 whole-genome sequences, covering the entire native range of the species. We used a combination of recent methods to quantify both positive and balancing selection in northern populations, including FST outlier methods, machine learning and ancestral recombination graphs. We naively scanned for genes of interest and assessed the overlap between multiple tests. Strikingly, we identified many genes involved in behaviour, suggesting that the recent successful colonization of northern environments may have been linked to behavioural shifts as well as physiological adaptation. Using a candidate genes strategy, we determined that genes involved in response to cold or behaviour displayed more frequently signals of selection, while controlling for local recombination rate, gene clustering and gene length. In addition, we found signatures of balancing selection at immune genes in all investigated genetic groups, but also at genes involved in neuronal and anatomical development.
Collapse
Affiliation(s)
- Yann Bourgeois
- New York University Abu Dhabi, Abu Dhabi, United Arab Emirates
| | | |
Collapse
|
19
|
Capblancq T, Mavárez J, Rioux D, Després L. Speciation with gene flow: Evidence from a complex of alpine butterflies ( Coenonympha, Satyridae). Ecol Evol 2019; 9:6444-6457. [PMID: 31236234 PMCID: PMC6580291 DOI: 10.1002/ece3.5220] [Citation(s) in RCA: 12] [Impact Index Per Article: 2.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/04/2019] [Revised: 04/14/2019] [Accepted: 04/14/2019] [Indexed: 12/19/2022] Open
Abstract
Until complete reproductive isolation is achieved, the extent of differentiation between two diverging lineages is the result of a dynamic equilibrium between genetic isolation and mixing. This is especially true for hybrid taxa, for which the degree of isolation in regard to their parental species is decisive in their capacity to rise as a new and stable entity. In this work, we explored the past and current patterns of hybridization and divergence within a complex of closely related butterflies in the genus Coenonympha in which two alpine species, C. darwiniana and C. macromma, have been shown to result from hybridization between the also alpine C. gardetta and the lowland C. arcania. By testing alternative scenarios of divergence among species, we show that gene flow has been uninterrupted throughout the speciation process, although leading to different degrees of current genetic isolation between species in contact zones depending on the pair considered. Nonetheless, at broader geographic scale, analyses reveal a clear genetic differentiation between hybrid lineages and their parental species, pointing out to an advanced stage of the hybrid speciation process. Finally, the positive correlation observed between ecological divergence and genetic isolation among these butterflies suggests a potential role for ecological drivers during their speciation processes.
Collapse
Affiliation(s)
- Thibaut Capblancq
- Laboratoire d’Écologie AlpineUMR5553CNRS‐Université Grenoble AlpesGrenobleFrance
| | - Jesús Mavárez
- Laboratoire d’Écologie AlpineUMR5553CNRS‐Université Grenoble AlpesGrenobleFrance
| | - Delphine Rioux
- Laboratoire d’Écologie AlpineUMR5553CNRS‐Université Grenoble AlpesGrenobleFrance
| | - Laurence Després
- Laboratoire d’Écologie AlpineUMR5553CNRS‐Université Grenoble AlpesGrenobleFrance
| |
Collapse
|