scholarly journals Demographic histories and genome-wide patterns of divergence in incipient species of shorebirds

2019 ◽  
Author(s):  
Xuejing Wang ◽  
Kathryn H. Maher ◽  
Nan Zhang ◽  
Pingjia Que ◽  
Chenqing Zheng ◽  
...  

AbstractUnderstanding how incipient species are maintained with gene flow is a fundamental question in evolutionary biology. Whole genome sequencing of multiple individuals holds great potential to illustrate patterns of genomic differentiation as well as the associated evolutionary histories. Kentish (Charadrius alexandrinus) and the white-faced (C. dealbatus) plovers, which differ in their phenotype, ecology and behaviour, are two incipient species and parapatrically distributed in East Asia. Previous studies show evidence of genetic diversification with gene flow between the two plovers. Under this scenario, it is of great importance to explore the patterns of divergence at the genomic level and to determine whether specific regions are involved in reproductive isolation and local adaptation. Here we present the first population genomic analysis of the two incipient species based on the de novo Kentish plover reference genome and resequenced populations. We show that the two plover lineages are distinct in both nuclear and mitochondrial genomes. Using model-based coalescence analysis, we found that population sizes of Kentish plover increased whereas white-faced plovers declined during the Last Glaciation Period. Moreover, the two plovers diverged allopatrically, with gene flow occurring after secondary contact. This has resulted in low levels of genome-wide differentiation, although we found evidence of a few highly differentiated genomic regions in both the autosomes and the Z-chromosome. This study illustrates that incipient shorebird species with gene flow after secondary contact can exhibit discrete divergence at specific genomic regions and provides basis to further exploration on the genetic basis of relevant phenotypic traits.

2020 ◽  
Author(s):  
Steven M. Van Belleghem ◽  
Jared M. Cole ◽  
Gabriela Montejo-Kovacevich ◽  
Caroline N. Bacquet ◽  
W. Owen McMillan ◽  
...  

AbstractCharacterizing the genetic architecture of species boundaries remains a difficult task. Hybridizing species provide a powerful system to identify the factors that shape genomic variation and, ultimately, identify the regions of the genome that maintain species boundaries. Unfortunately, complex histories of isolation, admixture and selection can generate heterogenous genomic landscapes of divergence which make inferences about the regions that are responsible for species boundaries problematic. However, as the signal of admixture and selection on genomic loci varies with recombination rate, their relationship can be used to infer their relative importance during speciation. Here, we explore patterns of genomic divergence, admixture and recombination rate among hybridizing lineages across the Heliconius erato radiation. We focus on the incipient species, H. erato and H. himera, and distinguish the processes that drive genomic divergence across three contact zones where they frequently hybridize. Using demographic modeling and simulations, we infer that periods of isolation and selection have been major causes of genome-wide correlation patterns between recombination rate and divergence between these incipient species. Upon secondary contact, we found surprisingly highly asymmetrical introgression between the species pair, with a paucity of H. erato alleles introgressing into the H. himera genomes. We suggest that this signal may result from a current polygenic species boundary between the hybridizing lineages. These results contribute to a growing appreciation for the importance of polygenic architectures of species boundaries and pervasive genome-wide selection during the early stages of speciation with gene flow.


Animals ◽  
2020 ◽  
Vol 10 (3) ◽  
pp. 493
Author(s):  
Salvatore Mastrangelo ◽  
Filippo Cendron ◽  
Gianluca Sottile ◽  
Giovanni Niero ◽  
Baldassare Portolano ◽  
...  

