How sea level change mediates genetic divergence in coastal species across regions with varying tectonic and sediment processes

Biogeography of Greater Antillean freshwater fishes, with a review of competing hypotheses

In biogeography, vicariance and long-distance dispersal are often characterised as competing scenarios. However, they are related concepts, both relying on collective geological, ecological, and phylogenetic evidence. This is illustrated by freshwater fishes, which may immigrate to islands either when freshwater connections are temporarily present and later severed (vicariance), or by unusual means when ocean gaps are crossed (long-distance dispersal). Marine barriers have a strong filtering effect on freshwater fishes, limiting immigrants to those most capable of oceanic dispersal. The roles of vicariance and dispersal are debated for freshwater fishes of the Greater Antilles. We review three active hypotheses [Cretaceous vicariance, Greater Antilles-Aves Ridge (GAARlandia), long-distance dispersal] and propose long-distance dispersal to be an appropriate model due to limited support for freshwater fish use of landspans. Greater Antillean freshwater fishes have six potential source bioregions (defined from faunal similarity): Northern Gulf of México, Western Gulf of México, Maya Terrane, Chortís Block, Eastern Panamá, and Northern South America. Faunas of the Greater Antilles are composed of taxa immigrating from many of these bioregions, but there is strong compositional disharmony between island and mainland fish faunas (>90% of Antillean species are cyprinodontiforms, compared to <10% in Northern Gulf of México and Northern South America, and ≤50% elsewhere), consistent with a hypothesis of long-distance dispersal. Ancestral-area reconstruction analysis indicates there were 16 or 17 immigration events over the last 51 million years, 14 or 15 of these by cyprinodontiforms. Published divergence estimates and evidence available for each immigration event suggests they occurred at different times and by different pathways, possibly with rafts of vegetation discharged from rivers or washed to sea during storms. If so, ocean currents likely provide critical pathways for immigration when flowing from one landmass to another. On the other hand, currents create dispersal barriers when flowing perpendicularly between landmasses. In addition to high salinity tolerance, cyprinodontiforms collectively display a variety of adaptations that could enhance their ability to live with rafts (small body size, viviparity, low metabolism, amphibiousness, diapause, self-fertilisation). These adaptations likely also helped immigrants establish island populations after arrival and to persist long term thereafter. Cichlids may have used a pseudo bridge (Nicaragua Rise) to reach the Greater Antilles. Gars (Lepisosteidae) may have crossed the Straits of Florida to Cuba, a relatively short crossing that is not a barrier to gene flow for several cyprinodontiform immigrants. Indeed, widespread distributions of Quaternary migrants (Cyprinodon, Gambusia, Kryptolebias), within the Greater Antilles and among neighbouring bioregions, imply that long-distance dispersal is not necessarily inhibitory for well-adapted species, even though it appears to be virtually impossible for all other freshwater fishes.

Genetic structure and historic demography of endangered unarmoured threespine stickleback at southern latitudes signals a potential new management approach

Habitat loss, flood control infrastructure, and drought have left most of southern California and northern Baja California's native freshwater fish near extinction, including the endangered unarmoured threespine stickleback (Gasterosteus aculeatus williamsoni). This subspecies, an unusual morph lacking the typical lateral bony plates of the G. aculeatus complex, occurs at arid southern latitudes in the eastern Pacific Ocean and survives in only three inland locations. Managers have lacked molecular data to answer basic questions about the ancestry and genetic distinctiveness of unarmoured populations. These data could be used to prioritize conservation efforts. We sampled G. aculeatus from 36 localities and used microsatellites and whole genome data to place unarmoured populations within the broader evolutionary context of G. aculeatus across southern California/northern Baja California. We identified three genetic groups with none consisting solely of unarmoured populations. Unlike G. aculeatus at northern latitudes, where Pleistocene glaciation has produced similar historical demographic profiles across populations, we found markedly different demographics depending on sampling location, with inland unarmoured populations showing steeper population declines and lower heterozygosity compared to low armoured populations in coastal lagoons. One exception involved the only high elevation population in the region, where the demography and alleles of unarmoured fish were similar to low armoured populations near the coast, exposing one of several cases of artificial translocation. Our results suggest that the current "management-by-phenotype" approach, based on lateral plates, is incidentally protecting the most imperilled populations; however, redirecting efforts toward evolutionary units, regardless of phenotype, may more effectively preserve adaptive potential.

