A phylogenomic resolution for the taxonomy of Aegean green lizards

Lacerta pamphylica and Lacerta trilineata are two currently recognized green lizard species with a historically problematic taxonomy. In cases of tangled phylogenies, next-generation sequencing and double-digest restriction-site-associated DNA protocols can provide a wealth of genomic data and resolve difficult taxonomic issues. Here, we generated genome-wide SNPs and mitochondrial sequences, and applied molecular species delimitation approaches to provide a stable taxonomy for the Aegean green lizards. Mitochondrial gene trees, genetic cluster delimitation and population structure analyses converged into recognizing the populations of (a) L. pamphylica, (b) east Aegean islands, Anatolia and Thrace (diplochondrodes lineage), (c) central Aegean islands (citrovittata), and (d) remaining Balkan populations and islands (trilineata), as separate clusters. Phylogenomic analyses revealed a split into two major clades, east and west of the Aegean Barrier, unambiguously showing a sister–clade relationship between pamphylica and diplochondrodes, rendering L. trilineata paraphyletic. Species delimitation models were tested in a Bayesian framework using the genomic SNPs: lumping all populations into a single ‘species’ had the lowest likelihood but the current taxonomy was also outperformed by all other models. All lines of evidence support the Pamphylian green lizard as a valid species; thus, east Aegean L. trilineata should also be considered a distinct species under the name Lacerta diplochondrodes. Finally, evidence from the mitochondrial and nuclear genomes is overwhelmingly in favour of recognizing the morphologically distinct Cycladian green lizards as a distinct species. We propose their elevation to full species under the name Lacerta citrovittata. All remaining insular and continental populations of the Balkan Peninsula represent the species L. trilineata.     

Metapopulation dynamics and foraging plasticity in a highly vagile seabird,the southern rockhopper penguin

Population connectivity is driven by individual dispersal potential and modulated by natal philopatry. In seabirds, high vagility facilitates dispersal yet philopatry is also common, with foraging area overlap often correlated with population connectivity. We assess the interplay between these processes by studying past and current connectivity and foraging niche overlap among southern rockhopper penguin colonies of the coast of southern South America using genomic and stable isotope analyses. We found two distinct genetic clusters and detected low admixture between northern and southern colonies. Stable isotope analysis indicated niche variability between colonies, with Malvinas/Falklands colonies encompassing the species entire isotopic foraging niche, while the remaining colonies had smaller, nonoverlapping niches. A recently founded colony in continental Patagonia differed in isotopic niche width and position with Malvinas/Falklands colonies, its genetically identified founder population, suggesting the exploitation of novel foraging areas and/or prey items. Additionally, dispersing individuals found dead across the Patagonian shore in an unusual mortality event were also assigned to the northern cluster, suggesting northern individuals reach southern localities, but do not breed in these colonies. Facilitated by variability in foraging strategies, and especially during unfavorable conditions, the number of dispersing individuals may increase and enhance the probability of founding new colonies. Metapopulation demographic dynamics in seabirds should account for interannual variability in dispersal behavior and pay special attention to extreme climatic events, classically related to negative effects on population trends.     

Canadian polar bear population structure using genome‐wide markers

Predicting the consequences of environmental changes, including human‐mediated climate change on species, requires that we quantify range‐wide patterns of genetic diversity and identify the ecological, environmental, and historical factors that have contributed to it. Here, we generate baseline data on polar bear population structure across most Canadian subpopulations (n = 358) using 13,488 genome‐wide single nucleotide polymorphisms (SNPs) identified with double‐digest restriction site‐associated DNA sequencing (ddRAD). Our ddRAD dataset showed three genetic clusters in the sampled Canadian range, congruent with previous studies based on microsatellites across the same regions; however, due to a lack of sampling in Norwegian Bay, we were unable to confirm the existence of a unique cluster in that subpopulation. These data on the genetic structure of polar bears using SNPs provide a detailed baseline against which future shifts in population structure can be assessed, and opportunities to develop new noninvasive tools for monitoring polar bears across their range.     

Spatiotemporally explicit demographic modelling supports a joint effect of historical barriers to dispersal and contemporary landscape composition on structuring genomic variation in a red‐listed grasshopper

