Latitudinal discontinuity in thermal conditions along the nearshore of central-northern Chile

Oceanographical-driven dispersal and environmental variation explain genetic structure in an upwelling coastal ecosystem

The seascape comprises multiple environmental variables that interact with species biology to determine patterns of spatial genetic variation. The environment imposes spatially variable selective forces together with homogenizing and diverging drivers that facilitate or restrict dispersal, which is a complex, time-dependent process. Understanding how the seascape influences spatial patterns of genetic variation remains elusive, particularly in coastal upwelling systems. Here, we combine genome-wide SNP data, Lagrangian larval dispersal simulated over a hydrodynamic model, and ocean environmental information to quantify the relative contribution of ocean circulation and environmental heterogeneity as drivers of the spatial genetic structure of two congeneric intertidal limpets, Scurria scurra and S. araucana, along the central coast of Chile. We find that a genetic break observed in both limpet species coincides with a break in connectivity shown by the Lagrangian dispersal, suggesting that mean ocean circulation is an important seascape feature, in particular for S. scurra. For S. araucana, environmental variation appears as a better predictor of genetic structure than ocean circulation. Overall, our study shows broad patterns of seascape forcing on genetic diversity and contributes to our understanding of the complex ecological and evolutionary interactions along coastal upwelling systems.

Physiological resilience of intertidal chitons in a persistent upwelling coastal region

Current climate projections for mid-latitude regions globally indicate an intensification of wind-driven coastal upwelling due to warming conditions. The dynamics of mid-latitude coastal upwelling are marked by environmental variability across temporal scales, which affect key physiological processes in marine calcifying organisms and can impact their large-scale distribution patterns. In this context, marine invertebrates often exhibit phenotypic plasticity, enabling them to adapt to environmental change. In this study, we examined the physiological performance (i.e., metabolism, Thermal Performance Curves, and biomass and calcification rates) of individuals of the intertidal mollusk Chiton granosus, a chiton found from northern Peru to Cape Horn (5° to 55°S). Our spatial study design indicated a pattern of contrasting conditions among locations. The Talcaruca site, characterized by persistent upwelling and serving as a biogeographic break, exhibited lower pH and carbonate saturation states, along with higher pCO2, compared to the sites located to the north and south of this location (Huasco and Los Molles, respectively). In agreement with the spatial pattern in carbonate system parameters, long-term temperature records showed lower temperatures that changed faster over synoptic scales (1-15 days) at Talcaruca, in contrast to the more stable conditions at the sites outside the break. Physiological performance traits from individuals from the Talcaruca population exhibited higher values and more significant variability, along with significantly broader and greater warming tolerance than chitons from the Huasco and Los Molles populations. Moreover, marked changes in local abundance patterns over three years suggested population-level responses to the challenging environmental conditions at the biogeographic break. Thus, C. granosus from the Talcaruca upwelling zone represents a local population with wide tolerance ranges that may be capable of withstanding future upwelling intensification on the Southern Eastern Pacific coast and likely serving as a source of propagules for less adapted populations.

Hydrographic shifts in coastal waters reflect climate-driven changes in hydrological regimes across Northwestern Patagonia

Climate-driven changes in freshwater inputs have been shown to affect the structure and function of coastal ecosystems. We evaluated changes in the influence of river runoff on coastal systems of Northwestern Patagonia (NWP) over recent decades (1993-2021) by combined analysis of long-term streamflow time series, hydrological simulation, satellite-derived and reanalysis data on sea surface conditions (temperature, turbidity, and salinity). Significant decreases in minimum streamflow across a zone spanning six major river basins were evident at weekly, monthly, and seasonal scales. These changes have been most pronounced in mixed-regime northern basins (e.g., Puelo River) but appear to be progressing southward to rivers characterised by a nival regime. In the adjacent two-layer inner sea, reduced freshwater input corresponds with a shallower halocline and increased surface temperatures across northern Patagonia. Our results underscore the rapidly evolving influence of rivers on adjacent estuarine and coastal waters in NWP. We highlight the need for cross-ecosystem observation, forecasting, mitigation and adaptation strategies in a changing climate, together with corresponding adaptive basin management of systems that supply runoff to the coastal marine waters.

Impact of persistent barrier to gene flow and catastrophic events on red algae evolutionary history along the Chilean coast

Historical vicariance events, linked to the existence of stable physical barriers to gene flow, generate concordant genetic breaks in co-distributed species while stochastic processes (e.g., costal uplift) could cause species-specific genetic breaks as a result of local strong demographic bottlenecks or extinction. In Chile, previous studies show that the area of the 30°S-33°S could correspond to a stable barrier to gene flow that have affected the genetic structure of various algae and marine invertebrates. Here we sequenced two organellar genes (COI and rbcL) in four taxonomically accepted co-distributed red seaweeds species characterized by a low dispersal potential: Mazzaella laminarioides, M. membranacea, Asterfilopsis disciplinalis, and Ahnfeltiopsis vermicularis. Our results revealed the existence of ten strongly differentiated linages in the taxa studied. Strong genetic breaks, concordant in both space and time (divergence estimated to have occurred some 2.9-12.4 million years ago), were observed between taxa distributed across the 33°S. Conversely, in the Central/South part of the Chilean coast, the localization of the genetic breaks/sub-structure observed varied widely (36°S, 38°S, 39°S, and 40°S). These results suggest that a major historical vicariance event has modeled the genetic structure of several Chilean marine organisms in the north of the Chilean coast during the mid-Miocene, while more recent stochastic events and genetic drift could be the driving forces of genetic divergence/structuration in the central-southern part of the coast.

Modeling the dynamics of harmful algal bloom events in two bays from the northern Chilean upwelling system

