Amanda J. Zellmer

Applications of next-generation sequencing to phylogeography and phylogenetics

This is a time of unprecedented transition in DNA sequencing technologies. Next-generation sequencing (NGS) clearly holds promise for fast and cost-effective generation of multilocus sequence data for phylogeography and phylogenetics. However, the focus on non-model organisms, in addition to uncertainty about which sample preparation methods and analyses are appropriate for different research questions and evolutionary timescales, have contributed to a lag in the application of NGS to these fields. Here, we outline some of the major obstacles specific to the application of NGS to phylogeography and phylogenetics, including the focus on non-model organisms, the necessity of obtaining orthologous loci in a cost-effective manner, and the predominate use of gene trees in these fields. We describe the most promising methods of sample preparation that address these challenges. Methods that reduce the genome by restriction digest and manual size selection are most appropriate for studies at the intraspecific level, whereas methods that target specific genomic regions (i.e., target enrichment or sequence capture) have wider applicability from the population level to deep-level phylogenomics. Additionally, we give an overview of how to analyze NGS data to arrive at data sets applicable to the standard toolkit of phylogeography and phylogenetics, including initial data processing to alignment and genotype calling (both SNPs and loci involving many SNPs). Even though whole-genome sequencing is likely to become affordable rather soon, because phylogeography and phylogenetics rely on analysis of hundreds of individuals in many cases, methods that reduce the genome to a subset of loci should remain more cost-effective for some time to come.

Does niche divergence accompany allopatric divergence in Aphelocoma jays as predicted under ecological speciation? Insights from tests with niche models.

The role of ecology in the origin of species has been the subject of long‐standing interest to evolutionary biologists. New sources of spatially explicit ecological data allow for large‐scale tests of whether speciation is associated with niche divergence or whether closely related species tend to be similar ecologically (niche conservatism). Because of the confounding effects of spatial autocorrelation of environmental variables, we generate null expectations for niche divergence for both an ecological‐niche modeling and a multivariate approach to address the question: do allopatrically distributed taxa occupy similar niches? In a classic system for the study of niche evolution—the Aphelocoma jays—we show that there is little evidence for niche divergence among Mexican Jay (A. ultramarina) lineages in the process of speciation, contrary to previous results. In contrast, Aphelocoma species that exist in partial sympatry in some regions show evidence for niche divergence. Our approach is widely applicable to the many cases of allopatric lineages in the beginning stages of speciation. These results do not support an ecological speciation model for Mexican Jay lineages because, in most cases, the allopatric environments they occupy are not significantly more divergent than expected under a null model.

Disentangling the effects of historic vs. contemporary landscape structure on population genetic divergence

Increasing habitat fragmentation poses an immediate threat to population viability, as gene flow patterns are changed in these altered landscapes. Patterns of genetic divergence can potentially reveal the impact of these shifts in landscape connectivity. However, divergence patterns not only carry the signature of altered contemporary landscapes, but also historical ones. When considered separately, both recent and historical landscape structure appear to significantly affect connectivity among 51 wood frog (Rana sylvatica) populations. However, by controlling for correlations among landscape structure from multiple time periods, we show that patterns of genetic divergence reflect recent landscape structure as opposed to landscape structure prior to European settlement of the region (before 1850s). At the same time, within‐population genetic diversities remain high and a genetic signature of population bottlenecks is lacking. Together, these results suggest that metapopulation processes – not drift‐induced divergence associated with strong demographic bottlenecks following habitat loss – underlie the strikingly rapid consequences of temporally shifting landscape structure on these amphibians. We discuss the implications of these results in the context of understanding the role of population demography in the adaptive variation observed in wood frog populations.

Deep Phylogeographic Structure and Environmental Differentiation in the Carnivorous Plant Sarracenia alata