Through the development of the high-throughput genotyping arrays, molecular markers and genes related to phenotypic traits have been identified in livestock species. In poultry, plumage color is an important qualitative trait that can be used as phenotypic marker for breed identification. In order to assess sources of genetic variation related to the Polverara chicken breed plumage colour (black vs. white), we carried out a genome-wide association study (GWAS) and a genome-wide fixation index (FST) scan to uncover the genomic regions involved. A total of 37 animals (17 white and 20 black) were genotyped with the Affymetrix 600 K Chicken single nucleotide polymorphism (SNP) Array. The combination of results from GWAS and FST revealed a total of 40 significant markers distributed on GGA 01, 03, 08, 12 and 21, and located within or near known genes. In addition to the well-known TYR, other candidate genes have been identified in this study, such as GRM5, RAB38 and NOTCH2. All these genes could explain the difference between the two Polverara breeds. Therefore, this study provides the basis for further investigation of the genetic mechanisms involved in plumage color in chicken.


2021 ◽  
Author(s):  
Jason Bertram

Resolving the role of natural selection is a basic objective of evolutionary biology. It is generally difficult to detect the influence of selection because ubiquitous non-selective stochastic change in allele frequencies (genetic drift) degrades evidence of selection. As a result, selection scans typically only identify genomic regions that have undergone episodes of intense selection. Yet it seems likely such episodes are the exception; the norm is more likely to involve subtle, concurrent selective changes at a large number of loci. We develop a new theoretical approach that uncovers a previously undocumented genome-wide signature of selection in the collective divergence of allele frequencies over time. Applying our approach to temporally-resolved allele frequency measurements from laboratory and wild Drosophila populations, we quantify the selective contribution to allele frequency divergence and find that selection has substantial effects on much of the genome. We further quantify the magnitude of the total selection coefficient (a measure of the combined effects of direct and linked selection) at a typical polymorphic locus, and find this to be large (of order 1%) even though most mutations are not directly under selection. We find that selective allele frequency divergence is substantial at intermediate allele frequencies, which we argue is most parsimoniously explained by positive --- not purifying --- selection. Thus, in these populations most mutations are far from evolving neutrally in the short term (tens of generations), including mutations with neutral fitness effects, and the result cannot be explained simply as a purging of deleterious mutations.


2017 ◽  
Vol 37 (11) ◽  
Author(s):  
Hua Chen ◽  
Kassa Semagn ◽  
Muhammad Iqbal ◽  
Neshat Pazooki Moakhar ◽  
Teketel Haile ◽  
...  

2015 ◽  
Author(s):  
Rob W Ness ◽  
Andrew D Morgan ◽  
Radhakrishnan B Vasanthakrishnan ◽  
Nick Colegrave ◽  
Peter D Keightley

Describing the process of spontaneous mutation is fundamental for understanding the genetic basis of disease, the threat posed by declining population size in conservation biology, and in much evolutionary biology. However, directly studying spontaneous mutation is difficult because of the rarity of de novo mutations. Mutation accumulation (MA) experiments overcome this by allowing mutations to build up over many generations in the near absence of natural selection. In this study, we sequenced the genomes of 85 MA lines derived from six genetically diverse wild strains of the green algaChlamydomonas reinhardtii. We identified 6,843 spontaneous mutations, more than any other study of spontaneous mutation. We observed seven-fold variation in the mutation rate among strains and that mutator genotypes arose, increasing the mutation rate dramatically in some replicates. We also found evidence for fine-scale heterogeneity in the mutation rate, driven largely by the sequence flanking mutated sites, and by clusters of multiple mutations at closely linked sites. There was little evidence, however, for mutation rate heterogeneity between chromosomes or over large genomic regions of 200Kbp. Using logistic regression, we generated a predictive model of the mutability of sites based on their genomic properties, including local GC content, gene expression level and local sequence context. Our model accurately predicted the average mutation rate and natural levels of genetic diversity of sites across the genome. Notably, trinucleotides vary 17-fold in rate between the most mutable and least mutable sites. Our results uncover a rich heterogeneity in the process of spontaneous mutation both among individuals and across the genome.


2021 ◽  
Author(s):  
Tyler K Chafin ◽  
Binod Regmi ◽  
Marlis R. Douglas ◽  
David R. Edds ◽  
Karma Wangchuk ◽  
...  

