Seascape genomics of common dolphins (Delphinus delphis) reveals adaptive diversity linked to regional and local oceanography

Background

High levels of standing genomic variation in wide-ranging marine species may enhance prospects for their long-term persistence. Patterns of connectivity and adaptation in such species are often thought to be influenced by spatial factors, environmental heterogeneity, and oceanographic and geomorphological features. Population-level studies that analytically integrate genome-wide data with environmental information (i.e., seascape genomics) have the potential to inform the spatial distribution of adaptive diversity in wide-ranging marine species, such as many marine mammals. We assessed genotype-environment associations (GEAs) in 214 common dolphins (Delphinus delphis) along > 3000 km of the southern coast of Australia.

Results

We identified 747 candidate adaptive SNPs out of a filtered panel of 17,327 SNPs, and five putatively locally-adapted populations with high levels of standing genomic variation were disclosed along environmentally heterogeneous coasts. Current velocity, sea surface temperature, salinity, and primary productivity were the key environmental variables associated with genomic variation. These environmental variables are in turn related to three main oceanographic phenomena that are likely affecting the dispersal of common dolphins: (1) regional oceanographic circulation, (2) localised and seasonal upwellings, and (3) seasonal on-shelf circulation in protected coastal habitats. Signals of selection at exonic gene regions suggest that adaptive divergence is related to important metabolic traits.

Conclusion

To the best of our knowledge, this represents the first seascape genomics study for common dolphins (genus Delphinus). Information from the associations between populations and their environment can assist population management in forecasting the adaptive capacity of common dolphins to climate change and other anthropogenic impacts.

Spatial and temporal variations in female size at maturity of a Southern Rock Lobster (Jasus edwardsii) population: A likely response to climate change

The size at which sexual maturity is reached is a key population parameter used to guide the setting of minimum legal size limits in fisheries. Understanding spatial and temporal variations in size at maturity is fundamental to management because the relationship between size at maturity and minimum legal size limits affects the fraction of the mature population biomass that is harvested, and resulting egg production, larval settlement and recruitment. This study measured the size at maturity of female Southern Rock Lobster (Jasus edwardsii) across South Australia between 1991 and 2015 in relation to known oceanographic characteristics, surface and subsurface temperature data, and relative changes in lobster abundance. There was pronounced north to south spatial variation in estimates of size at maturity. Larger average size at maturity was recorded in warmer north-western areas of the fishery relative to the cooler waters of the south-east. Estimates of size at maturity also differed over 25 years across the fishery. However, the nature of temporal responses varied spatially, and were more consistent with variations in surface and subsurface water temperature at local-scales than changes in lobster density. In the well-mixed waters of the north-western, western and south-eastern parts of the fishery, relatively high rates of increase in sea-surface temperature and size at maturity were recorded since 1991, indicating that size at maturity may be responding to ocean warming associated with global climate change. In more central parts of the fishery, contrasting temporal signals in sea-surface temperature (positive) and bottom temperature (negative) indicated increases in upwelling strength over the study period, and formation of a bottom cold pool below a warm surface layer, with corresponding decreases in size at maturity recorded. The spatio-temporal changes in size at maturity measured in this study highlight the need for oceanographic information to be integrated into future stock assessment models to enhance harvest strategy development, allow timely adaptive management decisions and increase the resilience of fisheries to the impacts of climate change.

Surface ocean carbon dioxide variability in South Pacific boundary currents and Subantarctic waters

To improve estimates of the long-term response of the marine carbon system to climate change a better understanding of the seasonal and interannual variability is needed. We use high-frequency multi-year data at three locations identified as climate change hotspots: two sites located close to South Pacific boundary currents and one in the Subantarctic Zone (SAZ). We investigate and identify the main drivers involved in the seasonal an interannual (2012-2016) variability of the carbon system. The seasonal variability at boundary current sites is temporally different and highly controlled by sea surface temperature. Advection processes also play a significant role on the monthly changes of the carbon system at the western boundary current site. The interannual variability at these sites most likely responds to long-term variability in oceanic circulation ultimately related to climatic indices such as the El Niño Southern Oscillation, the Pacific Decadal Oscillation and the Southern Annular Mode (SAM). In the SAZ, advection and entrainment processes drive most of the seasonality, augmented by the action of biological processes in spring. Given the relevance of advection and entrainment processes at SAZ, the interannual variability is most probably modulated by changes in the regional winds linked to the variability of the SAM.