We collected ~29 kb of sequence data using Roche 454 pyrosequencing in order to estimate the timing and pattern of diversification in the carnivorous pitcher plant Sarracenia alata. Utilizing modified protocols for reduced representation library construction, we generated sequence data from 86 individuals across 10 populations from throughout the range of the species. We identified 76 high-quality and high-coverage loci (containing over 500 SNPs) using the bioinformatics pipeline PRGmatic. Results from a Bayesian clustering analysis indicate that populations are highly structured, and are similar in pattern to the topology of a population tree estimated using *BEAST. The pattern of diversification within Sarracenia alata implies that riverine barriers are the primary factor promoting population diversification, with divergence across the Mississippi River occurring more than 60,000 generations before present. Further, significant patterns of niche divergence and the identification of several outlier loci suggest that selection may contribute to population divergence. Our results demonstrate the feasibility of using next-generation sequencing to investigate intraspecific genetic variation in nonmodel species.

Historical Species Distribution Models Predict Species Limits in Western Plethodon Salamanders

Allopatry is commonly used to predict boundaries in species delimitation investigations under the assumption that currently allopatric distributions are indicative of reproductive isolation; however, species ranges are known to change over time. Incorporating a temporal perspective of geographic distributions should improve species delimitation; to explore this, we investigate three species of western Plethodon salamanders that have shifted their ranges since the end of the Pleistocene. We generate species distribution models (SDM) of the current range, hindcast these models onto a climatic model 21 Ka, and use three molecular approaches to delimit species in an integrated fashion. In contrast to expectations based on the current distribution, we detect no independent lineages in species with allopatric and patchy distributions (Plethodon vandykei and Plethodon larselli). The SDMs indicate that probable habitat is more expansive than their current range, especially during the last glacial maximum (LGM) (21 Ka). However, with a contiguous distribution, two independent lineages were detected in Plethodon idahoensis, possibly due to isolation in multiple glacial refugia. Results indicate that historical SDMs are a better predictor of species boundaries than current distributions, and strongly imply that researchers should incorporate SDM and hindcasting into their investigations and the development of species hypotheses.

Genomic data reveal ancient microendemism in forest scorpions across the California Floristic Province

The California Floristic Province (CFP) in western North America is a globally significant biodiversity hotspot. Elucidating patterns of endemism and the historical drivers of this diversity has been an important challenge of comparative phylogeography for over two decades. We generated phylogenomic data using ddRADseq to examine genetic structure in Uroctonus forest scorpions, an ecologically restricted and dispersal‐limited organism widely distributed across the CFP north to the Columbia River. We coupled our genetic data with species distribution models (SDMs) to determine climatically suitable areas for Uroctonus both now and during the Last Glacial Maximum. Based on our analyses, Uroctonus is composed of two major genetic groups that likely diverged over 2 million years ago. Each of these groups itself contains numerous genetic groups that reveal a pattern of vicariance and microendemism across the CFP. Migration rates among these populations are low. SDMs suggest forest scorpion habitat has remained relatively stable over the last 21 000 years, consistent with the genetic data. Our results suggest tectonic plate rafting, mountain uplift, river drainage formation and climate‐induced habitat fragmentation have all likely played a role in the diversification of Uroctonus. The intricate pattern of genetic fragmentation revealed across a temporal continuum highlights the potential of low‐dispersing species to shed light on small‐scale patterns of biodiversity and the underlying processes that have generated this diversity in biodiversity hotspots.

Biogeographic barriers drive co-diversification within associated eukaryotes of the Sarracenia alata pitcher plant system

Understanding if the members of an ecological community have co-diversified is a central concern of evolutionary biology, as co-diversification suggests prolonged association and possible coevolution. By sampling associated species from an ecosystem, researchers can better understand how abiotic and biotic factors influence diversification in a region. In particular, studies of co-distributed species that interact ecologically can allow us to disentangle the effect of how historical processes have helped shape community level structure and interactions. Here we investigate the Sarracenia alata pitcher plant system, an ecological community where many species from disparate taxonomic groups live inside the fluid-filled pitcher leaves. Direct sequencing of the eukaryotes present in the pitcher plant fluid enables us to better understand how a host plant can shape and contribute to the genetic structure of its associated inquilines, and to ask whether genetic variation in the taxa are structured in a similar manner to the host plant. We used 454 amplicon-based metagenomics to demonstrate that the pattern of genetic diversity in many, but not all, of the eukaryotic community is similar to that of S. alata, providing evidence that associated eukaryotes share an evolutionary history with the host pitcher plant. Our work provides further evidence that a host plant can influence the evolution of its associated commensals.