Replicated evolutionary patterns are often attributed to recurrent emergence following parallel selective pressures. However, similar genetic patterns (e.g., 'genomic islands') can also emerge following extensive homogenization in secondary contact, as a by-product of heterogeneous introgression. For example, within Himalayan tributaries of the Ganges/Brahmaputra rivers, drainage-specific mtDNA clades of polyploid snowtrout (Cyprinidae: Schizothorax) are partitioned as co-occurring morphological 'ecotypes,' hypothesized to represent parallel divergence among adjacent streams. To evaluate this scenario, we utilized a reduced-representation genomic approach (N=35,319 de-novo and N=10,884 transcriptome-aligned SNPs) applied to high-altitude Nepali/Bhutanese snowtrout (N=48 each). We unambiguously quantified ploidy levels by first deriving genome-wide allelic depths followed by ploidy-aware Bayesian models that produced genotypes statistically consistent with diploid/tetraploid expectations. When genotyped SNPs were clustering within drainages, the convergence of eco-phenotypes was sustained. However, subsequent partitioned analyses of phylogeny and population admixture instead identified subsets of loci under selection which retained genealogical concordance with morphology, with apparent patterns of parallel ecotype emergence instead driven by widespread genomic homogenization. Here, prior isolation is effectively masked by admixture occurring in secondary contact. We note two salient factors:1) Polyploidy has promoted homogenization in tetraploid Himalayan snowtrout; and 2) Homogenization varies across Himalayan tributaries, presumably in lockstep with extent of anthropogenic modification.


PeerJ ◽  
2016 ◽  
Vol 4 ◽  
pp. e1910 ◽  
Author(s):  
Quentin Rougemont ◽  
Camille Roux ◽  
Samuel Neuenschwander ◽  
Jerome Goudet ◽  
Sophie Launey ◽  
...  

Inferring the history of isolation and gene flow during species divergence is a central question in evolutionary biology. The European river lamprey (Lampetra fluviatilis) and brook lamprey(L. planeri)show a low reproductive isolation but have highly distinct life histories, the former being parasitic-anadromous and the latter non-parasitic and freshwater resident. Here we used microsatellite data from six replicated population pairs to reconstruct their history of divergence using an approximate Bayesian computation framework combined with a random forest model. In most population pairs, scenarios of divergence with recent isolation were outcompeted by scenarios proposing ongoing gene flow, namely the Secondary Contact (SC) and Isolation with Migration (IM) models. The estimation of demographic parameters under the SC model indicated a time of secondary contact close to the time of speciation, explaining why SC and IM models could not be discriminated. In case of an ancient secondary contact, the historical signal of divergence is lost and neutral markers converge to the same equilibrium as under the less parameterized model allowing ongoing gene flow. Our results imply that models of secondary contacts should be systematically compared to models of divergence with gene flow; given the difficulty to discriminate among these models, we suggest that genome-wide data are needed to adequately reconstruct divergence history.


2018 ◽  
Author(s):  
Krzysztof M. Kozak ◽  
W. Owen McMillan ◽  
Mathieu Joron ◽  
Christopher D. Jiggins