Inferring the processes underlying spatial patterns of genomic variation is fundamental to understand how organisms interact with landscape heterogeneity and to identify the factors determining species distributional shifts. Here, we use genomic data (restriction site‐associated DNA sequencing) to test biologically informed models representing historical and contemporary demographic scenarios of population connectivity for the Iberian cross‐backed grasshopper Dociostaurus hispanicus, a species with a narrow distribution that currently forms highly fragmented populations. All models incorporated biological aspects of the focal taxon that could hypothetically impact its geographical patterns of genomic variation, including (a) spatial configuration of impassable barriers to dispersal defined by topographic landscapes not occupied by the species; (b) distributional shifts resulting from the interaction between the species bioclimatic envelope and Pleistocene glacial cycles; and (c) contemporary distribution of suitable habitats after extensive land clearing for agriculture. Spatiotemporally explicit simulations under different scenarios considering these aspects and statistical evaluation of competing models within an Approximate Bayesian Computation framework supported spatial configuration of topographic barriers to dispersal and human‐driven habitat fragmentation as the main factors explaining the geographical distribution of genomic variation in the species, with no apparent impact of hypothetical distributional shifts linked to Pleistocene climatic oscillations. Collectively, this study supports that both historical (i.e., topographic barriers) and contemporary (i.e., anthropogenic habitat fragmentation) aspects of landscape composition have shaped major axes of genomic variation in the studied species and emphasizes the potential of model‐based approaches to gain insights into the temporal scale at which different processes impact the demography of natural populations.     

Population genomic evidence for plant glacial survival in Scandinavia

Quaternary glaciations have played a major role in shaping the genetic diversity and distribution of plant species. Strong palaeoecological and genetic evidence supports a postglacial recolonization of most plant species to northern Europe from southern, eastern and even western glacial refugia. Although highly controversial, the existence of small in situ glacial refugia in northern Europe has recently gained molecular support. We used genomic analyses to examine the phylogeography of a species that is critical in this debate. Carex scirpoidea Michx subsp. scirpoidea is a dioecious, amphi‐Atlantic arctic–alpine sedge that is widely distributed in North America, but absent from most of Eurasia, apart from three extremely disjunct populations in Norway, all well within the limits of the Weichselian ice sheet. Range‐wide population sampling and variation at 5,307 single nucleotide polymorphisms show that the three Norwegian populations comprise unique evolutionary lineages divergent from Greenland with high between‐population divergence. The Norwegian populations have low within‐population genetic diversity consistent with having experienced genetic bottlenecks in glacial refugia, and host private alleles that probably accumulated in long‐term isolated populations. Demographic analyses support a single, pre‐Weichselian colonization into Norway from East Greenland, and subsequent divergence of the three populations in separate refugia. Other refugial areas are identified in North‐east Greenland, Minnesota/Michigan, Colorado and Alaska. Admixed populations in British Columbia and West Greenland indicate postglacial contact. Taken together, evidence from this study strongly indicates in situ glacial survival in Scandinavia.     

The population genomics of multiple tsetse fly (Glossina fuscipes fuscipes) admixture zones in Uganda

Understanding the mechanisms that enforce, maintain or reverse the process of speciation is an important challenge in evolutionary biology. This study investigates the patterns of divergence and discusses the processes that form and maintain divergent lineages of the tsetse fly Glossina fuscipes fuscipes in Uganda. We sampled 251 flies from 18 sites spanning known genetic lineages and the four admixture zones between them. We apply population genomics, hybrid zone and approximate Bayesian computation to the analysis of three types of genetic markers: 55,267 double‐digest restriction site‐associated DNA (ddRAD) SNPs to assess genome‐wide admixture, 16 microsatellites to provide continuity with published data and accurate biogeographic modelling, and a 491‐bp fragment of mitochondrial cytochrome oxidase I and II to infer maternal inheritance patterns. Admixture zones correspond with regions impacted by the reorganization of Uganda's river networks that occurred during the formation of the West African Rift system over the last several hundred thousand years. Because tsetse fly population distributions are defined by rivers, admixture zones likely represent both old and new regions of secondary contact. Our results indicate that older hybrid zones contain mostly parental types, while younger zones contain variable hybrid types resulting from multiple generations of interbreeding. These findings suggest that reproductive barriers are nearly complete in the older admixture zones, while nearly absent in the younger admixture zones. Findings are consistent with predictions of hybrid zone theory: Populations in zones of secondary contact transition rapidly from early to late stages of speciation or collapse all together.     

Genomic pedigree reconstruction identifies predictors of mating and reproductive success in an invasive vertebrate