The bays of Tongoy and Guanaqueros are located in the Humboldt Current system, where Argopecten purpuratus has been the subject of intense aquaculture development. These bays lie in one of the most productive marine ecosystems on Earth and are dominated by permanent coastal upwelling at Lengua de Vaca Point and Choros Point, one of the three upwelling centers on the Chilean coast. Significantly, this productive system experiences a high recurrence of harmful algal bloom (HAB) events. This paper examines 9-year (2010-2018) samples of three toxic microalgal species collected in different monitoring programs and research projects. During this period, nine HAB events were detected in Guanaqueros Bay and 14 in Tongoy Bay. Among these, three HAB events were produced simultaneously in both bays by Pseudo-nitzschia australis, and two events produced simultaneously were detected in one bay by Alexandrium spp. and the other by Dinophysis acuminata. Before El Niño 2015-16, there were more HAB events of longer duration by the three species. Since El Niño, the number and duration of events were reduced and only produced by P. australis. HAB events were simulated with the FVCOM model and a virtual particle tracker model to evaluate the dynamics of bays and their relationship with HAB events. The results showed retention in bays during the relaxation conditions of upwelling and low connectivity between bays, which explains why almost no simultaneous events were recorded.

Environmental DNA reveals temporal variation in mesophotic reefs of the Humboldt upwelling ecosystems of central Chile: Toward a baseline for biodiversity monitoring of unexplored marine habitats

Temperate mesophotic reef ecosystems (TMREs) are among the least known marine habitats. Information on their diversity and ecology is geographically and temporally scarce, especially in highly productive large upwelling ecosystems. Lack of information remains an obstacle to understanding the importance of TMREs as habitats, biodiversity reservoirs and their connections with better-studied shallow reefs. Here, we use environmental DNA (eDNA) from water samples to characterize the community composition of TMREs on the central Chilean coast, generating the first baseline for monitoring the biodiversity of these habitats. We analyzed samples from two depths (30 and 60 m) over four seasons (spring, summer, autumn, and winter) and at two locations approximately 16 km apart. We used a panel of three metabarcodes, two that target all eukaryotes (18S rRNA and mitochondrial COI) and one specifically targeting fishes (16S rRNA). All panels combined encompassed eDNA assigned to 42 phyla, 90 classes, 237 orders, and 402 families. The highest family richness was found for the phyla Arthropoda, Bacillariophyta, and Chordata. Overall, family richness was similar between depths but decreased during summer, a pattern consistent at both locations. Our results indicate that the structure (composition) of the mesophotic communities varied predominantly with seasons. We analyzed further the better-resolved fish assemblage and compared eDNA with other visual methods at the same locations and depths. We recovered eDNA from 19 genera of fish, six of these have also been observed on towed underwater videos, while 13 were unique to eDNA. We discuss the potential drivers of seasonal differences in community composition and richness. Our results suggest that eDNA can provide valuable insights for monitoring TMRE communities but highlight the necessity of completing reference DNA databases available for this region.

Heterozostera nigricaulis from the south-East Pacific coast of Chile: First insights into its physiology and growth

A short stretch (27°S and 30°S) along the coast of Chile is habitat for the seagrass Heterozostera nigricaulis. The seagrass is classified as endangered and grows only clonally, but there are no data on its physiology and growth. However, this information is important to gain insights into its acclimation potential and how disturbances may affect them. We therefore studied H. nigricaulis at 27° and 30°S, and determined their growth and physiology among seasons and depths over one year. Biomass was higher at 27° than at 30°S, and was always higher in summer than in autumn and winter. Increased photosynthesis supported growth in summer, and in winter carbonic anhydrase activity was in place to maintain these evergreen meadows. Our results suggest that these seagrass meadows are adapted to local conditions, which, together with their asexual reproduction, could make them more vulnerable to disturbance. Therefore, our results serve as a basis for future studies on seagrass growth dynamics, and are important for protection and management plans.

Spatio-temporal patterns of the crustacean demersal fishery discard from the south Humboldt Current System, based on scientific observer program (2014-2019)

This study summarises six years of spatio-temporal patterns of the discarded demersal community fauna recorded by onboard scientific observer program for both artisanal and industrial crustacean fisheries between 2014 and 2019, from mesophotic to aphotic depths (96 to 650 m) along the southern Humboldt Current System (28-38°S). In this period, one cold and two warm climatic events were observed during the austral summer 2014, 2015-2016 (ENSO Godzilla), and 2016-2017 (coastal ENSO), respectively. Satellite information showed that Chlorophyll-a concentration varied seasonally and latitudinally, associated with upwelling centres, while equatorial wind stress decreased southward of 36°S. Discards were composed of 108 species, dominated by finfish and molluscs. The Chilean hake Merluccius gayi was dominant and ubiquitous (occurrence, 95% of 9104 hauls), being the most vulnerable species of the bycatch. Three assemblages were identified: assemblage 1 (~200 m deep), dominated by flounders Hippoglossina macrops and lemon crabs Platymera gaudichaudii, assemblage 2 (~260 m deep), dominated by squat lobsters Pleuroncodes monodon and Cervimunida johni and assemblage 3 (~320 m depth), dominated by grenadiers Coelorinchus aconcagua and cardinalfish Epigonus crassicaudus. These assemblages were segregated by depth, and varied by year, and geographic zone. The latter represented changes in the width of the continental shelf, increasing southward of 36°S. Alpha-diversity indexes (richness, Shannon, Simpson, and Pielou) also varied with depth and latitude, with higher diversity in deeper continental waters (>300 m), between 2018-2019. Finally, at a spatial scale of tens of kilometres, and a monthly basis, interannual variations of biodiversity occurred in the demersal community. Surface sea temperature, chlorophyll-a, or wind stress did not correlate with discarded demersal fauna diversity of the crustacean fishery operating along central Chile.

Upwelling enhances seaweed nutrient quality, altering feeding behavior and growth rates in an intertidal sea urchin, Loxechinus albus

