scholarly journals Bioindicator snake shows genomic signatures of natural and anthropogenic barriers to gene flow

PLoS ONE ◽  
2021 ◽  
Vol 16 (10) ◽  
pp. e0259124
Author(s):  
Damian C. Lettoof ◽  
Vicki A. Thomson ◽  
Jari Cornelis ◽  
Philip W. Bateman ◽  
Fabien Aubret ◽  
...  

Urbanisation alters landscapes, introduces wildlife to novel stressors, and fragments habitats into remnant ‘islands’. Within these islands, isolated wildlife populations can experience genetic drift and subsequently suffer from inbreeding depression and reduced adaptive potential. The Western tiger snake (Notechis scutatus occidentalis) is a predator of wetlands in the Swan Coastal Plain, a unique bioregion that has suffered substantial degradation through the development of the city of Perth, Western Australia. Within the urban matrix, tiger snakes now only persist in a handful of wetlands where they are known to bioaccumulate a suite of contaminants, and have recently been suggested as a relevant bioindicator of ecosystem health. Here, we used genome-wide single nucleotide polymorphism (SNP) data to explore the contemporary population genomics of seven tiger snake populations across the urban matrix. Specifically, we used population genomic structure and diversity, effective population sizes (Ne), and heterozygosity-fitness correlations to assess fitness of each population with respect to urbanisation. We found that population genomic structure was strongest across the northern and southern sides of a major river system, with the northern cluster of populations exhibiting lower heterozygosities than the southern cluster, likely due to a lack of historical gene flow. We also observed an increasing signal of inbreeding and genetic drift with increasing geographic isolation due to urbanisation. Effective population sizes (Ne) at most sites were small (< 100), with Ne appearing to reflect the area of available habitat rather than the degree of adjacent urbanisation. This suggests that ecosystem management and restoration may be the best method to buffer the further loss of genetic diversity in urban wetlands. If tiger snake populations continue to decline in urban areas, our results provide a baseline measure of genomic diversity, as well as highlighting which ‘islands’ of habitat are most in need of management and protection.

2021 ◽  
Author(s):  
Zachary L Nikolakis ◽  
Richard Orton ◽  
Brian I Crother

Understanding the processes and mechanisms that promote lineage divergence is a central goal in evolutionary biology. For instance, studies investigating the spatial distribution of genomic variation often highlight biogeographic barriers underpinning geographic isolation, as well as patterns of isolation by environment and isolation by distance that can also lead to lineage divergence. However, the patterns and processes that shape genomic variation and drive lineage divergence may be taxa-specific, even across closely related taxa co-occurring within the same biogeographic region. Here, we use molecular data in the form of ultra-conserved elements (UCEs) to infer the evolutionary relationships and population genomic structure of the Eastern Pinesnake complex (Pituophis melanoleucus) – a polytypic wide-ranging species that occupies much of the Eastern Nearctic. In addition to inferring evolutionary relationships, population genomic structure, and gene flow, we also test relationships between genomic diversity and putative barriers to dispersal, environmental variation, and geographic distance. We present results that reveal shallow population genomic structure and ongoing gene flow, despite an extensive geographic range that transcends geographic features found to reduce gene flow among many taxa, including other squamate reptiles within the Eastern Nearctic. Further, our results indicate that the observed genomic diversity is spatially distributed as a pattern of isolation by distance and suggest that the current subspecific taxonomy do not adhere to independent lineages, but rather, show a significant amount of admixture across the entire P. melanoleucus range.


The Auk ◽  
2019 ◽  
Vol 136 (2) ◽  
Author(s):  
Joseph D Manthey ◽  
Stéphane Boissinot ◽  
Robert G Moyle

Abstract Evolutionary biologists have long used behavioral, ecological, and genetic data from contact zones between closely related species to study various phases of the speciation continuum. North America has several concentrations of avian contact zones, where multiple pairs of sister lineages meet, with or without hybridization. In a southern California contact zone, 2 species of woodpeckers, Nuttall’s Woodpecker (Dryobates nuttallii) and the Ladder-backed Woodpecker (D. scalaris), occasionally hybridize. We sampled these 2 species in a transect across this contact zone and included samples of their closest relative, the Downy Woodpecker (D. pubescens), to obtain large single nucleotide polymorphism panels using restriction-site associated DNA sequencing (RAD-seq). Furthermore, we used whole-genome resequencing data for 2 individuals per species to identify whether patterns of diversity inferred from RAD-seq were representative of whole-genome diversity. We found that these 3 woodpecker species are genomically distinct. Although low levels of gene flow occur between D. nuttallii and D. scalaris across the contact zone, there was no evidence for widespread genomic introgression between these 2 species. Overall patterns of genomic diversity from the RAD-seq and wholegenome datasets appear to be related to distributional range size and, by extension, are likely related to effective population sizes for each species.


