Erratum in J Hered. 2022 May 16;113(2):217.

The Atlantic spotted dolphin (Stenella frontalis) is endemic to tropical, subtropical, and warm temperate waters of the Atlantic Ocean. Throughout its distribution, both geographic distance and environmental variation may contribute to population structure of the species. In this study, we follow a seascape genetics approach to investigate population differentiation of Atlantic spotted dolphins based on a large worldwide dataset and the relationship with marine environmental variables. The results revealed that the Atlantic spotted dolphin exhibits population genetic structure across its distribution based on mitochondrial DNA control region (mtDNA-CR) data. Analyses based on the contemporary landscape suggested, at both the individual and population level, that the population genetic structure is consistent with the isolation-by-distance model. However, because geography and environmental matrices were correlated, and because in some, but not all analyses, we found a significant effect for the environment, we cannot rule out the addition contribution of environmental factors in structuring genetic variation. Future analyses based on nuclear data are needed to evaluate whether local processes, such as social structure and some level of philopatry within populations, may be contributing to the associations among genetic structure, geographic, and environmental distance.

Highly social species are successful because they cooperate in obligately integrated societies. We examined temporal genetic variation in the eusocial wasp Vespula maculifrons to gain a greater understanding of evolution in highly social taxa. First, we wished to test if effective population sizes of eusocial species were relatively low due to the reproductive division of labor that characterizes eusocial taxa. We thus estimated the effective population size of V. maculifrons by examining temporal changes in population allele frequencies. We sampled the genetic composition of a V. maculifrons population at 3 separate timepoints spanning a 13-year period. We found that effective population size ranged in the hundreds of individuals, which is similar to estimates in other, non-eusocial taxa. Second, we estimated levels of polyandry in V. maculifrons in different years to determine if queen mating system varied over time. We found no significant change in the number or skew of males mated to queens. In addition, mating skew was not significant within V. maculifrons colonies. Therefore, our data suggest that queen mate number may be subject to stabilizing selection in this taxon. Overall, our study provides novel insight into the selective processes operating in eusocial species by analyzing temporal genetic changes within populations.

Genetic changes underlying adaptation vary greatly in terms of complexity and, within the same species, genetic responses to similar selective pressures may or may not be the same. We examine both complex (supergene) and simple (SNP) genetic variants occurring in populations of rainbow trout (Oncorhynchus mykiss) independently isolated from ocean access and compared them to each other and to an anadromous below-barrier population representing their ancestral source to search for signatures of both parallel and nonparallel adaptation. All landlocked populations displayed an increased frequency of a large inversion on chromosome Omy05, while 3 of the 4 populations exhibited elevated frequencies of another inversion located on chromosome Omy20. In addition, we identified numerous regions outside these 2 inversions that also show significant shifts in allele frequencies consistent with adaptive evolution. However, there was little concordance among above-barrier populations in these specific genomic regions under selection. In part, the lack of concordance appears to arise from ancestral autopolyploidy in rainbow trout that provides duplicate genomic regions of similar functional composition for selection to act upon. Thus, while selection acting on landlocked populations universally favors the resident ecotype, outside of the major chromosomal inversions, the resulting genetic changes are largely distinct among populations. Our results indicate that selection on standing genetic variation is likely the primary mode of rapid adaptation, and that both supergene complexes and individual loci contribute to adaptive evolution, further highlighting the diversity of adaptive genomic variation involved in complex phenotypic evolution.

Ancient introgression can be an important source of genetic variation that shapes the evolution and diversification of many taxa. Here, we estimate the timing, direction, and extent of gene flow between two distantly related oak species in the same section (Quercus sect. Quercus). We estimated these demographic events using genotyping by sequencing data, which generated 25 702 single nucleotide polymorphisms for 24 individuals of California scrub oak (Quercus berberidifolia) and 23 individuals of Engelmann oak (Quercus engelmannii). We tested several scenarios involving gene flow between these species using the diffusion approximation-based population genetic inference framework and model-testing approach of the Python package DaDi. We found that the most likely demographic scenario includes a bottleneck in Q. engelmannii that coincides with asymmetric gene flow from Q. berberidifolia into Q. engelmannii. Given that the timing of this gene flow coincides with the advent of a Mediterranean-type climate in the California Floristic Province, we propose that changing precipitation patterns and seasonality may have favored the introgression of climate-associated genes from the endemic into the non-endemic California oak.

