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

The biogeographic and genomic signatures of dynamic river networks for terrestrial species in Amazonia

Amazonia contains Earth's largest freshwater basin, largest contiguous stretch of tropical forest, and most species-rich terrestrial biota on Earth. Rivers are key geographic features that drive diversification of the Amazonian biota, but they are also dynamic, which challenges their role as long-term barriers to dispersal and gene flow. The impacts of such river dynamics on organismal evolution have only recently been explored in detail. Here I examine biodiversity patterns and processes in Amazonia to elucidate how taxa diversify in the context of river network dynamics. I borrow the River Capture Hypothesis from ichthyology, and draw on evidence from speciation genomics, hybrid zones, and community assembly to demonstrate the effects of river network evolution on biodiversification. The idea is simple: populations of organisms whose dispersal is restricted by rivers become semi-isolated by rivers. Drift and selection against introgression drive divergence, but as rivers move, previously isolated populations come into secondary contact, facilitating lineage fusions or the migration of hybrid zones to other rivers. The basin's unique macroecological patterns and rich biota thus may have resulted from repeated divergences, lineage fusions, and range expansions around a network of non-stationary extrinsic barriers with variable results depending on the degree of intrinsic reproductive isolation that accumulates during this process. The evolutionary consequences of dynamic landscapes extend beyond Amazonia as "fission-fusion-fission" cycles modulate the diversification and spatial patterning of life on Earth in general.

The Intertidal North-South Split: Oceanographic Features and Life History Shape the Phylogeography of Chiton Acanthochitona rubrolineata

The genetic structure and demographic history of marine organisms are shaped by a variety of factors including biological and ecological characteristics, ocean currents, and the palaeogeological effects of sea-level fluctuations. Here we present a comprehensive method combining population genomics, laboratory experiments, and ocean modelling in 13 populations of the chiton Acanthochitona rubrolineata along the coast of China. Based on demographic and population genomic analyses, significant divergence was observed between the Northern and Southern population groups, which are separated by the Yangtze River Estuary. The numerical circulation model simulation showed that gene flow and population connectivity were strongly influenced by ocean currents and the larval dispersal ability of chiton A. rubrolineata. These data thus clearly revealed the presence of two separately evolving lineages in chiton-A. rubrolineata northern and A. rubrolineata southern. Our study highlights that a robust understanding of organisms in the intertidal zone requires a comprehensive consideration of factors that influence gene flow and genetic structure, including the life-history traits, coastal currents, geographic isolation, and habitat suitability. The life history of marine organisms, together with local oceanographic features, could ultimately drive the population divergence and lead to speciation. These findings provide a guideline for future analyses of non-model and potentially threatened species and will aid in the conservation of biodiversity.

Phylogenomics of Bivalvia Using Ultraconserved Elements Reveal New Topologies for Pteriomorphia and Imparidentia

Despite significant advances in phylogenetics over the past decades, the deep relationships within Bivalvia (phylum Mollusca) remain inconclusive. Previous efforts based on morphology or several genes have failed to resolve many key nodes in the phylogeny of Bivalvia. Advances have been made recently using transcriptome data, but the phylogenetic relationships within Bivalvia historically lacked consensus, especially within Pteriomorphia and Imparidentia. Here, we inferred the relationships of key lineages within Bivalvia using matrices generated from specifically designed ultraconserved elements (UCEs) with 16 available genomic resources and 85 newly sequenced specimens from 55 families. Our new probes (Bivalve UCE 2k v.1) for target sequencing captured an average of 849 UCEs with 1085 bp in mean length from in vitro experiments. Our results introduced novel schemes from 6 major clades (Protobranchina, Pteriomorphia, Palaeoheterodonta, Archiheterodonta, Anomalodesmata, and Imparidentia), though some inner nodes were poorly resolved, such as paraphyletic Heterodonta in some topologies potentially due to insufficient taxon sampling. The resolution increased when analyzing specific matrices for Pteriomorphia and Imparidentia. We recovered 3 Pteriomorphia topologies different from previously published trees, with the strongest support for ((Ostreida + (Arcida + Mytilida)) + (Pectinida + (Limida + Pectinida))). Limida were nested within Pectinida, warranting further studies. For Imparidentia, our results strongly supported the new hypothesis of (Galeommatida + (Adapedonta + Cardiida)), while the possible non-monophyly of Lucinida was inferred but poorly supported. Overall, our results provide important insights into the phylogeny of Bivalvia and show that target enrichment sequencing of UCEs can be broadly applied to study both deep and shallow phylogenetic relationships.