Upwelling systems deliver nutrient-rich water into coastal ecosystems, influencing primary productivity and potentially altering seaweed-herbivore interactions. Upwelling bottom-up effects on distinct trophic levels are well-known. However, their influence on seaweed biomolecules and on algae-herbivore interactions and growth are less known. The aim of this study was threefold: i) to compare physical-chemical characteristics and nutrient levels in the water of upwelling (U) and downwelling (DU) zones, ii) to quantify their influence on the content of protein and carbohydrates in seaweed tissues of representative U and DU locations, and iii) to experimentally assess their effect on the feeding behavior and growth of a prominent intertidal herbivore, the sea urchin Loxechinus albus. Waters from U zones showed lower temperatures and pH, and higher phosphate concentrations than those from downwelling zones. Similarly, the tissue of seaweeds from a U location had significantly more proteins and carbohydrates than those from a DU location. The origin location of the sea urchins had a significant influence on consumption and growth rates: in general, those coming from a site with U conditions consumed and grew more than those coming from DU conditions. The quality of the algae was a significant factor on consumption rates, although in the case of preference trials, this factor interacted with sea urchin origin location. Our results show that the availability and quality of the food in upwelling zones has an influence on herbivore-seaweed direct interactions. However, these interactions and the growth of the sea urchins were also related to the coastal site and conditions from which the sea urchins came from. These results are relevant considering the expected impact of climate change on the world's oceans, and the importance of U zones as thermal (cold water) refuges for marine ectotherms.

Long term environmental variability modulates the epigenetics of maternal traits of kelp crabs in the coast of Chile

The methylation of DNA is an environmentally inducible epigenetic mechanism reflecting the short-term ecological and environmental background of populations. Marine invertebrate populations, which spread along a latitudinal cline, are particularly suitable for profiling DNA methylation, due to the heterogenous environmental conditions experienced. We used the MSAP (Methylation Sensitive Amplified Polymorphism) technique to investigate the natural variation in DNA methylation of different female's tissues (muscle, gonads, and gills) and early-stage eggs from five populations of the kelp crab Taliepus dentatus, distributed along a latitudinal cline in the coast of Chile. We assessed whether, (1) the distribution of DNA methylation profiles can be associated with the temporal variability of long term (18 years) climatologies (sea surface temperature, turbidity and productivity) and (2) the epigenetic diversity of eggs is related to the population-level phenotypic variability of several maternal investment traits (egg volume, egg weight, egg lipids and fecundity). The DNA methylation of eggs correlated positively and negatively with the long term variability in productivity and sea surface temperature, respectively. Furthermore, the diversity of DNA methylation of eggs correlated positively with the population-level phenotypic variability of several maternal investment traits, suggesting a key role of epigenetic mechanisms in generating phenotypic variability at population level for this species. We provide evidence of a strong link between the temporal variability of long term climatologies with the epigenetic profiles of key early ontogenetic traits associated with the maternal investment of kelp crabs. These modulating mechanisms can hence contribute early to phenotypic variability at population levels in response to local and past environmental fluctuation.

Variations in phenotypic plasticity in a cosmopolitan copepod species across latitudinal hydrographic gradients

Studies assessing latitudinal variations in habitat conditions and phenotypic plasticity among populations yield evidence of the mechanisms governing differentiation in the potential to adapt to current/future habitat changes. The cosmopolitan copepod species Acartia tonsa thrives across ocean clines delimiting Seasonal (30–40° S) and Permanent (10–30° S) Upwelling coastal provinces established during the middle–late Pliocene (3.6–1.8 Ma) alongshore the South East Pacific (SEP), nowadays exhibiting contrasting variability features related to several ocean drivers (temperature, salinity, pH, and food availability). Latitudinal variation across the range of environmental conditions of the coastal provinces can contribute toward shaping divergent A. tonsa’s phenotypes, for example, through specific patterns of phenotypic plasticity in morphological and physiological traits and tolerance to environmental drivers. With the aim of contributing to the understanding of these adaptive processes in a relatively little studied oceanic region, here we compared the expression of parental (i.e., adult size, egg production, and ingestion rate) and offspring (i.e., egg size) traits in relation to variation in environmental habitat conditions across different cohorts of two distant (> 15° latitude) A. tonsa populations inhabiting estuarine and upwelling habitats located in the Seasonal and Permanent Upwelling province, respectively. Mean conditions and ranges of variability in the habitat conditions and phenotypic plasticity of parental and offspring traits within and among cohorts of A. tonsa populations varied significantly across the different examined regions (i.e., Seasonal vs. Permanent). We also found significant differences in the coupling of habitat variability and trait expression, suggesting that the differences in trait expressions might be related to habitat variability. The phenotypic divergence was translated to cohort-related patterns of trait trade-offs regulating reproduction and tolerance of egg production efficiency that can jointly determine the level of plasticity, genetic structure, or local adaptation. The current findings provide novel evidence of how divergent phenotypes might sustain A. tonsa populations across variable coastal provinces of the SEP.

Effect of UV-radiation on the physiology of the invasive green seaweed Codium fragile and its associated bacteria

Invasive species such as seaweeds often have a broad tolerance, allowing them to colonize novel habitats. During invasion, also new epibacteria can be formed on seaweeds, which have important chemo-ecological effects. Since UV-radiation (UVR) is one of the main factors affecting seaweeds and their epibacteria, we tested its effect on intertidal and subtidal thalli of the invasive seaweed Codium fragile from three sites and monitored photosynthesis, antioxidant activity and epibacteria. Exposure to UV-radiation resulted in photoinhibition with a subsequent low recovery in subtidal thalli from 23°S compared to 27°S and 30°S, which both showed a higher and almost complete recovery. However, a high antioxidant activity was present in all thalli, permitting to explain its relatively high tolerance to new environments. UV-radiation modified the composition of the epibacteria community by reducing its diversity and evenness. Our results showed that C. fragile responds plastic to variable UV-radiation (depending on site and water depth), which contributes to its high invasion potential.

Reef Fish Diversity Across the Temperate South Pacific Ocean