Towards a unified framework to study causality in Earth-life systems

There is considerable interest in better understanding how earth processes shape the generation and distribution of life on Earth. This question, at its heart, is one of causation. In this article I propose that at a regional level, earth processes can be thought of as behaving somewhat deterministically and may have an organized effect on the diversification and distribution of species. However, the study of how landscape features shape biology is challenged by pseudocongruent or collinear variables. I demonstrate that causal structures can be used to depict the cause–effect relationships between earth processes and biological patterns using recent examples from the literature about speciation and species richness in montane settings. This application shows that causal diagrams can be used to better decipher the details of causal relationships by motivating new hypotheses. Additionally, the abstraction of this knowledge into structural equation metamodels can be used to formulate theory about relationships within Earth–life systems more broadly. Causal structures are a natural point of collaboration between biologists and Earth scientists, and their use can mitigate against the risk of misassigning causality within studies. My goal is that by applying causal theory through application of causal structures, we can build a systems‐level understanding of what landscape features or earth processes most shape the distribution and diversification of species, what types of organisms are most affected, and why.

The hidden landscape: Evidence that sea-level change shaped the present population genomic patterns of marginal marine species

Oftentimes, to understand the genetic relatedness and diversity of today's populations requires considering the ancient landscape on which those populations evolved. Nowhere is this clearer than along Earth's coastline, which has been in its present-day configuration for only about 6.5% of the past 800,000 years (Dolby et al., 2020; Miller et al., 2005). During ice ages when glaciers expanded in the Northern Hemisphere, they stored enough of the planet's water to drop global sea level by ~120 m below present levels ("lowstand", Figure 1a), and there have been at least eight of these 100,000-year cycles preceding today. When glaciers melted, ocean water reflooded shorelines, shifting and re-forming marginal marine habitats globally and shaping the relatedness of populations (Dolby et al., 2016). In a From the Cover article in this issue of Molecular Ecology, Stiller et al. (2020) integrate population genomic analysis of leafy seadragons in southern Australia with estimates of available seabed area to reveal that the expansion of habitat that accompanied this reflooding led to strong demographic expansions. With statistical models, they also show that western populations were eliminated and then recolonized because the continental shelf there is narrow, leaving little available habitat when sea level was low (Figure 1b). Their results document the dynamic and interrelated nature of a hidden, changing landscape and the evolution of species inhabiting it.

Using ultraconserved elements to track the influence of sea-level change on leafy seadragon populations

During the Last Glacial Maximum (LGM), global sea levels were 120-130 m lower than today, resulting in the emergence of most continental shelves and extirpation of subtidal organisms from these areas. During the interglacial periods, rapid inundation of shelf regions created a dynamic environment for coastal organisms, such as the charismatic leafy seadragon (Phycodurus eques, Syngnathidae), a brooder with low dispersal ability inhabiting kelp beds in temperate Australia. Reconstructions of the palaeoshoreline revealed that the increase of shallow areas since the LGM was not uniform across the species' range and we investigated the effects of these asymmetries on genetic diversity and structuring. Using targeted capture of 857 variable ultraconserved elements (UCEs, 2,845 single nucleotide polymorphisms) in 68 individuals, we found that the regionally different shelf topographies were paralleled by contrasting population genetic patterns. In the west, populations may not have persisted through sea-level lows because shallow seabed was very limited. Shallow genetic structure, weak expansion signals and a westward cline in genetic diversity indicate a postglacial recolonization of the western part of the range from a more eastern location following sea-level rise. In the east, shallow seabed persisted during the LGM and increased considerably after the flooding of large bays, which resulted in strong demographic expansions, deeper genetic structure and higher genetic diversity. This study suggests that postglacial flooding with rising sea levels produced locally variable signatures in colonizing populations.