2010 ◽  
Vol 365 (1543) ◽  
pp. 1127-1138 ◽  
Author(s):  
Paquita E. A. Hoeck ◽  
Jennifer L. Bollmer ◽  
Patricia G. Parker ◽  
Lukas F. Keller

Small and isolated island populations provide ideal systems to study the effects of limited population size, genetic drift and gene flow on genetic diversity. We assessed genetic diversity within and differentiation among 19 mockingbird populations on 15 Galápagos islands, covering all four endemic species, using 16 microsatellite loci. We tested for signs of drift and gene flow, and used historic specimens to assess genetic change over the last century and to estimate effective population sizes. Within-population genetic diversity and effective population sizes varied substantially among island populations and correlated strongly with island size, suggesting that island size serves as a good predictor for effective population size. Genetic differentiation among populations was pronounced and increased with geographical distance. A century of genetic drift did not change genetic diversity on an archipelago-wide scale, but genetic drift led to loss of genetic diversity in small populations, especially in one of the two remaining populations of the endangered Floreana mockingbird. Unlike in other Galápagos bird species such as the Darwin's finches, gene flow among mockingbird populations was low. The clear pattern of genetically distinct populations reflects the effects of genetic drift and suggests that Galápagos mockingbirds are evolving in relative isolation.


2019 ◽  
Author(s):  
Arun Sethuraman ◽  
Melissa Lynch

AbstractUnsampled or extinct ‘ghost’ populations leave signatures on the genomes of individuals from extant, sampled populations, especially if they have exchanged genes with them over evolutionary time. This gene flow from ‘ghost’ populations can introduce biases when estimating evolutionary history from genomic data, often leading to data misinterpretation and ambiguous results. Here we assess these biases while accounting, or not accounting for gene flow from ‘ghost’ populations under the Isolation with Migration (IM) model. We perform extensive simulations under five scenarios with no gene flow (Scenario A), to extensive gene flow to- and from- an unsampled ‘ghost’ population (Scenarios B, C, D, and E). Estimates of evolutionary history across all scenarios A-E (effective population sizes, divergence times, and migration rates) indicate consistent a) under-estimation of divergence times between sampled populations, (b) over-estimation of effective population sizes of sampled populations, and (c) under-estimation of migration rates between sampled populations, with increased gene flow from the unsampled ‘ghost’ population. Without accounting for an unsampled ‘ghost’, summary statistics like FST are under-estimated, and π is over-estimated with increased gene flow from the‘ghost’. To show this persistent issue in empirical data, we use a 355 locus dataset from African Hunter-Gatherer populations and discuss similar biases in estimating evolutionary history while not accounting for unsampled ‘ghosts’. Considering the large effects of gene flow from these ‘ghosts’, we propose a multi-pronged approach to account for the presence of unsampled ‘ghost’ populations in population genomics studies to reduce erroneous inferences.


2019 ◽  
Author(s):  
Olga M. Pérez-Carrascal ◽  
Yves Terrat ◽  
Alessandra Giani ◽  
Nathalie Fortin ◽  
Charles W. Greer ◽  
...  