Three subspecies of Northern Bahamian Rock Iguanas, Cyclura cychlura, are currently recognized: C. c. cychlura, restricted to Andros Island, and C. c. figginsi and C. c. inornata, native to the Exuma Island chain. Populations on Andros are genetically distinct from Exuma Island populations, yet genetic divergence among populations in the Exumas is inconsistent with the 2 currently recognized subspecies from those islands. The potential consequences of this discrepancy might include the recognition of a single subspecies throughout the Exumas rather than 2. That inference also ignores evidence that populations of C. cychlura are potentially adaptively divergent. We compared patterns of population relatedness in a three-tiered host-parasite system: C. cychlura iguanas, their ticks (genus Amblyomma, preferentially parasitizing these reptiles), and Rickettsia spp. endosymbionts (within tick ectoparasites). Our results indicate that while C. c. cychlura on Andros is consistently supported as a separate clade, patterns of relatedness among populations of C. c. figginsi and C. c. inornata within the Exuma Island chain are more complex. The distribution of the hosts, different tick species, and Rickettsia spp., supports the evolutionary independence of C. c. inornata. Further, these patterns are also consistent with two independent evolutionarily significant units within C. c. figginsi. Our findings suggest coevolutionary relationships between the reptile hosts, their ectoparasites, and rickettsial organisms, suggesting local adaptation. This work also speaks to the limitations of using neutral molecular markers from a single focal taxon as the sole currency for recognizing evolutionary novelty in populations of endangered species.

Salmonids have emerged as important study systems for investigating molecular processes underlying parallel evolution given their tremendous life history variation. Kokanee, the resident form of anadromous sockeye salmon (Oncorhynchus nerka), have evolved multiple times across the species' pan-Pacific distribution, exhibiting multiple reproductive ecotypes including those that spawn in streams, on lake-shores, and at lake depths >50 m. The latter has only been detected in 5 locations in Japan and British Columbia, Canada. Here, we investigated the multiple origins hypothesis for deep-spawning kokanee, using 9721 single nucleotide polymorphisms distributed across the genome analyzed for the vast majority of known populations in Japan (Saiko Lake) and Canada (Anderson, Seton, East Barrière Lakes) relative to stream-spawning populations in both regions. We detected 397 outlier loci, none of which were robustly identified in paired-ecotype comparisons in Japan and Canada independently. Bayesian clustering and principal components analyses based on neutral loci revealed 6 distinct clusters, largely associated with geography or translocation history, rather than ecotype. Moreover, a high level of divergence between Canadian and Japanese populations, and between deep- and stream-spawning populations regionally, suggests the deep-spawning ecotype independently evolved on the 2 continents. On a finer level, Japanese kokanee populations exhibited low estimates of heterozygosity, significant levels of inbreeding, and reduced effective population sizes relative to Canadian populations, likely associated with transplantation history. Along with preliminary evidence for hybridization between deep- and stream-spawning ecotypes in Saiko Lake, these findings should be considered within the context of on-going kokanee fisheries management in Japan.

The southeastern United States is home to some of the richest biodiversity in the world. Over the last 200 years, however, rapid industrialization and urbanization have threatened many natural areas, including freshwater habitats. River impoundments have also rapidly altered freshwater habitats, often resulting in species extirpation or extinction. The Coosa River in Alabama experienced one of the largest faunal declines in modern history after impoundment, making it an ideal system for studying how invertebrate species are affected by reservoir creation. One such species, the Rough Hornsnail, Pleurocera foremani, is an endangered freshwater snail in the family Pleuroceridae. We sampled all known localities of P. foremani and used 2bRAD-seq to measure genetic diversity. We assessed riverscape genomic patterns across the current range of P. foremani and measured gene flow within and between impoundments. We also investigated the degree to which P. foremani displays an isolation by distance pattern and conforms to broad hypotheses that have been put forth for population genetics of riverine species like the Mighty Headwater Hypothesis that predicts greater genetic diversity in headwater reaches compared with mainstem populations. Like most other freshwater species, a pattern of isolation by distance was observed in P. foremani. We also found that Coosa River dams are a barrier to gene flow, and genetic fragmentation of P. foremani is likely to increase. However, gene flow appeared common within reservoirs and tributaries. Additionally, we found that spatial genetic structure of P. foremani deviates from what is expected under the Mighty Headwaters Hypothesis, adding to a growing body of research suggesting that the majority of genetic diversity in low-dispersing gastropods is found in mainstem populations.

As genomic-scale data sets become economically feasible for most organisms, a key question for conservation biology is whether the increased resolution offered by new genomic approaches justifies repeating earlier studies based on traditional markers, rather than investing those same time and monetary resources in less-known species. Genomic studies offer clear advantages when the objective is to identify adaptive loci that may be critical to conservation policy-makers. However, the answer is far less certain for the population and landscape studies based on neutral loci that dominate the conservation genetics research agenda. We used Restriction-site Associated DNA sequencing (RADseq) to revisit earlier molecular studies of the IUCN Critically Endangered Magdalena River turtle (Podocnemis lewyana), documenting the conservation insights gained by increasing the number of neutral markers by several orders of magnitude. Earlier research indicated that P. lewyana has the lowest genetic diversity known for any chelonian, and little or no population differentiation among independent rivers. In contrast, the RADseq data revealed discrete population structure with isolation-by-distance within river segments and identified precise population breaks clearly delineating management units. It also confirmed that the species does not have extremely low heterozygosity and that effective population sizes are probably sufficient to maintain long-term evolutionary potential. Contrary to earlier inferences from more limited population genetic markers, our genomic data suggest that management strategies should shift from active genetic rescue to more passive protection without extreme interventions. We conclude with a list of examples of conservation studies in other vertebrates indicating that for many systems a genomic update is worth the investment.