Mitochondrial DNA patterns describe the evolutionary history of the bonnethead shark Sphyrna tiburo (Linneus 1758) complex in the western Atlantic Ocean

The apparent lack of physical barriers in the marine realm has created the conception that many groups have a constant gene flow. However, changes in ocean circulation patterns, glacial cycles, temperature, and salinity gradients are responsible for vicariant events in many fish species, including sharks. The bonnethead shark, Sphyrna tiburo, is an endangered small coastal shark species. Although considerable efforts have recently been undertaken, little remains known about the possible biogeographic scenario that can explain its actual distribution within the western Atlantic (WA). Here, we used 599 mitochondrial sequences to assess the phylogeographic structure and implement Bayesian phylogenetic analyses to obtain divergence times and reconstruct the ancestral geographic range. This allowed us to infer processes responsible for the diversification of S. tiburo into major divergent lineages. Our results indicated that S. tiburo in the WA represents three independent lineages, with Brazilian samples differentiated into a distinct genetic cluster. The posterior probability of ancestral range analysis indicated that the species likely originated in the northern region (Carolina Province and the southern Gulf of Mexico), where it colonized southward through the uplifting of the Central American Isthmus (CAI). The Northern and Caribbean genetic clusters appear to have arisen from the intensification of the Loop Current, which currently flows northward passing the Yucatan Peninsula, Gulf of Mexico, and east Florida. Following initial colonization, the Northeastern Brazil group differentiated from the Caribbean region due to the sediment and freshwater discharge of the Amazon-Orinoco Plume. Thus, the evolutionary history of the S. tiburo complex can be explained by a combination of dispersal and vicariance events that occurred over the last ~5 million years (MY). We established and confirmed the species and population limits, demonstrating that the Amazon-Orinoco Plume constitutes a significant dispersal barrier for coastal sharks. Finally, we discuss some recommendations for the conservation of the bonnethead shark.

An assemblage-level comparison of genetic diversity and population genetic structure between island and mainland ant populations

Island biotas provide unparalleled opportunities to examine evolutionary processes. Founder effects and bottlenecks, e.g., typically decrease genetic diversity in island populations, while selection for reduced dispersal can increase population structure. Given that support for these generalities mostly comes from single-species analyses, assemblage-level comparisons are needed to clarify how (i) colonization affects the gene pools of interacting insular organisms, and (ii) patterns of genetic differentiation vary within assemblages of organisms. Here, we use genome-wide sequence data from ultraconserved elements (UCEs) to compare the genetic diversity and population structure of mainland and island populations of nine ant species in coastal southern California. As expected, island populations (from Santa Cruz Island) had lower expected heterozygosity and Watterson's theta compared to mainland populations (from the Lompoc Valley). Island populations, however, exhibited smaller genetic distances among samples, indicating less population subdivision. Within the focal assemblage, pairwise Fst values revealed pronounced interspecific variation in mainland-island differentiation, which increases with gyne body size. Our results reveal population differences across an assemblage of interacting species and illuminate general patterns of insularization in ants. Compared to single-species studies, our analysis of nine conspecific population pairs from the same island-mainland system offers a powerful approach to studying fundamental evolutionary processes.

Range-wide population genomics of common seadragons shows secondary contact over a former barrier and insights on illegal capture