Patterns of species richness and their structuring forces at multiple scales provide a critical context for research efforts focusing on ecology, evolution, and conservation. Diversity gradients have been demonstrated in tropical reef fish, but corresponding patterns and mechanisms remain poorly understood in temperate regions. We conducted hierarchical (spatially nested) sampling of temperate reef fish faunas across > 140 degrees of longitude in the eastern and western South Pacific Ocean. Our sampling efforts spanned five distinct provinces: the Southeast Australian Shelf (SAS), Northern and Southern New Zealand (N-SNZ), Juan Fernandez and Desventuradas Islands (JFD), and the Warm Temperate Southeastern Pacific (WTPA). We evaluated (i) spatial variation in patterns of species richness and abundance (using Chao 1 index), and distribution of functional diversity (using several functional attributes: max body size, trophic groups, feeding guilds, trophic level, habitat use, gregariousness, and activity patterns) and (ii) scale-dependencies in these patterns. Species richness declined from west to east across the temperate South Pacific, but this pattern was detectable only across larger spatial scales. A functional redundancy index was significantly higher in the western South Australian Shelf at multiple scales, revealing that species contribute in equivalent ways to an ecosystem function such that one species may substitute for another. We also detected that patterns of variation in functional diversity differed from patterns of variation in species richness, and were also dependent on the spatial scale of analysis. Lastly, we identified that species’ traits are not equally distributed among reef fish assemblages, where some provinces are characterized by a distinct functional component within their reef fish assemblages. Planktivorous and schooling species, for instance, dominated the assemblages in the eastern Pacific, which is characterized by higher primary productivity and steep bathymetric slopes favoring these traits. Demersal and pairing behavior traits dominated the reef fish assemblages in western Pacific provinces (SAS, SNZ). We conclude that combining the identifies and species’ traits allow us to disentangle historical, biogeographic and environmental factors that structure reef fish fauna.

Climate change in the coastal ocean: shifts in pelagic productivity and regionally diverging dynamics of coastal ecosystems

Climate change has led to intensification and poleward migration of the Southeastern Pacific Anticyclone, forcing diverging regions of increasing, equatorward and decreasing, poleward coastal phytoplankton productivity along the Humboldt Upwelling Ecosystem, and a transition zone around 31° S. Using a 20-year dataset of barnacle larval recruitment and adult abundances, we show that striking increases in larval arrival have occurred since 1999 in the region of higher productivity, while slower but significantly negative trends dominate poleward of 30° S, where years of recruitment failure are now common. Rapid increases in benthic adults result from fast recruitment-stock feedbacks following increased recruitment. Slower population declines in the decreased productivity region may result from aging but still reproducing adults that provide temporary insurance against population collapses. Thus, in this region of the ocean where surface waters have been cooling down, climate change is transforming coastal pelagic and benthic ecosystems through altering primary productivity, which seems to propagate up the food web at rates modulated by stock-recruitment feedbacks and storage effects. Slower effects of downward productivity warn us that poleward stocks may be closer to collapse than current abundances may suggest.

Testing the Test: A Comparative Study of Marine Microbial Corrosion under Laboratory and Field Conditions

Microbially influenced corrosion (MIC) is an aggressive type of corrosion that occurs in aquatic environments and is sparked by the development of a complex biological matrix over a metal surface. In marine environments, MIC is exacerbated by the frequent variability in environmental conditions and the typically high diversity of microbial communities; hence, local and in situ studies are crucial to improve our understanding of biofilm composition, biological interactions among its members, MIC characteristics, and corrosivity. Typically, material performance and anticorrosion strategies are evaluated under controlled laboratory conditions, where natural fluctuations and gradients (e.g., light, temperature, and microbial composition) are not effectively replicated. To determine whether MIC development and material deterioration observed in the laboratory are comparable to those that occur under service conditions (i.e., field conditions), we used two testing setups, in the lab and in the field. Stainless steel (SS) AISI 316L coupons were exposed to southeastern Pacific seawater for 70 days using (i) acrylic tanks in a running seawater laboratory and (ii) an offshore mooring system with experimental frames immersed at two depths (5 and 15 m). Results of electrochemical evaluation, together with those of microbial community analyses and micrographs of formed biofilms, demonstrated that the laboratory setup provides critical information on the early biofilm development process (days), but the information gathered does not predict deterioration or biofouling of SS surfaces exposed to natural conditions in the field. Our results highlight the need to conduct further research efforts to understand how laboratory experiments may better reproduce field conditions where applications are to be deployed, as well as to improve our understanding of the role of eukaryotes and the flux of nutrients and oxygen in marine MIC events.

Seasonality of primary productivity affects coastal species more than its magnitude

While the importance of extreme conditions is recognised, patterns in species' abundances are often interpreted through average environmental conditions within their distributional range. For marine species with pelagic larvae, temperature and phytoplankton concentration are key variables. Along the south coast of South Africa, conspicuous spatial patterns in recruitment rates and the abundances of different mussel species exist, with focal areas characterized by large populations. We studied 15 years of sea surface temperature (SST) and chlorophyll-a (chl-a) satellite data, using spectral analyses to partition their temporal variability over ecologically relevant time periods, including seasonal (101 to 365 days) and intra-seasonal cycles (20 to 100 days). Adult cover and mussel recruitment were measured at 10 sites along the south coast and regression models showed that about 70% of the variability in recruitment and adult cover was explained by seasonal variability in chl-a, while mean annual chl-a and SST only explained 30% of the recruitment, with no significant effect for adult cover. SST and chl-a at two upwelling centres showed less predictable seasonal cycles during the second half of the study period with a significant cooling trend during austral autumn, coinciding with one of the mussel reproductive peaks. This likely reflects recent changes in the Agulhas Current, the world's largest western boundary current, which affects coastal ecosystems by driving upwelling. Similar mechanisms probably operate in other marine systems with the potential to affect the distribution patterns of key ecosystem engineers. We propose that variability in the characteristic timescales of environmental fluctuations can explain the spatial patterns of abundance of foundational species by affecting larval recruitment. This is especially important in a context of global and pervasive climate change, as shifts in the periodicity of environmental fluctuations appear to reflect large scale climatic teleconnections driven by anthropogenic forcing.

Effect of environmental history on the habitat-forming kelp Macrocystis pyrifera responses to ocean acidification and warming: a physiological and molecular approach