ABSTRACTHow frequent is gene flow between species? The pattern of evolution is typically portrayed as a phylogenetic tree, implying that speciation is a series of splits between lineages. Yet gene flow between good species is increasingly recognized as an important mechanism in the diversification of radiations, often spreading adaptive traits and leading to a complex pattern of phylogenetic incongruence. This process has thus far been studied in cases involving few species, or geographically restricted to spaces like islands, but not on the scale of a continental radiation. Previous studies have documented gene flow, adaptive introgression and hybrid speciation in a small subsection of the charismatic Neotropical butterflies Heliconius. Using genome-wide resequencing of 40 out of 45 species in the genus we demonstrate for the first time that admixture has played a role throughout the evolution of Heliconius and the sister genus Eueides. Modelling of phylogenetic networks based on 6848 orthologous autosomal genes (Maximum Pseudo-Likelihood Networks) or 5,483,419 high quality SNPs (Ancestral Recombination Graph) uncovers nine new cases of interspecific gene flow at up to half of the genome. However, f4 statistics of admixture show that the extent of the process has varied between subgenera. Evidence for introgression is found at all five loci controlling the colour and shape of the mimetic wing patterns, including in the putative hybrid species H. hecalesia, characterised by an unusual hindwing. Due to hybridization and incomplete coalescence during rapid speciation, individual gene trees show rampant discordance. Although reduced gene flow and faster coalescence are expected at the Z chromosome, we discover high levels of conflict between the 416 sex-linked loci. Despite this discordant pattern, both concatenation and multispecies coalescent approaches yield surprisingly consistent and fully supported genome-wide phylogenies. We conclude that the imposition of the bifurcating tree model without testing for interspecific gene flow may distort our perception of adaptive radiations and thus the ability to study trait evolution in a comparative framework.


Author(s):  
Fanny E. Hartmann ◽  
Tiziana Vonlanthen ◽  
Nikhil Kumar Singh ◽  
Megan McDonald ◽  
Andrew Milgate ◽  
...  

AbstractConvergent evolution leads to identical phenotypic traits in different species or populations. Convergence can be driven by standing variation allowing selection to favor identical alleles in parallel or the same mutations can arise independently. However, the molecular basis of such convergent adaptation remains often poorly resolved. Pesticide resistance in agricultural ecosystems is a hallmark of convergence in phenotypic traits. Here, we analyze the major fungal pathogen Zymoseptoria tritici causing serious losses on wheat and with parallel fungicide resistance emergence across continents. We sampled three population pairs each from a different continent spanning periods early and late in the application of fungicides. To identify causal loci for resistance, we combined knowledge from molecular genetics work and performed genome-wide association studies (GWAS) on a global set of isolates. We discovered yet unknown factors in azole resistance including membrane stability functions. We found strong support for the ‘hotspot’ model of resistance evolution with parallel changes in a small set of loci but additional loci showed more population-specific allele frequency changes. Genome-wide scans of selection showed that half of all known resistance loci were overlapping a selective sweep region. Hence, the application of fungicides was one of the major selective agents acting on the pathogen over the past decades. Furthermore, loci identified through GWAS showed the highest overlap with selective sweep regions underlining the importance to map phenotypic trait variation in evolving populations. Our population genomic analyses showed that both de novo mutations and gene flow likely contributed to the parallel emergence of resistance.


Author(s):  
Matthew L Bendall ◽  
Keylie M Gibson ◽  
Margaret C Steiner ◽  
Uzma Rentia ◽  
Marcos Pérez-Losada ◽  
...  

Abstract Deep sequencing of viral populations using next generation sequencing (NGS) offers opportunities to understand and investigate evolution, transmission dynamics, and population genetics. Currently, the standard practice for processing NGS data to study viral populations is to summarize all the observed sequences from a sample as a single consensus sequence, thus discarding valuable information about the intra-host viral molecular epidemiology. Furthermore, existing analytical pipelines may only analyze genomic regions involved in drug resistance, thus are not suited for full viral genome analysis. Here we present HAPHPIPE, a HAplotype and PHylodynamics PIPEline for genome-wide assembly of viral consensus sequences and haplotypes. The HAPHPIPE protocol includes modules for quality trimming, error correction, de novo assembly, alignment, and haplotype reconstruction. The resulting consensus sequences, haplotypes, and alignments can be further analyzed using a variety of phylogenetic and population genetic software. HAPHPIPE is designed to provide users with a single pipeline to rapidly analyze sequences from viral populations generated from NGS platforms and provide quality output properly formatted for downstream evolutionary analyses.


Sign in / Sign up

Export Citation Format

Share Document