BACKGROUND

Common seadragons (Phyllopteryx taeniolatus, Syngnathidae) are an emblem of the diverse endemic fauna of Australia's southern rocky reefs, the newly recognized "Great Southern Reef." A lack of assessments spanning this global biodiversity hotspot in its entirety is currently hampering an understanding of the factors that have contributed to its diversity. The common seadragon has a wide range across Australia's entire temperate south and includes a geogenetic break over a former land bridge, which has called its status as a single species into question. As a popular aquarium display that sells for high prices, common seadragons are also vulnerable to illegal capture.

RESULTS

Here, we provide range-wide nuclear sequences (986 variable Ultraconserved Elements) for 198 individuals and mitochondrial genomes for 140 individuals to assess species status, identify genetic units and their diversity, and trace the source of two poached individuals. Using published data of the other two seadragon species, we found that lineages of common seadragons have diverged relatively recently (< 0.63 Ma). Within common seadragons, we found pronounced genetic structure, falling into three major groups in the western, central, and eastern parts of the range. While populations across the Bassian Isthmus were divergent, there is also evidence for secondary contact since the passage opened. We found a strong cline of genetic diversity from the range center tapering symmetrically towards the range peripheries. Based on their genetic similarities, the poached individuals were inferred to have originated from around Albany in southwestern Australia.

CONCLUSIONS

We conclude that common seadragons constitute a single species with strong geographic structure but coherence through gene flow. The low genetic diversity on the east and west coasts is concerning given that these areas are projected to face fast climate change. Our results suggest that in addition to their life history, geological events and demographic expansions have all played a role in shaping populations in the temperate south. These insights are an important step towards understanding the historical determinants of the diversity of species endemic to the Great Southern Reef.

The genetic basis of the leafy seadragon's unique camouflage morphology and avenues for its efficient conservation derived from habitat modeling

The leafy seadragon certainly is among evolution's most "beautiful and wonderful" species aptly named for its extraordinary camouflage mimicking its coastal seaweed habitat. However, limited information is known about the genetic basis of its phenotypes and conspicuous camouflage. Here, we revealed genomic signatures of rapid evolution and positive selection in core genes related to its camouflage, which allowed us to predict population dynamics for this species. Comparative genomic analysis revealed that seadragons have the smallest olfactory repertoires among all ray-finned fishes, suggesting adaptations to the highly specialized habitat. Other positively selected and rapidly evolving genes that serve in bone development and coloration are highly expressed in the leaf-like appendages, supporting a recent adaptive shift in camouflage appendage formation. Knock-out of bmp6 results in dysplastic intermuscular bones with a significantly reduced number in zebrafish, implying its important function in bone formation. Global climate change-induced loss of seagrass beds now severely threatens the continued existence of this enigmatic species. The leafy seadragon has a historically small population size likely due to its specific habitat requirements that further exacerbate its vulnerability to climate change. Therefore, taking climate change-induced range shifts into account while developing future protection strategies.

Historic and contemporary biogeographic perspectives on range-wide spatial genetic structure in a widespread seagrass

Historical and contemporary processes drive spatial patterns of genetic diversity. These include climate-driven range shifts and gene flow mediated by biogeographical influences on dispersal. Assessments that integrate these drivers are uncommon, but critical for testing biogeographic hypotheses. Here, we characterize intraspecific genetic diversity and spatial structure across the entire distribution of a temperate seagrass to test marine biogeographic concepts for southern Australia. Predictive modeling was used to contrast the current Posidonia australis distribution to its historical distribution during the Last Glacial Maximum (LGM). Spatial genetic structure was estimated for 44 sampled meadows from across the geographical range of the species using nine microsatellite loci. Historical and contemporary distributions were similar, with the exception of the Bass Strait. Genetic clustering was consistent with the three currently recognized biogeographic provinces and largely consistent with the finer-scale IMCRA bioregions. Discrepancies were found within the Flindersian province and southwest IMCRA bioregion, while two regions of admixture coincided with transitional IMCRA bioregions. Clonal diversity was highly variable but positively associated with latitude. Genetic differentiation among meadows was significantly associated with oceanographic distance. Our approach suggests how shared seascape drivers have influenced the capacity of P. australis to effectively track sea level changes associated with natural climate cycles over millennia, and in particular, the recolonization of meadows across the Continental Shelf following the LGM. Genetic structure associated with IMCRA bioregions reflects the presence of stable biogeographic barriers, such as oceanic upwellings. This study highlights the importance of biogeography to infer the role of historical drivers in shaping extant diversity and structure.

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.