The capacity of marine organisms to adapt and/or acclimate to climate change might differ among distinct populations, depending on their local environmental history and phenotypic plasticity. Kelp forests create some of the most productive habitats in the world, but globally, many populations have been negatively impacted by multiple anthropogenic stressors. Here, we compare the physiological and molecular responses to ocean acidification (OA) and warming (OW) of two populations of the giant kelp Macrocystis pyrifera from distinct upwelling conditions (weak vs strong). Using laboratory mesocosm experiments, we found that juvenile Macrocystis sporophyte responses to OW and OA did not differ among populations: elevated temperature reduced growth while OA had no effect on growth and photosynthesis. However, we observed higher growth rates and NO₃^- assimilation, and enhanced expression of metabolic-genes involved in the NO₃^- and CO₂ assimilation in individuals from the strong upwelling site. Our results suggest that despite no inter-population differences in response to OA and OW, intrinsic differences among populations might be related to their natural variability in CO₂, NO₃^- and seawater temperatures driven by coastal upwelling. Further work including additional populations and fluctuating climate change conditions rather than static values are needed to precisely determine how natural variability in environmental conditions might influence a species' response to climate change.

Geographical variation of multiplex ecological networks in marine intertidal communities

Understanding the drivers of geographical variation in species distributions, and the resulting community structure, constitutes one of the grandest challenges in ecology. Geographical patterns of species richness and composition have been relatively well studied. Less is known about how the entire set of trophic and non-trophic ecological interactions, and the complex networks that they create by gluing species together in complex communities, change across geographical extents. Here, we compiled data of species composition and three types of ecological interactions occurring between species in rocky intertidal communities across a large spatial extent (~970 km of shoreline) of central Chile, and analyzed the geographical variability in these multiplex networks (i.e., comprising several interaction types) of ecological interactions. We calculated nine network summary statistics common across interaction types, and additional network attributes specific to each of the different types of interactions. We then investigated potential environmental drivers of this multivariate network organization. These included variation in sea surface temperature and coastal upwelling, the main drivers of productivity in nearshore waters. Our results suggest that structural properties of multiplex ecological networks are affected by local species richness and modulated by factors influencing productivity and environmental predictability. Our results show that non-trophic negative interactions are more sensitive to spatially structured temporal environmental variation than feeding relationships, with non-trophic positive interactions being the least labile to it. We also show that environmental effects are partly mediated through changes in species richness and partly through direct influences on species interactions, probably associated to changes in environmental predictability and to bottom-up nutrient availability. Our findings highlight the need for a comprehensive picture of ecological interactions and their geographical variability if we are to predict potential effects of environmental changes on ecological communities.

Signatures of local adaptation in the spatial genetic structure of the ascidian Pyura chilensis along the southeast Pacific coast

The highly heterogeneous Humboldt Current System (HCS) and the 30°S transition zone on the southeast Pacific coast, represent an ideal scenario to test the influence of the environment on the spatial genomic structure in marine near-shore benthic organisms. In this study, we used seascape genomic tools to evaluate the genetic structure of the commercially important ascidian Pyura chilensis, a species that exhibits a low larval transport potential but high anthropogenic dispersal. A recent study in this species recorded significant genetic differentiation across a transition zone around 30°S in putatively adaptive SNPs, but not in neutral ones, suggesting an important role of environmental heterogeneity in driving genetic structure. Here, we aim to understand genomic-oceanographic associations in P. chilensis along the Southeastern Pacific coast using two combined seascape genomic approaches. Using 149 individuals from five locations along the HCS, a total of 2,902 SNPs were obtained by Genotyping-By-Sequencing, of which 29–585 were putatively adaptive loci, depending on the method used for detection. In adaptive loci, spatial genetic structure was better correlated with environmental differences along the study area (mainly to Sea Surface Temperature, upwelling-associated variables and productivity) than to the geographic distance between sites. Additionally, results consistently showed the presence of two groups, located north and south of 30°S, which suggest that local adaptation processes seem to allow the maintenance of genomic differentiation and the spatial genomic structure of the species across the 30°S biogeographic transition zone of the Humboldt Current System, overriding the homogenizing effects of gene flow.

Context-Dependence in parasite effects on keyhole limpets

Parasites alter the reproductive performance of their hosts, limit their growth, and thereby modify the energy budget of these hosts. Experimental studies and theoretical models suggest that the outcome of the host-parasite interactions could be determined by ecological factors such as food availability levels in the local habitats. Nutrient inputs may affect the host's food resource availability with positive or negative effects on parasite infection rates and tolerance of infection, however this has not been specifically evaluated in natural systems. In this study, we evaluate the effects of parasitism by Proctoeces humboldti on body size, gonadosomatic index (GSI), and metabolic rate (oxygen consumption) of their second intermediate host Fissurella crassa limpets, under contrasting natural conditions of productivity (upwelling center vs upwelling shadow sites). Our results evidenced that parasitized limpets collected from the intertidal habitat influenced by coastal upwelling site showed greater shell length, muscular foot biomass and GSI as compared to non-parasitized limpets collected in the same site, and compared to parasitized and non-parasitized limpets collected from the sites under the influence of upwelling shadow conditions. Oxygen consumption was lower in parasitized limpets collected from the upwelling-influenced site than in the other groups, independent of age, suggesting reduced metabolic stress in infected individuals inhabiting these productive sites. Our results suggest that increased productivity in upwelling sites could mitigate the conflict for resources in the P. humboldti - F. crassa system, influencing where such interaction is found in the continuum between parasitism and mutualism. Since parasitism is ubiquitous in natural systems, and play important roles in ecological and evolutionary processes, it is important to analyze host-parasite interaction across a variety of ecological conditions, especially in biological conservation.

Connections Between Benthic Populations and Local Strandings of the Southern Bull Kelp Durvillaea Antarctica Along the Continental Coast of Chile1

