1
|
Zhang Z, Liu G, Li M. Incomplete lineage sorting and gene flow within Allium (Amayllidaceae). Mol Phylogenet Evol 2024; 195:108054. [PMID: 38471599 DOI: 10.1016/j.ympev.2024.108054] [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: 11/30/2023] [Revised: 02/01/2024] [Accepted: 03/07/2024] [Indexed: 03/14/2024]
Abstract
The phylogeny and systematics of the genus Allium have been studied with a variety of diverse data types, including an increasing amount of molecular data. However, strong phylogenetic discordance and high levels of uncertainty have prevented the identification of a consistent phylogeny. The difficulty in establishing phylogenetic consensus and evidence for genealogical discordance make Allium a compelling test case to assess the relative contribution of incomplete lineage sorting (ILS), gene flow and gene tree estimation error on phylogenetic reconstruction. In this study, we obtained 75 transcriptomes of 38 Allium species across 10 subgenera. Whole plastid genome, single copy genes and consensus CDS were generated to estimate phylogenetic trees both using coalescence and concatenation methods. Multiple approaches including coalescence simulation, quartet sampling, reticulate network inference, sequence simulation, theta of ILS and reticulation index were carried out across the CDS gene trees to investigate the degrees of ILS, gene flow and gene tree estimation error. Afterward, a regression analysis was used to test the relative contributions of each of these forms of uncertainty to the final phylogeny. Despite extensive topological discordance among gene trees, we found a fully supported species tree that agrees with the most of well-accepted relationships and establishes monophyly of the genus Allium. We presented clear evidence for substantial ILS across the phylogeny of Allium. Further, we identified two ancient hybridization events for the formation of the second evolutionary line and subg. Butomissa as well as several introgression events between recently diverged species. Our regression analysis revealed that gene tree inference error and gene flow were the two most dominant factors explaining for the overall gene tree variation, with the difficulty in disentangling the effects of ILS and gene tree estimation error due to a positive correlation between them. Based on our efforts to mitigate the methodological errors in reconstructing trees, we believed ILS and gene flow are two principal reasons for the oft-reported phylogenetic heterogeneity of Allium. This study presents a strongly-supported and well-resolved phylogenetic backbone for the sampled Allium species, and exemplifies how to untangle heterogeneity in phylogenetic signal and reconstruct the true evolutionary history of the target taxa.
Collapse
Affiliation(s)
- ZengZhu Zhang
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, People's Republic of China
| | - Gang Liu
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, People's Republic of China
| | - Minjie Li
- State Key Laboratory of Herbage Improvement and Grassland Agro-ecosystems, College of Ecology, Lanzhou University, Lanzhou 730000, People's Republic of China.
| |
Collapse
|
2
|
Liu Y, Xiao W, Wang F, Wang Y, Dong Y, Nie W, Tan C, An S, Chang E, Jiang Z, Wang J, Jia Z. Adaptive divergence, historical population dynamics, and simulation of suitable distributions for Picea Meyeri and P. Mongolica at the whole-genome level. BMC PLANT BIOLOGY 2024; 24:479. [PMID: 38816690 PMCID: PMC11137980 DOI: 10.1186/s12870-024-05166-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/17/2024] [Accepted: 05/17/2024] [Indexed: 06/01/2024]
Abstract
The taxonomic classification of Picea meyeri and P. mongolica has long been controversial. To investigate the genetic relatedness, evolutionary history, and population history dynamics of these species, genotyping-by-sequencing (GBS) technology was utilized to acquire whole-genome single nucleotide polymorphism (SNP) markers, which were subsequently used to assess population structure, population dynamics, and adaptive differentiation. Phylogenetic and population structural analyses at the genomic level indicated that although the ancestor of P. mongolica was a hybrid of P. meyeri and P. koraiensis, P. mongolica is an independent Picea species. Additionally, P. mongolica is more closely related to P. meyeri than to P. koraiensis, which is consistent with its geographic distribution. There were up to eight instances of interspecific and intraspecific gene flow between P. meyeri and P. mongolica. The P. meyeri and P. mongolica effective population sizes generally decreased, and Maxent modeling revealed that from the Last Glacial Maximum (LGM) to the present, their habitat areas decreased initially and then increased. However, under future climate scenarios, the habitat areas of both species were projected to decrease, especially under high-emission scenarios, which would place P. mongolica at risk of extinction and in urgent need of protection. Local adaptation has promoted differentiation between P. meyeri and P. mongolica. Genotype‒environment association analysis revealed 96,543 SNPs associated with environmental factors, mainly related to plant adaptations to moisture and temperature. Selective sweeps revealed that the selected genes among P. meyeri, P. mongolica and P. koraiensis are primarily associated in vascular plants with flowering, fruit development, and stress resistance. This research enhances our understanding of Picea species classification and provides a basis for future genetic improvement and species conservation efforts.
Collapse
Affiliation(s)
- Yifu Liu
- Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, 100091, China
| | - Wenfa Xiao
- Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, 100091, China
| | - Fude Wang
- Heilongjiang Forestry Research Institute, Harbin, 150080, China
| | - Ya Wang
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, Beijing, 100091, China
| | - Yao Dong
- Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, 100091, China
| | - Wen Nie
- Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, 100091, China
| | - Cancan Tan
- Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, 100091, China
| | - Sanping An
- Research Institute of Forestry of Xiaolong Mountain, Gansu Provincial Key Laboratory of Secondary Forest Cultivation, Tianshui, 741022, China
| | - Ermei Chang
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, Beijing, 100091, China
| | - Zeping Jiang
- Key Laboratory of Forest Ecology and Environment of National Forestry and Grassland Administration, Ecology and Nature Conservation Institute, Chinese Academy of Forestry, Beijing, 100091, China
| | - Junhui Wang
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, Beijing, 100091, China.
| | - Zirui Jia
- State Key Laboratory of Tree Genetics and Breeding, Research Institute of Forestry, Chinese Academy of Forestry, Beijing, 100091, China.
| |
Collapse
|
3
|
Lin H, Li W, Zhao Y. A Diverging Species within the Stewartia gemmata (Theaceae) Complex Revealed by RAD-Seq Data. PLANTS (BASEL, SWITZERLAND) 2024; 13:1296. [PMID: 38794366 PMCID: PMC11124813 DOI: 10.3390/plants13101296] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/18/2024] [Revised: 05/04/2024] [Accepted: 05/07/2024] [Indexed: 05/26/2024]
Abstract
Informed species delimitation is crucial in diverse biological fields; however, it can be problematic for species complexes. Showing a peripatric distribution pattern, Stewartia gemmata and S. acutisepala (the S. gemmata complex) provide us with an opportunity to study species boundaries among taxa undergoing nascent speciation. Here, we generated genomic data from representative individuals across the natural distribution ranges of the S. gemmata complex using restriction site-associated DNA sequencing (RAD-seq). Based on the DNA sequence of assembled loci containing 41,436 single-nucleotide polymorphisms (SNPs) and invariant sites, the phylogenetic analysis suggested strong monophyly of both the S. gemmata complex and S. acutisepala, and the latter was nested within the former. Among S. gemmata individuals, the one sampled from Mt. Tianmu (Zhejiang) showed the closest evolutionary affinity with S. acutisepala (which is endemic to southern Zhejiang). Estimated from 2996 high-quality SNPs, the genetic divergence between S. gemmata and S. acutisepala was relatively low (an Fst of 0.073 on a per-site basis). Nevertheless, we observed a proportion of genomic regions showing relatively high genetic differentiation on a windowed basis. Up to 1037 genomic bins showed an Fst value greater than 0.25, accounting for 8.31% of the total. After SNPs subject to linkage disequilibrium were pruned, the principal component analysis (PCA) showed that S. acutisepala diverged from S. gemmata along the first and the second PCs to some extent. By applying phylogenomic analysis, the present study determines that S. acutisepala is a variety of S. gemmata and is diverging from S. gemmata, providing empirical insights into the nascent speciation within a species complex.
Collapse
Affiliation(s)
- Hanyang Lin
- Zhejiang Provincial Key Laboratory of Plant Evolutionary Ecology and Conservation, School of Life Sciences, Taizhou University, Taizhou 318000, China;
| | - Wenhao Li
- Systematic & Evolutionary Botany and Biodiversity Group, MOE Key Laboratory of Biosystems Homeostasis & Protection, College of Life Sciences, Zhejiang University, Hangzhou 310058, China;
| | - Yunpeng Zhao
- Systematic & Evolutionary Botany and Biodiversity Group, MOE Key Laboratory of Biosystems Homeostasis & Protection, College of Life Sciences, Zhejiang University, Hangzhou 310058, China;
| |
Collapse
|
4
|
Springer AL, Gompert Z. Considerable genetic diversity and structure despite narrow endemism and limited ecological specialization in the Hayden's ringlet, Coenonympha haydenii. Mol Ecol 2024; 33:e17310. [PMID: 38441401 DOI: 10.1111/mec.17310] [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/03/2023] [Revised: 11/26/2023] [Accepted: 02/15/2024] [Indexed: 03/26/2024]
Abstract
Understanding the processes that underlie the development of population genetic structure is central to the study of evolution. Patterns of genetic structure, in turn, can reveal signatures of isolation by distance (IBD), barriers to gene flow, or even the genesis of speciation. However, it is unclear how severe range restriction might impact the processes that dominate the development of genetic structure. In narrow endemic species, is population structure likely to be adaptive in nature, or rather the result of genetic drift? In this study, we investigated patterns of genetic diversity and structure in the narrow endemic Hayden's ringlet butterfly. Specifically, we asked to what degree genetic structure in the Hayden's ringlet can be explained by IBD, isolation by resistance (IBR) (in the form of geographic or ecological barriers to migration between populations), and isolation by environment (in the form of differences in host plant availability and preference). We employed a genotyping-by-sequencing (GBS) approach coupled with host preference assays, Bayesian modelling, and population genomic analyses to answer these questions. Our results suggest that despite their restricted range, levels of genetic diversity in the Hayden's ringlet are comparable to those seen in more widespread butterfly species. Hayden's ringlets showed a strong preference for feeding on grasses relative to sedges, but neither larval preference nor potential host availability at sampling sites correlated with genetic structure. We conclude that geography, in the form of IBR and simple IBD, was the major driver of contemporary patterns of differentiation in this narrow endemic species.
Collapse
Affiliation(s)
- Amy L Springer
- Department of Biology, Utah State University, Logan, Utah, USA
| | - Zachariah Gompert
- Department of Biology, Utah State University, Logan, Utah, USA
- Ecology Center, Utah State University, Logan, Utah, USA
| |
Collapse
|
5
|
Martins ARP, Warren NB, McMillan WO, Barrett RDH. Spatiotemporal dynamics in butterfly hybrid zones. INSECT SCIENCE 2024; 31:328-353. [PMID: 37596954 DOI: 10.1111/1744-7917.13262] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 05/02/2023] [Revised: 07/13/2023] [Accepted: 07/21/2023] [Indexed: 08/21/2023]
Abstract
Evaluating whether hybrid zones are stable or mobile can provide novel insights for evolution and conservation biology. Butterflies exhibit high sensitivity to environmental changes and represent an important model system for the study of hybrid zone origins and maintenance. Here, we review the literature exploring butterfly hybrid zones, with a special focus on their spatiotemporal dynamics and the potential mechanisms that could lead to their movement or stability. We then compare different lines of evidence used to investigate hybrid zone dynamics and discuss the strengths and weaknesses of each approach. Our goal with this review is to reveal general conditions associated with the stability or mobility of butterfly hybrid zones by synthesizing evidence obtained using different types of data sampled across multiple regions and spatial scales. Finally, we discuss spatiotemporal dynamics in the context of a speciation/divergence continuum, the relevance of hybrid zones for conservation biology, and recommend key topics for future investigation.
Collapse
Affiliation(s)
- Ananda R Pereira Martins
- Redpath Museum, McGill University, 859 Sherbrooke Street West, Montreal, Quebec, Canada
- Smithsonian Tropical Research Institute, Gamboa, Panama City, Panama
| | - Natalie B Warren
- Redpath Museum, McGill University, 859 Sherbrooke Street West, Montreal, Quebec, Canada
| | - W Owen McMillan
- Smithsonian Tropical Research Institute, Gamboa, Panama City, Panama
| | - Rowan D H Barrett
- Redpath Museum, McGill University, 859 Sherbrooke Street West, Montreal, Quebec, Canada
| |
Collapse
|
6
|
Sheinman M, Arndt PF, Massip F. Modeling the mosaic structure of bacterial genomes to infer their evolutionary history. Proc Natl Acad Sci U S A 2024; 121:e2313367121. [PMID: 38517978 PMCID: PMC10990148 DOI: 10.1073/pnas.2313367121] [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: 08/03/2023] [Accepted: 01/30/2024] [Indexed: 03/24/2024] Open
Abstract
The chronology and phylogeny of bacterial evolution are difficult to reconstruct due to a scarce fossil record. The analysis of bacterial genomes remains challenging because of large sequence divergence, the plasticity of bacterial genomes due to frequent gene loss, horizontal gene transfer, and differences in selective pressure from one locus to another. Therefore, taking advantage of the rich and rapidly accumulating genomic data requires accurate modeling of genome evolution. An important technical consideration is that loci with high effective mutation rates may diverge beyond the detection limit of the alignment algorithms used, biasing the genome-wide divergence estimates toward smaller divergences. In this article, we propose a novel method to gain insight into bacterial evolution based on statistical properties of genome comparisons. We find that the length distribution of sequence matches is shaped by the effective mutation rates of different loci, by the horizontal transfers, and by the aligner sensitivity. Based on these inputs, we build a model and show that it accounts for the empirically observed distributions, taking the Enterobacteriaceae family as an example. Our method allows to distinguish segments of vertical and horizontal origins and to estimate the time divergence and exchange rate between any pair of taxa from genome-wide alignments. Based on the estimated time divergences, we construct a time-calibrated phylogenetic tree to demonstrate the accuracy of the method.
Collapse
Affiliation(s)
- Michael Sheinman
- Institute for Advanced Studies, Sevastopol State University, Sevastopol299053, Crimea
| | - Peter F. Arndt
- Department of Computational Molecular Biology, Max Planck Institute for Molecular Genetics, Berlin12163, Germany
| | - Florian Massip
- Department U900, Centre for Computational Biology, Mines Paris, PSL University, Paris75006, France
- Department U900, Institut Curie, Université Paris Sciences et Lettres, Paris75005, France
- INSERM, U900, Paris75005, France
| |
Collapse
|
7
|
Groh JS, Coop G. The temporal and genomic scale of selection following hybridization. Proc Natl Acad Sci U S A 2024; 121:e2309168121. [PMID: 38489387 PMCID: PMC10962946 DOI: 10.1073/pnas.2309168121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/31/2023] [Accepted: 01/30/2024] [Indexed: 03/17/2024] Open
Abstract
Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the evolutionary dynamics within hybrid populations that underlie these patterns have been lacking. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of ancestry variation at varying spatial genomic scales through time. Here, we develop methods based on the Discrete Wavelet Transform to study the genomic scale of local ancestry variation and its association with recombination rates and show that these methods capture temporal dynamics of drift and genome-wide selection after hybridization. We apply these methods to published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio) and to inferred Neanderthal introgression in modern humans. Across systems, upward of 20% of variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. Signatures of selection at fine genomic scales suggest selection over longer time scales; however, we suggest that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from contiguous segments of genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available and can help shed light on generalities of the genomic consequences of interspecific hybridization.
Collapse
Affiliation(s)
- Jeffrey S. Groh
- Department of Evolution and Ecology and Center for Population Biology, University of California, Davis, CA95616
| | - Graham Coop
- Department of Evolution and Ecology and Center for Population Biology, University of California, Davis, CA95616
| |
Collapse
|
8
|
Innes PA, Goebl AM, Smith CCR, Rosenberger K, Kane NC. Gene expression and alternative splicing contribute to adaptive divergence of ecotypes. Heredity (Edinb) 2024; 132:120-132. [PMID: 38071268 PMCID: PMC10924094 DOI: 10.1038/s41437-023-00665-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/28/2023] [Revised: 11/23/2023] [Accepted: 11/24/2023] [Indexed: 03/10/2024] Open
Abstract
Regulation of gene expression is a critical link between genotype and phenotype explaining substantial heritable variation within species. However, we are only beginning to understand the ways that specific gene regulatory mechanisms contribute to adaptive divergence of populations. In plants, the post-transcriptional regulatory mechanism of alternative splicing (AS) plays an important role in both development and abiotic stress response, making it a compelling potential target of natural selection. AS allows organisms to generate multiple different transcripts/proteins from a single gene and thus may provide a source of evolutionary novelty. Here, we examine whether variation in alternative splicing and gene expression levels might contribute to adaptation and incipient speciation of dune-adapted prairie sunflowers in Great Sand Dunes National Park, Colorado, USA. We conducted a common garden experiment to assess transcriptomic variation among ecotypes and analyzed differential expression, differential splicing, and gene coexpression. We show that individual genes are strongly differentiated for both transcript level and alternative isoform proportions, even when grown in a common environment, and that gene coexpression networks are disrupted between ecotypes. Furthermore, we examined how genome-wide patterns of sequence divergence correspond to divergence in transcript levels and isoform proportions and find evidence for both cis and trans-regulation. Together, our results emphasize that alternative splicing has been an underappreciated mechanism providing source material for natural selection at short evolutionary time scales.
Collapse
Affiliation(s)
- Peter A Innes
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA.
| | - April M Goebl
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
- Research and Conservation Department, Denver Botanic Gardens, Denver, CO, USA
| | - Chris C R Smith
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
- Institute of Ecology and Evolution, University of Oregon, Eugene, OR, USA
| | - Kaylee Rosenberger
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
| | - Nolan C Kane
- Ecology and Evolutionary Biology Department, University of Colorado, Boulder, CO, USA
| |
Collapse
|
9
|
Johnson JA, Novak B, Athrey G, Sharo AG, Chase T, Toepfer J. Phylogenomics of the extinct Heath Hen provides support for sex-biased introgression among extant prairie grouse. Mol Phylogenet Evol 2023; 189:107927. [PMID: 37714443 DOI: 10.1016/j.ympev.2023.107927] [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: 01/11/2023] [Revised: 06/20/2023] [Accepted: 09/12/2023] [Indexed: 09/17/2023]
Abstract
Rapid divergence and subsequent reoccurring patterns of gene flow can complicate our ability to discern phylogenetic relationships among closely related species. To what degree such patterns may differ across the genome can provide an opportunity to extrapolate better how life history constraints may influence species boundaries. By exploring differences between autosomal and Z (or X) chromosomal-derived phylogenetic patterns, we can better identify factors that may limit introgression despite patterns of incomplete lineage sorting among closely related taxa. Here, using a whole-genome resequencing approach coupled with an exhaustive sampling of subspecies within the recently divergent prairie grouse complex (genus: Tympanuchus), including the extinct Heath Hen (T. cupido cupido), we show that their phylogenomic history differs depending on autosomal or Z-chromosome partitioned SNPs. Because the Heath Hen was allopatric relative to the other prairie grouse taxa, its phylogenetic signature should not be influenced by gene flow. In contrast, all the other extant prairie grouse taxa, except Attwater's Prairie-chicken (T. c. attwateri), possess overlapping contemporary geographic distributions and have been known to hybridize. After excluding samples that were likely translocated prairie grouse from the Midwest to the eastern coastal states or their resulting hybrids with mainland Heath Hens, species tree analyses based on autosomal SNPs consistently identified a paraphyletic relationship with regard to the Heath Hen with Lesser Prairie-chicken (T. pallidicinctus) sister to Greater Prairie-chicken (T. c. pinnatus) regardless of genic or intergenic partitions. In contrast, species trees based on the Z-chromosome were consistent with Heath Hen sister to a clade that included its conspecifics, Greater and Attwater's Prairie-chickens (T. c. attwateri). These results were further explained by historic gene flow, as shown with an excess of autosomal SNPs shared between Lesser and Greater Prairie-chickens but not with the Z-chromosome. Phylogenetic placement of Sharp-tailed Grouse (T. phasianellus), however, did not differ among analyses and was sister to a clade that included all other prairie grouse despite low levels of autosomal gene flow with Greater Prairie-chicken. These results, along with strong sexual selection (i.e., male hybrid behavioral isolation) and a lek breeding system (i.e., high variance in male mating success), are consistent with a pattern of female-biased introgression between prairie grouse taxa with overlapping geographic distributions. Additional study is warranted to explore how genomic components associated with the Z-chromosome influence the phenotype and thereby impact species limits among prairie grouse taxa despite ongoing contemporary gene flow.
Collapse
Affiliation(s)
- Jeff A Johnson
- The Peregrine Fund, Boise, ID 83709, USA; Department of Biological Sciences, University of North Texas, Denton, TX 76203, USA.
| | - Ben Novak
- Revive & Restore, Sausalito, CA 94965, USA
| | - Giridhar Athrey
- Department of Poultry Science & Faculty of Ecology and Evolutionary Biology, Texas A&M University, College Station, TX 77843, USA
| | - Andrew G Sharo
- Department of Ecology and Evolutionary Biology, University of California, Santa Cruz, Santa Cruz, CA 95064, USA
| | - Tom Chase
- Village and Wilderness, Martha's Vineyard, MA 02557, USA
| | - John Toepfer
- George Miksch Sutton Avian Research Center, Bartlesville, OK 74005, USA
| |
Collapse
|
10
|
Kayhani K, Barreto FS. Disproportionate role of nuclear-encoded proteins in organismal and mitochondrial thermal performance in a copepod. J Exp Biol 2023; 226:jeb246085. [PMID: 37947077 DOI: 10.1242/jeb.246085] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/08/2023] [Accepted: 10/30/2023] [Indexed: 11/12/2023]
Abstract
Determining the mechanisms by which organisms evolve thermal tolerance is crucial to predicting how populations may respond to changes in local temperature regimes. Although evidence of relationships between mitochondrial background and thermal adaptation have been found, the presence of both nuclear-encoded and mitochondrial DNA (mtDNA)-encoded proteins warrants experiments aimed at parsing out the relative role of each genome in thermal adaptation. We investigated the relative role of mtDNA-encoded products in thermal tolerance between two divergent populations of Tigriopus californicus using first-generation (F1) hybrids that vary in maternally inherited mtDNA but are heterozygous for population-specific alleles across nuclear loci. We tested two measures of thermal tolerance, (1) survivorship to acute thermal stress and (2) thermal stability of mitochondrial performance in Complex I-fueled ATP synthesis, both across a range of increasing temperatures. We found that the southern population (San Diego, CA, USA) outperformed the northern population (Strawberry Hill, OR, USA) in survivorship, and that both reciprocal F1 hybrid crosses had intermediate survival. Mitochondria from the San Diego population displayed greater stability in ATP synthesis with increasing temperatures compared with those from Strawberry Hill. Interestingly, hybrids from both cross directions had synthesis profiles that were very similar to that of Strawberry Hill. Taken together, these results suggest that the relative role of the mtDNA in these phenotypes is negligible compared with that of elements encoded by nuclear DNA in this system.
Collapse
Affiliation(s)
- Kamron Kayhani
- Department of Integrative Biology, Oregon State University, Corvallis, OR 97331, USA
| | - Felipe S Barreto
- Department of Integrative Biology, Oregon State University, Corvallis, OR 97331, USA
| |
Collapse
|
11
|
Deville D, Kawai K, Fujita H, Umino T. Genetic divergences and hybridization within the Sebastes inermis complex. PeerJ 2023; 11:e16391. [PMID: 38025733 PMCID: PMC10656903 DOI: 10.7717/peerj.16391] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2023] [Accepted: 10/11/2023] [Indexed: 12/01/2023] Open
Abstract
The Sebastes inermis complex includes three sympatric species (Sebastes cheni, viz Sebastes inermis, and Sebastes ventricosus) with clear ecomorphological differences, albeit incomplete reproductive isolation. The presence of putative morphological hybrids (PMH) with plausibly higher fitness than the parent species indicates the need to confirm whether hybridization occurs within the complex. In this sense, we assessed the dynamics of genetic divergence and hybridization within the species complex using a panel of 10 microsatellite loci, and sequences of the mitochondrial control region (D-loop) and the intron-free rhodopsin (RH1) gene. The analyses revealed the presence of three distinct genetic clusters, large genetic distances using D-loop sequences, and distinctive mutations within the RH1 gene. These results are consistent with the descriptions of the three species. Two microsatellite loci had signatures of divergent selection, indicating that they are linked to genomic regions that are crucial for speciation. Furthermore, nonsynonymous mutations within the RH1 gene detected in S. cheni and "Kumano" (a PMH) suggest dissimilar adaptations related to visual perception in dim-light environments. The presence of individuals with admixed ancestry between two species confirmed hybridization. The presence of nonsynonymous mutations within the RH1 gene and the admixed ancestry of the "Kumano" morphotype highlight the potential role of hybridization in generating novelties within the species complex. We discuss possible outcomes of hybridization within the species complex, considering hybrid fitness and assortative mating. Overall, our findings indicate that the genetic divergence of each species is maintained in the presence of hybridization, as expected in a scenario of speciation-with-gene-flow.
Collapse
Affiliation(s)
- Diego Deville
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashihiroshima, Hiroshima, Japón
| | - Kentaro Kawai
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashihiroshima, Hiroshima, Japón
| | - Hiroki Fujita
- Seto Marine Biological Laboratory, Field Science Education and Research Center, Kyoto University, Shirahama, Wakayama, Japan
| | - Tetsuya Umino
- Graduate School of Integrated Sciences for Life, Hiroshima University, Higashihiroshima, Hiroshima, Japón
| |
Collapse
|
12
|
Kajtoch Ł, Kolasa M, Mazur MA, Ścibior R, Zając K, Kubisz D. Limited congruence in phylogeographic patterns observed for riverine predacious beetles sharing distribution along the mountain rivers. Sci Rep 2023; 13:17883. [PMID: 37857828 PMCID: PMC10587157 DOI: 10.1038/s41598-023-44922-w] [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: 06/10/2023] [Accepted: 10/13/2023] [Indexed: 10/21/2023] Open
Abstract
Riverine predacious beetles (RPB) (Carabidae, Staphylinidae) are highly diverse and numerous elements of riverine ecosystems. Their historical and contemporary distribution and diversity are highly dependent on natural flow regimes and topography of watercourses. Despite broad knowledge of their ecology, data on population genetic diversity and connectivity are lacking. This study aimed to fill this gap in order to solve two principal hypotheses assuming (i) congruence of phylogeographic patterns observed for RPB indicating that they share a common history and the ecological adaptations to the dynamic environment, (ii) genetic structuration of populations according to river basins. The Carpathian populations of four ground beetles and three rove beetles were examined using cytochrome oxidase and arginine kinase sequencing. There are substantial differences in RPB demographic history and current genetic diversity. Star-like phylogeny of Bembidion and complex haplotype networks of Paederus/Paederidus, with some haplotypes being drainage-specific and others found in distant populations, indicate a general lack of isolation by distance. Signs of recent demographic expansion were detected for most RPB with the latest population collapse for some rove beetles. To some extent, migration of examined species has to be limited by watersheds. Observed phylogeographic patterns are essential for correctly understanding RPB meta-population functioning.
Collapse
Affiliation(s)
- Łukasz Kajtoch
- Institute of Systematics and Evolution of Animals, Polish Academy of Sciences, Sławkowska 17, 31-016, Kraków, Poland
| | - Michał Kolasa
- Institute of Systematics and Evolution of Animals, Polish Academy of Sciences, Sławkowska 17, 31-016, Kraków, Poland
- Institute of Environmental Sciences, Faculty of Biology, Jagiellonian University, Gronostajowa 7, 30-387, Kraków, Poland
| | - Miłosz A Mazur
- Institute of Biology, University of Opole, Oleska 22, 45-050, Opole, Poland
| | - Radosław Ścibior
- Department of Zoology and Animal Ecology, University of Life Sciences in Lublin, Akademicka 13, 20-950, Lublin, Poland.
| | | | - Daniel Kubisz
- Institute of Systematics and Evolution of Animals, Polish Academy of Sciences, Sławkowska 17, 31-016, Kraków, Poland
| |
Collapse
|
13
|
Bock DG, Cai Z, Elphinstone C, González-Segovia E, Hirabayashi K, Huang K, Keais GL, Kim A, Owens GL, Rieseberg LH. Genomics of plant speciation. PLANT COMMUNICATIONS 2023; 4:100599. [PMID: 37050879 PMCID: PMC10504567 DOI: 10.1016/j.xplc.2023.100599] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 01/26/2023] [Revised: 03/21/2023] [Accepted: 04/06/2023] [Indexed: 06/19/2023]
Abstract
Studies of plants have been instrumental for revealing how new species originate. For several decades, botanical research has complemented and, in some cases, challenged concepts on speciation developed via the study of other organisms while also revealing additional ways in which species can form. Now, the ability to sequence genomes at an unprecedented pace and scale has allowed biologists to settle decades-long debates and tackle other emerging challenges in speciation research. Here, we review these recent genome-enabled developments in plant speciation. We discuss complications related to identification of reproductive isolation (RI) loci using analyses of the landscape of genomic divergence and highlight the important role that structural variants have in speciation, as increasingly revealed by new sequencing technologies. Further, we review how genomics has advanced what we know of some routes to new species formation, like hybridization or whole-genome duplication, while casting doubt on others, like population bottlenecks and genetic drift. While genomics can fast-track identification of genes and mutations that confer RI, we emphasize that follow-up molecular and field experiments remain critical. Nonetheless, genomics has clarified the outsized role of ancient variants rather than new mutations, particularly early during speciation. We conclude by highlighting promising avenues of future study. These include expanding what we know so far about the role of epigenetic and structural changes during speciation, broadening the scope and taxonomic breadth of plant speciation genomics studies, and synthesizing information from extensive genomic data that have already been generated by the plant speciation community.
Collapse
Affiliation(s)
- Dan G Bock
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Zhe Cai
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Cassandra Elphinstone
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Eric González-Segovia
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | | | - Kaichi Huang
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Graeme L Keais
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Amy Kim
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada
| | - Gregory L Owens
- Department of Biology, University of Victoria, Victoria, BC, Canada
| | - Loren H Rieseberg
- Department of Botany and Biodiversity Research Centre, University of British Columbia, Vancouver, BC, Canada.
| |
Collapse
|
14
|
Wessinger CA, Katzer AM, Hime PM, Rausher MD, Kelly JK, Hileman LC. A few essential genetic loci distinguish Penstemon species with flowers adapted to pollination by bees or hummingbirds. PLoS Biol 2023; 21:e3002294. [PMID: 37769035 PMCID: PMC10538765 DOI: 10.1371/journal.pbio.3002294] [Citation(s) in RCA: 5] [Impact Index Per Article: 5.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/10/2023] [Accepted: 08/11/2023] [Indexed: 09/30/2023] Open
Abstract
In the formation of species, adaptation by natural selection generates distinct combinations of traits that function well together. The maintenance of adaptive trait combinations in the face of gene flow depends on the strength and nature of selection acting on the underlying genetic loci. Floral pollination syndromes exemplify the evolution of trait combinations adaptive for particular pollinators. The North American wildflower genus Penstemon displays remarkable floral syndrome convergence, with at least 20 separate lineages that have evolved from ancestral bee pollination syndrome (wide blue-purple flowers that present a landing platform for bees and small amounts of nectar) to hummingbird pollination syndrome (bright red narrowly tubular flowers offering copious nectar). Related taxa that differ in floral syndrome offer an attractive opportunity to examine the genomic basis of complex trait divergence. In this study, we characterized genomic divergence among 229 individuals from a Penstemon species complex that includes both bee and hummingbird floral syndromes. Field plants are easily classified into species based on phenotypic differences and hybrids displaying intermediate floral syndromes are rare. Despite unambiguous phenotypic differences, genome-wide differentiation between species is minimal. Hummingbird-adapted populations are more genetically similar to nearby bee-adapted populations than to geographically distant hummingbird-adapted populations, in terms of genome-wide dXY. However, a small number of genetic loci are strongly differentiated between species. These approximately 20 "species-diagnostic loci," which appear to have nearly fixed differences between pollination syndromes, are sprinkled throughout the genome in high recombination regions. Several map closely to previously established floral trait quantitative trait loci (QTLs). The striking difference between the diagnostic loci and the genome as whole suggests strong selection to maintain distinct combinations of traits, but with sufficient gene flow to homogenize the genomic background. A surprisingly small number of alleles confer phenotypic differences that form the basis of species identity in this species complex.
Collapse
Affiliation(s)
- Carolyn A. Wessinger
- Department of Biological Sciences, University of South Carolina, Columbia, South Carolina, United States of America
| | - Amanda M. Katzer
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas, United States of America
| | - Paul M. Hime
- Biodiversity Institute and Natural History Museum, University of Kansas, Lawrence, Kansas, United States of America
| | - Mark D. Rausher
- Department of Biology, Duke University, Durham, North Carolina, United States of America
| | - John K. Kelly
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas, United States of America
| | - Lena C. Hileman
- Department of Ecology and Evolutionary Biology, University of Kansas, Lawrence, Kansas, United States of America
| |
Collapse
|
15
|
Chen X, Avia K, Forler A, Remoué C, Venon A, Rousselet A, Lucas G, Kwarteng AO, Rover R, Le Guilloux M, Belcram H, Combes V, Corti H, Olverà-Vazquez S, Falque M, Alins G, Kirisits T, Ursu TM, Roman A, Volk GM, Bazot S, Cornille A. Ecological and evolutionary drivers of phenotypic and genetic variation in the European crabapple [Malus sylvestris (L.) Mill.], a wild relative of the cultivated apple. ANNALS OF BOTANY 2023; 131:1025-1037. [PMID: 37148364 PMCID: PMC10332392 DOI: 10.1093/aob/mcad061] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/30/2022] [Accepted: 05/03/2023] [Indexed: 05/08/2023]
Abstract
BACKGROUND AND AIMS Studying the relationship between phenotypic and genetic variation in populations distributed across environmental gradients can help us to understand the ecological and evolutionary processes involved in population divergence. We investigated the patterns of genetic and phenotypic diversity in the European crabapple, Malus sylvestris, a wild relative of the cultivated apple (Malus domestica) that occurs naturally across Europe in areas subjected to different climatic conditions, to test for divergence among populations. METHODS Growth rates and traits related to carbon uptake in seedlings collected across Europe were measured in controlled conditions and associated with the genetic status of the seedlings, which was assessed using 13 microsatellite loci and the Bayesian clustering method. Isolation-by-distance, isolation-by-climate and isolation-by-adaptation patterns, which can explain genetic and phenotypic differentiation among M. sylvestris populations, were also tested. KEY RESULTS A total of 11.6 % of seedlings were introgressed by M. domestica, indicating that crop-wild gene flow is ongoing in Europe. The remaining seedlings (88.4 %) belonged to seven M. sylvestris populations. Significant phenotypic trait variation among M. sylvestris populations was observed. We did not observe significant isolation by adaptation; however, the significant association between genetic variation and the climate during the Last Glacial Maximum suggests that there has been local adaptation of M. sylvestris to past climates. CONCLUSIONS This study provides insight into the phenotypic and genetic differentiation among populations of a wild relative of the cultivated apple. This might help us to make better use of its diversity and provide options for mitigating the impact of climate change on the cultivated apple through breeding.
Collapse
Affiliation(s)
- X Chen
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - K Avia
- Université de Strasbourg, INRAE, SVQV UMR-A 1131, F-68000 Colmar, France
| | - A Forler
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - C Remoué
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - A Venon
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - A Rousselet
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - G Lucas
- Université Paris-Saclay, CEA, CNRS, Institute for Integrative Biology of the Cell (I2BC), Gif-sur-Yvette 91198, France
| | - A O Kwarteng
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - R Rover
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - M Le Guilloux
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - H Belcram
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - V Combes
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - H Corti
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - S Olverà-Vazquez
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - M Falque
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| | - G Alins
- Institut de Recerca i Tecnologia Agroalimentàries, IRTA-Fruit Production, PCiTAL, Parc 21 de Gardeny, edifici Fruitcentre, 25003 Lleida, Spain
| | - T Kirisits
- Institute of Forest Entomology, Forest Pathology and Forest Protection (IFFF), Department of Forest and Soil Sciences, University of Natural Resources and Life Sciences, Vienna (BOKU), Peter-Jordan-Straße 82 (Franz Schwackhöfer-Haus), A-1190 Vienna, Austria
| | - T M Ursu
- NIRDBS, Institute of Biological Research Cluj-Napoca, 48 Republicii St., Cluj-Napoca, Romania
| | - A Roman
- NIRDBS, Institute of Biological Research Cluj-Napoca, 48 Republicii St., Cluj-Napoca, Romania
| | - G M Volk
- USDA-ARS National Laboratory for Genetic Resources Preservation, 1111 South Mason Street, Fort Collins, CO 80521, USA
| | - S Bazot
- Ecologie Systématique et Evolution, CNRS, AgroParisTech, Ecologie Systématique Evolution, Université Paris‐Saclay, Orsay, France
| | - A Cornille
- Université Paris Saclay, INRAE, CNRS, AgroParisTech, GQE – Le Moulon, 91190 Gif-sur-Yvette, France
| |
Collapse
|
16
|
Crawford NG, McGreevy TJ, Mullen SP, Schneider CJ. The genetic basis of conspicuous coloration in the Guadeloupean anole: Evolution by sexual and ecological selection. Ecol Evol 2023; 13:e10266. [PMID: 37435022 PMCID: PMC10330958 DOI: 10.1002/ece3.10266] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/11/2023] [Revised: 06/05/2023] [Accepted: 06/06/2023] [Indexed: 07/13/2023] Open
Abstract
Understanding how natural selection acts on the genome and contributes to the process of speciation is a primary aim of the study of evolution. Here we used natural variation in two subspecies of the Guadeloupean anole (Anolis marmoratus ssp.), from the island of Guadeloupe in the Lesser Antilles, to explore the genomic basis of adaptation and speciation in Anolis lizards. These subspecies inhabit distinct ecological environments and display marked differences in adult male color and pattern. We sequenced the complete genomes of 20 anoles, 10 from each subspecies, at 1.4× coverage. We used genome-wide scans of population differentiation, allele frequency spectrum, and linkage disequilibrium to characterize the genomic architecture within and between the subspecies. While most of the genome was undifferentiated, we observed five large divergent regions. Within these regions we identified blocks, 5 kb pairs in length, enriched for fixed single nucleotide polymorphisms. These blocks encompass 97 genes, two of which are candidate pigmentation genes. One is melanophilin (mlph), which helps transport melanosomes within melanocytes. The other is a cluster of differentiation 36 (cd36), which regulates carotenoid pigment sequestration. We used high-pressure liquid chromatography to confirm that carotenoid pigments are significantly more abundant in the conspicuous orange-pigmented skin of male A. m. marmoratus suggesting that cd36 may be regulating pigment deposition in this tissue. We identified for the first time a carotenoid gene that is a potential target of divergent sexual selection and may be contributing to the early stages of speciation in Anolis lizards.
Collapse
Affiliation(s)
| | - Thomas J. McGreevy
- Department of BiologyBoston UniversityBostonMassachusettsUSA
- Department of Natural Resources ScienceUniversity of Rhode IslandKingstonRhode IslandUSA
| | - Sean P. Mullen
- Department of BiologyBoston UniversityBostonMassachusettsUSA
| | | |
Collapse
|
17
|
Turbek SP, Funk WC, Ruegg KC. Where to draw the line? Expanding the delineation of conservation units to highly mobile taxa. J Hered 2023; 114:300-311. [PMID: 36815497 DOI: 10.1093/jhered/esad011] [Citation(s) in RCA: 8] [Impact Index Per Article: 8.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/27/2022] [Accepted: 02/21/2023] [Indexed: 02/24/2023] Open
Abstract
Conservation units (CUs) are an essential tool for maximizing evolutionary potential and prioritizing areas across a species' range for protection when implementing conservation and management measures. However, current workflows for identifying CUs on the basis of neutral and adaptive genomic variation largely ignore information contained in patterns of isolation by distance (IBD), frequently the primary signal of population structure in highly mobile taxa, such as birds, bats, and marine organisms with pelagic larval stages. While individuals located on either end of a species' distribution may exhibit clear genetic, phenotypic, and ecological differences, IBD produces subtle changes in allele frequencies across space, making it difficult to draw clear boundaries for conservation purposes in the absence of discrete population structure. Here, we highlight potential pitfalls that arise when applying common methods for delineating CUs to continuously distributed organisms and review existing methods for detecting subtle breakpoints in patterns of IBD that can indicate barriers to gene flow in highly mobile taxa. In addition, we propose a new framework for identifying CUs in all organisms, including those characterized by continuous genomic differentiation, and suggest several possible ways to harness the information contained in patterns of IBD to guide conservation and management decisions.
Collapse
Affiliation(s)
- Sheela P Turbek
- Department of Biology, Colorado State University, Fort Collins, CO, United States
| | - W Chris Funk
- Department of Biology, Colorado State University, Fort Collins, CO, United States
| | - Kristen C Ruegg
- Department of Biology, Colorado State University, Fort Collins, CO, United States
| |
Collapse
|
18
|
Mittermayer F, Helmerson C, Duvetorp M, Johannesson K, Panova M. The molecular background of the aspartate aminotransferase polymorphism in Littorina snails maintained by strong selection on small spatial scales. Gene 2023:147517. [PMID: 37257792 DOI: 10.1016/j.gene.2023.147517] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/17/2023] [Revised: 05/14/2023] [Accepted: 05/24/2023] [Indexed: 06/02/2023]
Abstract
Allozymes present several classical examples of divergent selection, including the variation in the cytosolic aspartate aminotransferase (AAT) in the intertidal snails Littorina saxatilis. AAT is a part of the asparate-malate shuttle, in the interidal molluscs involved in the anaerobic respiration during desiccation. Previous allozyme studies reported the sharp gradient in the frequencies of the AAT100and the AAT120 alleles between the low and high shores in the Northern Europe and the differences in their enzymatic activity, supporting the role of AAT in adaptation to desiccation. However, the populations in the Iberian Peninsula showed the opposite allele cline. Using the mRNA sequencing and the genome pool-seq analyses we characterize DNA sequences of the different AAT alleles, report the amino acid replacements behind the allozyme variation and show that same allozyme alleles in Northern and Southern populations have different protein sequences. Gene phylogeny reveals that the AAT100 and the northern AAT120 alleles represent the old polymorphism, shared among the closely related species of Littorina, while the southern AAT120 allele is more recently derived from AAT100. Further, we show that the Aat gene is expressed at constitutive level in different genotypes and conditions, supporting the role of structural variation in regulation of enzyme activity. Finally, we report the location and the structure of the gene in the L. saxatilis genome and the presence of two additional non-functional gene copies. Altogether, we provide a missing link between the classical allozyme studies and the genome scans and bring together the results produced over decades of the genetic research.
Collapse
Affiliation(s)
- Felix Mittermayer
- Research Division Marine Ecology, Research Unit Marine Evolutionary Ecology, GEOMAR Helmholtz Centre for Ocean Research Kiel, Düsternbrooker Weg 20, 24105 Kiel, Germany; Department of Marine Sciences, University of Gothenburg, Tjärnö, 45296 Strömstad, Sweden
| | - Cecilia Helmerson
- Centre for Ecological and Evolutionary Synthesis, Department of Biosciences, University of Oslo; Department of Marine Sciences, University of Gothenburg, Tjärnö, 45296 Strömstad, Sweden
| | - Mårten Duvetorp
- Department of Marine Sciences, University of Gothenburg, Tjärnö, 45296 Strömstad, Sweden
| | - Kerstin Johannesson
- Department of Marine Sciences, University of Gothenburg, Tjärnö, 45296 Strömstad, Sweden
| | - Marina Panova
- Department of Marine Sciences, University of Gothenburg, Tjärnö, 45296 Strömstad, Sweden.
| |
Collapse
|
19
|
Groh J, Coop G. The temporal and genomic scale of selection following hybridization. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.05.25.542345. [PMID: 37337589 PMCID: PMC10276902 DOI: 10.1101/2023.05.25.542345] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/21/2023]
Abstract
Genomic evidence supports an important role for selection in shaping patterns of introgression along the genome, but frameworks for understanding the dynamics underlying these patterns within hybrid populations have been lacking. Here, we develop methods based on the Wavelet Transform to understand the spatial genomic scale of local ancestry variation and its association with recombination rates. We present theory and use simulations to show how wavelet-based decompositions of ancestry variance along the genome and the correlation between ancestry and recombination reflect the joint effects of recombination, genetic drift, and genome-wide selection against introgressed alleles. Due to the clock-like effect of recombination in hybrids breaking up parental haplotypes, drift and selection produce predictable patterns of local ancestry variation at varying spatial genomic scales through time. Using wavelet approaches to identify the genomic scale of variance in ancestry and its correlates, we show that these methods can detect temporally localized effects of drift and selection. We apply these methods to previously published datasets from hybrid populations of swordtail fish (Xiphophorus) and baboons (Papio), and to inferred Neanderthal introgression in modern humans. Across systems, we find that upwards of 20% of the variation in local ancestry at the broadest genomic scales can be attributed to systematic selection against introgressed alleles, consistent with strong selection acting on early-generation hybrids. We also see signals of selection at fine genomic scales and much longer time scales. However, we show that our ability to confidently infer selection at fine scales is likely limited by inherent biases in current methods for estimating local ancestry from genomic similarity. Wavelet approaches will become widely applicable as genomic data from systems with introgression become increasingly available, and can help shed light on generalities of the genomic consequences of interspecific hybridization.
Collapse
Affiliation(s)
- Jeffrey Groh
- Department of Evolution and Ecology, and Center for Population Biology, University of California, Davis, CA 95616
| | - Graham Coop
- Department of Evolution and Ecology, and Center for Population Biology, University of California, Davis, CA 95616
| |
Collapse
|
20
|
Wu Y, Hipp AL, Fargo G, Stith N, Ricklefs RE. Improving species delimitation for effective conservation: a case study in the endemic maple-leaf oak (Quercus acerifolia). THE NEW PHYTOLOGIST 2023; 238:1278-1293. [PMID: 36707920 DOI: 10.1111/nph.18777] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 08/24/2022] [Accepted: 01/22/2023] [Indexed: 06/18/2023]
Abstract
Species delimitation is challenging in lineages that exhibit both high plasticity and introgression. This challenge can be compounded by collection biases, which may downweight specimens morphologically intermediate between traditional species. Additionally, mismatch between named species and observable phenotypes can compromise species conservation. We studied the species boundaries of Quercus acerifolia, a tree endemic to Arkansas, U.S. We performed morphometric analyses of leaves and acorns from 527 field and 138 herbarium samples of Q. acerifolia and its close relatives, Q. shumardii and Q. rubra. We employed two novel approaches: sampling ex situ collections to detect phenotypic plasticity caused by environmental variation and comparing random field samples with historical herbarium samples to identify collection biases that might undermine species delimitation. To provide genetic evidence, we also performed molecular analyses on genome-wide SNPs. Quercus acerifolia shows distinctive morphological, ecological, and genomic characteristics, rejecting the hypothesis that Q. acerifolia is a phenotypic variant of Q. shumardii. We found mismatches between traditional taxonomy and phenotypic clusters. We detected underrepresentation of morphological intermediates in herbarium collections, which may bias species discovery and recognition. Rare species conservation requires considering and addressing taxonomic problems related to phenotypic plasticity, mismatch between taxonomy and morphological clusters, and collection biases.
Collapse
Affiliation(s)
- Yingtong Wu
- Biology Department, University of Missouri - St Louis, St Louis, MO, 63121, USA
- Whitney R. Harris World Ecology Center, University of Missouri - St Louis, St Louis, MO, 63121, USA
| | - Andrew L Hipp
- The Morton Arboretum, 4100 Illinois Route 53, Lisle, IL, 60532, USA
- The Field Museum, Integrative Research Center, 1400S Lake Shore Dr., Chicago, IL, 60605, USA
| | - Gregory Fargo
- Biology Department, University of Missouri - St Louis, St Louis, MO, 63121, USA
| | - Nora Stith
- Biology Department, University of Missouri - St Louis, St Louis, MO, 63121, USA
| | - Robert E Ricklefs
- Biology Department, University of Missouri - St Louis, St Louis, MO, 63121, USA
| |
Collapse
|
21
|
Evans A, de Kort H, Brys R, Duffy KJ, Jersáková J, Kull T, Selosse MA, Tsiftsis S, Minasiewicz J, Jacquemyn H. Historical biogeography and local adaptation explain population genetic structure in a widespread terrestrial orchid. ANNALS OF BOTANY 2023; 131:623-634. [PMID: 36680796 PMCID: PMC10147325 DOI: 10.1093/aob/mcad010] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/02/2022] [Accepted: 01/11/2023] [Indexed: 05/20/2023]
Abstract
BACKGROUND AND AIMS Historical changes in environmental conditions and colonization-extinction dynamics have a direct impact on the genetic structure of plant populations. However, understanding how past environmental conditions influenced the evolution of species with high gene flow is challenging when signals for genetic isolation and adaptation are swamped by gene flow. We investigated the spatial distribution and genetic structure of the widespread terrestrial orchid Epipactis helleborine to identify glacial refugia, characterize postglacial population dynamics and assess its adaptive potential. METHODS Ecological niche modelling was used to locate possible glacial refugia and postglacial recolonization opportunities of E. helleborine. A large single-nucleotide polymorphism (SNP) dataset obtained through genotyping by sequencing was used to define population genetic diversity and structure and to identify sources of postglacial gene flow. Outlier analyses were used to elucidate how adaptation to the local environment contributed to population divergence. KEY RESULTS The distribution of climatically suitable areas was restricted during the Last Glacial Maximum to the Mediterranean, south-western Europe and small areas in the Alps and Carpathians. Within-population genetic diversity was high in E. helleborine (mean expected heterozygosity, 0.373 ± 0.006; observed heterozygosity, 0.571 ± 0.012; allelic richness, 1.387 ± 0.007). Italy and central Europe are likely to have acted as important genetic sources during postglacial recolonization. Adaptive SNPs were associated with temperature, elevation and precipitation. CONCLUSIONS Forests in the Mediterranean and Carpathians are likely to have acted as glacial refugia for Epipactis helleborine. Postglacial migration northwards and to higher elevations resulted in the dispersal and diversification of E. helleborine in central Europe and Italy, and to geographical isolation and divergent adaptation in Greek and Italian populations. Distinguishing adaptive from neutral genetic diversity allowed us to conclude that E. helleborine has a high adaptive potential to climate change and demonstrates that signals of adaptation and historical isolation can be identified even in species with high gene flow.
Collapse
Affiliation(s)
- Alexandra Evans
- Department of Biology, Plant Conservation and Population Biology, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Hanne de Kort
- Department of Biology, Plant Conservation and Population Biology, Katholieke Universiteit Leuven, Leuven, Belgium
| | - Rein Brys
- Research Institute for Forest and Nature, Geraardsbergen, Belgium
| | - Karl J Duffy
- Department of Biology, University of Naples Federico II, Complesso Monte Sant’Angelo, Naples 80126, Italy
| | - Jana Jersáková
- Department of Biology of Ecosystems, Faculty of Science, University of South Bohemia, České Budějovice, Czech Republic
| | - Tiiu Kull
- Institute of Agricultural and Environmental Sciences, Estonian University of Life Sciences, 5 Kreutzwaldi, 51014 Tartu, Estonia
| | - Marc-André Selosse
- Institut Systématique Evolution Biodiversité, Muséum national d’Histoire naturelle, CNRS, Sorbonne Université, Paris, France
- Department of Plant Taxonomy and Nature Conservation, University of Gdansk, Wita Stwosza 59, 80-308 Gdansk, Poland
| | - Spyros Tsiftsis
- Department of Forest and Natural Environment Sciences, International Hellenic University, GR-66132, Drama, Greece
| | - Julita Minasiewicz
- Department of Plant Taxonomy and Nature Conservation, University of Gdansk, Wita Stwosza 59, 80-308 Gdansk, Poland
| | - Hans Jacquemyn
- Department of Biology, Plant Conservation and Population Biology, Katholieke Universiteit Leuven, Leuven, Belgium
| |
Collapse
|
22
|
Zhou C, Xia S, Wen Q, Song Y, Jia Q, Wang T, Liu L, Ouyang T. Genetic structure of an endangered species Ormosia henryi in southern China, and implications for conservation. BMC PLANT BIOLOGY 2023; 23:220. [PMID: 37098472 PMCID: PMC10131447 DOI: 10.1186/s12870-023-04231-w] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 01/06/2023] [Accepted: 04/15/2023] [Indexed: 06/19/2023]
Abstract
BACKGROUND The evergreen broadleaved forest (EBLF) is an iconic vegetation type of East Asia, and it contributes fundamentally to biodiversity-based ecosystem functioning and services. However, the native habitat of EBLFs keeps on decreasing due to anthropogenic activities. Ormosia henryi is a valuable rare woody species in EBLFs that is particularly sensitive to habitat loss. In this study, ten natural populations of O. henryi in southern China were sampled, and then genotyping by sequencing (GBS) was applied to elucidate the standing genetic variation and population structure of this endangered species. RESULTS In ten O. henryi populations, 64,158 high-quality SNPs were generated by GBS. Based on these markers, a relatively low level of genetic diversity was found with the expected heterozygosity (He) ranging from 0.2371 to 0.2901. Pairwise FST between populations varied from 0.0213 to 0.1652, indicating a moderate level of genetic differentiation. However, contemporary gene flow between populations were rare. Assignment test and principal component analysis (PCA) both supported that O. henryi populations in southern China could be divided into four genetic groups, and prominent genetic admixture was found in those populations located in southern Jiangxi Province. Mantel tests and multiple matrix regression with randomization (MMRR) analyses suggested that isolation by distance (IBD) could be the possible reason for describing the current population genetic structure. In addition, the effective population size (Ne) of O. henryi was extremely small, and showed a continuous declining trend since the Last Glacial Period. CONCLUSIONS Our results indicate that the endangered status of O. henryi is seriously underestimated. Artificial conservation measures should be applied as soon as possible to prevent O. henryi from the fate of extinction. Further studies are needed to elucidate the mechanism that leading to the continuous loss of genetic diversity in O. henryi and help to develop a better conservation strategy.
Collapse
Affiliation(s)
- Chengchuan Zhou
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China
| | - Shiqi Xia
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China
| | - Qiang Wen
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China
| | - Ying Song
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China
| | - Quanquan Jia
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China
| | - Tian Wang
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China
| | - Liting Liu
- Identification and Evaluation Center for Forest Germplasm Resources in Jiangxi Province, Jiangxi Academy of Forestry, Nanchang, China.
| | - Tianlin Ouyang
- Jiangxi Provincial Forestry Science and Technology Experiment Center, Ganzhou, China.
| |
Collapse
|
23
|
Maas DL, Prost S, de Leeuw CA, Bi K, Smith LL, Purwanto P, Aji LP, Tapilatu RF, Gillespie RG, Becking LE. Sponge diversification in marine lakes: Implications for phylogeography and population genomic studies on sponges. Ecol Evol 2023; 13:e9945. [PMID: 37066063 PMCID: PMC10099488 DOI: 10.1002/ece3.9945] [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: 12/06/2022] [Revised: 03/06/2023] [Accepted: 03/08/2023] [Indexed: 04/18/2023] Open
Abstract
The relative influence of geography, currents, and environment on gene flow within sessile marine species remains an open question. Detecting subtle genetic differentiation at small scales is challenging in benthic populations due to large effective population sizes, general lack of resolution in genetic markers, and because barriers to dispersal often remain elusive. Marine lakes can circumvent confounding factors by providing discrete and replicated ecosystems. Using high-resolution double digest restriction-site-associated DNA sequencing (4826 Single Nucleotide Polymorphisms, SNPs), we genotyped populations of the sponge Suberites diversicolor (n = 125) to test the relative importance of spatial scales (1-1400 km), local environmental conditions, and permeability of seascape barriers in shaping population genomic structure. With the SNP dataset, we show strong intralineage population structure, even at scales <10 km (average F ST = 0.63), which was not detected previously using single markers. Most variation was explained by differentiation between populations (AMOVA: 48.8%) with signatures of population size declines and bottlenecks per lake. Although the populations were strongly structured, we did not detect significant effects of geographic distance, local environments, or degree of connection to the sea on population structure, suggesting mechanisms such as founder events with subsequent priority effects may be at play. We show that the inclusion of morphologically cryptic lineages that can be detected with the COI marker can reduce the obtained SNP set by around 90%. Future work on sponge genomics should confirm that only one lineage is included. Our results call for a reassessment of poorly dispersing benthic organisms that were previously assumed to be highly connected based on low-resolution markers.
Collapse
Affiliation(s)
- Diede L. Maas
- Marine Animal EcologyWageningen University & ResearchWageningenThe Netherlands
| | - Stefan Prost
- LOEWE Centre for Translational Biodiversity GenomicsSenckenberg Natural History MuseumFrankfurt am MainGermany
- South African National Biodiversity InstituteNational Zoological Gardens of South AfricaPretoriaSouth Africa
| | | | - Ke Bi
- Museum of Vertebrate ZoologyUniversity of California BerkeleyBerkeleyCaliforniaUSA
- Computational Genomics Resource Laboratory, California Institute for Quantitative BiosciencesUniversity of California BerkeleyBerkeleyCaliforniaUSA
| | - Lydia L. Smith
- Museum of Vertebrate ZoologyUniversity of California BerkeleyBerkeleyCaliforniaUSA
| | | | - Ludi P. Aji
- Marine Animal EcologyWageningen University & ResearchWageningenThe Netherlands
- Research Centre for Oceanography, Indonesian Institute of SciencesLembaga Ilmu Pengetahuan IndonesiaJakartaIndonesia
| | - Ricardo F. Tapilatu
- Marine Science and Fisheries Departments and Research Center of Pacific Marine ResourcesState University of PapuaManokwariIndonesia
| | - Rosemary G. Gillespie
- Department of Environmental Science, Policy and ManagementUniversity of California BerkeleyBerkeleyCaliforniaUSA
| | - Leontine E. Becking
- Department of Environmental Science, Policy and ManagementUniversity of California BerkeleyBerkeleyCaliforniaUSA
- Aquaculture and Fisheries, Naturalis Biodiversity CenterWageningen University & ResearchWageningenThe Netherlands
| |
Collapse
|
24
|
Insights into the differentiation and adaptation within Circaeasteraceae from Circaeaster agrestis genome sequencing and resequencing. iScience 2023; 26:106159. [PMID: 36895650 PMCID: PMC9988679 DOI: 10.1016/j.isci.2023.106159] [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/20/2022] [Revised: 01/26/2023] [Accepted: 02/03/2023] [Indexed: 02/11/2023] Open
Abstract
Circaeaster agrestis and Kingdonia uniflora are sister species that reproduce sexually and mainly asexually respectively, providing a good system for comparative genome evolution between taxa with different reproductive models. Comparative genome analyses revealed the two species have similar genome size, but C. agrestis encodes many more genes. The gene families specific to C. agrestis show significant enrichment of genes associated with defense response, while those gene families specific to K. uniflora are enriched in genes regulating root system development. Collinearity analyses revealed C. agrestis experienced two rounds of whole-genome duplication. Fst outlier test across 25 C. agrestis populations uncovered a close inter-relationship between abiotic stress and genetic variability. Genetic feature comparisons showed K. uniflora presents much higher genome heterozygosity, transposable element load, linkage disequilibrium degree, and πN/πS ratio. This study provides new insights into understanding the genetic differentiation and adaptation within ancient lineages characterized by multiple reproductive models.
Collapse
|
25
|
Guo JF, Zhao W, Andersson B, Mao JF, Wang XR. Genomic clines across the species boundary between a hybrid pine and its progenitor in the eastern Tibetan Plateau. PLANT COMMUNICATIONS 2023:100574. [PMID: 36906801 PMCID: PMC10363505 DOI: 10.1016/j.xplc.2023.100574] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Subscribe] [Scholar Register] [Received: 10/18/2022] [Revised: 02/09/2023] [Accepted: 03/08/2023] [Indexed: 06/18/2023]
Abstract
Most species have clearly defined distribution ranges and ecological niches. The genetic and ecological causes of species differentiation and the mechanisms that maintain species boundaries between newly evolved taxa and their progenitors are, however, less clearly defined. This study investigated the genetic structure and clines in Pinus densata, a pine of hybrid origin on the southeastern Tibetan Plateau, to gain an understanding of the contemporary dynamics of species barriers. We analyzed genetic diversity in a range-wide collection of P. densata and representative populations of its progenitors, Pinus tabuliformis and Pinus yunnanensis, using exome capture sequencing. We detected four distinct genetic groups within P. densata that reflect its migration history and major gene-flow barriers across the landscape. The demographies of these genetic groups in the Pleistocene were associated with regional glaciation histories. Interestingly, population sizes rebounded rapidly during interglacial periods, suggesting persistence and resilience of the species during the Quaternary ice age. In the contact zone between P. densata and P. yunnanensis, 3.36% of the analyzed loci (57 849) showed exceptional patterns of introgression, suggesting their potential roles in either adaptive introgression or reproductive isolation. These outliers showed strong clines along critical climate gradients and enrichment in a number of biological processes relevant to high-altitude adaptation. This indicates that ecological selection played an important role in generating genomic heterogeneity and a genetic barrier across a zone of species transition. Our study highlights the forces that operate to maintain species boundaries and promote speciation in the Qinghai-Tibetan Plateau and other mountain systems.
Collapse
Affiliation(s)
- Jing-Fang Guo
- National Engineering Research Center of Tree Breeding and Ecological Restoration; State Key Laboratory of Tree Genetics and Breeding; Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education; College of Biological Sciences and Technology, Beijing Forestry University, Beijing 100083, China
| | - Wei Zhao
- Department of Ecology and Environmental Science, Umeå Plant Science Centre, Umeå University, 90187 Umeå, Sweden
| | - Bea Andersson
- Department of Ecology and Environmental Science, Umeå Plant Science Centre, Umeå University, 90187 Umeå, Sweden
| | - Jian-Feng Mao
- National Engineering Research Center of Tree Breeding and Ecological Restoration; State Key Laboratory of Tree Genetics and Breeding; Key Laboratory of Genetics and Breeding in Forest Trees and Ornamental Plants, Ministry of Education; College of Biological Sciences and Technology, Beijing Forestry University, Beijing 100083, China
| | - Xiao-Ru Wang
- Department of Ecology and Environmental Science, Umeå Plant Science Centre, Umeå University, 90187 Umeå, Sweden.
| |
Collapse
|
26
|
Wersebe MJ, Sherman RE, Jeyasingh PD, Weider LJ. The roles of recombination and selection in shaping genomic divergence in an incipient ecological species complex. Mol Ecol 2023; 32:1478-1496. [PMID: 35119153 DOI: 10.1111/mec.16383] [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: 09/06/2021] [Revised: 01/16/2022] [Accepted: 01/28/2022] [Indexed: 11/24/2022]
Abstract
Speciation genomic studies have revealed that genomes of diverging lineages are shaped jointly by the actions of gene flow and selection. These evolutionary forces acting in concert with processes such as recombination and genome features such as gene density shape a mosaic landscape of divergence. We investigated the roles of recombination and gene density in shaping the patterns of differentiation and divergence between the cyclically parthenogenetic ecological sister-taxa, Daphnia pulicaria and Daphnia pulex. First, we assembled a phased chromosome-scale genome assembly using trio-binning for D. pulicaria and constructed a genetic map using an F2-intercross panel to understand sex-specific recombination rate heterogeneity. Finally, we used a ddRADseq data set with broad geographic sampling of D. pulicaria, D. pulex, and their hybrids to understand the patterns of genome-scale divergence and demographic parameters. Our study provides the first sex-specific estimates of recombination rates for a cyclical parthenogen, and unlike other eukaryotic species, we observed male-biased heterochiasmy in D. pulicaria, which may be related to this somewhat unique breeding mode. Additionally, regions of high gene density and recombination are generally more divergent than regions of suppressed recombination. Outlier analysis indicated that divergent genomic regions are probably driven by selection on D. pulicaria, the derived lineage colonizing a novel lake habitat. Together, our study supports a scenario of selection acting on genes related to local adaptation shaping genome-wide patterns of differentiation despite high local recombination rates in this species complex. Finally, we discuss the limitations of our data in light of demographic uncertainty.
Collapse
Affiliation(s)
- Matthew J Wersebe
- Department of Biology, Program in Ecology and Evolutionary Biology, University of Oklahoma, Norman, Oklahoma, USA
| | - Ryan E Sherman
- Department of Integrative Biology, Oklahoma State University, Stillwater, Oklahoma, USA
| | - Punidan D Jeyasingh
- Department of Integrative Biology, Oklahoma State University, Stillwater, Oklahoma, USA
| | - Lawrence J Weider
- Department of Biology, Program in Ecology and Evolutionary Biology, University of Oklahoma, Norman, Oklahoma, USA
| |
Collapse
|
27
|
Pampoulie C, Berg PR, Jentoft S. Hidden but revealed: After years of genetic studies behavioural monitoring combined with genomics uncover new insight into the population dynamics of Atlantic cod in Icelandic waters. Evol Appl 2023; 16:223-233. [PMID: 36793686 PMCID: PMC9923494 DOI: 10.1111/eva.13471] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/28/2021] [Revised: 05/29/2022] [Accepted: 08/18/2022] [Indexed: 11/27/2022] Open
Abstract
Stock structure is of paramount importance for sustainable management of exploited resources. In that context, genetic markers have been used for more than two decades to resolve spatial structure of marine exploited resources and to fully fathom stock dynamics and interactions. While genetic markers such as allozymes and RFLP dominated the debate in the early era of genetics, technology advances have provided scientists with new tools every decade to better assess stock discrimination and interactions (i.e. gene flow). Here, we provide a review of genetic studies performed to understand stock structure of Atlantic cod in Icelandic waters, from the early allozyme approaches to the genomic work currently carried out. We further highlight the importance of the generation of a chromosome-anchored genome assembly together with whole-genome population data, which drastically changed our perception of the possible management units to consider. After nearly 60 years of genetic investigation of Atlantic cod structure in Icelandic waters, genetic (and later genomic) data combined with behavioural monitoring using Data Storage Tags shifted the attention from geographical population structures to behavioural ecotypes. This review also demonstrates the need for future research to further disentangle the impact of these ecotypes (and gene flow among them) on the population structure of Atlantic cod in Icelandic waters. It also highlights the importance of whole-genome data to unravel unexpected within-species diversity related to chromosomal inversions and associated supergenes, which are important to consider for future development of sustainable management programmes of the species within the North Atlantic.
Collapse
Affiliation(s)
| | - Paul Ragnar Berg
- Norwegian Institute for Water Research Oslo Norway.,Department of Natural Sciences, Centre for Coastal Research (CCR) University of Agder Kristiansand Norway
| | - Sissel Jentoft
- Centre for Ecological and Evolutionary Synthesis Oslo Norway
| |
Collapse
|
28
|
Nocchi G, Wang J, Yang L, Ding J, Gao Y, Buggs RJA, Wang N. Genomic signals of local adaptation and hybridization in Asian white birch. Mol Ecol 2023; 32:595-612. [PMID: 36394364 DOI: 10.1111/mec.16788] [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: 07/04/2022] [Revised: 11/03/2022] [Accepted: 11/08/2022] [Indexed: 11/19/2022]
Abstract
Disentangling the numerous processes that affect patterns of genome-wide diversity in widespread tree species has important implications for taxonomy, conservation, and forestry. Here, we investigate the population genomic structure of Asian white birch (Betula platyphylla) in China and seek to explain it in terms of hybridization, demography and adaptation. We generate whole genome sequence data from 83 individuals across the species range in China. Combining this with an existing data set for 79 European and Russian white birches, we show a clear distinction between B. pendula and B. platyphylla, which have sometimes been lumped taxonomically. Genomic diversity of B. platyphylla in north-western China and Central Russia is affected greatly by hybridization with B. pendula. Excluding these hybridized populations, B. platyphylla in China has a linear distribution from north-eastern to south-western China, along the edge of the inland mountainous region. Within this distribution, three genetic clusters are found, which we model as long diverged with subsequent episodes of gene flow. Patterns of covariation between allele frequencies and environmental variables in B. platyphylla suggest the role of natural selection in the distribution of diversity at 7609 SNPs of which 3767 were significantly differentiated among the genetic clusters. The putative adaptive SNPs are distributed throughout the genome and span 1633 genic regions. Of these genic regions, 87 were previously identified as candidates for selective sweeps in Eurasian B. pendula. We use the 7609 environmentally associated SNPs to estimate the risk of nonadaptedness for each sequenced B. platyphylla individual under a scenario of future climate change, highlighting areas where populations may be under future threat from rising temperatures.
Collapse
Affiliation(s)
- Gabriele Nocchi
- School of Biological and Behavioural Sciences, Queen Mary University of London, London, UK.,Royal Botanic Gardens Kew, Richmond, Surrey, UK
| | - Jing Wang
- Key Laboratory for Bio-resources and Eco-environment, College of Life Science, Sichuan University, Chengdu, China
| | - Long Yang
- Agricultural Big-Data Research Centre and College of Plant Protection, Shandong Agricultural University, Tai'an, China
| | - Junyi Ding
- State Forestry and Grassland Administration Key Laboratory of Silviculture in Downstream Areas of the Yellow River, College of Forestry, Shandong Agricultural University, Tai'an, China.,Mountain Tai Forest Ecosystem Research Station of State Forestry and Grassland Administration, College of Forestry, Shandong Agricultural University, Tai'an, China
| | - Ying Gao
- Agricultural Big-Data Research Centre and College of Plant Protection, Shandong Agricultural University, Tai'an, China
| | - Richard J A Buggs
- School of Biological and Behavioural Sciences, Queen Mary University of London, London, UK.,Royal Botanic Gardens Kew, Richmond, Surrey, UK
| | - Nian Wang
- State Forestry and Grassland Administration Key Laboratory of Silviculture in Downstream Areas of the Yellow River, College of Forestry, Shandong Agricultural University, Tai'an, China.,Mountain Tai Forest Ecosystem Research Station of State Forestry and Grassland Administration, College of Forestry, Shandong Agricultural University, Tai'an, China.,State Key Laboratory of Crop Biology, Shandong Agricultural University, Tai'an, China
| |
Collapse
|
29
|
Ottenburghs J, Honka J, Heikkinen ME, Madsen J, Müskens GJDM, Ellegren H. Highly differentiated loci resolve phylogenetic relationships in the Bean Goose complex. BMC Ecol Evol 2023; 23:2. [PMID: 36658479 PMCID: PMC9854053 DOI: 10.1186/s12862-023-02103-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: 08/26/2022] [Accepted: 01/11/2023] [Indexed: 01/20/2023] Open
Abstract
BACKGROUND Reconstructing phylogenetic relationships with genomic data remains a challenging endeavor. Numerous phylogenomic studies have reported incongruent gene trees when analyzing different genomic regions, complicating the search for a 'true' species tree. Some authors have argued that genomic regions of increased divergence (i.e. differentiation islands) reflect the species tree, although other studies have shown that these regions might produce misleading topologies due to species-specific selective sweeps or ancient introgression events. In this study, we tested the extent to which highly differentiated loci can resolve phylogenetic relationships in the Bean Goose complex, a group of goose taxa that includes the Taiga Bean Goose (Anser fabalis), the Tundra Bean Goose (Anser serrirostris) and the Pink-footed Goose (Anser brachyrhynchus). RESULTS First, we show that a random selection of genomic loci-which mainly samples the undifferentiated regions of the genome-results in an unresolved species complex with a monophyletic A. brachyrhynchus embedded within a paraphyletic cluster of A. fabalis and A. serrirostris. Next, phylogenetic analyses of differentiation islands converged upon a topology of three monophyletic clades in which A. brachyrhynchus is sister to A. fabalis, and A. serrirostris is sister to the clade uniting these two species. Close inspection of the locus trees within the differentiated regions revealed that this topology was consistently supported over other phylogenetic arrangements. As it seems unlikely that selection or introgression events have impacted all differentiation islands in the same way, we are convinced that this topology reflects the 'true' species tree. Additional analyses, based on D-statistics, revealed extensive introgression between A. fabalis and A. serrirostris, which partly explains the failure to resolve the species complex with a random selection of genomic loci. Recent introgression between these taxa has probably erased the phylogenetic branching pattern across a large section of the genome, whereas differentiation islands were unaffected by the homogenizing gene flow and maintained the phylogenetic patterns that reflect the species tree. CONCLUSIONS The evolution of the Bean Goose complex can be depicted as a simple bifurcating tree, but this would ignore the impact of introgressive hybridization. Hence, we advocate that the evolutionary relationships between these taxa are best represented as a phylogenetic network.
Collapse
Affiliation(s)
- Jente Ottenburghs
- Department of Evolutionary Biology, University of Uppsala, Uppsala, Sweden.
| | - Johanna Honka
- grid.10858.340000 0001 0941 4873Ecology and Genetics Research Unit, University of Oulu, PO Box 3000, 90014 Oulu, Finland
| | - Marja E. Heikkinen
- grid.10858.340000 0001 0941 4873Ecology and Genetics Research Unit, University of Oulu, PO Box 3000, 90014 Oulu, Finland
| | - Jesper Madsen
- grid.7048.b0000 0001 1956 2722Department of Ecoscience, Aarhus University, C. F. Møllers Allé 8, 8000 Aarhus C, Denmark
| | - Gerhard J. D. M. Müskens
- grid.4818.50000 0001 0791 5666Team Animal Ecology, Wageningen Environmental Research, Wageningen University & Research, Droevendaalsesteeg 3-3A, 6708 PB Wageningen, The Netherlands
| | - Hans Ellegren
- grid.8993.b0000 0004 1936 9457Department of Evolutionary Biology, University of Uppsala, Uppsala, Sweden
| |
Collapse
|
30
|
Strickland K, Räsänen K, Kristjánsson BK, Phillips JS, Einarsson A, Snorradóttir RG, Bartrons M, Jónsson ZO. Genome-phenotype-environment associations identify signatures of selection in a panmictic population of threespine stickleback. Mol Ecol 2023; 32:1708-1725. [PMID: 36627230 DOI: 10.1111/mec.16845] [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: 07/14/2022] [Revised: 12/01/2022] [Accepted: 12/13/2022] [Indexed: 01/12/2023]
Abstract
Adaptive genetic divergence occurs when selection imposed by the environment causes the genomic component of the phenotype to differentiate. However, genomic signatures of natural selection are usually identified without information on which trait is responding to selection by which selective agent(s). Here, we integrate whole-genome sequencing with phenomics and measures of putative selective agents to assess the extent of adaptive divergence in threespine stickleback occupying the highly heterogeneous lake Mývatn, NE Iceland. We find negligible genome wide divergence, yet multiple traits (body size, gill raker structure and defence traits) were divergent along known ecological gradients (temperature, predatory bird densities and water depth). SNP based heritability of all measured traits was high (h2 = 0.42-0.65), indicating adaptive potential for all traits. Environment-association analyses further identified thousands of loci putatively involved in selection, related to genes linked to, for instance, neuron development and protein phosphorylation. Finally, we found that loci linked to water depth were concurrently associated with pelvic spine length variation - supporting the conclusion that divergence in pelvic spine length occurred in the face of gene flow. Our results suggest that whilst there is substantial genetic variation in the traits measured, phenotypic divergence of Mývatn stickleback is mostly weakly associated with environmental gradients, potentially as a result of substantial gene flow. Our study illustrates the value of integrative studies that combine genomic assays of multivariate trait variation with landscape genomics.
Collapse
Affiliation(s)
- Kasha Strickland
- Institute of Ecology and Evolution, School of Biological Sciences, University of Edinburgh, Edinburgh, UK.,Department of Aquaculture and Fish Biology, Hólar University, Sauðárkrókur, Iceland
| | - Katja Räsänen
- Department of Aquatic Ecology, EAWAG and Institute of Integrative Biology, ETH, Zurich, Switzerland.,Department of Biological and Environmental Science, University of Jyväskylä, Jyväskylä, Finland
| | | | - Joseph S Phillips
- Department of Aquaculture and Fish Biology, Hólar University, Sauðárkrókur, Iceland.,Department of Biology, Creighton University, Omaha, Nebraska, USA
| | | | - Ragna G Snorradóttir
- Department of Aquaculture and Fish Biology, Hólar University, Sauðárkrókur, Iceland
| | - Mireia Bartrons
- Aquatic Ecology Group, University of Vic (UVic-UCC), Catalonia, Spain
| | | |
Collapse
|
31
|
Ahmad HI, Khan FA, Khan MA, Imran S, Akhtar RW, Pandupuspitasari NS, Negara W, Chen J. Molecular Evolution of the Bactericidal/Permeability-Increasing Protein (BPIFA1) Regulating the Innate Immune Responses in Mammals. Genes (Basel) 2022; 14:genes14010015. [PMID: 36672756 PMCID: PMC9858190 DOI: 10.3390/genes14010015] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/21/2022] [Revised: 12/06/2022] [Accepted: 12/13/2022] [Indexed: 12/24/2022] Open
Abstract
Bactericidal/permeability-increasing protein, a primary factor of the innate immune system of mammals, participates in natural immune protection against invading bacteria. BPIFA1 actively contributes to host defense via multiple mechanisms, such as antibacterial, surfactant, airway surface liquid control, and immunomodulatory activities. However, the evolutionary history and selection forces on the BPIFA1 gene in mammals during adaptive evolution are poorly understood. This study examined the BPIFA1 gene of humans compared with that of other mammalian species to estimate the selective pressure derived by adaptive evolution. To assess whether or not positive selection occurred, we employed several different possibility tests (M1 vs. M2 and M7 vs. M8). The proportions of positively selected sites were significant, with a likelihood log value of 93.63 for the BPIFA1 protein. The Selecton server was used on the same dataset to reconfirm positive selection for specific sites by employing the Mechanistic-Empirical Combination model, thus providing additional evidence supporting the findings of positive selection. There was convincing evidence for positive selection signals in the BPIFA1 genes of mammalian species, which was more significant for selection signs and creating signals. We performed probability tests comparing various models based on dN/dS ratios to recognize specific codons under positive selection pressure. We identified positively selected sites in the LBP-BPI domain of BPIFA1 proteins in the mammalian genome, including a lipid-binding domain with a very high degree of selectivity for DPPC. BPIFA1 activates the upper airway's innate immune system in response to numerous genetic signals in the mammalian genome. These findings highlight evolutionary advancements in immunoregulatory effects that play a significant role in the antibacterial and antiviral defenses of mammalian species.
Collapse
Affiliation(s)
- Hafiz Ishfaq Ahmad
- Department of Animal Breeding and Genetics, Faculty of Veterinary and Animal Sciences, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
- Correspondence: (H.I.A.); (J.C.)
| | - Faheem Ahmed Khan
- Laboratory of Molecular Biology and Genomics, Faculty of Science and Technology, University of Central Punjab, Lahore 54000, Pakistan
- Research Center for Animal Husbandry, National Research and Innovation Agency, South Tangerang 15314, Indonesia
| | - Musarrat Abbas Khan
- Department of Animal Breeding and Genetics, Faculty of Veterinary and Animal Sciences, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
| | - Safdar Imran
- Department of Animal Breeding and Genetics, Faculty of Veterinary and Animal Sciences, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
| | - Rana Waseem Akhtar
- Department of Animal Breeding and Genetics, Faculty of Veterinary and Animal Sciences, The Islamia University of Bahawalpur, Bahawalpur 63100, Pakistan
| | - Nuruliarizki Shinta Pandupuspitasari
- Laboratory of Animal Nutrition and Feed Science, Animal Science Department, Faculty of Animal and Agricultural Sciences, Universitas Diponegoro, Semarang 50275, Indonesia
| | - Windu Negara
- Research Center for Animal Husbandry, National Research and Innovation Agency, South Tangerang 15314, Indonesia
| | - Jinping Chen
- Guangdong Key Laboratory of Animal Conservation and Resource Utilization, Guangdong Public Laboratory of Wild Animal Conservation and Utilization, Institute of Zoology, Guangdong Academy of Sciences, Guangzhou 510260, China
- Correspondence: (H.I.A.); (J.C.)
| |
Collapse
|
32
|
Zhang J, Zhang S, Zheng Z, Lu Z, Yang Y. Genomic divergence between two sister Ostrya species through linked selection and recombination. Ecol Evol 2022; 12:e9611. [PMID: 36540075 PMCID: PMC9754895 DOI: 10.1002/ece3.9611] [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: 09/06/2022] [Revised: 11/16/2022] [Accepted: 11/22/2022] [Indexed: 12/23/2022] Open
Abstract
Studying the evolution of genomic divergence between lineages is a topical issue in evolutionary biology. However, the evolutionary forces that shape the heterogeneous divergence of the genomic landscape are still poorly understood. Here, two wind-pollinated sister-species (Ostrya japonica and O. chinensis) are used to explore what these potential forces might be. A total of 40 individuals from 16 populations across their main distribution areas in China were sampled for genome-wide resequencing. Population demography analyses revealed that these two sister-species diverged at 3.06-4.43 Mya. Both population contraction and increased gene flow were detected during glacial periods, suggesting secondary contact at those times. All three parameters (D XY, π, and ρ) decreased in those regions showing high levels of differentiation (F ST). These findings indicate that linked selection and recombination played a key role in the genomic heterogeneous differentiation between the two Ostrya species. Genotype-environment association analyses showed that precipitation was the most important ecological factor for speciation. Such environmentally related genes and positive selection genes may have contributed to local adaptation and the maintenance of species boundaries.
Collapse
Affiliation(s)
- Jin Zhang
- State Key Laboratory of Grassland Agro‐Ecosystems, College of EcologyLanzhou UniversityLanzhouChina
| | - Shangzhe Zhang
- State Key Laboratory of Grassland Agro‐Ecosystems, College of EcologyLanzhou UniversityLanzhouChina
| | - Zeyu Zheng
- State Key Laboratory of Grassland Agro‐Ecosystems, College of EcologyLanzhou UniversityLanzhouChina
| | - Zhiqiang Lu
- CAS Key Laboratory of Tropical Forest Ecology, Xishuangbanna Tropical Botanical GardenChinese Academy of SciencesMenglaChina
| | - Yongzhi Yang
- State Key Laboratory of Grassland Agro‐Ecosystems, College of EcologyLanzhou UniversityLanzhouChina
| |
Collapse
|
33
|
Anderson SAS, López-Fernández H, Weir JT. Ecology and the origin of non-ephemeral species. Am Nat 2022; 201:619-638. [PMID: 37130236 DOI: 10.1086/723763] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/13/2022]
Abstract
AbstractResearch over the past three decades has shown that ecology-based extrinsic reproductive barriers can rapidly arise to generate incipient species-but such barriers can also rapidly dissolve when environments change, resulting in incipient species collapse. Understanding the evolution of unconditional, "intrinsic" reproductive barriers is therefore important for understanding the longer-term buildup of biodiversity. In this article, we consider ecology's role in the evolution of intrinsic reproductive isolation. We suggest that this topic has fallen into a gap between disciplines: while evolutionary ecologists have traditionally focused on the rapid evolution of extrinsic isolation between co-occurring ecotypes, speciation geneticists studying intrinsic isolation in other taxa have devoted little attention to the ecological context in which it evolves. We argue that for evolutionary ecology to close this gap, the field will have to expand its focus beyond rapid adaptation and its traditional model systems. Synthesizing data from several subfields, we present circumstantial evidence for and against different forms of ecological adaptation as promoters of intrinsic isolation and discuss alternative forces that may be significant. We conclude by outlining complementary approaches that can better address the role of ecology in the evolution of nonephemeral reproductive barriers and, by extension, less ephemeral species.
Collapse
|
34
|
Xi L, Sun Y, Xu T, Wang Z, Chiu MY, Plouviez S, Jollivet D, Qiu J. Phylogenetic divergence and population genetics of the hydrothermal vent annelid genus
Hesiolyra
along the East Pacific Rise: Reappraisal using multi‐locus data. DIVERS DISTRIB 2022. [DOI: 10.1111/ddi.13653] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/03/2022] Open
Affiliation(s)
- Leyi Xi
- Department of Biology Hong Kong Baptist University Hong Kong China
| | - Yanan Sun
- Department of Biology Hong Kong Baptist University Hong Kong China
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) Guangzhou China
- Department of Ocean Science The Hong Kong University of Science and Technology Hong Kong China
| | - Ting Xu
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) Guangzhou China
- Department of Ocean Science The Hong Kong University of Science and Technology Hong Kong China
| | - Zhi Wang
- State Key Laboratory of Marine Environmental Science, College of Ocean and Earth Sciences Xiamen University Xiamen China
| | - Man Ying Chiu
- Department of Biology Hong Kong Baptist University Hong Kong China
| | - Sophie Plouviez
- Department of Biology University of Louisiana at Lafayette Lafayette Louisiana USA
| | - Didier Jollivet
- Sorbonne Université‐CNRS, UMR 7144 Adaptation et Diversité en Milieu Marin, Equipe DyDiv Station Biologique de Roscoff Roscoff France
| | - Jian‐Wen Qiu
- Department of Biology Hong Kong Baptist University Hong Kong China
- Southern Marine Science and Engineering Guangdong Laboratory (Guangzhou) Guangzhou China
| |
Collapse
|
35
|
Yusuf LH, Tyukmaeva V, Hoikkala A, Ritchie MG. Divergence and introgression among the virilis group of Drosophila. Evol Lett 2022; 6:537-551. [PMID: 36579165 PMCID: PMC9783487 DOI: 10.1002/evl3.301] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/10/2022] [Revised: 09/23/2022] [Accepted: 10/12/2022] [Indexed: 12/03/2022] Open
Abstract
Speciation with gene flow is now widely regarded as common. However, the frequency of introgression between recently diverged species and the evolutionary consequences of gene flow are still poorly understood. The virilis group of Drosophila contains 12 species that are geographically widespread and show varying levels of prezygotic and postzygotic isolation. Here, we use de novo genome assemblies and whole-genome sequencing data to resolve phylogenetic relationships and describe patterns of introgression and divergence across the group. We suggest that the virilis group consists of three, rather than the traditional two, subgroups. Some genes undergoing rapid sequence divergence across the group were involved in chemical communication and desiccation tolerance, and may be related to the evolution of sexual isolation and adaptation. We found evidence of pervasive phylogenetic discordance caused by ancient introgression events between distant lineages within the group, and more recent gene flow between closely related species. When assessing patterns of genome-wide divergence in species pairs across the group, we found no consistent genomic evidence of a disproportionate role for the X chromosome as has been found in other systems. Our results show how ancient and recent introgressions confuse phylogenetic reconstruction, but may play an important role during early radiation of a group.
Collapse
Affiliation(s)
- Leeban H. Yusuf
- Centre for Biological Diversity, School of BiologyUniversity of St AndrewsSt AndrewsKY16 9THUnited Kingdom
| | - Venera Tyukmaeva
- Centre for Biological Diversity, School of BiologyUniversity of St AndrewsSt AndrewsKY16 9THUnited Kingdom,Department of Evolution, Ecology and BehaviourUniversity of LiverpoolLiverpoolL69 7ZBUnited Kingdom
| | - Anneli Hoikkala
- Department of Biological and Environmental ScienceUniversity of JyväskyläJyväskylä40014Finland
| | - Michael G. Ritchie
- Centre for Biological Diversity, School of BiologyUniversity of St AndrewsSt AndrewsKY16 9THUnited Kingdom
| |
Collapse
|
36
|
Wang L, Liu S, Yang Y, Meng Z, Zhuang Z. Linked selection, differential introgression and recombination rate variation promote heterogeneous divergence in a pair of yellow croakers. Mol Ecol 2022; 31:5729-5744. [PMID: 36111361 PMCID: PMC9828471 DOI: 10.1111/mec.16693] [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: 04/24/2021] [Accepted: 09/06/2022] [Indexed: 01/13/2023]
Abstract
Understanding the mechanisms underlying heterogeneous genomic divergence is of particular interest in evolutionary biology. Highly differentiated genomic regions, known as genomic islands, often evolve between diverging lineages. These genomic islands may be related to selection promoting adaptation or reproductive isolation. Based on whole genome assembly and genome-wide RAD sequencing in a pair of yellow croakers (genus: Larimichthys), we investigated the evolutionary processes shaping genomic landscapes of divergence. Demographic modelling indicated that the two species diverged following a secondary contact scenario, where differential introgression and linked selection were suggested to be involved in heterogeneous genomic divergence. We identified reduced recombination rate in genomic islands and a relatively good conservation of both genetic diversity and recombination landscapes between species, which highlight the roles of linked selection and recombination rate variation in promoting heterogeneous divergence in the common ancestral lineage of the two species. In addition, we found a positive correlation between differentiation (FST ) and absolute sequence divergence (Dxy ), and elevated Dxy in genomic islands, indicating that the genomic landscape of divergence was not shaped by linked selection alone. Restricted gene flow in highly differentiated regions has probably remodelled the landscape of heterogeneous genomic divergence. This study highlights that highly differentiated genomic regions can also arise from a combination of linked selection and differential gene flow in interaction with varying recombination rates.
Collapse
Affiliation(s)
- Le Wang
- State Key Laboratory of Biocontrol, Institute of Aquatic Economic Animals and the Guangdong Province Key Laboratory for Aquatic Economic Animals, School of Life SciencesSun Yat‐sen UniversityGuangzhouChina,Molecular Population Genetics Group, Temasek Life Sciences Laboratory, 1 Research LinkNational University of SingaporeSingapore CitySingapore
| | - Shufang Liu
- Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences & Function Laboratory for Marine Biology and BiotechnologyQingdao National Laboratory for Marine Science and TechnologyQingdaoChina
| | - Yang Yang
- State Key Laboratory of Biocontrol, Institute of Aquatic Economic Animals and the Guangdong Province Key Laboratory for Aquatic Economic Animals, School of Life SciencesSun Yat‐sen UniversityGuangzhouChina
| | - Zining Meng
- State Key Laboratory of Biocontrol, Institute of Aquatic Economic Animals and the Guangdong Province Key Laboratory for Aquatic Economic Animals, School of Life SciencesSun Yat‐sen UniversityGuangzhouChina,Southern Laboratory of Ocean Science and EngineeringZhuhaiChina
| | - Zhimeng Zhuang
- Yellow Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences & Function Laboratory for Marine Biology and BiotechnologyQingdao National Laboratory for Marine Science and TechnologyQingdaoChina
| |
Collapse
|
37
|
Spies I, Tarpey C, Kristiansen T, Fisher M, Rohan S, Hauser L. Genomic differentiation in Pacific cod using
P
ool‐
S
eq. Evol Appl 2022; 15:1907-1924. [PMID: 36426128 PMCID: PMC9679252 DOI: 10.1111/eva.13488] [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: 03/15/2022] [Revised: 09/05/2022] [Accepted: 09/21/2022] [Indexed: 11/28/2022] Open
Abstract
Patterns of genetic differentiation across the genome can provide insight into selective forces driving adaptation. We used pooled whole genome sequencing, gene annotation, and environmental covariates to evaluate patterns of genomic differentiation and to investigate mechanisms responsible for divergence among proximate Pacific cod (Gadus macrocephalus) populations from the Bering Sea and Aleutian Islands and more distant Washington Coast cod. Samples were taken from eight spawning locations, three of which were replicated to estimate consistency in allele frequency estimation. A kernel smoothing moving weighted average of relative divergence (FST) identified 11 genomic islands of differentiation between the Aleutian Islands and Bering Sea samples. In some islands of differentiation, there was also elevated absolute divergence (dXY) and evidence for selection, despite proximity and potential for gene flow. Similar levels of absolute divergence (dXY) but roughly double the relative divergence (FST) were observed between the distant Bering Sea and Washington Coast samples. Islands of differentiation were much smaller than the four large inversions among Atlantic cod ecotypes. Islands of differentiation between the Bering Sea and Aleutian Island were associated with SNPs from five vision system genes, which can be associated with feeding, predator avoidance, orientation, and socialization. We hypothesize that islands of differentiation between Pacific cod from the Bering Sea and Aleutian Islands provide evidence for adaptive differentiation despite gene flow in this commercially important marine species.
Collapse
Affiliation(s)
- Ingrid Spies
- Resource Ecology and Fisheries Management Division Alaska Fisheries Science Center Seattle Washington USA
| | - Carolyn Tarpey
- School of Aquatic and Fishery Sciences University of Washington Seattle Washington USA
| | | | - Mary Fisher
- School of Aquatic and Fishery Sciences University of Washington Seattle Washington USA
| | - Sean Rohan
- Resource Assessment and Conservation Engineering Division Alaska Fisheries Science Center Seattle Washington USA
| | - Lorenz Hauser
- Resource Ecology and Fisheries Management Division Alaska Fisheries Science Center Seattle Washington USA
| |
Collapse
|
38
|
Hollenbeck CM, Portnoy DS, Garcia de la Serrana D, Magnesen T, Matejusova I, Johnston IA. Temperature-associated selection linked to putative chromosomal inversions in king scallop ( Pecten maximus). Proc Biol Sci 2022; 289:20221573. [PMID: 36196545 PMCID: PMC9532988 DOI: 10.1098/rspb.2022.1573] [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: 11/27/2022] Open
Abstract
The genomic landscape of divergence—the distribution of differences among populations or species across the genome—is increasingly characterized to understand the role that microevolutionary forces such as natural selection and recombination play in causing and maintaining genetic divergence. This line of inquiry has also revealed chromosome structure variation to be an important factor shaping the landscape of adaptive genetic variation. Owing to a high prevalence of chromosome structure variation and the strong pressure for local adaptation necessitated by their sessile nature, bivalve molluscs are an ideal taxon for exploring the relationship between chromosome structure variation and local adaptation. Here, we report a population genomic survey of king scallop (Pecten maximus) across its natural range in the northeastern Atlantic Ocean, using a recent chromosome-level genome assembly. We report the presence of at least three large (12–22 Mb), putative chromosomal inversions associated with sea surface temperature and whose frequencies are in contrast to neutral population structure. These results highlight a potentially large role for recombination-suppressing chromosomal inversions in local adaptation and suggest a hypothesis to explain the maintenance of differences in reproductive timing found at relatively small spatial scales across king scallop populations.
Collapse
Affiliation(s)
- Christopher M Hollenbeck
- Department of Life Sciences, Texas A&M University Corpus Christi, 6300 Ocean Drive, Corpus Christi, TX 78412, USA.,Texas A&M AgriLife Research, College Station, TX, USA
| | - David S Portnoy
- Department of Life Sciences, Texas A&M University Corpus Christi, 6300 Ocean Drive, Corpus Christi, TX 78412, USA
| | - Daniel Garcia de la Serrana
- Department of Cell Biology, Physiology and Immunology, Faculty of Biology, University of Barcelona, Barcelona, Spain
| | - Thorolf Magnesen
- Department of Biological Sciences, University of Bergen, Thormøhlensgt 53B, Bergen, Norway
| | - Iveta Matejusova
- Marine Science Scotland, Marine Laboratory, 375 Victoria Road, Aberdeen AB11 9DB, UK
| | - Ian A Johnston
- Scottish Oceans Institute, School of Biology, University of St Andrews, St Andrews, Fife KY16 8LB, UK.,Xelect Ltd, Horizon House, Abbey Walk, St Andrews KY16 9LB, UK
| |
Collapse
|
39
|
Nikolakis ZL, Adams RH, Wade KJ, Lund AJ, Carlton EJ, Castoe TA, Pollock DD. Prospects for genomic surveillance for selection in schistosome parasites. FRONTIERS IN EPIDEMIOLOGY 2022; 2:932021. [PMID: 38455290 PMCID: PMC10910990 DOI: 10.3389/fepid.2022.932021] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/29/2022] [Accepted: 09/12/2022] [Indexed: 03/09/2024]
Abstract
Schistosomiasis is a neglected tropical disease caused by multiple parasitic Schistosoma species, and which impacts over 200 million people globally, mainly in low- and middle-income countries. Genomic surveillance to detect evidence for natural selection in schistosome populations represents an emerging and promising approach to identify and interpret schistosome responses to ongoing control efforts or other environmental factors. Here we review how genomic variation is used to detect selection, how these approaches have been applied to schistosomes, and how future studies to detect selection may be improved. We discuss the theory of genomic analyses to detect selection, identify experimental designs for such analyses, and review studies that have applied these approaches to schistosomes. We then consider the biological characteristics of schistosomes that are expected to respond to selection, particularly those that may be impacted by control programs. Examples include drug resistance, host specificity, and life history traits, and we review our current understanding of specific genes that underlie them in schistosomes. We also discuss how inherent features of schistosome reproduction and demography pose substantial challenges for effective identification of these traits and their genomic bases. We conclude by discussing how genomic surveillance for selection should be designed to improve understanding of schistosome biology, and how the parasite changes in response to selection.
Collapse
Affiliation(s)
- Zachary L. Nikolakis
- Department of Biology, University of Texas at Arlington, Arlington, TX, United States
| | - Richard H. Adams
- Department of Biological and Environmental Sciences, Georgia College and State University, Milledgeville, GA, United States
| | - Kristen J. Wade
- Department of Biochemistry and Molecular Genetics, University of Colorado School of Medicine, Aurora, CO, United States
| | - Andrea J. Lund
- Department of Environmental and Occupational Health, Colorado School of Public Health, University of Colorado, Anschutz, Aurora, CO, United States
| | - Elizabeth J. Carlton
- Department of Environmental and Occupational Health, Colorado School of Public Health, University of Colorado, Anschutz, Aurora, CO, United States
| | - Todd A. Castoe
- Department of Biology, University of Texas at Arlington, Arlington, TX, United States
| | - David D. Pollock
- Department of Biochemistry and Molecular Genetics, University of Colorado School of Medicine, Aurora, CO, United States
| |
Collapse
|
40
|
Proulx SR, Teotónio H. Selection on modifiers of genetic architecture under migration load. PLoS Genet 2022; 18:e1010350. [PMID: 36070315 PMCID: PMC9484686 DOI: 10.1371/journal.pgen.1010350] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/13/2022] [Revised: 09/19/2022] [Accepted: 07/20/2022] [Indexed: 11/22/2022] Open
Abstract
Gene flow between populations adapting to differing local environmental conditions might be costly because individuals can disperse to habitats where their survival is low or because they can reproduce with locally maladapted individuals. The amount by which the mean relative population fitness is kept below one creates an opportunity for modifiers of the genetic architecture to spread due to selection. Prior work that separately considered modifiers changing dispersal, recombination rates, or altering dominance or epistasis, has typically focused on the direction of selection rather than its absolute magnitude. We here develop methods to determine the strength of selection on modifiers of the genetic architecture, including modifiers of the dispersal rate, in populations that have previously evolved local adaptation. We consider scenarios with up to five loci contributing to local adaptation and derive a new model for the deterministic spread of modifiers. We find that selection for modifiers of epistasis and dominance is stronger than selection for decreased recombination, and that selection for partial reductions in recombination are extremely weak, regardless of the number of loci contributing to local adaptation. The spread of modifiers that reduce dispersal depends on the number of loci, epistasis and extent of local adaptation in the ancestral population. We identify a novel effect, that modifiers of dominance are more strongly selected when they are unlinked to the locus that they modify. These findings help explain population differentiation and reproductive isolation and provide a benchmark to compare selection on modifiers under finite population sizes and demographic stochasticity. When populations of a species are spread over different habitats the populations can adapt to their local conditions, provided dispersal between habitats is low enough. Natural selection allows the populations to maintain local adaptation, but dispersal and gene flow create a cost called the migration load. The migration load measures how much fitness is lost because of dispersal between different habitats, and also creates an opportunity for selection to act on the arrangement and interaction between genes that are involved in local adaptation. Modifier genes can spread in these linked populations and cause functional, local adaptation genes, to become more closely linked on a chromosome, or change the way that these genes are expressed so that the locally adapted gene copy becomes dominant. We modeled this process and found that selection on modifiers that create tighter linkage between locally adapted genes is generally weak, and modifiers that cause gene interactions are more strongly selected. Even after these gene interactions have begun to evolve, further selection for increased gene interaction is still strong. Our results show that populations are more likely to adapt to local conditions by evolving new gene interactions than by evolving tightly linked gene clusters.
Collapse
Affiliation(s)
- Stephen R. Proulx
- Department of Ecology, Evolution, and Marine Biology, UC Santa Barbara, Santa Barbara, California, United States of America
- * E-mail:
| | | |
Collapse
|
41
|
Thomas-Bulle C, Bertrand D, Nagarajan N, Copley RR, Corre E, Hourdez S, Bonnivard É, Claridge-Chang A, Jollivet D. Genomic patterns of divergence in the early and late steps of speciation of the deep-sea vent thermophilic worms of the genus Alvinella. BMC Ecol Evol 2022; 22:106. [PMID: 36057769 PMCID: PMC9441076 DOI: 10.1186/s12862-022-02057-y] [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: 02/28/2022] [Accepted: 08/05/2022] [Indexed: 11/10/2022] Open
Abstract
Background The transient and fragmented nature of the deep-sea hydrothermal environment made of ridge subduction, plate collision and the emergence of new rifts is currently acting to separate of vent populations, promoting local adaptation and contributing to bursts of speciation and species specialization. The tube-dwelling worms Alvinella pompejana called the Pompeii worm and its sister species A. caudata live syntopically on the hottest part of deep-sea hydrothermal chimneys along the East Pacific Rise. They are exposed to extreme thermal and chemical gradients, which vary greatly in space and time, and thus represent ideal candidates for understanding the evolutionary mechanisms at play in the vent fauna evolution. Results We explored genomic patterns of divergence in the early and late stages of speciation of these emblematic worms using transcriptome assemblies and the first draft genome to better understand the relative role of geographic isolation and habitat preference in their genome evolution. Analyses were conducted on allopatric populations of Alvinella pompejana (early stage of separation) and between A. pompejana and its syntopic species Alvinella caudata (late stage of speciation). We first identified divergent genomic regions and targets of selection as well as their position in the genome over collections of orthologous genes and, then, described the speciation dynamics by documenting the annotation of the most divergent and/or positively selected genes involved in the isolation process. Gene mapping clearly indicated that divergent genes associated with the early stage of speciation, although accounting for nearly 30% of genes, are highly scattered in the genome without any island of divergence and not involved in gamete recognition or mito-nuclear incompatibilities. By contrast, genomes of A. pompejana and A. caudata are clearly separated with nearly all genes (96%) exhibiting high divergence. This congealing effect however seems to be linked to habitat specialization and still allows positive selection on genes involved in gamete recognition, as a possible long-duration process of species reinforcement.
Conclusion Our analyses highlight the non-negligible role of natural selection on both the early and late stages of speciation in the iconic thermophilic worms living on the walls of deep-sea hydrothermal chimneys. They shed light on the evolution of gene divergence during the process of speciation and species specialization over a very long period of time. Supplementary Information The online version contains supplementary material available at 10.1186/s12862-022-02057-y.
Collapse
|
42
|
Wang Z, Cai X, Jiang X, Xia Q, Li L, Lu B. Sympatric genetic divergence between early- and late-season weedy rice populations. THE NEW PHYTOLOGIST 2022; 235:2066-2080. [PMID: 35637631 PMCID: PMC9544748 DOI: 10.1111/nph.18288] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 12/26/2021] [Accepted: 05/13/2022] [Indexed: 06/15/2023]
Abstract
Sympatric genetic divergence is the most appealing and controversial pattern in the theory of ecological speciation. Examples that support sympatric genetic divergence in plant species are extremely rare. Solid evidence of sympatric genetic divergence will provide deep insights for revealing the underlying mechanisms of ecological speciation. We analysed the total genomic DNA sequences of 120 weedy rice (WR; Oryza sativa f. spontanea) plants, representing three WR population pairs separately from three early- and late-season rice fields, in comparison with those of the co-occurring rice cultivars and other rice materials. We detected substantial genetic divergence within the pairs of the sympatric early- and late-season WR populations, although genetic divergence was unevenly distributed across the genomes. Restricted gene flow was determined between the sympatric WR populations, resulting in their distinct genetic structures. We also detected relatively low genetic diversity that was likely to be associated with stronger selection in early-season WR populations. Our findings provide strong evidence for sympatric genetic divergence between the WR populations in the same fields but in different seasons. We conclude that temporal isolation plays an important role in creating genetic divergence between sympatric populations/species in plants.
Collapse
Affiliation(s)
- Zhi Wang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Department of Ecology and Evolutionary BiologyFudan UniversitySonghu Road 2005Shanghai200438China
| | - Xingxing Cai
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Department of Ecology and Evolutionary BiologyFudan UniversitySonghu Road 2005Shanghai200438China
| | - Xiao‐Qi Jiang
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Department of Ecology and Evolutionary BiologyFudan UniversitySonghu Road 2005Shanghai200438China
| | - Qi‐Yu Xia
- Key Laboratory of Biology and Genetic Resources of Tropical Crops, Hainan Key Laboratory for Biosafety Monitoring and Molecular Breeding in Off‐Season Reproduction RegionsInstitute of Tropical Bioscience and Biotechnology, CATASHaikou571101China
| | - Lin‐Feng Li
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Department of Ecology and Evolutionary BiologyFudan UniversitySonghu Road 2005Shanghai200438China
| | - Bao‐Rong Lu
- Ministry of Education Key Laboratory for Biodiversity Science and Ecological Engineering, Department of Ecology and Evolutionary BiologyFudan UniversitySonghu Road 2005Shanghai200438China
| |
Collapse
|
43
|
Daoui O, Bennaid H, Kbaich MA, Mhaidi I, Aderdour N, Rhinane H, Bouhout S, Akarid K, Lemrani M. Environmental, Climatic, and Parasite Molecular Factors Impacting the Incidence of Cutaneous Leishmaniasis Due to Leishmania tropica in Three Moroccan Foci. Microorganisms 2022; 10:microorganisms10091712. [PMID: 36144314 PMCID: PMC9506065 DOI: 10.3390/microorganisms10091712] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/28/2022] [Revised: 08/18/2022] [Accepted: 08/22/2022] [Indexed: 11/16/2022] Open
Abstract
Cutaneous leishmaniasis (CL) occurring due to Leishmania tropica is a public health problem in Morocco. The distribution and incidence of this form of leishmaniasis have increased in an unusual way in the last decade, and the control measures put in place are struggling to slow down the epidemic. This study was designed to assess the impact of climatic and environmental factors on CL in L. tropica foci. The data collected included CL incidence and climatic and environmental factors across three Moroccan foci (Foum Jemaa, Imintanout, and Ouazzane) from 2000 to 2019. Statistical analyses were performed using the linear regression model. An association was found between the occurrence of CL in Imintanout and temperature and humidity (r2 = 0.6076, df = (1.18), p-value = 3.09 × 10−5; r2 = 0.6306, df = (1.18), p-value = 1.77 × 10−5). As a second objective of our study, we investigated the population structure of L.tropica in these three foci, using the nuclear marker internal transcribed spacer 1 (ITS1). Our results showed a low-to-medium level of geographic differentiation among the L.tropica populations using pairwise differentiation. Molecular diversity indices showed a high genetic diversity in Foum Jemaa and Imintanout; indeed, 29 polymorphic sites were identified, leading to the definition of 13 haplotypes. Tajima’s D and Fu’s F test statistics in all populations were not statistically significant, and consistent with a population at drift–mutation equilibrium. Further analysis, including additional DNA markers and a larger sample size, could provide a more complete perspective of L. tropica’s population structure in these three regions. In addition, further research is needed to better understand the impact of climatic conditions on the transmission cycle of Leishmania, allowing both for the development of effective control measures, and for the development of a predictive model for this parasitosis.
Collapse
Affiliation(s)
- Othmane Daoui
- Laboratory of Parasitology and Vector-Borne-Diseases, Institut Pasteur du Maroc, Casablanca 20250, Morocco
- Molecular Genetics and Immunophysiopathology Research Team, Health and Environment Laboratory, Aïn Chock Faculty of Sciences, Hassan II University of Casablanca, Casablanca 20000, Morocco
| | - Hamza Bennaid
- Laboratory of Parasitology and Vector-Borne-Diseases, Institut Pasteur du Maroc, Casablanca 20250, Morocco
- Information Retrieval and Data Analytics Laboratory, National School of Computer Science and System Analysis (Ensias), Rabat 10112, Morocco
| | - Mouad Ait Kbaich
- Laboratory of Parasitology and Vector-Borne-Diseases, Institut Pasteur du Maroc, Casablanca 20250, Morocco
| | - Idris Mhaidi
- Laboratory of Parasitology and Vector-Borne-Diseases, Institut Pasteur du Maroc, Casablanca 20250, Morocco
| | - Nacer Aderdour
- Geosciences Laboratory, Faculty of Sciences Ain Chock, Hassan II University of Casablanca, Casablanca 20000, Morocco
| | - Hassan Rhinane
- Geosciences Laboratory, Faculty of Sciences Ain Chock, Hassan II University of Casablanca, Casablanca 20000, Morocco
| | - Souad Bouhout
- Directorate of Epidemiology, Division of Infectious Diseases, Service of Parasitic Diseases, Ministry of Health of Morocco, Rabat 10020, Morocco
| | - Khadija Akarid
- Molecular Genetics and Immunophysiopathology Research Team, Health and Environment Laboratory, Aïn Chock Faculty of Sciences, Hassan II University of Casablanca, Casablanca 20000, Morocco
| | - Meryem Lemrani
- Laboratory of Parasitology and Vector-Borne-Diseases, Institut Pasteur du Maroc, Casablanca 20250, Morocco
- Correspondence:
| |
Collapse
|
44
|
Introgression at the emerging secondary contact zone of magpie Pica pica subspecies (Aves: Corvidae): integrating data on nuclear and mitochondrial markers, vocalizations, and field observations. ORG DIVERS EVOL 2022. [DOI: 10.1007/s13127-022-00568-6] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/26/2022]
Abstract
Abstract
Zones of secondary contact provide a good opportunity to investigate the origin and dynamics of reproductive isolation between related populations. We analyzed genetic and phenotypic patterns and gene flow between two subspecies of the Eurasian magpie Pica pica s.l. which recently came into contact after presumably long periods of isolation. We describe the distribution of subspecies in a young contact zone at Argun’ river basin in southern Siberia where populations occur in parapatry and an older hybrid population in eastern Mongolia. Based on genome-wide SNP data, we analyzed patterns and strength of gene flow between the subspecies. Our results indicate occasional hybridization with backcrossing and asymmetric introgression along a wide range in Transbaikalia and locally in eastern Mongolia. Males of P. p. jankowskii apparently exhibit higher dispersal ability towards the west compared to P. p. leucoptera (towards the east). The former occasionally migrates to eastern Mongolia and Transbaikalia where introgression of nuclear, but not mitochondrial DNA was evident. Bioacoustic investigations showed differences between the subspecies in speed and structure of vocalization. We discovered intermediate calls of hybrid magpies and bilingual birds alternating calls that are typical for the two taxa. Furthermore, we found dramatically decreased reproductive success in hybridogeneous populations. By complementing our results with established phylogeographic patterns of P. pica s.l. based on a mitochondrial marker sequence, and considering indications of sterility of hybrids in the contact zone, we propose to elevate the two corresponding subspecies to species level: P. pica for the western form and P. serica for the eastern form.
Collapse
|
45
|
Wang G, Lai H, Bi S, Guo D, Zhao X, Chen X, Liu S, Liu X, Su Y, Yi H, Li G. ddRAD‐Seq
reveals evolutionary insights into population differentiation and the cryptic phylogeography of
Hyporhamphus intermedius
in Mainland China. Ecol Evol 2022; 12:e9053. [PMID: 35813915 PMCID: PMC9251877 DOI: 10.1002/ece3.9053] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/08/2021] [Revised: 05/28/2022] [Accepted: 06/08/2022] [Indexed: 11/12/2022] Open
Abstract
Species differentiation and local adaptation in heterogeneous environments have attracted much attention, although little is known about the mechanisms involved. Hyporhamphus intermedius is an anadromous, brackish‐water halfbeak that is widely distributed in coastal areas and hyperdiverse freshwater systems in China, making it an interesting model for research on phylogeography and local adaptation. Here, 156 individuals were sampled at eight sites from heterogeneous aquatic habitats to examine environmental and genetic contributions to phenotypic divergence. Using double‐digest restriction‐site‐associated DNA sequencing (ddRAD‐Seq) in the specimens from the different watersheds, 5498 single nucleotide polymorphisms (SNPs) were found among populations, with obvious population differentiation. We find that present‐day Mainland China populations are structured into distinct genetic clusters stretching from southern and northern ancestries, mirroring geography. Following a transplant event in Plateau Lakes, there were virtually no variations of genetic diversity occurred in two populations, despite the fact two main splits were unveiled in the demographic history. Additionally, dorsal, and anal fin traits varied widely between the southern group and the others, which highlighted previously unrecognized lineages. We then explore genotype–phenotype‐environment associations and predict candidate loci. Subgroup ranges appeared to correspond to geographic regions with heterogeneous hydrological factors, indicating that these features are likely important drivers of diversification. Accordingly, we conclude that genetic and phenotypic polymorphism and a moderate amount of genetic differentiation occurred, which might be ascribed to population subdivision, and the impact of abiotic factors.
Collapse
Affiliation(s)
- Gongpei Wang
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- State Key Laboratory of Ophthalmology Zhongshan Ophthalmic Center Sun Yat‐Sen University Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Han Lai
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Sheng Bi
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Dingli Guo
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Xiaopin Zhao
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Xiaoli Chen
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Shuang Liu
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Xuange Liu
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Yuqin Su
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Huadong Yi
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| | - Guifeng Li
- Guangdong Province Key Laboratory for Aquatic Economic Animals State Key Laboratory of Biocontrol School of Life Sciences Sun Yat‐Sen University Southern Marine Science and Engineering Guangdong Laboratory (Zhuhai) Guangzhou China
- Guangdong Provincial Engineering Technology Research Center for Healthy Breeding of Important Economic Fish Guangzhou China
| |
Collapse
|
46
|
Xiong T, Li X, Yago M, Mallet J. Admixture of evolutionary rates across a butterfly hybrid zone. eLife 2022; 11:e78135. [PMID: 35703474 PMCID: PMC9246367 DOI: 10.7554/elife.78135] [Citation(s) in RCA: 4] [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: 02/23/2022] [Accepted: 06/14/2022] [Indexed: 12/26/2022] Open
Abstract
Hybridization is a major evolutionary force that can erode genetic differentiation between species, whereas reproductive isolation maintains such differentiation. In studying a hybrid zone between the swallowtail butterflies Papilio syfanius and Papilio maackii (Lepidoptera: Papilionidae), we made the unexpected discovery that genomic substitution rates are unequal between the parental species. This phenomenon creates a novel process in hybridization, where genomic regions most affected by gene flow evolve at similar rates between species, while genomic regions with strong reproductive isolation evolve at species-specific rates. Thus, hybridization mixes evolutionary rates in a way similar to its effect on genetic ancestry. Using coalescent theory, we show that the rate-mixing process provides distinct information about levels of gene flow across different parts of genomes, and the degree of rate-mixing can be predicted quantitatively from relative sequence divergence ([Formula: see text]) between the hybridizing species at equilibrium. Overall, we demonstrate that reproductive isolation maintains not only genomic differentiation, but also the rate at which differentiation accumulates. Thus, asymmetric rates of evolution provide an additional signature of loci involved in reproductive isolation.
Collapse
Affiliation(s)
- Tianzhu Xiong
- Department of Organismic and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| | - Xueyan Li
- Kunming Institute of Zoology, Chinese Academy of SciencesKunmingChina
| | - Masaya Yago
- The University Museum, The University of TokyoTokyoJapan
| | - James Mallet
- Department of Organismic and Evolutionary Biology, Harvard UniversityCambridgeUnited States
| |
Collapse
|
47
|
Hou Z, Li A. Genomic Differentiation and Demographic Histories of Two Closely Related Salicaceae Species. FRONTIERS IN PLANT SCIENCE 2022; 13:911467. [PMID: 35747877 PMCID: PMC9210983 DOI: 10.3389/fpls.2022.911467] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 04/02/2022] [Accepted: 05/06/2022] [Indexed: 06/15/2023]
Abstract
Populus alba (P. alba) and Populus davidiana (P. davidiana) are important plant species for answering a variety of issues on species evolution due to their wide distribution and ability to adapt to a variety of environments and climates. Even though P. alba and P. davidiana belong to ecologically and economically important forest trees in the Northern Hemisphere, little is known about their genomic landscape and genome divergence during speciation. We re-sequenced 20 and 19 members of P. davidiana and P. alba, respectively, and found that the Dxy value between P. alba and P. davidiana was 0.2658, whereas the F ST values were 0.2988, indicating that the genetic divergence was fairly clear. Populus davidiana and P. alba diverged from the ancestor in the middle Pleistocene, c. 0.80 Ma (95% HPD: 0.79-0.81 Ma). The population sizes of P. davidiana increased ~20,000 years ago after a considerable long-term decline following divergence. However, after differentiation, the effective population size of P. alba expanded slightly before experiencing a long-term bottleneck effect. According to the expectation of allopatric speciation, we found a significant number of genomic differentiation sites in both species' speciation events, and the majority of these genomic differentiation regions can be attributed to neutral evolutionary processes. Nevertheless, the regions with extreme divergence exist in abundance, indicating that natural selection has had an impact. Positive selection can be found in highly differentiated regions, while long-term balancing selection traits can be easily observed in low differentiated regions. According to these findings, climate differences over the Quaternary, as well as variance in linked selection and recombination, all contributed significantly to genomic divergence during allopatric speciation of the two aspens.
Collapse
|
48
|
Li Q, Lindtke D, Rodríguez-Ramírez C, Kakioka R, Takahashi H, Toyoda A, Kitano J, Ehrlich RL, Chang Mell J, Yeaman S. Local Adaptation and the Evolution of Genome Architecture in Threespine Stickleback. Genome Biol Evol 2022; 14:6589818. [PMID: 35594844 PMCID: PMC9178229 DOI: 10.1093/gbe/evac075] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Accepted: 05/16/2022] [Indexed: 12/11/2022] Open
Abstract
Theory predicts that local adaptation should favor the evolution of a concentrated genetic architecture, where the alleles driving adaptive divergence are tightly clustered on chromosomes. Adaptation to marine versus freshwater environments in threespine stickleback has resulted in an architecture that seems consistent with this prediction: divergence among populations is mainly driven by a few genomic regions harboring multiple quantitative trait loci for environmentally adapted traits, as well as candidate genes with well-established phenotypic effects. One theory for the evolution of these "genomic islands" is that rearrangements remodel the genome to bring causal loci into tight proximity, but this has not been studied explicitly. We tested this theory using synteny analysis to identify micro- and macro-rearrangements in the stickleback genome and assess their potential involvement in the evolution of genomic islands. To identify rearrangements, we conducted a de novo assembly of the closely related tubesnout (Aulorhyncus flavidus) genome and compared this to the genomes of threespine stickleback and two other closely related species. We found that small rearrangements, within-chromosome duplications, and lineage-specific genes (LSGs) were enriched around genomic islands, and that all three chromosomes harboring large genomic islands have experienced macro-rearrangements. We also found that duplicates and micro-rearrangements are 9.9× and 2.9× more likely to involve genes differentially expressed between marine and freshwater genotypes. While not conclusive, these results are consistent with the explanation that strong divergent selection on candidate genes drove the recruitment of rearrangements to yield clusters of locally adaptive loci.
Collapse
Affiliation(s)
- Qiushi Li
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Canada T2N 1N4
| | - Dorothea Lindtke
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Canada T2N 1N4
| | - Carlos Rodríguez-Ramírez
- Division of Evolutionary Ecology, Institute of Ecology and Evolution, University of Bern, Bern, Switzerland
| | - Ryo Kakioka
- Tropical Biosphere Research Center, University of the Ryukyus, Nishihara, Nakagami-gun, Okinawa 903-0213, Japan
| | - Hiroshi Takahashi
- National Fisheries University, 2-7-1 Nagata-honmachi, Shimonoseki, Yamaguchi 759-6595, Japan
| | - Atsushi Toyoda
- Comparative Genomics Laboratory, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
| | - Jun Kitano
- Ecological Genetics Laboratory, National Institute of Genetics, Mishima, Shizuoka 411-8540, Japan
| | - Rachel L Ehrlich
- Department of Microbiology & Immunology, Drexel University College of Medicine, Philadelphia 19102, PA, USA
| | - Joshua Chang Mell
- Department of Microbiology & Immunology, Drexel University College of Medicine, Philadelphia 19102, PA, USA
| | - Sam Yeaman
- Department of Biological Sciences, University of Calgary, 2500 University Drive NW, Calgary, Canada T2N 1N4
| |
Collapse
|
49
|
Klimova A, Ruiz Mondragón KY, Molina Freaner F, Aguirre-Planter E, Eguiarte LE. Genomic Analyses of Wild and Cultivated Bacanora Agave (Agave angustifolia var. pacifica) Reveal Inbreeding, Few Signs of Cultivation History and Shallow Population Structure. PLANTS 2022; 11:plants11111426. [PMID: 35684199 PMCID: PMC9183054 DOI: 10.3390/plants11111426] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 04/13/2022] [Revised: 05/04/2022] [Accepted: 05/05/2022] [Indexed: 01/12/2023]
Abstract
Due to the recent increase in demand for agave-based beverages, many wild agave populations have experienced rapid decline and fragmentation, whereas cultivated plants are now managed at monocultural plantations, in some cases involving clonal propagation. We examined the relative effect of migration, genetic drift, natural selection and human activities on the genetic repertoire of Agave angustifolia var. pacifica, an agave used for bacanora (an alcoholic spirit similar to tequila) production in northwestern Mexico. We sampled 34 wild and cultivated sites and used over eleven thousand genome-wide SNPs. We found shallow genetic structure among wild samples, although we detected differentiation between coastal and inland sites. Surprisingly, no differentiation was found between cultivated and wild populations. Moreover, we detected moderate inbreeding (FIS ~ 0.13) and similar levels of genomic diversity in wild and cultivated agaves. Nevertheless, the cultivated plants had almost no private alleles and presented evidence of clonality. The overall low genetic structure in A. angustifolia var. pacifica is apparently the result of high dispersibility promoted by pollinators and the possibility of clonal reproduction. Incipient cultivation history and reliance on wild seeds and plants are probably responsible for the observed patterns of high genetic connectivity and considerable diversity in cultivated samples.
Collapse
Affiliation(s)
- Anastasia Klimova
- Departamento de Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Circuito Exterior s/n Annex to the Botanical Garden, Mexico City 04510, Mexico; (K.Y.R.M.); (E.A.-P.)
- Correspondence: (A.K.); (L.E.E.); Tel.: +052-6121716113 (A.K.); +052-5556229006 (L.E.E.)
| | - Karen Y. Ruiz Mondragón
- Departamento de Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Circuito Exterior s/n Annex to the Botanical Garden, Mexico City 04510, Mexico; (K.Y.R.M.); (E.A.-P.)
| | - Francisco Molina Freaner
- Departamento de Ecología de la Biodiversidad, Instituto de Ecología, Universidad Nacional Autónoma de México Hermosillo, Sonora 83250, Mexico;
| | - Erika Aguirre-Planter
- Departamento de Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Circuito Exterior s/n Annex to the Botanical Garden, Mexico City 04510, Mexico; (K.Y.R.M.); (E.A.-P.)
| | - Luis E. Eguiarte
- Departamento de Ecología Evolutiva, Instituto de Ecología, Universidad Nacional Autónoma de México, Ciudad Universitaria, Circuito Exterior s/n Annex to the Botanical Garden, Mexico City 04510, Mexico; (K.Y.R.M.); (E.A.-P.)
- Correspondence: (A.K.); (L.E.E.); Tel.: +052-6121716113 (A.K.); +052-5556229006 (L.E.E.)
| |
Collapse
|
50
|
Abstract
SignificancePhylogenies are the basis of many ecological and evolutionary studies. However, zokor phylogeny and speciation patterns are heavily debated. This study disentangled the phylogeny and speciation of zokors genomically. Six species of the Eospalax were separated into high-altitude E. baileyi and E. smithi and the rest four low-altitude species by recent Qinghai-Tibet Plateau uplift 3.6 million y ago. E. rothschildi and E. smithi speciated south of the Qinling-Huaihe Line, where refuges were supplied during glaciation. Introgression and incomplete lineage sorting led to the complex phylogeny of zokors. Genomic islands were formed due to ancient polymorphisms and divergence hitchhiking. This study concluded that climatic, geological, and tectonic events shaped the phylogeny and speciation of zokors in China.
Collapse
|