Temporal Stability of Genetic Structure in a Mesopelagic Copepod

Life in the Midwater: The Ecology of Deep Pelagic Animals

The water column of the deep ocean is dark, cold, low in food, and under crushing pressures, yet it is full of diverse life. Due to its enormous volume, this mesopelagic zone is home to some of the most abundant animals on the planet. Rather than struggling to survive, they thrive-owing to a broad set of adaptations for feeding, behavior, and physiology. Our understanding of these adaptations is constrained by the tools available for exploring the deep sea, but this tool kit is expanding along with technological advances. Each time we apply a new method to the depths, we gain surprising insights about genetics, ecology, behavior, physiology, diversity, and the dynamics of change. These discoveries show structure within the seemingly uniform habitat, limits to the seemingly inexhaustible resources, and vulnerability in the seemingly impervious environment. To understand midwater ecology, we need to reimagine the rules that govern terrestrial ecosystems. By spending more time at depth-with whatever tools are available-we can fill the knowledge gaps and better link ecology to the environment throughout the water column.

Oceanic dispersal barriers in a holoplanktonic gastropod

Pteropods, a group of holoplanktonic gastropods, are regarded as bioindicators of the effects of ocean acidification on open ocean ecosystems, because their thin aragonitic shells are susceptible to dissolution. While there have been recent efforts to address their capacity for physiological acclimation, it is also important to gain predictive understanding of their ability to adapt to future ocean conditions. However, little is known about the levels of genetic variation and large-scale population structuring of pteropods, key characteristics enabling local adaptation. We examined the spatial distribution of genetic diversity in the mitochondrial cytochrome c oxidase I (COI) and nuclear 28S gene fragments, as well as shell shape variation, across a latitudinal transect in the Atlantic Ocean (35°N-36°S) for the pteropod Limacina bulimoides. We observed high levels of genetic variability (COI π = 0.034, 28S π = 0.0021) and strong spatial structuring (COI ΦST  = 0.230, 28S ΦST  = 0.255) across this transect. Based on the congruence of mitochondrial and nuclear differentiation, as well as differences in shell shape, we identified a primary dispersal barrier in the southern Atlantic subtropical gyre (15-18°S). This barrier is maintained despite the presence of expatriates, a gyral current system, and in the absence of any distinct oceanographic gradients in this region, suggesting that reproductive isolation between these populations must be strong. A secondary dispersal barrier supported only by 28S pairwise ΦST comparisons was identified in the equatorial upwelling region (between 15°N and 4°S), which is concordant with barriers observed in other zooplankton species. Both oceanic dispersal barriers were congruent with regions of low abundance reported for a similar basin-scale transect that was sampled 2 years later. Our finding supports the hypothesis that low abundance indicates areas of suboptimal habitat that result in barriers to gene flow in widely distributed zooplankton species. Such species may in fact consist of several populations or (sub)species that are adapted to local environmental conditions, limiting their potential for adaptive responses to ocean changes. Future analyses of genome-wide diversity in pteropods could provide further insight into the strength, formation and maintenance of oceanic dispersal barriers.

Genetic diversity and novel lineages in the cosmopolitan copepod Pleuromamma abdominalis in the Southeast Pacific

Across boundary currents, zooplankton are subject to strong oceanographic gradients and hence strong selective pressures. How such gradients interact with the speciation process of pelagic organisms is still poorly understood in the open ocean realm. Here we report on genetic diversity within the pelagic copepod Pleuromamma abdominalis in the poorly known Southeast Pacific region, with samples spanning an ocean gradient from coastal upwelling to the oligotrophic South Pacific Subtropical Gyre. We assessed variation in fragments of the mitochondrial (mt) genes cytochrome c oxidase subunit I (COI) and Cytochrome b as well as in the nuclear internal transcribed spacer (ITS) region and 28 S rRNA. Phylogenetic analyses revealed the presence of 8 divergent lineages occurring across the gradient with genetic distances in the range of 0.036 and 0.44 (mt genes), and GMYC species delimitation methods support their inference as distinct (undescribed) species. Genetic lineages occurring across the zonal gradient showed strong genetic structuring, with the presence of at least two new lineages within the coastal upwelling zone, revealing an unexpectedly high level of endemism within the Humboldt Current System. Multivariate analyses found strong correlation between genetic variation and surface chlorophyll-a and salinity, suggesting an important role for hydrographic gradients in maintaining genetic diversity. However, the presence of cryptic lineages within the upwelling zone cannot be easily accounted for by environmental heterogeneity and poses challenging questions for understanding the speciation process for oceanic zooplankton.

Assessing species boundaries in the open sea: an integrative taxonomic approach to the pteropod genus Diacavolinia