Species and genetic diversity relationships in benthic macroinvertebrate communities along a salinity gradient

BACKGROUND

Species- and genetic diversity can change in parallel, resulting in a species-genetic diversity correlation (SGDC) and raising the question if the same drivers influence both biological levels of diversity. The SGDC can be either positive or negative, depending on whether the species diversity and the genetic diversity of the measured species respond in the same or opposite way to drivers. Using a traditional species diversity approach together with ultra-conserved elements and high throughput sequencing, we evaluated the SGDCs in benthic macrofauna communities in the Baltic Sea, a geologically young brackish water sea characterised by its steep salinity gradient and low species richness. Assessing SGDCs from six focal marine invertebrate species from different taxonomic groups and with differing life histories and ecological functions on both a spatial and temporal scale gives a more comprehensive insight into the community dynamics of this young ecosystem and the extrinsic factors that might drive the SGDCs.

RESULTS

No significant correlations between species diversity and genetic diversity were found for any of the focal species. However, both negative and positive trends of SGDCs for the individual focal species were observed. When examining the environmental drivers, no common trends between the species were found, even when restricting the analysis to specific taxonomic classes. Additionally, there were no common environmental factors driving the diversity relationships for species sharing the same SGDC trend (positive or negative). Local population dynamics, together with the invasion history of the individual species and their unique adaptation to the distinctive environment of the Baltic Sea, are expected to be of major influence on the outcome of the SGDCs.

CONCLUSIONS

The present results highlight the importance of assessing SGDCs using multiple species, not just a single indicator species. This emphasises a need to pay attention to the ecology and life history of the focal species. This study also provides insight into the large differences in both patterns and drivers of genetic diversity, which is important when including genetic biodiversity in conservation plans. We conclude that the effects of environmental and biological factors and processes that affects diversity patterns at both the community and genetic levels are likely species dependent, even in an environment such as the Baltic Sea with strong environmental gradients.

Paleogene emergence and evolutionary history of the Amazonian fossorial fish genus Tarumania (Teleostei: Tarumaniidae)

Tarumania walkerae is a rare fossorial freshwater fish species from the lower Rio Negro, Central Amazonia, composing the monotypic and recently described family Tarumaniidae. The family has been proposed as the sister group of Erythrinidae by both morphological and molecular studies despite distinct arrangements of the superfamily Erythrinoidea within Characiformes. Recent phylogenomic studies and time-calibrated analyses of characoid fishes have not included specimens of Tarumania in their analyses. We obtained genomic data for T. walkerae and constructed a phylogeny based on 1795 nuclear loci with 488,434 characters of ultraconserved elements (UCEs) for 108 terminals including specimens of all 22 characiform families. The phylogeny confirms the placement of Tarumaniidae as sister to Erythrinidae but differs from the morphological hypothesis in the placement of the two latter families as sister to the clade with Hemiodontidae, Cynodontidae, Serrasalmidae, Parodontidae, Anostomidae, Prochilodontidae, Chilodontidae, and Curimatidae. The phylogeny calibrated with five characoid fossils indicates that Erythrinoidea diverged from their relatives during the Late Cretaceous circa 90 Ma (108–72 Ma), and that Tarumania diverged from the most recent common ancestor of Erythrinidae during the Paleogene circa 48 Ma (66–32 Ma). The occurrence of the erythrinoid-like †Tiupampichthys in the Late Cretaceous–Paleogene formations of the El Molino Basin of Bolivia supports our hypothesis for the emergence of the modern Erythrinidae and Tarumaniidae during the Paleogene.