AbstractMicrocystis is a genus of freshwater cyanobacteria which causes harmful blooms in ecosystems worldwide. Some Microcystis strains produce harmful toxins such as microcystin, impacting drinking water quality. Microcystis colony morphology, rather than genetic similarity, is often used to classify Microcystis into morphospecies. However, colony morphology is a plastic trait which can change depending on environmental and laboratory culture conditions, and is thus an inadequate criterion for species delineation. Furthermore, Microcystis populations are thought to disperse globally and constitute a homogeneous gene pool. However, this assertion is based on relatively incomplete characterization of Microcystis genomic diversity. To better understand these issues, we performed a population genomic analysis of 33 newly sequenced genomes (of which 19 were resequenced to check for mutation in culture) mainly from Canada and Brazil. We identified eight Microcystis clusters of genomic similarity, only four of which correspond to named morphospecies and monophyletic groups. Notably, M. aeruginosa is paraphyletic, distributed across four genomic clusters, suggesting it is not a coherent species. Most monophyletic groups are specific to a unique geographic location, suggesting biogeographic structure over relatively short evolutionary time scales. Higher homologous recombination rates within than between clusters further suggest that monophyletic groups might adhere to a Biological Species-like concept, in which barriers to gene flow maintain species distinctness. However, certain genes – including some involved in microcystin and micropeptin biosynthesis – are recombined between monophyletic groups in the same geographic location, suggesting local adaptation. Together, our results show the importance of using genomic criteria for Microcystis species delimitation and suggest the existence of locally adapted lineages and genes.ImportanceThe genus Microcystis is responsible for harmful and often toxic cyanobacterial blooms across the world, yet it is unclear how and if the genus should be divided into ecologically and genomically distinct species. To resolve the controversy and uncertainty surrounding Microcystis species, we performed a population genomic analysis of Microcystis genome from public databases, along with new isolates from Canada and Brazil. We inferred that significant genetic substructure exists within Microcystis, with several species being maintained by barriers to gene flow. Thus, Microcystis appears to be among a growing number of bacteria that adhere to a Biological Species-like Concept (BSC). Barriers to gene flow are permeable, however, and we find evidence for relatively frequent cross-species horizontal gene transfer (HGT) of genes that may be involved in local adaptation. Distinct clades of Microcystis (putative species) tend to have distinct profiles of toxin biosynthesis genes, and yet toxin genes are also subject to cross-species HGT and local adaptation. Our results thus pave the way for more informed classification, monitoring and understanding of harmful Microcystis blooms.


2020 ◽  
Author(s):  
Thomas L Schmidt ◽  
T. Swan ◽  
Jessica Chung ◽  
Stephan Karl ◽  
Samuel Demok ◽  
...  

AbstractPopulation genomic approaches can characterise dispersal across a single generation through to many generations in the past, bridging the gap between individual movement and intergenerational gene flow. These approaches are particularly useful when investigating dispersal in recently altered systems, where they provide a way of inferring long-distance dispersal between newly established populations and their interactions with existing populations. Human-mediated biological invasions represent such altered systems which can be investigated with appropriate study designs and analyses. Here we apply temporally-restricted sampling and a range of population genomic approaches to investigate dispersal in a 2004 invasion of Aedes albopictus (the Asian tiger mosquito) in the Torres Strait Islands (TSI) of Australia. We sampled mosquitoes from 13 TSI villages simultaneously and genotyped 373 mosquitoes at genome-wide single nucleotide polymorphisms (SNPs): 331 from the TSI, 36 from Papua New Guinea (PNG), and 4 incursive mosquitoes detected in uninvaded regions. Within villages, spatial genetic structure varied substantially but overall displayed isolation by distance and a neighbourhood size of 232–577. Close kin dyads revealed recent movement between islands 31–203 km apart, and deep learning inferences showed incursive Ae. albopictus had travelled to uninvaded regions from both adjacent and non-adjacent islands. Private alleles and a coancestry matrix indicated direct gene flow from PNG into nearby islands. Outlier analyses also detected four linked alleles introgressed from PNG, with the alleles surrounding 12 resistance-associated cytochrome P450 genes. By treating dispersal as both an intergenerational process and a set of discrete events, we describe a highly interconnected invasive system.


2014 ◽  
Author(s):  
Jonathan Puritz ◽  
Christopher M. Hollenbeck ◽  
John R. Gold