Genome-wide signals of drift and local adaptation during rapid lineage divergence in a songbird

The formation of independent evolutionary lineages involves neutral and selective factors, and understanding their relative roles in population divergence is a fundamental goal of speciation research. Correlations between allele frequencies and environmental variability can reveal the role of selection, yet the relative contribution of drift can be difficult to establish. Recently diversified systems such as that of the Oregon junco (Aves: Emberizidae) of western North America provide ideal scenarios to apply genetic-environment association analyses (GEA) while controlling for population structure. Genome-wide SNP loci analyses revealed marked genetic structure consisting of differentiated populations in isolated, dry southern mountain ranges, and more admixed recently expanded populations in humid northern latitudes. We used correlations between genomic and environmental variance to test for three specific modes of evolutionary divergence: (i) drift in geographic isolation, (ii) differentiation along continuous selective gradients, and (iii) isolation by adaptation. We found evidence of strong drift in southern mountains, but also signals of local adaptation in several populations, driven by temperature, precipitation, elevation and vegetation, especially when controlling for population history. We identified numerous variants under selection scattered across the genome, suggesting that local adaptation can promote rapid differentiation over short periods when acting over multiple independent loci.

Microgeographic morphological variation across larval wood frog populations associated with environment despite gene flow

Abstract Gene flow has historically been thought to constrain local adaptation; yet, recent research suggests that populations can diverge despite exchanging genes. Here I use a common garden experiment to assess the combined effects of gene flow and natural selection on morphological variation of 16 wood frog (Rana sylvatica) populations, a species known to experience divergent selection pressures in open‐ and closed‐canopy ponds across relatively small geographic scales. Wood frog tadpoles from different ponds showed significant morphological variation associated with canopy type with a trade‐off between tail length and body depth consistent with previous research. In contrast, neutral genetic differentiation of nine microsatellite loci as measured by Josts D was not associated with canopy type, indicating no pattern of isolation by environment. Genetic structure analyses indicated some substructure across the 16 ponds (K = 4); however, three out of four assigned clusters included both open‐ and closed‐canopy ponds. Together, these results suggest that morphological divergence among these wood frog populations is occurring despite gene flow and that selection within these environments is strong. Furthermore, morphological variation among ponds differed across two sampling periods during larval development, demonstrating the importance of evaluating phenotypic divergence over multiple time periods and at a time relevant to the processes being studied.

Culturally inclusive STEM education

Brazil’s government is reopening bidding rounds for deep-sea oil and gas exploration after 4 years of economic and political turmoil. According to its ambitious 4-year plan (1), through which the government expects to profit from licenses and production royalties, Brazil will lease hundreds of offshore areas for exploration in depths below 200 m. Allowing such exploration will substantially expand the offshore industry in regions that are of biological and ecological relevance for deep-sea conservation (2). In deep-sea basins within Brazil’s exclusive economic zone, where more than 70% of the current offshore oil production is concentrated, there is already substantial overlap between leased areas and vulnerable marine ecosystems, including cold-water corals and submarine canyons (3). Other vulnerable ecosystems, such as cold seeps, are poorly reported even in basins where there is biological and geophysical evidence for their presence (4, 5). The spatial overlap and depth distribution of pockmarks within oil fields suggest that seeps may be common in areas currently offered on bidding rounds. As a result, the planned expansion of offshore leasing areas will increase the impacts of the offshore oil industry from the Amazon to the Edited by Jennifer Sills Deep-sea drilling threatens cold seep ecosystems.