Floating seaweeds are important dispersal vectors in marine ecosystems. However, the relationship between benthic populations and stranded seaweeds has received little attention. After detachment, a fraction of floating specimens returns to the shore, resulting in strandings that fluctuate in space and time. It has been hypothesized that the availability of stranded seaweeds is related to their benthic abundance on adjacent coasts. Using the large fucoid Durvillaea antarctica, we tested whether stranded biomasses are higher at sites with dense adjacent benthic populations. Benthic abundance of D. antarctica along the continental coast of Chile was estimated using three approximations: (i) availability of potentially suitable habitat (PSH), (ii) categorical visual abundance estimates in the field, and (iii) abundance measurements in the intertidal zone. Higher PSH for D. antarctica was observed between 31° S-32° S and 40° S-42° S than between 33° S and 39° S. Lowest benthic biomasses were estimated for the northern latitudes (31° S-32° S). Regression models showed that the association between stranded biomass and PSH was highest when only the extent of rocky shore 10 km to the south of each beach was included, suggesting relatively short-distance dispersal and asymmetrical transport of floating kelps, which is further supported by low proportions of rafts with Lepas spp. (indicator of rafting). The results indicate that stranded biomasses are mostly subsidized by nearby benthic populations, which can partly explain the low genetic connectivity among populations in the study region. Future studies should also incorporate other local factors (e.g., winds, currents, wave-exposure) that influence stranding dynamics.

Antagonistic interplay between pH and food resources affects copepod traits and performance in a year-round upwelling system

Linking pH/pCO2 natural variation to phenotypic traits and performance of foundational species provides essential information for assessing and predicting the impact of ocean acidification (OA) on marine ecosystems. Yet, evidence of such linkage for copepods, the most abundant metazoans in the oceans, remains scarce, particularly for naturally corrosive Eastern Boundary Upwelling systems (EBUs). This study assessed the relationship between pH levels and traits (body and egg size) and performance (ingestion rate (IR) and egg reproduction rate (EPR)) of the numerically dominant neritic copepod Acartia tonsa, in a year-round upwelling system of the northern (23° S) Humboldt EBUs. The study revealed decreases in chlorophyll (Chl) ingestion rate, egg production rate and egg size with decreasing pH as well as egg production efficiency, but the opposite for copepod body size. Further, ingestion rate increased hyperbolically with Chl, and saturated at ~1 µg Chl. L-1. Food resources categorized as high (H, >1 µg L-1) and low (L, <1 µg L-1) levels, and pH-values categorized as equivalent to present day (≤400 µatm pCO2, pH > 7.89) and future (>400 µatm pCO2, pH < 7.89) were used to compare our observations to values globally employed to experimentally test copepod sensitivity to OA. A comparison (PERMANOVA) test with Chl/pH (2*2) design showed that partially overlapping OA levels expected for the year 2100 in other ocean regions, low-pH conditions in this system negatively impacted traits and performance associated with copepod fitness. However, interacting antagonistically with pH, food resource (Chl) maintained copepod production in spite of low pH levels. Thus, the deleterious effects of ocean acidification are modulated by resource availability in this system.

Upwelling Bays: How Coastal Upwelling Controls Circulation, Habitat, and Productivity in Bays

Bays in coastal upwelling regions are physically driven and biochemically fueled by their interaction with open coastal waters. Wind-driven flow over the shelf imposes a circulation in the bay, which is also influenced by local wind stress and thermal bay-ocean density differences. Three types of bays are recognized based on the degree of exposure to coastal currents and winds (wide-open bays, square bays, and elongated bays), and the characteristic circulation and stratification patterns of each type are described. Retention of upwelled waters in bays allows for dense phytoplankton blooms that support productive bay ecosystems. Retention is also important for the accumulation of larvae, which accounts for high recruitment in bays. In addition, bays are coupled to the shelf ecosystem through export of plankton-rich waters during relaxation events. Ocean acidification and deoxygenation are a concern in bays because local extrema can develop beneath strong stratification.

Long-term persistence of the floating bull kelp Durvillaea antarctica from the South-East Pacific: Potential contribution to local and transoceanic connectivity

Current knowledge about the performance of floating seaweeds as dispersal vectors comes mostly from mid latitudes (30°-40°), but phylogeographic studies suggest that long-distance dispersal (LDD) is more common at high latitudes (50°-60°). To test this hypothesis, long-term field experiments with floating southern bull kelp Durvillaea antarctica were conducted along a latitudinal gradient (30°S, 37°S and 54°S) in austral winter and summer. Floating time exceeded 200d in winter at the high latitudes but in summer it dropped to 90d, being still higher than at low latitudes (<45d). Biomass variations were due to loss of buoyant fronds. Reproductive activity diminished during long floating times. Physiological changes included mainly a reduction in photosynthetic (Fv/Fm and pigments) rather than in defence variables (phlorotannins and antioxidant activity). The observed long floating persistence and long-term acclimation responses at 54°S support the hypothesis of LDD by kelp rafts at high latitudes.

Embryonic and larval traits of the temperate damselfish Chromis crusma reveal important similarities with other Pomacentridae throughout the family's thermal range

Embryonic development and larval morphology of Chromis crusma was described from five nests sampled between 21 and 25 m depth in central Chile (33°S). From each nest, a set of c. 100 randomly selected eggs were hand-collected and transported in seawater to the laboratory. Subsets of c. 30 eggs per nest were maintained in 50 ml glass containers at a constant ambient temperature of c. 12°C (range 11.5-12.9°C). Egg length (L) and width (W) and larval notochordal length (L_N ) were measured from photographs. Geometric morphometric analyses were performed in newly hatched and 1 week old larvae to quantify shape changes. Ellipsoid eggs had an average (mean ± SE) size of 1.12 ± 0.05 mm L and 0.67 ± 0.02 mm W, with volume being similar throughout 15 developmental stages (i.e., ellipsoid-shaped; 0.27 mm³ ). Planktonic larvae hatched between 5 and 11 days at 12°C and had a mean L_N of 3.13 ± 0.25 mm, a yolk sack volume of 0.03 mm³ and an oil droplet volume of 0.005 mm³ . Morphological traits at hatching included: (a) lack of paired fins and jaws; (b) single medial fin fold; (c) lack of eye pigmentation; (d) yolk sac present near anterior tip; (e) melanophores distributed along ventral surface with one pair over the forehead. In order to generate an up-to-date summary of developmental traits within Pomacentridae, we reviewed literature on egg development (e.g., shape and number of oil droplets), hatching and larval traits (e.g., morphology, pigmentation patterns). Thirty-two publications accounting for 35 species were selected, where eggs, embryonic development, hatching and larval traits were found for 26, 21, 24 and 34 species, respectively. In order to evaluate potential phylogenetic and environmental relationships within the early stages of Pomacentridae, cluster analyses (Bray Curtis similarity, group average) were also performed on egg and larval traits of 22 species divided by subfamily (Stegastinae, Chrominae, Abudefdufinae, Pomacentrinae) and thermal ranges (i.e., low: 16.5°C (range: 12-21°C), medium: 24.5°C (range:21-28°C) and high: 27°C (range: 26-28°C)), suggesting that early developmental patterns can be segregated by both temperature and phylogenetic relationships.