Restriction-site associated DNA sequencing (RADseq) has become a powerful and useful approach for population genomics. Currently, no software exists that utilizes both paired-end reads from RADseq data to efficiently produce population-informative variant calls, especially for organisms with large effective population sizes and high levels of genetic polymorphism but for which no genomic resources exist. dDocent is an analysis pipeline with a user-friendly, command-line interface designed to process individually barcoded RADseq data (with double cut sites) into informative SNPs/Indels for population-level analyses. The pipeline, written in BASH, uses data reduction techniques and other stand-alone software packages to perform quality trimming and adapter removal, de novo assembly of RAD loci, read mapping, SNP and Indel calling, and baseline data filtering. Double-digest RAD data from population pairings of three different marine fishes were used to compare dDocent with Stacks, the first generally available, widely used pipeline for analysis of RADseq data. dDocent consistently identified more SNPs shared across greater numbers of individuals and with higher levels of coverage. This is most likely due to the fact that dDocent quality trims instead of filtering and incorporates both forward and reverse reads in assembly, mapping, and SNP calling, thus enabling use of reads with Indel polymorphisms. The pipeline and a comprehensive user guide can be found at (http://dDocent.wordpress.com).


2019 ◽  
Vol 110 (5) ◽  
pp. 587-600
Author(s):  
A Millie Burrell ◽  
Jeffrey H R Goddard ◽  
Paul J Greer ◽  
Ryan J Williams ◽  
Alan E Pepper

Abstract Globally, a small number of plants have adapted to terrestrial outcroppings of serpentine geology, which are characterized by soils with low levels of essential mineral nutrients (N, P, K, Ca, Mo) and toxic levels of heavy metals (Ni, Cr, Co). Paradoxically, many of these plants are restricted to this harsh environment. Caulanthus ampexlicaulis var. barbarae (Brassicaceae) is a rare annual plant that is strictly endemic to a small set of isolated serpentine outcrops in the coastal mountains of central California. The goals of the work presented here were to 1) determine the patterns of genetic connectivity among all known populations of C. ampexlicaulis var. barbarae, and 2) estimate contemporary effective population sizes (Ne), to inform ongoing genomic analyses of the evolutionary history of this taxon, and to provide a foundation upon which to model its future evolutionary potential and long-term viability in a changing environment. Eleven populations of this taxon were sampled, and population-genetic parameters were estimated using 11 nuclear microsatellite markers. Contemporary effective population sizes were estimated using multiple methods and found to be strikingly small (typically Ne &lt; 10). Further, our data showed that a substantial component of genetic connectivity of this taxon is not at equilibrium, and instead showed sporadic gene flow. Several lines of evidence indicate that gene flow between isolated populations is maintained through long-distance seed dispersal (e.g., &gt;1 km), possibly via zoochory.


2022 ◽  
Vol 13 (1) ◽  
Author(s):  
José Martín Pujolar ◽  
Mozes P. K. Blom ◽  
Andrew Hart Reeve ◽  
Jonathan D. Kennedy ◽  
Petter Zahl Marki ◽  
...  

AbstractTropical mountains harbor exceptional concentrations of Earth’s biodiversity. In topographically complex landscapes, montane species typically inhabit multiple mountainous regions, but are absent in intervening lowland environments. Here we report a comparative analysis of genome-wide DNA polymorphism data for population pairs from eighteen Indo-Pacific bird species from the Moluccan islands of Buru and Seram and from across the island of New Guinea. We test how barrier strength and relative elevational distribution predict population differentiation, rates of historical gene flow, and changes in effective population sizes through time. We find population differentiation to be consistently and positively correlated with barrier strength and a species’ altitudinal floor. Additionally, we find that Pleistocene climate oscillations have had a dramatic influence on the demographics of all species but were most pronounced in regions of smaller geographic area. Surprisingly, even the most divergent taxon pairs at the highest elevations experience gene flow across barriers, implying that dispersal between montane regions is important for the formation of montane assemblages.


2020 ◽  
Author(s):  
Eeva Jansson ◽  
Francois Besnier ◽  
Ketil Malde ◽  
Carl André ◽  
Geir Dahle ◽  
...  

Abstract Background Marine fish populations are often characterized by high levels of gene flow and correspondingly low genetic divergence. This presents a challenge to define management units. Goldsinny wrasse ( Ctenolabrus rupestris ) is a heavily exploited species due to its importance as a cleaner-fish in commercial salmonid aquaculture. However, at the present, the population genetic structure of this species is still largely unresolved. Here, full-genome sequencing was used to produce the first genomic reference for this species, to study population-genomic divergence among four geographically distinct populations, and, to identify informative SNP markers for future studies. Results After construction of a de novo assembly, the genome was estimated to be highly polymorphic and of ~600Mbp in size. 33 235 genome wide SNPs were thereafter selected to assess genomic diversity and differentiation among four populations collected from Scandinavia, Scotland, and Spain. Global F ST among these populations was 0.015–0.092. Approximately 4% of the investigated loci were identified as putative global outliers, and ~1% within Scandinavia. SNPs showing large divergence ( F ST >0.15) were picked as candidate diagnostic markers for population assignment. 173 of the most diagnostic SNPs between the two Scandinavian populations were validated by genotyping 47 individuals from each end of the species’ Scandinavian distribution range. 69 of these SNPs were significantly ( p <0.05) differentiated (mean F ST_173_loci = 0.065, F ST_69_loci = 0.140). Using these validated SNPs, individuals were assigned with high probability (≥ 94%) to their populations of origin. Conclusions Goldsinny wrasse displays a highly polymorphic genome, and substantial population genomic structure. Diversifying selection likely affects population structuring globally and within Scandinavia. The diagnostic loci identified now provide a promising and cost-efficient tool to investigate goldsinny wrasse populations further.


Sign in / Sign up

Export Citation Format

Share Document