Baseline characterization of aerobic hydrocarbon degrading microbial communities in deep-sea sediments of the Great Australian Bight, Australia

Exploratory drilling for deep-sea oil and gas resources is planned for the Great Australian Bight (GAB). There is scant knowledge of the region's benthic ecosystems and no baseline information of the region's indigenous oil degrading bacteria. To address this knowledge gap, we used next generation sequencing (NGS) of three marker genes (alkB, c23o and pmoA) to detect and characterize the microbial communities capable of aerobic hydrocarbon degradation. Unique, highly novel microbial communities capable of degrading hydrocarbons occur in surface sediments at depths between 200 and 2800 m. Clustering at 97% demonstrated differences in community structure with depth, changing most markedly between 400 and 1000 m depth on the continental slope, and identified putative functional 'ecotypes' related to depth. Observed differences in community structure showed strong correlations with temperature, other physicochemical properties of the overlying water column and are further modulated by differences in sediment grain size. This study provides important baseline data on hydrocarbon degrading microbial communities prior to the start of petroleum resource extraction. Our data will inform future ecological monitoring of the GAB deep-sea ecosystem.

Migration dynamics of juvenile southern bluefin tuna

Large scale migrations are a key component of the life history of many marine species. We quantified the annual migration cycle of juvenile southern bluefin tuna (Thunnus maccoyii; SBT) and spatiotemporal variability in this cycle, based on a multi-decadal electronic tagging dataset. Behaviour-switching models allowed for the identification of cohesive areas of residency and classified the temporal sequence of movements within a migration cycle from austral summer foraging grounds in the Great Australian Bight (GAB) to winter foraging grounds in the Indian Ocean and Tasman Sea and back to the GAB. Although specific regions within the Indian Ocean were frequented, individuals did not always return to the same area in consecutive years. Outward migrations from the GAB were typically longer than return migrations back to the GAB. The timing of individual arrivals to the GAB, which may be driven by seasonality in prey availability, was more cohesive than the timing of departures from the GAB, which may be subject to the physiological condition of SBT. A valuable fishery for SBT operates in the GAB, as do a number of scientific research programs designed to monitor SBT for management purposes; thus, understanding SBT migration to and from the area is of high importance to a number of stakeholders.

Macrofaunal Patterns in and around du Couedic and Bonney Submarine Canyons, South Australia

Two South Australian canyons, one shelf-incising (du Couedic) and one slope-limited (Bonney) were compared for macrofaunal patterns on the shelf and slope that spanned three water masses. It was hypothesized that community structure would (H1) significantly differ by water mass, (H2) show significant regional differences and (H3) differ significantly between interior and exterior of each canyon. Five hundred and thirty-one species of macrofauna ≥1 mm were captured at 27 stations situated in depth stratified transects inside and outside the canyons from 100 to1500 m depth. The macrofauna showed a positive relationship to depth in abundance, biomass, species richness and community composition while taxonomic distinctness and evenness remained high at all depths. Biotic variation on the shelf was best defined by variation in bottom water primary production while sediment characteristics and bottom water oxygen, temperature and nutrients defined biotic variation at greater depth. Community structure differed significantly (p<0.01) among the three water masses (shelf-flowing South Australian current, upper slope Flinders current and lower slope Antarctic Intermediate Water) (H1). Although community differences between the du Couedic and Bonney regions were marginally above significance at p = 0.05 (H2), over half of the species captured were unique to each region. This supports the evidence from fish and megafaunal distributions that the du Couedic and Bonney areas are in different bioregions. Overall, the canyon interiors were not significantly different in community composition from the exterior (H3). However, both canyons had higher abundance and/or biomass, increased species dominance, different species composition and coarser sediments near the canyon heads compared to outside the canyons at the same depth (500 m), suggestive of heightened currents within the canyons that influence community composition there. At 1000–1500 m, the canyon interiors were depauperate, typical of V-shaped canyons elsewhere. The large number of species captured, given the relatively low sampling effort and focus on the larger macrofauna, support previous studies that identify the South Australian coast as a high biodiversity area.

Cetacean strandings in South Australia (1881–2008)