Leafy and weedy seadragon genomes connect genic and repetitive DNA features to the extravagant biology of syngnathid fishes

Seadragons are a remarkable lineage of teleost fishes in the family Syngnathidae, renowned for having evolved male pregnancy. Comprising three known species, seadragons are widely recognized and admired for their fantastical body forms and coloration, and their specific habitat requirements have made them flagship representatives for marine conservation and natural history interests. Until recently, a gap has been the lack of significant genomic resources for seadragons. We have produced gene-annotated, chromosome-scale genome models for the leafy and weedy seadragon to advance investigations of evolutionary innovation and elaboration of morphological traits in seadragons as well as their pipefish and seahorse relatives. We identified several interesting features specific to seadragon genomes, including divergent noncoding regions near a developmental gene important for integumentary outgrowth, a high genome-wide density of repetitive DNA, and recent expansions of transposable elements and a vesicular trafficking gene family. Surprisingly, comparative analyses leveraging the seadragon genomes and additional syngnathid and outgroup genomes revealed striking, syngnathid-specific losses in the family of fibroblast growth factors (FGFs), which likely involve reorganization of highly conserved gene regulatory networks in ways that have not previously been documented in natural populations. The resources presented here serve as important tools for future evolutionary studies of developmental processes in syngnathids and hold value for conservation of the extravagant seadragons and their relatives.

Pleistocene expansion, anthropogenic pressure and ocean currents: Disentangling the past and ongoing evolutionary history of Patella aspera Röding, 1798 in the archipelago of Madeira

AIMS

Rising sea-level following the Last Glacial Maximum lead to fragmentation of coastal limpet populations between islands of the Archipelago of Madeira. This fragmentation is reinforced by recent heavy exploitation reducing effective population size on Madeira Island. We use the limpet P. aspera to understand how the role of processes at different time scales (i.e. changes in the sea level and overexploitation) can influence the genetic composition of an extant species, relating these processes to reproductive phenology and seasonal shifts in ocean currents.

LOCATION

Madeira Island, Porto Santo and Desertas (Archipelago of Madeira, NE Atlantic Ocean).

TAXON

The limpet Patella aspera.

METHODS

Twelve microsatellite genetic markers were used. A power analysis was used to evaluate the power of the microsatellite markers to detect a signal of population differentiation. Long-term past migrations were assessed using a Bayesian Markov Montecarlo approach in the software MIGRATE-n to estimate mutation-scaled migration rates (M = m/μ; m, probability of a lineage immigrating per generation; μ, mutation rate). Two scenarios were evaluated using an Approximate Bayesian Computation (ABC) in the software DIYABC 2.1 (i) Scenario 1: considered a population scenario from a reduced Ne at time t₃ to a higher Ne at time t₂; and (ii) Scenario 2 considering a reduction of Ne from a time t₃ to a time t₂.

RESULTS

Colonization of the archipelago by Portuguese settlers six centuries ago probably led to an important decrease in the genetic diversity of the species (N_e). Contemporary gene flow strongly support a pattern of high asymmetric connectivity explained by the reproductive phenology of the species and spatio-temporal seasonal changes in the ocean currents. Spatio-temporal reconstructions using Bayesian methods, including coalescent and Approximate Bayesian Computation (ABC) approaches, suggest changes in the migration patterns from highly symmetric to highly asymmetric connectivity with subtle population differentiation as consequence of post-glacial maximum sea level rise during the Holocene.

MAIN CONCLUSIONS

Our results suggest that anthropogenic activity could have had serious effects on the genetic diversity of heavily exploited littoral species since the end of the Pleistocene, probably accelerating in recent years.

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.