Phylogeographic analyses of a migratory freshwater fish (Megalobrama terminalis) reveal a shallow genetic structure and pronounced effects of sea-level changes

Different types of fish taxa generally seem to present diverse phylogeographic structures and show different responses to environmental changes. In southern China, however, a large number of phylogeographic studies have been employed for small and/or benthic fish species, while phylogeographic patterns of larger and migratory species are not well understood. In this study, we chose Megalobrama terminalis, an endemic cyprinid from Southern China with a median size and relatively high migratory potential as a candidate to explore the phylogeographic structure and to seek the relevant driving forces using a multilocus approach (two mitochondrial markers and three nuclear markers). Our results show that M. terminalis exhibits three mtDNA genetic groups, each presenting genetic structure to the local geography, but the differentiation level was much weaker than that of small and/or benthic fish species reported in earlier studies. Nuclear loci did not observe pronounced genetic group subdivision, but did examine noteworthy genetic differentiation between the Hainan Island population and the mainland populations. Divergence time estimation and demographic analyses suggest that sea-level changes associated with the Late Pleistocene glacial cycles have exerted pronounced effects on the divergence of the three groups and the expansion of M. terminalis populations. Above all, our study provides important knowledge that can be used to improve monitoring programs and develop a conservation strategy for this endemic organism.

Aquatic Adaptation and Depleted Diversity: A Deep Dive into the Genomes of the Sea Otter and Giant Otter

Despite its recent invasion into the marine realm, the sea otter (Enhydra lutris) has evolved a suite of adaptations for life in cold coastal waters, including limb modifications and dense insulating fur. This uniquely dense coat led to the near-extinction of sea otters during the 18th-20th century fur trade and an extreme population bottleneck. We used the de novo genome of the southern sea otter (E. l. nereis) to reconstruct its evolutionary history, identify genes influencing aquatic adaptation, and detect signals of population bottlenecks. We compared the genome of the southern sea otter with the tropical freshwater-living giant otter (Pteronura brasiliensis) to assess common and divergent genomic trends between otter species, and with the closely related northern sea otter (E. l. kenyoni) to uncover population-level trends. We found signals of positive selection in genes related to aquatic adaptations, particularly limb development and polygenic selection on genes related to hair follicle development. We found extensive pseudogenization of olfactory receptor genes in both the sea otter and giant otter lineages, consistent with patterns of sensory gene loss in other aquatic mammals. At the population level, the southern sea otter and the northern sea otter showed extremely low genomic diversity, signals of recent inbreeding, and demographic histories marked by population declines. These declines may predate the fur trade and appear to have resulted in an increase in putatively deleterious variants that could impact the future recovery of the sea otter.

Living on the edge: Exploring the role of coastal refugia in the Alexander Archipelago of Alaska

Although islands are of long-standing interest to biologists, only a handful of studies have investigated the role of climatic history in shaping evolutionary diversification in high-latitude archipelagos. In this study of the Alexander Archipelago (AA) of Southeast Alaska, we address the impact of glacial cycles on geographic genetic structure for three mammals co-distributed along the North Pacific Coast. We examined variation in mitochondrial and nuclear loci for long-tailed voles (Microtus longicaudus), northwestern deermice (Peromyscus keeni), and dusky shrews (Sorex monticola), and then tested hypotheses derived from Species Distribution Models, reconstructions of paleoshorelines, and island area and isolation. In all three species, we identified paleoendemic clades that likely originated in coastal refugia, a finding consistent with other paleoendemic lineages identified in the region such as ermine. Although there is spatial concordance at the regional level for endemism, finer scale spatial and temporal patterns are less clearly defined. Demographic expansion across the region for these distinctive clades is also evident and highlights the dynamic history of Late Quaternary contraction and expansion that characterizes high-latitude species.