Long-term monitoring of cetacean strandings is essential for good management. This study updates previous summaries for South Australia by adding up to 20 years of comprehensive data, including results of necropsy examinations. A total of 1078 records were examined. Thirty-one species were recorded: 9 (7% of records) mysticetes, 22 (88%) odontocetes and the rest (5%) unidentified. The number of species new to South Australia did not reach an asymptote, with potential for at least five additional species. Small cetaceans were more frequently recorded after 1990, possibly due to increased reporting effort. Stranding records increased markedly after 1970. Records for all species occurred year-round. Beaked whales stranded primarily during January–April, baleen whales during July–January and common dolphins during February–May. Geographic hotspots were identified and related to upwelling and reporting effort. A necropsy program since 1990 resulted in 315 of 856 records being assigned to a circumstance of death, with anthropogenic circumstances accounting for 42% of these. Known Entanglement (21%, 66 of 315) and Probable Entanglement (12%, 37 of 315) were the most recorded anthropogenic circumstances of death. Future research correlating strandings with oceanographic/climatic conditions may help to explain the documented patterns but first the effects of reporting effort need to be accounted for.

Steady as he goes: at-sea movement of adult male Australian sea lions in a dynamic marine environment

The southern coastline of Australia forms part of the worlds' only northern boundary current system. The Bonney Upwelling occurs every austral summer along the south-eastern South Australian coastline, a region that hosts over 80% of the worlds population of an endangered endemic otariid, the Australian sea lion. We present the first data on the movement characteristics and foraging behaviour of adult male Australian sea lions across their South Australian range. Synthesizing telemetric, oceanographic and isotopic datasets collected from seven individuals enabled us to characterise individual foraging behaviour over an approximate two year time period. Data suggested seasonal variability in stable carbon and nitrogen isotopes that could not be otherwise explained by changes in animal movement patterns. Similarly, animals did not change their foraging patterns despite fine-scale spatial and temporal variability of the upwelling event. Individual males tended to return to the same colony at which they were tagged and utilized the same at-sea regions for foraging irrespective of oceanographic conditions or time of year. Our study contrasts current general assumptions that male otariid life history strategies should result in greater dispersal, with adult male Australian sea lions displaying central place foraging behaviour similar to males of other otariid species in the region.

Biophysical factors affecting the distribution of demersal fish around the head of a Submarine Canyon off the Bonney Coast, South Australia

We sampled the demersal fish community of the Bonney Canyon, South Australia at depths (100–1,500 m) and locations that are poorly known. Seventy-eight species of demersal fish were obtained from 12 depth-stratified trawls along, and to either side, of the central canyon axis. Distributional patterns in species richness and biomass were highly correlated. Three fish assemblage groupings, characterised by small suites of species with narrow depth distributions, were identified on the shelf, upper slope and mid slope. The assemblage groupings were largely explained by depth (ρw = 0.78). Compared to the depth gradient, canyon-related effects are weak or occur at spatial or temporal scales not sampled in this study. A conceptual physical model displayed features consistent with the depth zonational patterns in fish, and also indicated that canyon upwelling can occur. The depth zonation of the fish assemblage was associated with the depth distribution of water masses in the area. Notably, the mid-slope community (1,000 m) coincided with a layer of Antarctic Intermediate Water, the upper slope community (500 m) resided within the core of the Flinders Current, and the shelf community was located in a well-mixed layer of surface water (<450 m depth).

Spawning dynamics and biomass estimates of an anchovy Engraulis australis population in contrasting gulf and shelf environments

The spawning biomass of Australian anchovy Engraulis australis in gulf and shelf waters of South Australia was compared using the daily egg production method (DEPM). The total survey area was 128 700 km2 with recorded spawning areas in gulf and shelf waters of 4898 and 44 618 km2, respectively. High egg densities in the warm, shallow gulf waters were produced by small, young (<1 year old) E. australis that spawned relatively small batches of eggs (c. 855) approximately every 3 days. In cooler, deeper shelf waters, where larger, older E. australis are found, lower egg densities occurred despite individuals producing much larger batches of eggs (c. 15,572) approximately every 7 days. In shelf waters, the highest densities were recorded at inshore sampling stations. Spawning appeared to peak between 0000 and 0100 hours. Females were more abundant than males in samples from both gulf and shelf waters with sex ratios of 0.61 and 0.56, respectively. The spawning biomass of E. australis in shelf waters was 101 522 t, whereas the estimate for gulf waters was 25 374 t. Due to the differences in mean size of the spawning females, however, c. 6x10(9)E. australis were present in each region. The results support the hypothesis that variability in habitat conditions may directly influence E. australis reproduction. A large reserve of young fish in the relatively stable gulf environment may increase the resilience of the E. australis population in South Australia to unfavourable interannual changes in offshore environmental conditions.