Genetic and morphological divergence at a biogeographic break in the beach-dwelling brooder Excirolana hirsuticauda Menzies (Crustacea, Peracarida)

BACKGROUND

There is a biogeographic break located at 30°S in the southeast Pacific, in a coastal area of strong environmental discontinuities. Several marine benthic taxa with restricted dispersal have a coincident phylogeographic break at 30°S, indicating that genetic structure is moulded by life history traits that limit gene flow and thereby promote divergence and speciation. In order to evaluate intraspecific divergence at this biogeographic break, we investigated the genetic and morphological variation of the directly developing beach isopod Excirolana hirsuticauda along 1900 km of the southeast Pacific coast, across 30°S.

RESULTS

The COI sequences and microsatellite data both identified a strong discontinuity between populations of E. hirsuticauda to the north and south of 30°S, and a second weaker phylogeographic break at approximately 35°S. The three genetic groups were evidenced by different past demographic and genetic diversity signatures, and were also clearly distinguished with microsatellite data clustering. The COI sequences established that the genetic divergence of E. hirsuticauda at 30°S started earlier than divergence at 35°. Additionally, the three groups have different past demographic signatures, with probable demographic expansion occurring earlier in the southern group (south of 35°S), associated with Pleistocene interglacial periods. Interestingly, body length, multivariate morphometric analyses, and the morphology of a fertilization-related morphological character in males, the appendix masculina, reinforced the three genetic groups detected with genetic data.

CONCLUSIONS

The degree of divergence of COI sequences, microsatellite data, and morphology was concordant and showed two geographic areas in which divergence was promoted at differing historical periods. Variation in the appendix masculina of males has probably promoted reproductive isolation. This variation together with gene flow restrictions promoted by life history traits, small body size, oceanographic discontinuities and sandy-beach habitat continuity, likely influenced species divergence at 30°S in the southeast Pacific coast. The degree of genetic and morphological differentiation of populations to the north and south of 30°S suggests that E. hirsuticauda harbours intraspecific divergence consistent with reproductive isolation and an advanced stage of speciation. The speciation process within E. hirsuticauda has been shaped by both restrictions to gene flow and a prezygotic reproductive barrier.

Oceanographic upwelling conditions influence signaling pathways involved in muscle growth of intertidal fish

Few studies have addressed the impact of the upwelling oceanographical conditions on biological processes, such as growth, using a molecular and physiological approach. Upwelling conditions are characterized by low temperature seawater and high nutrient availability, which represents an ideal opportunity to understand how habitat modulates animal performance at different levels of biological complexity. We aimed to assess intraspecific variations in weight, oxygen consumption, protein content, and key signaling pathways involved in muscle-growth (protein kinase B (AKT) and extracellular signal-regulated kinase (ERK)) under experimental trials considering high/low seawater temperatures with full/restricted food rations. For this purpose, we studied Girella laevifrons, one of the most abundant fish species inhabiting rocky intertidal zones along the Eastern South Pacific coasts. Using fish obtained from upwelling (U) and non-upwelling (NU) zones, we reported that U animals displayed higher growth performance during both contrasting trials, with a weight gain (~3 g), lower oxygen consumption (~12%), and higher protein contents (~20%). Only ERK showed significant differences during the trials (~2-fold downregulation between NU and U fish). We also found that U fish increased protein ubiquitination in high water temperature and restricted food ration in contrast to NU fish. Our results help to elucidate how upwelling conditions may influence fish growth at physiological and molecular levels. Still, future analyses are necessary to improve the information regarding the impact of U and NU condition on animals, as well as the possible applications of this data in the aquaculture industry.

Historical isolation and contemporary gene flow drive population diversity of the brown alga Sargassum thunbergii along the coast of China

BACKGROUND

Long-term survival in isolated marginal seas of the China coast during the late Pleistocene ice ages is widely believed to be an important historical factor contributing to population genetic structure in coastal marine species. Whether or not contemporary factors (e.g. long-distance dispersal via coastal currents) continue to shape diversity gradients in marine organisms with high dispersal capability remains poorly understood. Our aim was to explore how historical and contemporary factors influenced the genetic diversity and distribution of the brown alga Sargassum thunbergii, which can drift on surface water, leading to long-distance dispersal.

RESULTS

We used 11 microsatellites and the plastid RuBisCo spacer to evaluate the genetic diversity of 22 Sargassum thunbergii populations sampled along the China coast. Population structure and differentiation was inferred based on genotype clustering and pairwise F ST and allele-frequency analyses. Integrated genetic analyses revealed two genetic clusters in S. thunbergii that dominated in the Yellow-Bohai Sea (YBS) and East China Sea (ECS) respectively. Higher levels of genetic diversity and variation were detected among populations in the YBS than in the ECS. Bayesian coalescent theory was used to estimate contemporary and historical gene flow. High levels of contemporary gene flow were detected from the YBS (north) to the ECS (south), whereas low levels of historical gene flow occurred between the two regions.

CONCLUSIONS

Our results suggest that the deep genetic divergence in S. thunbergii along the China coast may result from long-term geographic isolation during glacial periods. The dispersal of S. thunbergii driven by coastal currents may facilitate the admixture between southern and northern regimes. Our findings exemplify how both historical and contemporary forces are needed to understand phylogeographical patterns in coastal marine species with long-distance dispersal.

Upwelling-derived oceanographic conditions impact growth performance and growth-related gene expression in intertidal fish