To track changes in pelagic biodiversity in response to climate change, it is essential to accurately define species boundaries. Shelled pteropods are a group of holoplanktonic gastropods that have been proposed as bio-indicators because of their vulnerability to ocean acidification. A particularly suitable, yet challenging group for integrative taxonomy is the pteropod genus Diacavolinia, which has a circumglobal distribution and is the most species-rich pteropod genus, with 24 described species. We assessed species boundaries in this genus, with inferences based on geometric morphometric analyses of shell-shape variation, genetic (cytochrome c oxidase subunit I, 28S rDNA sequences) and geographic data. We found support for a total of 13 species worldwide, with observations of 706 museum and 263 freshly collected specimens across a global collection of material, including holo‐ and paratype specimens for 14 species. In the Atlantic Ocean, two species are well supported, in contrast to the eight currently described, and in the Indo‐Pacific we found a maximum of 11 species, partially merging 13 of the described species. Distributions of these revised species are congruent with well-known biogeographic provinces. Combining varied datasets in an integrative framework may be suitable for many diverse taxa and is an important first step to predicting species-specific responses to global change.

Beyond the snapshot: Landscape genetic analysis of time series data reveal responses of American black bears to landscape change

Landscape genetic studies typically focus on the evolutionary processes that give rise to spatial patterns that are quantified at a single point in time. Although landscape change is widely recognized as a strong driver of microevolutionary processes, few landscape genetic studies have directly evaluated the change in spatial genetic structure (SGS) over time with concurrent changes in landscape pattern. We introduce a novel approach to analyze landscape genetic data through time. We demonstrate this approach using genotyped samples (n = 569) from a large black bear (Ursus americanus) population in Michigan (USA) that were harvested during 3 years (2002, 2006, and 2010). We identified areas that were consistently occupied over this 9-year period and quantified temporal variation in SGS. Then, we evaluated alternative hypotheses about effects of changes in landscape features (e.g., deforestation or crop conversion) on fine-scale SGS among years using spatial autoregressive modeling and model selection. Relative measures of landscape change such as magnitude of landscape change (i.e., number of patches changing from suitable to unsuitable states or vice versa), and during later periods, measures of fragmentation (i.e., patch aggregation and cohesion) were associated with change in SGS. Our results stress the importance of conducting time series studies for the conservation and management of wildlife inhabiting rapidly changing landscapes.

Temperature gradient affects differentiation of gene expression and SNP allele frequencies in the dominant Lake Baikal zooplankton species

Local adaptation and phenotypic plasticity are main mechanisms of organisms' resilience in changing environments. Both are affected by gene flow and are expected to be weak in zooplankton populations inhabiting large continuous water bodies and strongly affected by currents. Lake Baikal, the deepest and one of the coldest lakes on Earth, experienced epilimnion temperature increase during the last 100 years, exposing Baikal's zooplankton to novel selective pressures. We obtained a partial transcriptome of Epischura baikalensis (Copepoda: Calanoida), the dominant component of Baikal's zooplankton, and estimated SNP allele frequencies and transcript abundances in samples from regions of Baikal that differ in multiyear average surface temperatures. The strongest signal in both SNP and transcript abundance differentiation is the SW-NE gradient along the 600+ km long axis of the lake, suggesting isolation by distance. SNP differentiation is stronger for nonsynonymous than synonymous SNPs and is paralleled by differential survival during a laboratory exposure to increased temperature, indicating directional selection operating on the temperature gradient. Transcript abundance, generally collinear with the SNP differentiation, shows samples from the warmest, less deep location clustering together with the southernmost samples. Differential expression is more frequent among transcripts orthologous to candidate thermal response genes previously identified in model arthropods, including genes encoding cytoskeleton proteins, heat-shock proteins, proteases, enzymes of central energy metabolism, lipid and antioxidant pathways. We conclude that the pivotal endemic zooplankton species in Lake Baikal exists under temperature-mediated selection and possesses both genetic variation and plasticity to respond to novel temperature-related environmental pressures.

Temporal genetic stability in natural populations of the waterflea Daphnia magna in response to strong selection pressure

Studies monitoring changes in genetic diversity and composition through time allow a unique understanding of evolutionary dynamics and persistence of natural populations. However, such studies are often limited to species with short generation times that can be propagated in the laboratory or few exceptional cases in the wild. Species that produce dormant stages provide powerful models for the reconstruction of evolutionary dynamics in the natural environment. A remaining open question is to what extent dormant egg banks are an unbiased representation of populations and hence of the species' evolutionary potential, especially in the presence of strong environmental selection. We address this key question using the water flea Daphnia magna, which produces dormant stages that accumulate in biological archives over time. We assess temporal genetic stability in three biological archives, previously used in resurrection ecology studies showing adaptive evolutionary responses to rapid environmental change. We show that neutral genetic diversity does not decline with the age of the population and it is maintained in the presence of strong selection. In addition, by comparing temporal genetic stability in hatched and unhatched populations from the same biological archive, we show that dormant egg banks can be consulted to obtain a reliable measure of genetic diversity over time, at least in the multidecadal time frame studied here. The stability of neutral genetic diversity through time is likely mediated by the buffering effect of the resting egg bank.