The persistence of an invasive species is influenced by its reproductive ecology, and a successful control program must operate on this premise. However, the reproductive ecology of invasive species may be enigmatic due to factors that also limit their management, such as cryptic coloration and behavior. We explored the mating and reproductive ecology of the invasive Brown Treesnake (BTS: Boiga irregularis) by reconstructing a multigenerational genomic pedigree based on 654 single nucleotide polymorphisms for a geographically closed population established in 2004 on Guam (N = 426). The pedigree allowed annual estimates of individual mating and reproductive success to be inferred for snakes in the study population over a 14‐year period. We then employed generalized linear mixed models to gauge how well phenotypic and genomic data could predict sex‐specific annual mating and reproductive success. Average snout–vent length (SVL), average body condition index (BCI), and trappability were significantly related to annual mating success for males, with average SVL also related to annual mating success for females. Male and female annual reproductive success was positively affected by SVL, BCI, and trappability. Surprisingly, the degree to which individuals were inbred had no effect on annual mating or reproductive success. When juxtaposed with current control methods, these results indicate that baited traps, a common interdiction tool, may target fecund BTS in some regards but not others. Our study emphasizes the importance of reproductive ecology as a focus for improving BTS control and promotes genomic pedigree reconstruction for such an endeavor in this invasive species and others.     

Combining allele frequency and tree‐based approaches improves phylogeographic inference from natural history collections

Model selection approaches in phylogeography have allowed researchers to evaluate the support for competing demographic histories, which provides a mode of inference and a measure of uncertainty in understanding climatic and spatial influences on intraspecific diversity. Here, to rank all models in the comparison set and determine what proportion of the total support the top‐ranked model garners, we conduct model selection using two analytical approaches—allele frequency‐based, implemented in fastsimcoal2 , and gene tree‐based, implemented in phrapl . We then expand this model selection framework by including an assessment of absolute fit of the models to the data. For this, we utilize DNA isolated from existing natural history collections that span the distribution of red alder (Alnus rubra) in the Pacific Northwest of North America to generate genomic data for the evaluation of 13 demographic scenarios. The quality of DNA recovered from herbarium specimen leaf tissue was assessed for its utility and effectiveness in demographic model selection, specifically in the two approaches mentioned. We present strong support for the use of herbarium tissue in the generation of genomic DNA, albeit with the inclusion of additional quality control checks prior to library preparation and analyses with multiple approaches that incorporate various data. Analyses with allele frequency spectra and gene trees predominantly support A. rubra having experienced an ancient vicariance event with intermittent and frequent gene flow between the disjunct populations. Additionally, the data consistently fit the most frequently selected model, corroborating the model selection techniques. Finally, these results suggest that the A. rubra disjunct populations do not represent separate species.     

Oxygen concentration drives local adaptation in the greenlip abalone (Haliotis laevigata)

Understanding adaptation has become one of the major biological questions especially in the light of rapid environmental changes induced by climate change. Ocean temperatures are rising which triggers massive changes in water chemistry and thereby alters the living environment of all marine organisms. Studying adaptation, however, can be tricky because spatial genetic patterns might also occur due to random effects, for example, genetic drift. Genetic drift is reduced in very large and well‐connected populations, such as in broadcast marine spawning organisms. Here, spatial genetic divergence is likely to be produced by selection. In this issue of Molecular Ecology, Sandoval‐Castillo et al. (2018) investigated patterns of spatial genetic divergence and their association with environmental factors in the greenlip abalone (Haliotis laevigata). This commercially important species of mollusc is a broadcast spawner with large population sizes, rendering genetic drift an unlikely factor in the genetic divergence of wild populations. Sandoval‐Castillo et al. (2018) used a ddRAD genomic approach to test for genetic divergence between sampled populations while also measuring different environmental factors, for example, water temperature and oxygen content. The majority of identified SNPs was putatively neutral and showed only low levels of genetic divergence between field sites. However, 323 candidate adaptive markers were identified that clearly separated the individuals into five different clusters. These genetic clusters correlated with environmental clusters mainly determined by water temperature and (correlated) oxygen concentration. Gene annotation of the candidate SNPs revealed a large proportion of loci being involved in biological processes influenced by oxygen availability. The study by Sandoval‐Castillo et al. (2018) in this issue of Molecular Ecology exemplifies the benefits of combining genomic studies with ecological data. It is a great starting point for more detailed (gene function, physiology) as well as broader (biodiversity) investigations that might help us to better understand adaptation and predict ecosystems' resilience and resistance to environmental disturbances. In addition, this information can be applied to implement optimal conservation regime policies and sustainable harvesting strategies, hopefully protecting biodiversity as well as commercial interests in marine life.