Growth is one of the main biological processes in aquatic organisms that is affected by environmental fluctuations such as upwelling (characterized by food-rich waters). In fish, growth is directly related with skeletal muscle increase; which represents the largest tissue of body mass. However, the effects of upwelling on growth, at the physiological and molecular level, are unknown. This study used Girella laevifrons (one of the most abundant intertidal fish in Eastern South Pacific) as a biological model, considering animals from upwelling (U) and non-upwelling (NU) areas. Here, we evaluated the effect of nutritional composition and food availability on growth performance and expression of key growth-related genes (insulin-kike growth factor 1 (igf1) and myosin heavy-chain (myhc)) and atrophy-related genes (muscle ring-finger 1 (murf1), F-box only protein 32 (atrogin-1) and BCL2/adenovirus E1B 19kDa-interacting protein 3 (bnip3)). We reported that, among zones, U fish displayed higher growth performance in response to nutritional composition, specifically between protein- and fiber-rich diets (~1g). We also found in NU fish that atrophy-related genes were upregulated with fiber-rich diet and during fasting (~2-fold at minimum respect U). In conclusion, our results suggest that the growth potential of upwelling fish may be a consequence of differential muscle gene expression. Our data provide a preliminary approach contributing on how upwelling influence fish growth at the physiological and molecular levels. Future studies are required to gain further knowledge about molecular differences between U and NU animals, as well as the possible applications of this knowledge in the aquaculture industry.

The variable routes of rafting: stranding dynamics of floating bull kelp Durvillaea antarctica (Fucales, Phaeophyceae) on beaches in the SE Pacific

Dispersal on floating seaweeds depends on availability, viability, and trajectories of the rafts. In the southern hemisphere, the bull kelp Durvillaea antarctica is one of the most common floating seaweeds, but phylogeographic studies had shown low connectivity between populations from continental Chile, which could be due to limitations in local supply and dispersal of floating kelps. To test this hypothesis, the spatiotemporal dynamics of kelp strandings were examined in four biogeographic districts along the Chilean coast (28°-42°S). We determined the biomass and demography of stranded individuals on 33 beaches for three subsequent years (2013, 2014, 2015) to examine whether rafting is restricted to certain districts and seasons (winter or summer). Stranded kelps were found on all beaches. Most kelps had only one stipe (one individual), although we also frequently found coalesced holdfasts with mature males and females, which would facilitate successful rafting dispersal, gamete release, and reproduction upon arrival. High biomasses of stranded kelps occurred in the northern-central (30°S-33°S) and southernmost districts (37°S-42°S), and lower biomasses in the northernmost (28°S-30°S) and southern-central districts (33°S-37°S). The highest percentages and sizes of epibionts (Lepas spp.), indicative of prolonged floating periods, were found on stranded kelps in the northernmost and southernmost districts. Based on these results, we conclude that rafting dispersal can vary regionally, being more common in the northernmost and southernmost districts, depending on intrinsic (seaweed biology) and extrinsic factors (shore morphology and oceanography) that affect local supply of kelps and regional hydrodynamics.

Short- and long-term acclimation patterns of the giant kelp Macrocystis pyrifera (Laminariales, Phaeophyceae) along a depth gradient

The giant kelp, Macrocystis pyrifera, is exposed to highly variable irradiance and temperature regimes across its geographic and vertical depth gradients. The objective of this study was to extend our understanding of algal acclimation strategies on different temporal scales to those varying abiotic conditions at various water depths. Different acclimation strategies to various water depths (0.2 and 4 m) between different sampling times (Jan/Feb and Aug/Sept 2012; long-term acclimation) and more rapid adjustments to different depths (0.2, 2 and 4 m; short-term acclimation) during 14 d of transplantation were found. Adjustments of variable Chl a fluorescence, pigment composition (Chl c, fucoxanthin), and the de-epoxidation state of the xanthophyll cycle pigments were responsible for the development of different physiological states with respect to various solar radiation and temperature climates. Interestingly, the results indicated that phlorotannins are important during long-term acclimation while antioxidants have a crucial role during short-term acclimation. Furthermore, the results suggested that modifications in total lipids and fatty acid compositions apparently also might play a role in depth acclimation. In Aug/Sept (austral winter), M. pyrifera responded to the transplantation from 4 m to 0.2 m depth with a rise in the degree of saturation and a switch from shorter- to longer-chain fatty acids. These changes seem to be essential for the readjustment of thylakoid membranes and might, thus, facilitate efficient photosynthesis under changing irradiances and temperatures. Further experiments are needed to disentangle the relative contribution of solar radiation, temperature and also other abiotic parameters in the observed physiological changes.

Herbivore-Alga Interaction Strength Influences Spatial Heterogeneity in a Kelp-Dominated Intertidal Community

There is a general consensus that marine herbivores can affect algal species composition and abundance, but little empirical work exists on the role of herbivores as modifiers of the spatial structure of resource assemblages. Here, we test the consumption/bulldozing effects of the molluscan grazer Enoplochiton niger and its influence on the spatial structure of a low intertidal community dominated by the bull kelp Durvillaea antarctica and the kelp Lessonia spicata. Through field experiments conducted at a rocky intertidal shore in north-central Chile (~30°-32°S), the edge of the grazer and algae geographic distributions, we estimated the strength and variability of consumptive effects of the grazer on different functional group of algae. We also used data from abundance field surveys to evaluate spatial co-occurrence patterns of the study species. Exclusion-enclosure experiments showed that E. niger maintained primary space available by preventing algal colonization, even of large brown algae species. The grazing activity of E. niger also reduced spatial heterogeneity of the ephemeral algal species, increasing bare space availability and variability through time in similar ways to those observed for the collective effect with other grazers. Overall, our result suggests that E. niger can be considered an important modifier of the spatial structure of the large brown algae-dominated community. Effects of E. niger on resource variability seem to be directly related to its foraging patterns, large body size, and population densities, which are all relevant factors for management and conservation of the large brown algae community. Our study thus highlights the importance of considering functional roles and identity of generalist consumers on spatial structure of the entire landscape.