Front. Mar. Sci. Frontiers in Marine Science Front. Mar. Sci. 2296-7745 Frontiers Media S.A. 10.3389/fmars.2023.1197468 Marine Science Original Research Changing nutritional seascapes of kelp forests Shalders Tanika C. 1 2 * Champion Curtis 1 2 Benkendorff Kirsten 1 Davis Tom 1 2 Wernberg Thomas 3 4 Morris Stephen 5 Coleman Melinda A. 1 2 3 1 Faculty of Science and Engineering, National Marine Science Centre, Southern Cross University, Coffs Harbour, NSW, Australia 2 Fisheries Research, New South Wales (NSW) Department of Primary Industries, National Marine Science Centre, Coffs Harbour, NSW, Australia 3 The University of Western Australia (UWA) Oceans Institute & School of Biological Sciences, University of Western Australia, Perth, WA, Australia 4 Flødevigen Research Station, Institute of Marine Research, His, Norway 5 New South Wales Department of Primary Industries, Wollongbar Primary Industries Institute, Wollongbar, NSW, Australia

Edited by: George Jackson, Loma Linda University, United States

Reviewed by: Lydia Ladah, Center for Scientific Research and Higher Education in Ensenada (CICESE), Mexico; Patti Virtue, University of Tasmania, Australia

*Correspondence: Tanika C. Shalders, tanika.shalders@gmail.com

21 11 2023 2023 10 1197468 31 03 2023 07 11 2023 Copyright © 2023 Shalders, Champion, Benkendorff, Davis, Wernberg, Morris and Coleman 2023 Shalders, Champion, Benkendorff, Davis, Wernberg, Morris and Coleman

This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

Ocean warming and short-term extreme events (e.g. marine heatwaves) are becoming more intense and frequent and have had major impacts on ecosystems. Seaweeds are foundational components of temperate reefs, providing nutrition for a diversity of species and underpinning temperate food webs. While the impacts of climate-driven environmental change on seaweed distribution, abundance and interactions are well studied, potential impacts on the provision of nutrients from seaweeds and their availability to consumers remain poorly understood. Here, we use metabolomic and lipidomic approaches to understand how the provision of nutrients from key seaweeds may be altered under climate change. We optimize a small-scale microextraction technique to test the effects of warming and marine heatwaves on the nutritional quality of two Australian habitat-forming seaweeds; Ecklonia and Sargassum. We then model changes in lipid availability in response to climate-driven ocean warming throughout Ecklonia’s eastern Australian distribution. Contrary to expectations, ocean warming and marine heatwaves had limited effects on seaweed nutritional quality, indicating resilience in the nutritional value of these species to climate-driven warming. However, nutritional quality varied significantly between species, which presents new implications for predicted species redistributions from the base of the food chain. Modelling revealed a projected net loss of 3.5% of lipids across the east coast of Australia by 2100 under RCP 8.5 in response to shifts in the redistribution of Ecklonia biomass. The climate-driven redistribution of seaweeds is likely to alter the availability of seaweed-derived nutrients to consumers, which may have broad implications for the transfer of energy within temperate marine food webs.

Ecklonia radiata marine heatwaves nutritional quality ocean warming Sargassum linearifolium species distributions fatty acids section-in-acceptance Marine Molecular Biology and Ecology

香京julia种子在线播放

    1. <form id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></form>
      <address id=HxFbUHhlv><nobr id=HxFbUHhlv><nobr id=HxFbUHhlv></nobr></nobr></address>

      Introduction

      Primary producers, such as seaweeds, are foundational components of temperate coasts providing habitat and nutrition that underpins diverse ecosystems (Christie et al., 2009; Bennett et al., 2015; Teagle et al., 2017). Seaweeds are threatened globally by ocean warming and marine heatwaves, which cause shifts in distributions and restructure temperate ecosystems (Vergés et al., 2014; Krumhansl et al., 2016; Smale, 2020; Michaud et al., 2022). Many studies suggest that continuous and discrete warming could have substantial effects on seaweed morphology (Geppi and Riera, 2022), physiology (Van Den Hoek, 1982; Wernberg et al., 2013a; Román et al., 2020), phenology (Coleman and Brawley, 2005; de Bettignies et al., 2018), abundance (Davis et al., 2021), distribution (Lima et al., 2007; Wernberg et al., 2011; Casado-Amezúa et al., 2019) and genetics (Coleman et al., 2020b; Gurgel et al., 2020). Such changes may have subsequent effects to seaweed consumers that perpetuate throughout the entire food chain (Parrish, 2009; Lozano-Montes et al., 2011; Jin et al., 2020). Despite their importance to temperate marine food webs, limited research has investigated the effects of climate change on the small bioactive compounds that contribute to seaweed nutritional quality (e.g. metabolites, but see Park et al., 2023).

      Seaweeds contribute to the health of marine life through the provision of essential nutrients, such as essential fatty acids (McCauley et al., 2015; Skrzypczyk et al., 2019), antioxidants, vitamins and minerals (MacArtain et al., 2007; Wells et al., 2017). Whilst these essential nutrients form a minor component of a seaweeds composition, they have a disproportional importance in nutrition due to intrinsic bioactive properties (Lopes et al., 2020). It remains poorly understood, however, whether climate-driven environmental change will impact these nutrients in seaweeds and thus their availability to seaweed consumers. Reductions to the nutritional quality of seaweeds (e.g. lower levels of essential fatty acids) are likely to affect the trophic transfer of nutrients within marine ecosystems (Parrish, 2009; Lozano-Montes et al., 2011; Jin et al., 2020), potentially compromising reproduction in marine invertebrates and fish (Litzow et al., 2006), with potential flow on effects to higher tropic levels ( Figure 1 ; Wanless et al., 2005). Given the implications for consumers and the subsequent trophic transfer of energy within marine food webs, understanding the impacts of climate change on the nutritional quality of seaweeds is an important research frontier. The application of molecular techniques, such as metabolomics and lipidomics, presents a new pathway to assess such impacts but have not yet been widely employed to understand the change in nutritional composition of seaweeds.

      A representative temperate food web and biosynthetic pathway illustrating: (A) the potential trophic transfer of nutrients derived from Ecklonia (green) and Sargassum (orange) to higher trophic levels under current and future scenarios. (B) The biosynthetic pathway of omega-3 fatty acids EPA and DHA illustrating potential future changes to temperate food webs. The bars indicate the absolute difference in percent content in each fatty acid by species (shown by colour). Greater total levels of omega-3 fatty acids in Ecklonia, relative to Sargassum, indicate a likely decrease in availability at low latitude temperate marine ecosystems under future ocean conditions in east Australia.

      Lipid- and metabolomics methods can be applied to quantify nutritionally valuable compounds in seaweed (e.g. omega-3 fatty acids, sugars and amino acids; Kumar et al., 2016; Britton et al., 2021; Choudhary et al., 2021). These compounds function as building blocks and end products of cellular processes in seaweeds and have beneficial properties for both seaweeds and their consumers. Differences in the composition of these compounds among species and in response to environmental change can provide insights into the health, nutritional quality and physiology of seaweeds (Britton et al., 2020; Ghaderiardakani et al., 2022; Park et al., 2023). The new application of molecular techniques to seaweeds is improving our understanding of biological responses to ocean warming and marine heatwaves, including in the concentration of key nutritional compounds. For example, Britton et al. (2020) found that exposure to marine heatwaves and future ocean conditions led to the saturation of fatty acids in the macroalga Phyllospora comosa. However, considerable scope remains to expand our understanding of biological responses to environmental change using molecular techniques. For instance, combining nutritional responses derived through molecular methods to model the condition and quality of marine species in response to climate change (Bolin et al., 2021) presents a valuable opportunity to understand likely future implications across marine food webs.

      Extreme events such as marine heatwaves have had major impacts on the cooccurring temperate seaweeds, Ecklonia radiata and Sargassum linearifolium, throughout their Australian distributions (Vergés et al., 2016; Straub et al., 2019; Gurgel et al., 2020; Wernberg, 2021). Combined with Phyllospora comosa, our study species, Ecklonia radiata and S. linearifolium are the dominant habitat forming species found on subtidal rocky reefs, throughout New South Wales, Australia (Layton et al., 2020). Within the region and throughout Australia’s Great Southern Reef, these species provide essential ecosystem services (including food and habitat) supporting diverse marine communities (Coleman and Wernberg, 2017; Wernberg et al., 2019). As primary producers, these seaweeds are at the base of the food chain supporting species such as abalone, sea urchins and many species of fish, which in turn are a vital food source for larger marine species ( Figure 1A ; Vanderklift and Wernberg, 2010; Provost et al., 2017; Miranda et al., 2019). A better understanding of how the nutritional quality of these primary producers will be affected by climate change is vital to anticipate and manage changes throughout the temperate Australian marine environment, including the ecological and economic values it supports.

      Here, we test the effects of ocean warming and MHWs on the nutritional quality of two of the most dominant seaweed species on Australian temperate reefs, E. radiata, and S. linearifolium. Specifically, we (1) tested the effects of warming and different types of MHWs on seaweed nutritional quality and the ability of these parameters to recover from impacts, (2) compared the nutritional quality of the two species’, and (3) estimated climate-driven shifts in Ecklonia lipid availability along Australia’s east coast by integrating our results with projections of the redistribution of biomass for this species under multiple climate change scenarios.

      Materials and methods Experimental setup and treatment conditions

      Individual whole plants of Ecklonia radiata and Sargassum linearifolium (hereafter ‘Ecklonia’ and ‘Sargassum’), were collected by SCUBA in Sydney, Australia (34°03’59.0”S 151°09’44.3”E) under NSW DPI permits (P01/0059(A)-3.0 and OUT17/7525). Carefully keeping the holdfasts intact, 150 whole juvenile plants (< 30cm total length) were collected per species from a shallow subtidal rocky reef between 3 to 6m depth. Within 4h of collection, the plants were transported (by air) to the National Marine Science Centre in Coffs Harbour, Australia and placed in 1000 L mesocosms with flowing seawater, weighed down by the holdfast. Following the methods described by Straub et al. (2021) plants were acclimated for 6d before the application of experimental temperature treatments. Four ocean temperature treatments were applied over 48d. The four treatments were i) ambient (22°C), ii) warming (24°C), iii) heatwave non-variable (HWnon) – remaining at 22°C for 20d before moving to 27°C for 20d and back to 22°C for 8d, and iv) heatwave variable (HWvar) – exposed to rapid fluctuations between 22 and 24°C for 8d (temperature variability), remaining at 22°C for 6d moving to a heatwave at 27°C for 20d and back to 22°C for 8d ( Figure 2 ). The magnitude of cumulative intensity was similar among treatments throughout the experiment (warming = 96 degree-days, HWnon = 85 degree-days, and HWvar = 97 degree-days cumulative intensity above the ambient treatment; Hobday et al., 2016).

      Experimental temperature treatment profiles: ambient (22°C), warming (24°C), HWnon (MHW from 22°C to 27°C and back to 22°C), and HWvar (initial fluctuations between 22-24°C, followed by a MHW to 27°C and back to 22°C). HWnon and HWvar treatments were offset by ± 0.2°C to avoid stacking of lines. Vertical dashed lines indicate sampling events (days 20, 33 and 48; modified from Straub et al., 2021).

      Twenty aerated mesocosms (80cm diameter × 45cm high, 230L water volume) received a continuous inflow (flowrate ≈ 2L min-1) of filtered seawater (50μm) from Charlesworth Bay (30°16’3.78”S, 153° 8’25.60”E). Three header tanks with heater/chillers units (Aquahort Ltd.) controlled experimental temperatures at 22, 24, and 27°C, respectively. At the start of the experimental period, five tanks were assigned per treatment, and seven individual Ecklonia and Sargassum plants were placed within each tank (see Supplementary Material ). Seaweeds acclimated at 22°C were randomly allocated to ambient, HWnon, and HWvar treatments (which start at 22°C) and seaweeds acclimated at 24°C were randomly allocated to a warming treatment tank (n = 14 plants/tank).

      Tissue samples were collected on day 20 (after temperature variation, before marine heatwave), day 33 (12d into the marine heatwave), and day 48 (one-week recovery post-heatwave). At each sampling, a 9 mm hole punch was used to collect samples from the meristem of Ecklonia and tweezers were used to remove approximately 20 fronds of Sargassum. Each time, three independent plants, per species (n = 3) were sampled (see Supplementary Material ). Samples were collected from three of the five independent tanks to generate experimental replicates and no plant was sampled twice, allowing the inclusion of ‘exposure time’ as an independent experimental factor. Following best practice, tissue samples were immediately snap-frozen in liquid nitrogen and frozen at – 80°C until analysis (Kumari, 2017).

      Lipid and fatty acid extraction and analysis

      Lipids were extracted from seaweed samples using a scaled-down approach modified from Folch et al. (1957) standard procedures, which accounted for the small sample weight obtained from the Ecklonia and Sargassum (0.07 ± 0.01 and 0.13 ± 0.07 g wet weight, respectively). Using analytical grade solvents (Sigma-Aldrich), samples were soaked in 2:1 chloroform:methanol, then polarized using 0.9% NaCl solution to obtain the lipophilic (chloroform) layer. The lipophilic layer was, transferred to a pre-weighed vial and dried using high-purity nitrogen gas to obtain the lipid extract weight. Lipid extracts were resuspended at 1mg mL-1 in dichloromethane:methanol (50/50 v/v). Total lipid concentration was expressed as a percentage of wet tissue weight for each extracted sample.

      Lipid extracts were analyzed using direct infusion mass spectrometry (DIMS) using a high-resolution SCIEX 5600 Triple-TOF (Simons et al., 2012). Samples were run in both positive and negative modes with a flow rate of 200µL min-1. MS/MS reactions were obtained for peaks in the range of 500–1250 Da. Obtained spectra were analyzed using Lipidview™ 1.3 software (SCIEX) for lipid class and fatty acid composition (Ejsing et al., 2006). Thirty-three lipid classes and 43 fatty acids were detected using both positive and negative detection modes. Because the negative mode detected more lipid classes and fatty acids, the subsequent analysis used negative mode data only to avoid overrepresenting lipids detected in both modes. Lipid classes and fatty acids less than 0.01% were excluded from statistical analyses, and the remaining relative proportion of 13 lipid classes and 19 fatty acids were analyzed separately. The relative compositions (total percent) summed for saturated, monounsaturated and polyunsaturated fatty acids, the ω-3 index and ω-6/ω-3 ratio were also analyzed. The fatty acid double bond isomers cannot be identified using this method of mass spectrometry which limited our ability to distinguish between some fatty acid families (e.g. C20:3n-3/6/9 and C22:5n-3/6). Consequently, the ω-3 index was calculated as the proportion of eicosapentaenoic acid (EPA) and docosahexaenoic acid (DHA) expressed as a percent of total fatty acids. The ω-6/ω-3 PUFA ratio was calculated as arachidonic acid (ARA) divided by the sum of EPA and DHA.

      Metabolomic extraction and analysis

      Tissue subsamples were freeze-dried and derivatized for metabolite analysis using a modified, scaled-down approach based on methods described by Mabuchi et al. (2018; adapted from Lisec et al., 2006; Pongsuwan et al., 2007). A mixed solution of 2.5:1 methanol/ultrapure water/chloroform was used to extract hydrophilic primary metabolites in 15.73 (± 2.00) and 13.19 (± 1.30) mg of freeze-dried and macerated Ecklonia and Sargassum, respectively, with 0.2 µg mg-1 sample (dry weight) of ribitol (Sigma-Aldrich) added as an internal standard. The mixture was stirred for 5 min and centrifuged for a further 5 min (16,000xg, 4°C) and freeze-dried (~16h). Oxime formation was carried out by adding methoxyamine hydrochloride solubilized with pyridine (20mg/mL, 50µL) to the freeze-dried sample (35°C for 80 min). A trimethylsilylation reaction at 40°C for 40 min was induced using N-methyl-N-(trimethylsilyl) trifluoroacetamide (100µL).

      Derivatized metabolite samples were analyzed by gas chromatography-mass spectrometry (GC-MS) using an Agilent 6890A Gas Chromatograph fitted with an Agilent 5973N Mass Selective Detector. The samples were injected with a split ratio of 50:1 using helium carrier gas at 32cm sec–1 constant flow and an inlet temperature of 280°C. Separating volatile compounds by using a Phenomenex ZB-5 column (60m x 0.25mm x 250µm) starting at 60°C with a 3 min hold, then ramped at 8°C min–1 to 330°C (44.75min) then held for 8min. The temperature of the mass spectrometer transfer line was 280°C. The mass spectra were recorded at 70 eV ionisation voltage over the mass range of 45 to 600amu. Peaks were identified by comparison of spectra with MS databases (WILEY 275 and NIST14), along with retention times and elution order. In instances where the compounds could not be positively identified, the compounds were assigned to a metabolite class when possible, based on examination of similarities in the major fragment ions corresponding to functional groups (see Supplementary Material ). The relative abundance of each metabolite was quantified by TIC peak area integration from the GC and expressed as a percentage of the total in each run.

      Statistical analysis Nutritional quality

      Linear mixed models were fitted to assess the effects of seaweed species, treatment, exposure time and their interactions on each nutritional quality response measured. Models used a gaussian distribution and included a random effect for tank ID to accommodate the nesting of tanks within treatments. Model residuals were inspected graphically to check for variance heterogeneity and non-normality. A natural log transformation of the data was required for several responses to meet the assumptions of normality and homogeneity of variance. An analysis of variance was derived from each model to test the effect of species, treatment, exposure time and their interaction. The models were then used to estimate the effects of species, treatment and time in various combinations as required. Post-hoc statistical comparisons of the estimates were made by Tukey’s test. The data analysis was conducted in the R environment (R Core Team, 2022) using the lme4 (Bates et al., 2015) and emmeans (Lenth, 2021) packages for modelling and estimation. Multivariate analysis of overall nutritional change confirmed species-specific differences in the compositional data (see Supplementary Material ).

      Climate change projections

      Projections of Ecklonia lipid availability off the eastern Australian coast (28.0 – 37.5°S) were created using models of future Ecklonia biomass developed by Davis et al. (2021). Briefly, Ecklonia biomass projections were constructed using generalized additive mixed effect models (GAMMs) for the period 2090–2100 under representative concentration pathway (RCP) scenarios 2.6, 6.0 and 8.5 (IPCC, 2014). These projections incorporate ecological interaction effects with the urchin Centrostephanus rodgersii and accounted for the effects of key environmental parameters (habitat availability, light and in situ temperatures) that drive Ecklonia occurrence and cover (see more detail in Davis et al., 2021). As the total extracted lipids per gram wet weight of Ecklonia did not differ between temperature treatments ( Table 1 ), the mean percentage of total Ecklonia lipids (1.4%) was used as a conversion factor to translate projected Ecklonia biomass values into lipid availability. The projected lipid availability maps were created in the R environment using the tidyverse suite of packages (Wickham et al., 2019). We focus our future projections of lipid availability on Ecklonia only as this is the dominant seaweed species occurring throughout the study extent and because comparable biomass projections are still underway for Sargassum.

      Summary table of 3-factor ANOVA of relative percentage of lipid content (% wet weight) and lipid class composition from Ecklonia and Sargassum tissue exposed to four different temperature profiles (ambient, HWnon, HWvar, and warming), sampled before (20 days), during (33 days) and after (48 days) a simulated marine heatwave.

      Lipid ANOVA
      Species Treatment Time Species:Treatment Treatment:Time Species:Time Species:Treatment:Time
      F (df) p F (df) p F (df) p F (df) p F (df) p F (df) p F (df) p
      Content 3.2(1) 0.080 1.0(3) 0.933 2.0(2) 0.152 0.6(3) 0.635 1.5(6) 0.196 0.4(2) 0.680 0.4(6) 0.867
      Class
      Phosphatidic acid (PA)a 241.4(1) <0.001 0.0(3) 0.997 1.0(2) 0.393 0.8(3) 0.499 1.8(6) 0.115 1.9(2) 0.162 3.0(6) 0.016
      Phosphatidylcholine (PC) 180.6(1,32) <0.001 0.4(3,12) 0.792 0.4(2,41) 0.655 1.1(3,32) 0.382 1.0(6) 0.433 0.8(2,36) 0.441 1.4(6,36) 0.261
      Phosphatidylethanolamine (PE)a 187.5(1,35) <0.001 0.9(3,12) 0.482 2.7(2,43) 0.079 0.7(3,35) 0.539 0.8(6,43) 0.583 1.0(2,40) 0.370 1.2(6,39) 0.347
      Phosphatidylglycerol (PG) 23.2(1) <0.001 1.4(3) 0.257 1.7(2) 0.192 0.4(3) 0.762 0.5(6) 0.776 0.7(2) 0.493 1.5(6) 0.187
      Phosphatidylinositol (PI)a 37.7(1,39) <0.001 1.1(3,17) 0.397 7.6(2,44) 0.001 0.5(3,39) 0.679 1.6(6,44) 0.183 1.0(2,42) 0.372 2.2(6,42) 0.065
      Phosphatidylserine (PS) 7.1(1) 0.011 1.0(3) 0.418 0.3(2) 0.717 1.4(3) 0.262 1.7(6) 0.139 0.1(2) 0.904 2.0(6) 0.080
      Cardiolipin (CL) 92.2(1) <0.001 1.5(3) 0.233 1.1(2) 0.334 2.9(3) 0.046 0.8(6) 0.565 3.6(2) 0.034 0.3(6) 0.954
      Diacylglycerol pyrophosphate (DGPP) 14.9(1) <0.001 0.8(3) 0.485 0.5(2) 0.632 0.4(3) 0.778 1.0(6) 0.456 0.6(2) 0.541 1.6(6) 0.172
      Monomethylphosphatidyl ethanolamine (MMPE) 4.3(1) 0.043 0.7(3) 0.549 1.6(2) 0.224 0.4(3) 0.758 0.7(6) 0.645 0.4(2) 0.699 1.2(6) 0.345
      1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE) 15.8(1) <0.001 1.0(3) 0.442 3.8(2) 0.033 7.0(3) 0.001 2.8(6) 0.024 9.4(2) <0.001 1.8(6) 0.122
      Phosphatidylinositol-4-phosphate (PIP)a 20.9(1,31) <0.001 1.2(3,11) 0.346 1.2(2,41) 0.320 3.1(3,31) 0.042 1.2(6,40) 0.343 0.4(2,36) 0.658 0.8(6,35) 0.591
      Cytidine diphosphate diacylglycerol (CDPDAG)a 267.0(1) <0.001 1.4(3) 0.244 0.9(2) 0.404 1.9(3) 0.149 2.0(6) 0.083 5.2(2) 0.009 1.9(6) 0.101
      Sulfatide galactosylceramide (SGalCer)a 3.7(1,27) 0.066 0.7(3,6) 0.600 1.0(2,41) 0.365 0.7(3,27) 0.586 0.7(6,40) 0.640 0.4(2,34) 0.652 2.0(6,34) 0.093

      Significant values in bold. aLog-transformed. Degrees of freedom (df) differ depending on whether a linear model or mixed linear model was used to analyse the data.

      Ambient and warming treatments were kept at a steady temperature throughout.

      Results Effects of ocean warming and marine heatwaves

      Overall, the ocean warming and marine heatwave treatments showed minimal impacts on the nutrient content of Ecklonia and Sargassum. However, univariate analysis of the lipid classes showed that the interaction between species and treatment were statistically significant for 1,2-Dimyristoyl-sn-glycero-3-phosphorylethanolamine (DMPE), cardiolipin (CL) and phosphatidic acid (PA; Table 1 ). Post hoc pairwise comparisons showed that the interactive effect of species and treatment on DMPE was driven by significantly lower levels of DMPE in Sargassum from the ambient treatment relative to the HWnon treatment (p < 0.05). This difference was not observed for Ecklonia. Pair-wise tests did not detect any significant differences in CL composition among treatments for either species (p > 0.05). A significant interaction was observed for PA between species, treatment and exposure time. The post-hoc pairwise comparisons showed that after 33 days (12 days into the heatwaves) phosphatidic acid was significantly lower in Ecklonia from the HWnon temperature treatment compared to the ambient treatment (1.8%, p < 0.05); this trend was not observed for Sargassum or after a period of recovery (at 48 days). The trends in HWnon phosphatidic acid concentration throughout the experiment aligned with the HWnon maximum quantum yield of the Ecklonia meristem (as reported in Straub et al., 2021, Figure 3 ). The concentration of phosphatidic acid was higher under ambient conditions at this sampling time relative to both heatwave treatments.

      Ecklonia (A) phosphatidic acid concentration, and (B) maximum quantum yield [FV/FM; data derived from Straub et al. (2021)]; compared for samples exposed to the non-variable heatwave ‘HWnon’ treatment over time (means ± SE). Horizontal dotted lines mark the start and end of HWnon, letters indicate significant differences of the Tukey’s post-hoc test.

      Ecklonia and Sargassum each showed minimal differences in fatty acids among treatments or among treatments over time ( Table 2 ). However, there was a significant interaction between species and treatment for the saturated fatty acids C15:0 ( Table 2 ). The post-hoc pairwise comparisons showed that the interaction between species and treatment for C15:0 was largely associated with a statistically significant (p < 0.05) increase under the warming treatment compared to the HWnon and HWvar treatments for Ecklonia which was not observed for Sargassum. Although statistically significant, the changes in C15:0 were less than 1% of the total fatty acid composition. A significant interaction between species, treatment and exposure time was detected for total saturated fatty acid (SFA) and the polyunsaturated fatty acids C20:3 and C20:4 ( Table 2 ). However, these three-way interactions were not supported by the detection of any two-way interactions or main effects and no meaningful pairwise differences were detected by the post-hoc tests.

      Summary table of 3-factor ANOVA of relative percentage composition of selected fatty acids from Ecklonia and Sargassum tissue exposed to four different temperature profiles (ambient, HWnon, HWvar, and warming), sampled before (20 days), during (33 days) and after (48 days) a simulated marine heatwave.

      Fatty acids ANOVA
      Species Treatment Time Treatment:Species Treatment:Time Species:Time Species:Treatment : Time
      F (df) p F (df) p F (df) p F (df) p F (df) p F (df) p F (df) p
      C12:0a 9.5(1) 0.004 1.3(3) 0.273 0.5(2) 0.632 0.8(3) 0.524 1.0(6) 0.468 0.5(2) 0.621 1.4(6) 0.236
      C14:0a 96.3(1) <0.001 1.7(3) 0.207 3.4(2) 0.042 0.8(3) 0.511 1.1(6) 0.369 0.2(2) 0.828 0.7(6) 0.667
      C15:0a 0.0(1) 0.944 1.7(3) 0.177 1.7(2) 0.192 7.6(3) <0.001 1.6(6) 0.160 3.7(2) 0.033 1.2(6) 0.310
      C16:0 151.4(1) <0.001 0.9(3) 0.472 3.1(2) 0.054 0.3(3) 0.844 2.2(6) 0.061 0.4(2) 0.655 2.0(6) 0.083
      C18:0a 0.4(1) 0.549 0.6(3) 0.604 0.9(2) 0.436 1.5(3) 0.227 1.2(6) 0.344 1.2(2) 0.302 0.3(6) 0.948
      C20:0a 59.4(1) <0.001 4.9(3) 0.005 1.0(2) 0.364 2.4(3) 0.076 1.3(6) 0.274 0.3(2) 0.715 1.3(6) 0.294
      SFA 76.2(1,32) <0.001 1.4(3,13) 0.283 0.9(2,40) 0.426 2.0(3,32) 0.140 1.7(6,40) 0.156 1.1(2,35) 0.333 3.0(6,35) 0.019
      C16:1 37.5(1) <0.001 2.0(3) 0.122 5.5(2) 0.007 2.1(3) 0.112 2.6(6) 0.028 1.1(2) 0.344 0.5(6) 0.835
      C18:1 0.8(1) 0.372 1.9(3) 0.145 1.2(2) 0.327 0.8(3) 0.510 0.9(6) 0.480 1.0(2) 0.376 1.3(6) 0.288
      C20:1a 2.5(1) 0.123 0.9(3) 0.451 1.6(2) 0.216 0.0(3) 0.992 1.9(6) 0.103 0.3(2) 0.748 1.9(6) 0.101
      C22:1a 362.4(1) <0.001 4.4(3) 0.008 2.7(2) 0.076 0.7(3) 0.587 1.2(6) 0.326 1.6(2) 0.222 2.3(6) 0.051
      MUFA 14.9(1) <0.001 1.5(3) 0.218 5.6(2) 0.006 0.7(3) 0.580 0.3(6) 0.927 0.1(2) 0.893 0.9(6) 0.512
      C18:2 6.3(1,35) 0.017 4.0(3,11) 0.039 0.2(2,44) 0.825 0.4(3,35) 0.727 2.0(6,43) 0.085 3.9(2,40) 0.030 0.6(6,40) 0.743
      C18:3a 91.3(1) <0.001 1.0(3) 0.415 10.7(2) <0.001 0.6(3) 0.614 0.4(6) 0.882 5.3(2) 0.009 1.5(6) 0.204
      C20:2a 30.2(1) <0.001 0.5(3) 0.713 0.7(2) 0.499 1.2(3) 0.321 0.7(6) 0.665 0.4(2) 0.675 1.7(6) 0.140
      C20:3a 126.0(1,31) <0.001 8.7(3,9) 0.802 3.0(2,42) 0.063 0.6(3,31) 0.626 2.1(6,41) 0.077 0.0(2,36) 0.983 3.2(6,36) 0.012
      C20:4 (ARA) 85.8(1,33) <0.001 1.6(3,12) 0.242 0.8(2,41) 0.455 0.4(3,33) 0.743 0.6(6,41) 0.711 0.1(2,37) 0.893 2.6(6,37) 0.031
      C20:5 (EPA) 25.4(3,12) <0.001 0.7(3,12) 0.554 7.4(2,43) 0.002 0.8(3,34) 0.520 1.0(6,42) 0.459 7.3(2,38) 0.002 1.7(6,38) 0.160
      C22:4a 87.1(1) <0.001 0.8(3) 0.508 1.3(2) 0.281 0.6(3) 0.594 1.2(6) 0.313 3.7(2) 0.033 1.1(6) 0.358
      C22:5 15.1(1) <0.001 0.3(3) 0.857 0.6(2) 0.577 0.4(3) 0.752 1.0(6) 0.460 2.0(2) 0.155 0.7(6) 0.628
      C22:6 (DHA) 34.1(1,33) <0.001 1.5(3,9) 0.276 0.9(2,43) 0.420 0.8(3,33) 0.531 1.8(6,42) 0.132 0.1(2,39) 0.877 1.3(6,38) 0.302
      PUFA 66.5(1,34) <0.001 0.0(3,13) 0.997 1.4(2,43) 0.257 0.4(3,35) 0.793 0.9(6,42) 0.530 0.7(2,39) 0.515 1.9(6,38) 0.104
      ω-3 index 18.2(1,34) <0.001 0.5(3,12) 0.663 5.3(2,42) 0.008 0.8(3,34) 0.530 1.2(6,42) 0.324 6.0(2,38) 0.005 1.8(6,38) 0.123
      ω-6/ω -3¤a 7.5(1) 0.009 0.2(3) 0.909 5.3(2) 0.009 0.3(3) 0.829 0.4(6) 0.892 3.8(2) 0.030 0.4(6) 0.874

      Significant values in bold. aLog-transformed. Degrees of freedom (df) differ depending on whether a linear model or mixed linear model was used to analyse the data.

      Ambient and warming treatments were kept at a steady temperature throughout. The ω-3 index is calculated as the amount of EPA and DHA expressed as the percent of total fatty acids. ¤The ω-6/3 ratio is calculated at C20:4/(C20:5 + C22:6).

      When grouped into broad metabolite classifications, alkaloids showed a significant interaction between species and treatments. The post-hoc pairwise comparisons showed that the interaction between species and treatment was largely associated with a statistically significant (p < 0.05) decrease in alkaloids under the HWvar (p < 0.05) compared to the ambient tanks for Ecklonia which was not observed for Sargassum. No other significant interactions in broad metabolite groups, including amino acids and sugars, were found within the seaweed species in response to treatments or with treatments and exposure time ( Table 3 ).

      Summary table of 3-factor ANOVA of relative percentage composition (± SE) of selected metabolites from Ecklonia and Sargassum tissue exposed to four different temperature profiles (ambient, HWnon, HWvar, and warming), sampled before (20 days), during (33 days) and after (48 days) a simulated marine heatwave.

      Metabolite group ANOVA
      Species Treatment Time Species:Treatment Treatment:Time Species:Time Species:Treatment : Time
      F (df) p F (df) p F (df) p F (df) p F (df) p F (df) p F (df) p
      Alcohola 59.1(1) <0.001 8.8(3) <0.001 25.8(2) <0.001 0.9(3) 0.460 1.3(6) 0.279 2.2(2) 0.126 2.1(6) 0.066
      Alkaloida 84.8(1,32) <0.001 3.4(3,10) 0.060 6.9(2,42) 0.003 3.1(3,32) 0.042 0.7(6,42) 0.661 5.9(2,40) 0.006 0.6(6,39) 0.761
      Amino acidsa 74.5(1) <0.001 1.8(3) 0.164 9.5(2) <0.001 0.8(3) 0.519 2.5(6) 0.037 6.1(2) 0.005 2.1(6) 0.067
      Aromatic alcohola 24.6(1,32) <0.001 0.6(3,9) 0.602 4.2(2,43) 0.021 0.1(3,32) 0.963 0.6(6,42) 0.727 1.7(2,41) 0.192 0.6(6,41) 0.733
      Aromatic alkaloida 4.8(1) 0.037 1.1(3) 0.400 9.1(2) <0.001 0.8(3) 0.501 0.4(6) 0.893 7.8(2) 0.001 0.8(6) 0.592
      Aromatic compounda 24.6(1) <0.001 0.6(3) 0.602 4.2(2) 0.021 0.1(3) 0.963 0.6(6) 0.727 1.7(2) 0.192 0.6(6) 0.733
      Complex sugara 0.5(1) 0.470 2.1(3) 0.121 2.2(2) 0.118 0.8(3) 0.508 0.8(6) 0.612 4.5(2) 0.016 0.6(6) 0.734
      Heterocyclea 58.9(1,30) <0.001 1.5(3,9) 0.291 1.7(2,41) 0.202 2.0(3,30) 0.132 0.9(6,41) 0.505 1.1(2,39) 0.339 0.9(6,38) 0.491
      Organic acida 9.9(1) 0.003 2.6(3) 0.064 3.0(2) 0.064 1.6(3) 0.212 3.3(6) 0.008 2.0(2) 0.144 1.2(6) 0.323
      Simple sugars 16.9(1,32) <0.001 1.3(3,11) 0.320 1.0(2,41) 0.384 1.3(3,32) 0.293 1.2(6,41) 0.344 2.0(2,39) 0.144 1.0(6,39) 0.438
      Sterola 62.9(1) 0.000 2.3(3) 0.087 2.0(2) 0.154 0.7(3) 0.583 0.6(6) 0.750 8.8(2) <0.001 0.9(6) 0.540
      Unidentifieda 4.5(1) 0.039 3.4(3) 0.026 8.7(2) <0.001 0.2(3) 0.920 0.8(6) 0.556 3.8(2) 0.030 0.7(6) 0.668

      Significant values in bold. aLog-transformed. Degrees of freedom (df) differ depending on whether a linear model or mixed linear model was used to analyse the data.

      Ambient and warming treatments were kept at a steady temperature throughout.

      Variation between species

      The total extracted lipids of Ecklonia and Sargassum comprised on average 1.400 ± 0.113 and 1.100 ± 0.114% g-1 of the seaweeds’ total wet weight respectively, with no significant differences detected between species, experimental treatments or sampling period ( Table 1 ). The percent composition of lipid classes showed that 92.3% (12 of 13) were significantly different between Ecklonia and Sargassum ( Figure 4 ; Table 1 ). The relative concentration of phosphatidylcholine (PC), phosphatidylethanolamine (PE), phosphatidylglycerol (PG), phosphatidylinositol (PI) and monomethylphosphatidyl ethanolamine (MMPE) were significantly higher in Ecklonia compared to Sargassum ( Table 1 ). While PA, phosphatidylserine (PS), CL, diacylglycerol pyrophosphate (DGPP), DMPE, phosphatidylinositol-4-phosphate (PIP) and cytidine diphosphate diacylglycerol (CDPDAG) were significantly higher in Sargassum in comparison to Ecklonia ( Table 1 ). The percent of sulfatide galactosylceramide (SGalCer) was not significantly different between the two seaweed species ( Table 1 ).

      Absolute difference (%, +SE of the difference) in (A) lipid class; and (B) fatty acid composition between Ecklonia (green) and Sargassum (orange), for the main factor of species. The bar color indicates the species with greater percent composition of each variable. Asterisks (*) indicate a significant difference between species.

      Univariate analysis of the fatty acids showed that 83% (20 of 24) were significantly different between Ecklonia and Sargassum ( Table 2 ). Sargassum was 7.4% higher in total saturated fatty acids than Ecklonia ( Figure 4 ; Table 2 ), driven by significantly higher levels of C12:0, C16:0 and C20:0 when compared to Ecklonia ( Table 2 ). Monounsaturated fatty acids were also 3.4% higher in Sargassum compared to Ecklonia ( Figure 4 ; Table 2 ), with C16:1 and C22:1 significantly higher in Sargassum than Ecklonia ( Table 2 ). Polyunsaturated fatty acids were found to be 10.8% higher in Ecklonia when compared to Sargassum ( Figure 4 ; Table 2 ). Omega-3 (C18:3, C20:5), omega-6 (C20:4, C22:4) and C20:3 fatty acids were significantly higher in Ecklonia when compared to Sargassum ( Table 2 ). Compared to Ecklonia, Sargassum was significantly higher in C18:2 (omega-6), C20:2, C22:5 and C22:6 (omega-3; Table 2 ). Importantly, arachidonic acid (ARA, 20:4), a crucial long-chain polyunsaturated fatty acid in the omega-6 pathway, was 6.2% higher in Ecklonia than Sargassum ( Figure 4 ). α–Linolenic (ALA, 18:3), important to the omega-3 pathway ( Figure 1 ), was also 3.9% higher in Ecklonia ( Figure 4 ). Omega-3 docosahexaenoic acid and omega-6 docosapentaenoic acid showed statistically significant changes however, with only ≤ 0.05% difference between the two species ( Figure 4 ). Omega-3 eicosapentaenoic acid (EPA, 20:5) was, however, 1.6% higher in Ecklonia when compared to Sargassum ( Figure 4 ). The ω-3 index was 1.5% higher in Ecklonia which shows the seaweed species has a higher omega-3 content than Sargassum ( Figure 4 ). Sargassum had a 10.8% higher ω-6:ω-3 ratio of 16.8 compared to 6.0 in Ecklonia ( Figure 4 ).

      Metabolites were grouped into broad classes to further assess functional differences between species (see Supplementary Material ). We found 91.7% (11 of 12) of metabolite classes significantly differed in content between species (see Supplementary Material ). Ecklonia and Sargassum metabolites consist largely of simple sugars (93.3 and 96.5%, respectively; Figure 5A ), of which Sargassum had significantly higher levels. The remaining metabolite composition of Sargassum comprised mostly of complex sugars (2.5%) and alcohols (0.5%; Figure 5B ). Ecklonia metabolite composition was more diverse containing complex sugar (2.4%), and significantly higher levels of heterocycles (1.3%), aromatic alcohols (1.2%), alcohols (0.9%) and alkaloids (0.5%) than Sargassum ( Figure 5B ; Table 3 ).

      Comparisons of metabolite concentration between Ecklonia (n = 36) and Sargassum (n = 34) grouped based on their chemical class; (A) simple sugars (mono-, di- and tri- saccharides) compared to other metabolites, and (B) the compounds included in ‘other metabolites’, for the main factor of species. All compounds are significantly different between species except for complex sugars. The different colours within the bars represent different compounds, as detailed in the keys within each panel.

      Redistribution of <italic>Ecklonia</italic> lipids under climate change

      Projections of Ecklonia lipid availability off eastern Australia (28.0 – 37.5°S) revealed lipid reductions in warmer, lower latitudes where Ecklonia will range contract and increases within cooler, poleward and deeper regions of the study extent by 2100 where Ecklonia biomass will increase under all climate scenarios ( Figure 6 ). Overall, total net lipid gains of 18.5% (4.9 g.m-2; Figure 6A ) and 29.4% (7.7 g.m-2; Figure 6B ) are projected under RCPs 2.6 and 6.0, respectively. In contrast, a 3.5% total net loss in Ecklonia derived lipids are projected under RCP 8.5 (-0.9 g.m-2; Figure 6C ). As the majority of current Ecklonia biomass and subsequent lipid availability occurs at shallow to intermediate reef depths (i.e. <20 m), projected losses are associated with these depths and areas of the coast (primarily lower latitudes) where seabed temperatures are above ~26°C ( Figure 6 ). Projected increases in lipid availability at high latitudes occur across the entire depth range under RCP 2.6, but are mostly restricted to deeper reefs under RCP 8.5 (>20 m; Figure 6 ).

      Projected changes to Ecklonia lipid availability (t) under climate change scenarios (A) RCP 2.6, (B) RCP 6.0 and (C) RCP 8.5 by 2100. Increases in lipid availability is represented by green and decreases by orange.

      Discussion

      Climate-driven environmental change is having profound impacts on the structure and function of temperate marine ecosystems (Wernberg et al., 2016a; Vergés et al., 2019; Holland et al., 2020). Although ocean temperature is one of the most crucial factors influencing seaweed distribution and physiological performance (Van Den Hoek, 1982; Van den Hoek et al., 1990; Krumhansl et al., 2016; Wernberg et al., 2016b; Veenhof et al., 2022), the impacts of temperature on seaweed nutrient composition and the implications for ecosystems have received limited attention (Stengel et al., 2011). In this study, we used molecular techniques to investigate the effects of warming and MHWs on the nutritional quality of two foundational Australian seaweeds, Ecklonia and Sargassum, through the optimization of a small-scale microextraction technique. We found that simulated ocean warming and MHWs had minimal impact on the nutritional properties of Ecklonia and Sargassum. However, nutritional quality differed considerably between the two species, suggesting that the climate-driven redistribution of these species could have implications for the provision of nutrients to marine food webs.

      Ocean warming and marine heatwaves

      Despite both ocean warming and MHWs having caused range contractions, local losses and genetic and physiological changes in Ecklonia and Sargassum in Australia (Vergés et al., 2016; Wernberg et al., 2016a; Straub et al., 2019; Coleman et al., 2020a; Gurgel et al., 2020), we observed minimal changes in the lipid and metabolite composition of these species in response to these ocean stressors. Our results were unexpected and contrary to other studies that observed changes in the fatty acid and metabolite composition of seaweeds in response to temperature changes in both field (Schmid et al., 2017; Gaubert et al., 2019) and laboratory settings (Renaud et al., 2002; Gosch et al., 2015). No biologically significant changes were detected in the fatty acids of Ecklonia or Sargassum in response to the warming and marine heatwave profiles. Increases in fatty acid saturation with increasing temperatures have been widely observed in marine algae (Renaud et al., 2002; Barbosa et al., 2020) and other marine plant species (Beca-Carretero et al., 2018). The unsaturation of fatty acids at low temperatures and saturation at high temperatures is an adaptative strategy algae use to maintain cell membrane fluidity and photosynthetic capacity with changing temperatures (Los and Murata, 2004; Schmid et al., 2021). Although this mechanism has commonly been observed in response to seasonal changes in temperature (Schmid et al., 2014; Schmid et al., 2017), this thermal adaptation may also assist seaweeds in coping with predicted ocean warming and MHWs (Britton et al., 2020). This response was not observed in our experiment for either Ecklonia or Sargassum, which could indicate that these species cannot adjust their fatty acid biosynthesis rapidly in response to unpredictable extreme temperature events. However, Britton et al. (2021) suggest that species that can also occupy intertidal or shallower areas such as Ecklonia and Sargassum may be better adapted to temperature fluctuations and potentially do not need to adjust fatty acid composition to maintain optimum membrane fluidity, or that the adjustments are small and hard to identify. Alternatively, our results could indicate a level of resilience, as the individuals sampled at each time point were survivors, which may have been better able to cope with warmer temperatures (e.g. thermally tolerant genotypes, Vranken et al., 2021; see Straub et al., 2021 for mortality throughout the experiment).

      The time frame over which seaweeds might display a temperature-driven change in fatty acid composition is poorly understood. Changes in seaweed lipids in response to elevated temperature have been observed after as little as three days in the giant kelp Macrocystis pyrifera (Schmid et al., 2020), and differences in seaweed fatty acids have been observed during and after a period of recovery from MHWs (Britton et al., 2020). However, we do not yet understand whether these temperature-driven changes in physiological performance and nutritional properties are species-specific, adaptations of distinct genetic populations or correlate to natural processes (e.g. growth and reproduction). Twelve days into the simulated non-variable heatwave we observed substantially lower levels of phosphatidic acid in Ecklonia when compared to the ambient temperature treatment. However, following a period of recovery from the heatwave, this difference was no longer detected, suggesting recovery. Levels of phosphatidic acid were initially also relatively low in the warming and variable heatwave treatments but increased over time. Phosphatidic acid is an important component of seaweed chloroplast membranes (Melo et al., 2015), which are important for photosynthesis (Wang and Benning, 2012). During the marine heatwave treatments, Ecklonia photosynthetic quantum yield was significantly lower than the warming and ambient treatments but these differences were no longer observed after the recovery period ( Figure 3 ; Straub et al., 2021). As only surviving plants were sampled at each time point, the performance improvement could be reflective of the responses of the better adapted survivors at 56 days compared to individuals sampled at 33 days. However, similar reductions in photosynthesis have been observed in other brown seaweeds (e.g. Phyllospora comosa) in response to MHWs, and following a period of recovery no lasting negative effects were detected (Britton et al., 2020). Therefore, the reduction in phosphatidic acid during the marine heatwave may be indicative of a stress response and connected to photosynthetic performance, however, further research is required to establish phosphatidic acid as a biomarker for photosynthetic stress responses in seaweeds.

      Simulated ocean warming and marine heatwave treatments showed minimal effects on the metabolite content of both seaweed species. In particular, there was no evidence for impacts of temperature on complex and simple sugars which indicates there were no major effects on the potential beneficial properties (e.g. anti-inflammatory), energy content or calorific value of these species (Bayu et al., 2021). For example, fucoidan sugars like fructose have several biological activities including anti-inflammatory and anti-tumour properties which were unaffected by the temperature changes (Deniaud-Bouët et al., 2017). Amino acids function as important building blocks for protein (Maschek and Baker, 2008) and were unaffected by simulated ocean warming and marine heatwave scenarios. This suggests that temperature treatments did not trigger stress-induced protein degradation (Kumar et al., 2016), which would be expected to increase the levels of some amino acids in the tissues (Araújo et al., 2011). Additional studies are required to determine if the stability in metabolite content is a sign of resilience in the species’ or attributed to the sampling of survivors at each time point. While the application of metabolomics to marine system biology is increasing globally, the main obstacle that remains is the accurate annotation and identification of metabolites (Kumar et al., 2016; Lawson et al., 2022). A change was detected in an unidentified metabolite in response to the simulated temperatures, however, without sufficient context on the function of these compounds, the potential impacts on Ecklonia metabolism and nutritional quality could not be interpreted. For the future application of metabolomics techniques to marine seaweeds, access to advanced technology and the establishment of large mass spectral libraries centred on model organisms are needed (Kumar et al., 2016; Lawson et al., 2022).

      Variation in nutritional quality between species

      The nutritional composition between Ecklonia and Sargassum varied considerably. All 68 nutritional properties detected were present in both species, however, 82% differed substantially in content between species. The magnitude of difference between species was larger among lipid classes and fatty acids than metabolites. The lipid content of both Ecklonia and Sargassum was low (1.4 and 1.1%, respectively), with no difference detected between species and were similar to levels previously measured (Verma et al., 2017; Skrzypczyk et al., 2019). Ecklonia and Sargassum had lower total lipid levels than other Australian seaweed species such as Cystophora torulosa and C. polycystidea (8.4 and 5.7%, respectively), but were higher than the commercial species Sargassum fusiforme (0.7%, Skrzypczyk et al., 2019). Macroalgae are typically low in total fatty acids (in comparison to other algae e.g. microalgae) but have a high proportion of valuable polyunsaturated fatty acids (Schmid et al., 2014). This was supported by our results, as Ecklonia and Sargassum fatty acid composition was comprised of 39.3 and 28.6% polyunsaturated fatty acids, respectively. Although both species contained relatively high proportions of polyunsaturated fatty acids, higher levels were found in Ecklonia, while Sargassum contained substantially higher levels of saturated and monounsaturated fatty acids. Concentrations of polyunsaturated fatty acids tend to be higher in species from cold water and saturated fatty acids higher in species from warmer water, which could explain why temperate-tropical adapted species Sargassum displays higher levels of saturation (Choudhary et al., 2021).

      The higher levels of polyunsaturated fatty acids in Ecklonia were driven by higher levels of omega-3 fatty acids (alpha-linolenic and docosapentaenoic) and arachidonic acid (an omega-6 fatty acid) when compared to Sargassum. Despite omega-3 fatty acids being primarily known for their benefits to human nutrition, these valuable fatty acids are also important to marine systems as an essential component of the cells of marine plants, animals and bacteria (Valentine and Valentine, 2004; Gladyshev et al., 2018). For example, numerous fish species have high eicosapentaenoic and docosahexaenoic acid content that is believed to support fast continuous swimming through its role in muscle cell metabolism (Gladyshev et al., 2018). Ecklonia displayed a higher omega-3 index than Sargassum and could be a better source of these fatty acids in temperate ecosystems in comparison. By contrast, Sargassum had a higher omega-6 to omega-3 ratio. Due to the opposing effects of omega-6 and omega-3 fatty acids (i.e. pro- and anti-inflammatory), a balanced consumption is recommended in human diets (Simopoulos, 2002). However, research into the importance of omega fatty acids for marine food chains is limited and the importance of this ratio for marine consumers is unclear. The species differences in nutritional composition are important in the context of nutritional seascapes because these two species are projected to change significantly and contrastingly under climate change.

      Seaweed-herbivore interactions are complex (Hay and Fenical, 1988; Poore et al., 2012; Wernberg et al., 2013b). Consequently, structural shifts in seaweed availability may not equate to a suitable dietary replacement for consumers. Palatability is a key factor influencing herbivory and is influenced by macronutrient composition, for example, high protein and fat content often lead to higher palatability. Several other palatability factors including secondary metabolites, morphology and physical stress can also influence the herbivory of seaweeds (Cheng et al., 2022). Research has shown abalone species prefer a diet of brown algae over red or green algae which may be driven by high levels of polyunsaturated acids, in particular, eicosapentaenoic and arachidonic fatty acids (Nelson et al., 2002; Guest et al., 2008; Holland et al., 2021). Secondary metabolites are typically active in low concentrations and can impact the palatability for herbivores (Hay, 1996). These compounds were found in higher content in Ecklonia, including aromatic alcohols, heterocycles and alkaloids and could be used in defence against predators and pathogens (Hay, 1996; Engel et al., 2002). Some seaweeds also actively deter herbivory through chemical defences, for example, brown seaweeds contain phlorotannins which can discourage consumption by herbivores (Toth and Pavia, 2007). Simple sugars can affect the sweetness of seaweeds (Patrick et al., 2013), which may influence the consumption preferences of herbivores. These sugars were higher in Sargassum than Ecklonia, however, both species contained high quantities. These palatability factors may impact the uptake of nutrients from seaweeds despite their nutritional quality, providing a novel insight into potential nutrient transfer in marine food webs through metabolomics.

      Redistribution of nutrients under climate change

      Seaweeds supply coastal ecosystems with nutrients directly through macroalgal detritus (Duggins et al., 1989) and herbivory (Dunton and Schell, 1987) that are transferred to higher trophic levels ( Figure 1 ; Vanderklift and Wernberg, 2010; Jin et al., 2020). Global changes in seaweed distribution and biomass have been driven by ocean warming (Krumhansl et al., 2016; Martínez et al., 2018), localized marine heatwave events (Wernberg et al., 2016a; Straub et al., 2019; Félix-Loaiza et al., 2022), and local anthropogenic stressors (Reed et al., 1994; Coleman et al., 2008; Strain et al., 2014), altering the nutritional seascapes of temperate environments. By 2100, the Australian distribution of Ecklonia is projected to contract by 49-71% under RCPs 2.6 and 6.0, respectively. Similarly, a 41-57% contraction in the distribution of Sargassum is projected under RCPs 2.6 and 6.0, respectively (Martínez et al., 2018). These estimates are based on correlative responses between each species and environmental variables (primarily SST) only and recent work incorporating key ecological interactions into future projections (e.g. urchin herbivory) of kelp biomass loss are more optimistic (Castro et al., 2020; Davis et al., 2021). For example, incorporating ecological interactions with urchins moderated projected net changes in Ecklonia biomass on the south-east coast of Australia to +18.5, +29.4% and -3.5% under RCPs 2.6, 6.0 and 8.5, respectively (Davis et al., 2021).

      Projections of change in Ecklonia lipid availability across Australia’s east coast were made by combining our quantified lipid content with predicted alterations in distribution and biomass (Davis et al., 2021). While we found total net gains in Ecklonia lipid availability across south-east Australia under RCPs 2.6 and 6.0 and a small magnitude of total loss under RCP 8.5, the changes in lipid availability are associated with considerable spatial variation due to the shifting distribution of Ecklonia biomass off eastern Australia. For example, in northern areas of the study extent where it is projected that Ecklonia will undergo a 55-276 km range contraction (Davis et al., 2021), a complete loss of Ecklonia derived lipid availability is likely while increases are expected in more poleward regions as biomass is projected to increase due to weakening urchin herbivory. Increasing water temperatures are also projected to drive Ecklonia redistributions from shallow reefs to cooler, deeper reefs in mid-latitude areas (Davis et al., 2021). These changes may have major implications for marine consumers as the availability of nutrients shifts across not only a latitudinal but also a vertical dimension, particularly if consumers do not experience similar shifts in distribution. However, local drivers such as eutrophication may buffer the negative physiological effects of ocean warming in some areas (Fernández et al., 2020). The results also suggest that the redistribution of Ecklonia and Sargassum will considerably alter the availability of nutrients to marine ecosystems in areas where patterns of species co-occurrence are affected and there is a transition from one species to the other. For example, the predicted shift from Ecklonia to Sargassum dominated reefs at lower latitudes (Vergés et al., 2014; Vergés et al., 2019) would signify increasing saturation in available fatty acids and loss of metabolite diversity, while the loss of canopy cover, in general, would lead to an overall loss of valuable essential fatty acids for secondary consumers ( Figure 1 ). Changes at the primary level could scale up through multiple trophic levels (e.g. primary, secondary, tertiary and finally humans) altering the trophic transfer of nutrients throughout entire marine systems (Jin et al., 2020). Combining lipidomic and metabolomic techniques with traditionally measured carbon, nitrogen and phosphate ratios will provide a greater depth to quantifying the nutritional changes in seaweeds (Peters et al., 2005). Moreover, to fully understand potential change in nutritional seascapes under climate change, we still require better understanding of food web dynamics as they relate to consumer and producer climate-redistributions and their interactions.

      Conclusion

      Lipidomic and metabolomic approaches harbor considerable value for quantifying ecosystem production and could be a powerful way to understand how nutritional seascapes may be altered under projected climate change. However, several knowledge gaps remain that must be addressed before such approaches can be widely adopted. First, future application of metabolomics techniques to marine seaweeds will rely on access to advanced technology and the establishment of large mass spectral libraries that allow the characterization of compounds and products that are currently uncharacterised (Kumar et al., 2016; Lawson et al., 2022). In addition, further research is required to establish links between seaweed compounds and biological responses (e.g. phosphatidic acid and photosynthetic stress). Nevertheless, through the optimization of a small-scale microextraction technique, we were able to demonstrate little influence of marine heatwaves and warming on the nutritional quality of two dominant habitat-forming seaweeds, but major differences between the species. These results could have major implications for climate-mediated shifts in species distributions and replacements which will fundamentally alter the nutritional seascapes of marine ecosystems under climate change.

      Data availability statement

      The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

      Author contributions

      This project was conceptualized by TCS, MAC, TW, CC and KB. The experiment was run by MAC. Lipidomic and metabolomic analysis was undertaken by TCS and KB. Statistical analysis was carried out by TCS, SM, TD and MAC. Writing the manuscript was led by TCS with assistance from all authors. All authors contributed to the article and approved the submitted version.

      Funding

      This work is part of the NSW Primary Industries Climate Change Research Strategy, funded by the NSW Climate Change Fund. This work was undertaken with funding from the Australian Research Council (DP190100058 awarded to MC and TW) and the Research Training Program Scholarship awarded to TS. The SCIEX equipment used for lipid analyses was funded by the Australian Research Council Linkage Infrastructure, Equipment and Facilities grant LE170100192.

      Acknowledgments

      We extend our thanks to Sandra Straub, Ziggy Marzinelli, Adriana Verges, Reina Veenhof, Peter Butcherine and Endurance Ewere for their assistance and Bryan Mead and Ashley Dowell from the Southern Cross Analytical Research Facility for technical support. The vectors used in Figure 1 are courtesy of the Integration and Application Network, the University of Maryland Center for Environmental Science (ian.umces.edu/symbols/) and PhyloPic (phylopic.org/). We thank the reviewers for their invaluable feedback, which improved our manuscript.

      Conflict of interest

      The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

      Publisher’s note

      All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

      Supplementary material

      The Supplementary Material for this article can be found online at: /articles/10.3389/fmars.2023.1197468/full#supplementary-material

      References Araújo W. L. Tohge T. Ishizaki K. Leaver C. J. Fernie A. R. (2011). Protein degradation – an alternative respiratory substrate for stressed plants. Trends Plant Sci. 16 (9), 489498. doi: 10.1016/j.tplants.2011.05.008 Barbosa M. Fernandes F. Pereira D. M. Azevedo I. C. Sousa-Pinto I. Andrade P. B. . (2020). Fatty acid patterns of the kelps Saccharina latissima, Saccorhiza polyschides and Laminaria ochroleuca: Influence of changing environmental conditions. Arabian. J. Chem. 13 (1), 4558. doi: 10.1016/j.arabjc.2017.01.015 Bates D. Mächler M. Bolker B. Walker S. (2015). Fitting linear mixed-effects models using lme4. J. Stat. Software. 67 (1), 148. doi: 10.18637/jss.v067.i01 Bayu A. Warsito M. F. Putra M. Y. Karnjanakom S. Guan G. (2021). Macroalgae-derived rare sugars: Applications and catalytic synthesis. Carbon. Resour. Conversion. 4, 150163. doi: 10.1016/j.crcon.2021.04.002 Beca-Carretero P. Guihéneuf F. Marín-Guirao L. Bernardeau-Esteller J. García-Muñoz R. Stengel D. B. . (2018). Effects of an experimental heat wave on fatty acid composition in two Mediterranean seagrass species. Mar. Pollut. Bull. 134, 2737. doi: 10.1016/j.marpolbul.2017.12.057 Bennett S. Wernberg T. Connell S. D. Hobday A. J. Johnson C. R. Poloczanska E. S. (2015). The ‘Great Southern Reef’: social, ecological and economic value of Australia’s neglected kelp forests. Mar. Freshw. Res. 67 (1), 4756. doi: 10.1071/MF15232 Bolin J. A. Schoeman D. S. Evans K. J. Cummins S. F. Scales K. L. (2021). Achieving sustainable and climate-resilient fisheries requires marine ecosystem forecasts to include fish condition. Fish. Fisheries. 118. doi: 10.1111/faf.12569 Britton D. Schmid M. Noisette F. Havenhand J. N. Paine E. R. McGraw C. M. . (2020). Adjustments in fatty acid composition is a mechanism that can explain resilience to marine heatwaves and future ocean conditions in the habitat-forming seaweed Phyllospora comosa (Labillardière) C.Agardh. Global Change Biol. 26 (6), 35123524. doi: 10.1111/gcb.15052 Britton D. Schmid M. Revill A. T. Virtue P. Nichols P. D. Hurd C. L. . (2021). Seasonal and site-specific variation in the nutritional quality of temperate seaweed assemblages: implications for grazing invertebrates and the commercial exploitation of seaweeds. J. Appl. Phycol. 33 (1), 603616. doi: 10.1007/s10811-020-02302-1 Casado-Amezúa P. Araújo R. Bárbara I. Bermejo R. Borja Á. Díez I. . (2019). Distributional shifts of canopy-forming seaweeds from the Atlantic coast of Southern Europe. Biodivers. Conserv. 28 (5), 11511172. doi: 10.1007/s10531-019-01716-9 Castro L. C. Cetina-Heredia P. Roughan M. Dworjanyn S. Thibaut L. Chamberlain M. A. . (2020). Combined mechanistic modelling predicts changes in species distribution and increased co-occurrence of a tropical urchin herbivore and a habitat-forming temperate kelp. Divers. Distrib. 26 (9), 12111226. doi: 10.1111/ddi.13073 Cheng A. Lim W. Y. Lim P.-E. Yang Amri A. Poong S.-W. Song S.-L. . (2022). Marine autotroph-herbivore synergies: Unravelling the roles of macroalgae in marine ecosystem dynamics. Biology 11 (8). doi: 10.3390/biology11081209 Choudhary B. Chauhan O. P. Mishra A. (2021). Edible seaweeds: A potential novel source of bioactive metabolites and nutraceuticals with human health benefits. Front. Mar. Sci. 8. doi: 10.3389/fmars.2021.740054 Christie H. Norderhaug K. M. Fredriksen S. (2009). Macrophytes as habitat for fauna. Mar. Ecol. Prog. Ser. 396, 221233. doi: 10.3354/meps08351 Coleman M. A. Brawley S. H. (2005). Spatial and temporal variability in dispersal and population genetic structure of a rockpool alga. Mar. Ecol. Prog. Ser. 300, 6377. doi: 10.3354/meps300063 Coleman M. A. Kelaher B. P. Steinberg P. D. Millar A. J. K. (2008). Absence of a large brown macroalga on urbanized rocky reefs around Sydney, Australia, and evidence for historical decline. J. Phycol. 44 (4), 897901. doi: 10.1111/j.1529-8817.2008.00541.x Coleman M. A. Minne A. J. P. Vranken S. Wernberg T. (2020a). Genetic tropicalisation following a marine heatwave. Sci. Rep. 10 (1), 12726. doi: 10.1038/s41598-020-69665-w Coleman M. A. Wernberg T. (2017). Forgotten underwater forests: The key role of fucoids on Australian temperate reefs. Ecol. Evol. 7 (20), 84068418. doi: 10.1002/ece3.3279 Coleman M. A. Wood G. Filbee-Dexter K. Minne A. J. P. Goold H. D. Vergés A. . (2020b). Restore or redefine: Future trajectories for restoration. Front. Mar. Sci. 7 (237). doi: 10.3389/fmars.2020.00237 Davis T. R. Champion C. Coleman M. A. (2021). Climate refugia for kelp within an ocean warming hotspot revealed by stacked species distribution modelling. Mar. Environ. Res. 166, 105267. doi: 10.1016/j.marenvres.2021.105267 de Bettignies T. Wernberg T. Gurgel C. F. D. (2018). Exploring the influence of temperature on aspects of the reproductive phenology of temperate seaweeds. Front. Mar. Sci. 5 (218). doi: 10.3389/fmars.2018.00218 Deniaud-Bouët E. Hardouin K. Potin P. Kloareg B. Hervé C. (2017). A review about brown algal cell walls and fucose-containing sulfated polysaccharides: Cell wall context, biomedical properties and key research challenges. Carbohydr. Polymers. 175, 395408. doi: 10.1016/j.carbpol.2017.07.082 Duggins D. O. Simenstad C. A. Estes J. A. (1989). Magnification of secondary production by kelp detritus in coastal marine ecosystems. Science 245 (4914), 170173. doi: 10.1126/science.245.4914.170 Dunton K. H. Schell D. M. (1987). Dependence of consumers on macroalgal (Laminaria solidungula) carbon in an arctic kelp community: δ13C evidence. Mar. Biol. 93 (4), 615625. doi: 10.1007/BF00392799 Ejsing C. S. Duchoslav E. Sampaio J. Simons K. Bonner R. Thiele C. . (2006). Automated identification and quantification of glycerophospholipid molecular species by multiple precursor ion scanning. Analytical. Chem. 78 (17), 62026214. doi: 10.1021/ac060545x Engel S. Jensen P. R. Fenical W. (2002). Chemical ecology of marine microbial defense. J. Chem. Ecol. 28 (10), 19711985. doi: 10.1023/A:1020793726898 Félix-Loaiza A. C. Rodríguez-Bravo L. M. Beas-Luna R. Lorda J. de la Cruz-González E. Malpica-Cruz L. (2022). Marine heatwaves facilitate invasive algae takeover as foundational kelp. Bot. Mar. 65 (5), 315–319. doi: 10.1515/bot-2022-0037 Fernández P. A. Gaitán-Espitia J. D. Leal P. P. Schmid M. Revill A. T. Hurd C. L. (2020). Nitrogen sufficiency enhances thermal tolerance in habitat-forming kelp: implications for acclimation under thermal stress. Sci. Rep. 10 (1), 3186. doi: 10.1038/s41598-020-60104-4 Folch J. Lees M. Sloane Stanley G. H. (1957). A simple method for the isolation and purification of total lipides from animal tissues. J. Biol. Chem. 226 (1), 497509. doi: 10.1016/S0021-9258(18)64849-5 Gaubert J. Payri C. E. Vieira C. Solanki H. Thomas O. P. (2019). High metabolic variation for seaweeds in response to environmental changes: a case study of the brown algae Lobophora in coral reefs. Sci. Rep. 9 (1), 993. doi: 10.1038/s41598-018-38177-z Geppi E. F. Riera R. (2022). Responses of intertidal seaweeds to warming: A 38- year time series shows differences of sizes. Estuarine. Coast. Shelf. Sci. 270, 107841. doi: 10.1016/j.ecss.2022.107841 Ghaderiardakani F. Langhans L. Kurbel V. B. Fenizia S. Wichard T. (2022). Metabolite profiling reveals insights into the species-dependent cold stress response of the green seaweed holobiont Ulva (Chlorophyta). Environ. Exp. Bot. 200, 104913. doi: 10.1016/j.envexpbot.2022.104913 Gladyshev M. I. Sushchik N. N. Tolomeev A. P. Dgebuadze Y. Y. (2018). Meta-analysis of factors associated with omega-3 fatty acid contents of wild fish. Rev. Fish. Biol. Fisheries. 28 (2), 277299. doi: 10.1007/s11160-017-9511-0 Gosch B. J. Lawton R. J. Paul N. A. de Nys R. Magnusson M. (2015). Environmental effects on growth and fatty acids in three isolates of Derbesia tenuissima (Bryopsidales, Chlorophyta). Algal. Res. 9, 8293. doi: 10.1016/j.algal.2015.02.022 Guest M. A. Nichols P. D. Frusher S. D. Hirst A. J. (2008). Evidence of abalone (Haliotis rubra) diet from combined fatty acid and stable isotope analyses. Mar. Biol. 153 (4), 579588. doi: 10.1007/s00227-007-0831-9 Gurgel C. F. D. Camacho O. Minne A. J. P. Wernberg T. Coleman M. A. (2020). Marine heatwave drives cryptic loss of genetic diversity in underwater forests. Curr. Biol. 30 (7), 11991206.e1192. doi: 10.1016/j.cub.2020.01.051 Hay M. E. (1996). Marine chemical ecology: what’s known and what’s next? J. Exp. Mar. Biol. Ecol. 200 (1), 103134. doi: 10.1016/S0022-0981(96)02659-7 Hay M. E. Fenical W. (1988). Marine plant-herbivore interactions: The ecology of chemical defense. Annu. Rev. Ecol. Syst. 19 (1), 111145. doi: 10.1146/annurev.es.19.110188.000551 Hobday A. J. Alexander L. V. Perkins S. E. Smale D. A. Straub S. C. Oliver E. C. . (2016). A hierarchical approach to defining marine heatwaves. Prog. Oceanogr. 141, 227238. doi: 10.1016/j.pocean.2015.12.014 Holland M. M. Smith J. A. Everett J. D. Vergés A. Suthers I. M. (2020). Latitudinal patterns in trophic structure of temperate reef-associated fishes and predicted consequences of climate change. Fish. Fisheries. 21 (6), 10921108. doi: 10.1111/faf.12488 Holland O. J. Young M. A. Sherman C. D. H. Tan M. H. Gorfine H. Matthews T. . (2021). Ocean warming threatens key trophic interactions supporting a commercial fishery in a climate change hotspot. Global Change Biol. 27 (24), 64986511. doi: 10.1111/gcb.15889 IPCC . (2014). Climate change 2014: synthesis report. Contribution of working groups I, II and III to the fifth assessment report f the Intergovernmental Panel on Climate Change (Geneva, Switzerland: IPCC). Jin P. Hutchins D. A. Gao K. (2020). The impacts of ocean acidification on marine food quality and its potential food chain consequences. Front. Mar. Sci. 7 (780). doi: 10.3389/fmars.2020.543979 Krumhansl K. A. Okamoto D. K. Rassweiler A. Novak M. Bolton J. J. Cavanaugh K. C. . (2016). Global patterns of kelp forest change over the past half-century. Proc. Natl. Acad. Sci. 113 (48), 1378513790. doi: 10.1073/pnas.1606102113 Kumar M. Kuzhiumparambil U. Pernice M. Jiang Z. Ralph P. J. (2016). Metabolomics: an emerging frontier of systems biology in marine macrophytes. Algal. Res. 16, 7692. doi: 10.1016/j.algal.2016.02.033 Kumari P. (2017). “Seaweed lipidomics in the era of ‘omics’ biology: A contemporary perspective,” in Systems biology of Marine Ecosystems. Eds. Kumar M. Ralph P. (Cham: Springer). Lawson C. A. Camp E. Davy S. K. Ferrier-Pagès C. Matthews J. Suggett D. J. (2022). “Informing coral reef conservation through metabolomic approaches,” in Coral reef conservation and restoration in the omics age, vol. Vol. 1 . Eds. van Oppen M. J. H. Aranda Lastra M. (Cham: Springer International Publishing) 179202. Layton C. Coleman M. A. Marzinelli E. M. Steinberg P. D. Swearer S. E. Vergés A. . (2020). Kelp forest restoration in Australia. Front. Mar. Sci. 7. doi: 10.3389/fmars.2020.00074 Lenth R. (2021) emmeans: Estimated marginal means, aka leastsquares means (Version R package 1.7.5). Available at: https://CRAN.Rproject.org/package=emmeans. Lima F. P. Ribeiro P. A. Queiroz N. Hawkins S. J. Santos A. M. (2007). Do distributional shifts of northern and southern species of algae match the warming pattern? Global Change Biol. 13 (12), 25922604. doi: 10.1111/j.1365-2486.2007.01451.x Lisec J. Schauer N. Kopka J. Willmitzer L. Fernie A. R. (2006). Gas chromatography mass spectrometry–based metabolite profiling in plants. Nat. Protoc. 1 (1), 387396. doi: 10.1038/nprot.2006.59 Litzow M. ,. A. Bailey K. M. Prahl F. G. Heintz R. (2006). Climate regime shifts and reorganization of fish communities: the essential fatty acid limitation hypothesis. Mar. Ecol. Prog. Ser. 315, 111. doi: 10.3354/meps315001 Lopes D. Melo T. Rey F. Meneses J. Monteiro F. L. Helguero L. A. . (2020). Valuing bioactive lipids from green, red and brown macroalgae from aquaculture, to foster functionality and biotechnological applications. Molecules 25 (17). doi: 10.3390/molecules25173883 Los D. A. Murata N. (2004). Membrane fluidity and its roles in the perception of environmental signals. Biochim. Biophys. Acta (BBA). - Biomembranes. 1666 (1), 142157. doi: 10.1016/j.bbamem.2004.08.002 Lozano-Montes H. M. Loneragan N. R. Babcock R. C. Jackson K. (2011). Using trophic flows and ecosystem structure to model the effects of fishing in the Jurien Bay Marine Park, temperate Western Australia. Mar. Freshw. Res. 62 (5), 421431. doi: 10.1071/MF09154 Mabuchi R. Adachi M. Kikutani H. Tanimoto S. (2018). Discriminant analysis of muscle tissue type in yellowtail Seriola quinqueradiata muscle based on metabolic component profiles. Food Sci. Technol. Res. 24 (5), 883891. doi: 10.3136/fstr.24.883 MacArtain P. Gill C. I. R. Brooks M. Campbell R. Rowland I. R. (2007). Nutritional value of edible seaweeds. Nutr. Rev. 65 (12), 535543. doi: 10.1111/j.1753-4887.2007.tb00278.x Martínez B. Radford B. Thomsen M. S. Connell S. D. Carreño F. Bradshaw C. J. A. . (2018). Distribution models predict large contractions of habitat-forming seaweeds in response to ocean warming. Divers. Distrib. 24 (10), 13501366. doi: 10.1111/ddi.12767 Maschek J. A. Baker B. J. (2008). “The chemistry of algal secondary metabolism,” in Algal Chemical Ecology. Ed. Amsler C. D. (Berlin, Heidelberg: Springer Berlin Heidelberg), 124. McCauley J. I. Meyer B. J. Winberg P. C. Ranson M. Skropeta D. (2015). Selecting Australian marine macroalgae based on the fatty acid composition and anti-inflammatory activity. J. Appl. Phycol. 27 (5), 21112121. doi: 10.1007/s10811-014-0465-5 Melo T. Alves E. Azevedo V. Martins A. S. Neves B. Domingues P. . (2015). Lipidomics as a new approach for the bioprospecting of marine macroalgae — Unraveling the polar lipid and fatty acid composition of Chondrus crispus . Algal. Res. 8, 181191. doi: 10.1016/j.algal.2015.02.016 Michaud K. M. Reed D. C. Miller R. J. (2022). The Blob marine heatwave transforms California kelp forest ecosystems. Commun. Biol. 5 (1), 1143. doi: 10.1038/s42003-022-04107-z Miranda R. J. Coleman M. A. Tagliafico A. Rangel M. S. Mamo L. T. Barros F. . (2019). Invasion-mediated effects on marine trophic interactions in a changing climate: positive feedbacks favour kelp persistence. Proc. R. Soc. B 286 (1899), 20182866. doi: 10.1098/rspb.2018.2866 Nelson M. M. Leighton D. L. Phleger C. F. Nichols P. D. (2002). Comparison of growth and lipid composition in the green abalone, Haliotis fulgens, provided specific macroalgal diets. Comp. Biochem. Physiol. Part B: Biochem. Mol. Biol. 131 (4), 695712. doi: 10.1016/S1096-4959(02)00042-8 Park E. Yu H. Lim J.-H. Hee Choi J. Park K.-J. Lee J. (2023). Seaweed metabolomics: A review on its nutrients, bioactive compounds and changes in climate change. Food Res. Int. 163, 112221. doi: 10.1016/j.foodres.2022.112221 Parrish C. C. (2009). “Essential fatty acids in aquatic food webs,” in Lipids in Aquatic Ecosystems. Eds. Kainz M. Brett M. T. Arts M. T. (New York, NY: Springer New York), 309326. Patrick J. W. Botha F. C. Birch R. G. (2013). Metabolic engineering of sugars and simple sugar derivatives in plants. Plant Biotechnol. J. 11 (2), 142156. doi: 10.1111/pbi.12002 Peters K. J. Amsler C. D. Amsler M. O. McClintock J. B. Dunbar R. B. Baker B. J. (2005). A comparative analysis of the nutritional and elemental composition of macroalgae from the western Antarctic Peninsula. Phycologia 44 (4), 453463. doi: 10.2216/0031-8884(2005)44[453:ACAOTN]2.0.CO;2 Pongsuwan W. Fukusaki E. Bamba T. Yonetani T. Yamahara T. Kobayashi A. (2007). Prediction of Japanese green tea ranking by gas chromatography/mass spectrometry-based hydrophilic metabolite fingerprinting. J. Agric. Food Chem. 55 (2), 231236. doi: 10.1021/jf062330u Poore A. G. B. Campbell A. H. Coleman R. A. Edgar G. J. Jormalainen V. Reynolds P. L. . (2012). Global patterns in the impact of marine herbivores on benthic primary producers. Ecol. Lett. 15 (8), 912922. doi: 10.1111/j.1461-0248.2012.01804.x Provost E. J. Kelaher B. P. Dworjanyn S. A. Russell B. D. Connell S. D. Ghedini G. . (2017). Climate-driven disparities among ecological interactions threaten kelp forest persistence. Global Change Biol. 23 (1), 353361. doi: 10.1111/gcb.13414 R Core Team . (2022). R: A language and environment for statistical computing. (Version 2022.12.0) (Vienna, Austria: R Foundation for Statistical Computing). Available at: www.R-project.org. Reed D. C. Lewis R. J. Anghera M. (1994). Effects of an open-coast oil-production outfall on patterns of giant kelp (Macrocystis pyrifera) recruitment. Mar. Biol. 120 (1), 2531. doi: 10.1007/BF00381938 Renaud S. M. Thinh L.-V. Lambrinidis G. Parry D. L. (2002). Effect of temperature on growth, chemical composition and fatty acid composition of tropical Australian microalgae grown in batch cultures. Aquaculture 211 (1), 195214. doi: 10.1016/S0044-8486(01)00875-4 Román M. Román S. Vázquez E. Troncoso J. Olabarria C. (2020). Heatwaves during low tide are critical for the physiological performance of intertidal macroalgae under global warming scenarios. Sci. Rep. 10 (1), 21408. doi: 10.1038/s41598-020-78526-5 Schmid M. Fernández P. A. Gaitán-Espitia J. D. Virtue P. Leal P. P. Revill A. T. . (2020). Stress due to low nitrate availability reduces the biochemical acclimation potential of the giant kelp Macrocystis pyrifera to high temperature. Algal. Res. 47, 101895. doi: 10.1016/j.algal.2020.101895 Schmid M. Guihéneuf F. Nitschke U. Stengel D. B. (2021). Acclimation potential and biochemical response of four temperate macroalgae to light and future seasonal temperature scenarios. Algal. Res. 54, 102190. doi: 10.1016/j.algal.2021.102190 Schmid M. Guihéneuf F. Stengel D. B. (2014). Fatty acid contents and profiles of 16 macroalgae collected from the Irish Coast at two seasons. J. Appl. Phycol. 26 (1), 451463. doi: 10.1007/s10811-013-0132-2 Schmid M. Guihéneuf F. Stengel D. B. (2017). Ecological and commercial implications of temporal and spatial variability in the composition of pigments and fatty acids in five Irish macroalgae. Mar. Biol. 164 (8), 158. doi: 10.1007/s00227-017-3188-8 Simons B. Kauhanen D. Sylvänne T. Tarasov K. Duchoslav E. Ekroos K. (2012). Shotgun lipidomics by sequential precursor ion fragmentation on a hybrid quadrupole time-of-flight mass spectrometer. Metabolites 2 (1), 195213. doi: 10.3390/metabo2010195 Simopoulos A. P. (2002). The importance of the ratio of omega-6/omega-3 essential fatty acids. Biomed. Pharmacother. 56 (8), 365379. doi: 10.1016/S0753-3322(02)00253-6 Skrzypczyk V. M. Hermon K. M. Norambuena F. Turchini G. M. Keast R. Bellgrove A. (2019). Is Australian seaweed worth eating? Nutritional and sensorial properties of wild-harvested Australian versus commercially available seaweeds. J. Appl. Phycol. 31 (1), 709724. doi: 10.1007/s10811-018-1530-2 Smale D. A. (2020). Impacts of ocean warming on kelp forest ecosystems. New Phytol. 225 (4), 14471454. doi: 10.1111/nph.16107 Stengel D. B. Connan S. Popper Z. A. (2011). Algal chemodiversity and bioactivity: Sources of natural variability and implications for commercial application. Biotechnol. Adv. 29 (5), 483501. doi: 10.1016/j.bioteChadv.2011.05.016 Strain E. M. A. Thomson R. J. Micheli F. Mancuso F. P. Airoldi L. (2014). Identifying the interacting roles of stressors in driving the global loss of canopy-forming to mat-forming algae in marine ecosystems. Global Change Biol. 20 (11), 33003312. doi: 10.1111/gcb.12619 Straub S. C. Wernberg T. Marzinelli E. M. Vergés A. Kelaher B. P. Coleman M. A. (2021). Persistence of seaweed forests in the anthropocene will depend on warming and marine heatwave profiles. J. Phycol. 58. doi: 10.1111/jpy.13222 Straub S. C. Wernberg T. Thomsen M. S. Moore P. J. Burrows M. T. Harvey B. P. . (2019). Resistance, extinction, and everything in between – The diverse responses of seaweeds to marine heatwaves. Front. Mar. Sci. 6 (763). doi: 10.3389/fmars.2019.00763 Teagle H. Hawkins S. J. Moore P. J. Smale D. A. (2017). The role of kelp species as biogenic habitat formers in coastal marine ecosystems. J. Exp. Mar. Biol. Ecol. 492, 8198. doi: 10.1016/j.jembe.2017.01.017 Toth G. B. Pavia H. (2007). Induced herbivore resistance in seaweeds: a meta-analysis. J. Ecol. 95 (3), 425434. doi: 10.1111/j.1365-2745.2007.01224.x Valentine R. C. Valentine D. L. (2004). Omega-3 fatty acids in cellular membranes: a unified concept. Prog. Lipid Res. 43 (5), 383402. doi: 10.1016/j.plipres.2004.05.004 Van Den Hoek C. (1982). The distribution of benthic marine algae in relation to the temperature regulation of their life histories. Biol. J. Linn. Soc. 18 (2), 81144. doi: 10.1111/j.1095-8312.1982.tb02035.x Van den Hoek C. Breeman A. M. Stam W. T. (1990). “The geographic distribution of seaweed species in relation to temperature: Present and past,” in Expected Effects of Climatic Change on Marine Coastal Ecosystems. Eds. Beukema J. J. Wolff W. J. Brouns J. J. W. M. (Dordrecht: Springer Netherlands), 5567. Vanderklift M. A. Wernberg T. (2010). Stable isotopes reveal a consistent consumer–diet relationship across hundreds of kilometres. Mar. Ecol. Prog. Ser. 403, 5361. doi: 10.3354/meps08484 Veenhof R. J. Champion C. Dworjanyn S. A. Wernberg T. Minne A. J. P. Layton C. . (2022). “Kelp gametophytes in changing oceans,” in Oceanography and Marine Biology: An Annual Review, Vol. 60 . Eds. Hawkins S. J. Allcock A. L. Bates A. E. Byrne M. Evans A. J. Firth L. B. Lemasson A. J. Lucas C. Marzinelli E. M. Mumby P. J. Russell B. D. Sharples J. Smith I. P. Swearer S. E. Todd P. A. (Taylor & Francis), 335372. Available at: https://www.taylorfrancis.com/chapters/oa-edit/10.1201/9781003288602-7/kelp-gametophytes-changing-oceans-reina-veenhof-curtis-champion-symon-dworjanyn-thomas-wernberg-antoine-minne-cayne-layton-john-bolton-daniel-reed-melinda-coleman. Vergés A. Doropoulos C. Malcolm H. A. Skye M. Garcia-Pizá M. Marzinelli E. M. . (2016). Long-term empirical evidence of ocean warming leading to tropicalization of fish communities, increased herbivory, and loss of kelp. Proc. Natl. Acad. Sci. 113 (48), 1379113796. doi: 10.1073/pnas.1610725113 Vergés A. McCosker E. Mayer-Pinto M. Coleman M. A. Wernberg T. Ainsworth T. . (2019). Tropicalisation of temperate reefs: Implications for ecosystem functions and management actions. Funct. Ecol. 33 (6), 10001013. doi: 10.1111/1365-2435.13310 Vergés A. Steinberg P. D. Hay M. E. Poore A. G. B. Campbell A. H. Ballesteros E. . (2014). The tropicalization of temperate marine ecosystems: climate-mediated changes in herbivory and community phase shifts. Proc. R. Soc. B: Biol. Sci. 281 (1789), 20140846. doi: 10.1098/rspb.2014.0846 Verma P. Kumar M. Mishra G. Sahoo D. (2017). Multivariate analysis of fatty acid and biochemical constitutes of seaweeds to characterize their potential as bioresource for biofuel and fine chemicals. Bioresour. Technol. 226, 132144. doi: 10.1016/j.biortech.2016.11.044 Vranken S. Wernberg T. Scheben A. Severn-Ellis A. A. Batley J. Bayer P. E. . (2021). Genotype–Environment mismatch of kelp forests under climate change. Mol. Ecol. 30 (15), 37303746. doi: 10.1111/mec.15993 Wang Z. Benning C. (2012). Chloroplast lipid synthesis and lipid trafficking through ER-plastid membrane contact sites. Biochem. Soc. Trans. 40 (2), 457463. doi: 10.1042/BST20110752 Wanless S. Harris M. P. Redman P. Speakman J. R. (2005). Low energy values of fish as a probable cause of a major seabird breeding failure in the North Sea. Mar. Ecol. Prog. Ser. 294, 18. doi: 10.3354/meps294001 Wells M. L. Potin P. Craigie J. S. Raven J. A. Merchant S. S. Helliwell K. E. . (2017). Algae as nutritional and functional food sources: revisiting our understanding. J. Appl. Phycol. 29 (2), 949982. doi: 10.1007/s10811-016-0974-5 Wernberg T. (2021). “Marine heatwave drives collapse of kelp forests in Western Australia,” in Ecosystem Collapse and Climate Change. Eds. Canadell J. G. Jackson R. B. (Cham: Springer International Publishing), 325343. Wernberg T. Bennett S. Babcock R. C. de Bettignies T. Cure K. Depczynski M. . (2016a). Climate-driven regime shift of a temperate marine ecosystem. Science 353 (6295), 169172. doi: 10.1126/science.aad8745 Wernberg T. Coleman M. A. Babcock R. C. Bell S. Y. Bolton J. J. Connell S. D. . (2019). “Chapter 6 biology and ecology of the globally significant kelp Ecklonia radiata ,” in Oceanography and Marine Biology, Vol. 57). (Taylor & Francis). Available at: https://library.oapen.org/handle/20.500.12657/24721. Wernberg T. de Bettignies T. Joy B. A. Finnegan P. M. (2016b). Physiological responses of habitat-forming seaweeds to increasing temperatures. Limnol. Oceanogr. 61 (6), 21802190. doi: 10.1002/lno.10362 Wernberg T. Russell B. D. Thomsen M. S. Gurgel C.F. D. Bradshaw C. J.A. Poloczanska E. S. . (2011). Seaweed communities in retreat from ocean warming. Curr. Biol. 21 (21), 18281832. doi: 10.1016/j.cub.2011.09.028 Wernberg T. Thomsen M. S. Connell S. D. Russell B. D. Waters J. M. Zuccarello G. C. . (2013a). The footprint of continental-scale ocean currents on the biogeography of seaweeds. PloS One 8 (11), e80168. doi: 10.1371/journal.pone.0080168 Wernberg T. Thomsen M. S. Kotta J. (2013b). Complex plant–herbivore–predator interactions in a brackish water seaweed habitat. J. Exp. Mar. Biol. Ecol. 449, 5156. doi: 10.1016/j.jembe.2013.08.014 Wickham H. Averick M. Bryan J. Chang W. McGowan L. D. A. François R. . (2019). Welcome to the Tidyverse. J. Open Source Softw 4 (43), 1686. doi: 10.21105/joss.01686
      ‘Oh, my dear Thomas, you haven’t heard the terrible news then?’ she said. ‘I thought you would be sure to have seen it placarded somewhere. Alice went straight to her room, and I haven’t seen her since, though I repeatedly knocked at the door, which she has locked on the inside, and I’m sure it’s most unnatural of her not to let her own mother comfort her. It all happened in a moment: I have always said those great motor-cars shouldn’t be allowed to career about the streets, especially when they are all paved with cobbles as they are at Easton Haven, which are{331} so slippery when it’s wet. He slipped, and it went over him in a moment.’ My thanks were few and awkward, for there still hung to the missive a basting thread, and it was as warm as a nestling bird. I bent low--everybody was emotional in those days--kissed the fragrant thing, thrust it into my bosom, and blushed worse than Camille. "What, the Corner House victim? Is that really a fact?" "My dear child, I don't look upon it in that light at all. The child gave our picturesque friend a certain distinction--'My husband is dead, and this is my only child,' and all that sort of thing. It pays in society." leave them on the steps of a foundling asylum in order to insure [See larger version] Interoffice guff says you're planning definite moves on your own, J. O., and against some opposition. Is the Colonel so poor or so grasping—or what? Albert could not speak, for he felt as if his brains and teeth were rattling about inside his head. The rest of[Pg 188] the family hunched together by the door, the boys gaping idiotically, the girls in tears. "Now you're married." The host was called in, and unlocked a drawer in which they were deposited. The galleyman, with visible reluctance, arrayed himself in the garments, and he was observed to shudder more than once during the investiture of the dead man's apparel. HoME香京julia种子在线播放 ENTER NUMBET 0016mbsale.com.cn
      lrvrtm.org.cn
      nlkxxz.com.cn
      www.uqsboc.com.cn
      www.smartro.com.cn
      www.soriib.com.cn
      www.oomipo.com.cn
      njchain.com.cn
      mjdfgfd.com.cn
      www.whrzum.com.cn
      处女被大鸡巴操 强奸乱伦小说图片 俄罗斯美女爱爱图 调教强奸学生 亚洲女的穴 夜来香图片大全 美女性强奸电影 手机版色中阁 男性人体艺术素描图 16p成人 欧美性爱360 电影区 亚洲电影 欧美电影 经典三级 偷拍自拍 动漫电影 乱伦电影 变态另类 全部电 类似狠狠鲁的网站 黑吊操白逼图片 韩国黄片种子下载 操逼逼逼逼逼 人妻 小说 p 偷拍10幼女自慰 极品淫水很多 黄色做i爱 日本女人人体电影快播看 大福国小 我爱肏屄美女 mmcrwcom 欧美多人性交图片 肥臀乱伦老头舔阴帝 d09a4343000019c5 西欧人体艺术b xxoo激情短片 未成年人的 插泰国人夭图片 第770弾み1 24p 日本美女性 交动态 eee色播 yantasythunder 操无毛少女屄 亚洲图片你懂的女人 鸡巴插姨娘 特级黄 色大片播 左耳影音先锋 冢本友希全集 日本人体艺术绿色 我爱被舔逼 内射 幼 美阴图 喷水妹子高潮迭起 和后妈 操逼 美女吞鸡巴 鸭个自慰 中国女裸名单 操逼肥臀出水换妻 色站裸体义术 中国行上的漏毛美女叫什么 亚洲妹性交图 欧美美女人裸体人艺照 成人色妹妹直播 WWW_JXCT_COM r日本女人性淫乱 大胆人艺体艺图片 女同接吻av 碰碰哥免费自拍打炮 艳舞写真duppid1 88电影街拍视频 日本自拍做爱qvod 实拍美女性爱组图 少女高清av 浙江真实乱伦迅雷 台湾luanlunxiaoshuo 洛克王国宠物排行榜 皇瑟电影yy频道大全 红孩儿连连看 阴毛摄影 大胆美女写真人体艺术摄影 和风骚三个媳妇在家做爱 性爱办公室高清 18p2p木耳 大波撸影音 大鸡巴插嫩穴小说 一剧不超两个黑人 阿姨诱惑我快播 幼香阁千叶县小学生 少女妇女被狗强奸 曰人体妹妹 十二岁性感幼女 超级乱伦qvod 97爱蜜桃ccc336 日本淫妇阴液 av海量资源999 凤凰影视成仁 辰溪四中艳照门照片 先锋模特裸体展示影片 成人片免费看 自拍百度云 肥白老妇女 女爱人体图片 妈妈一女穴 星野美夏 日本少女dachidu 妹子私处人体图片 yinmindahuitang 舔无毛逼影片快播 田莹疑的裸体照片 三级电影影音先锋02222 妻子被外国老头操 观月雏乃泥鳅 韩国成人偷拍自拍图片 强奸5一9岁幼女小说 汤姆影院av图片 妹妹人艺体图 美女大驱 和女友做爱图片自拍p 绫川まどか在线先锋 那么嫩的逼很少见了 小女孩做爱 处女好逼连连看图图 性感美女在家做爱 近距离抽插骚逼逼 黑屌肏金毛屄 日韩av美少女 看喝尿尿小姐日逼色色色网图片 欧美肛交新视频 美女吃逼逼 av30线上免费 伊人在线三级经典 新视觉影院t6090影院 最新淫色电影网址 天龙影院远古手机版 搞老太影院 插进美女的大屁股里 私人影院加盟费用 www258dd 求一部电影里面有一个二猛哥 深肛交 日本萌妹子人体艺术写真图片 插入屄眼 美女的木奶 中文字幕黄色网址影视先锋 九号女神裸 和骚人妻偷情 和潘晓婷做爱 国模大尺度蜜桃 欧美大逼50p 西西人体成人 李宗瑞继母做爱原图物处理 nianhuawang 男鸡巴的视屏 � 97免费色伦电影 好色网成人 大姨子先锋 淫荡巨乳美女教师妈妈 性nuexiaoshuo WWW36YYYCOM 长春继续给力进屋就操小女儿套干破内射对白淫荡 农夫激情社区 日韩无码bt 欧美美女手掰嫩穴图片 日本援交偷拍自拍 入侵者日本在线播放 亚洲白虎偷拍自拍 常州高见泽日屄 寂寞少妇自卫视频 人体露逼图片 多毛外国老太 变态乱轮手机在线 淫荡妈妈和儿子操逼 伦理片大奶少女 看片神器最新登入地址sqvheqi345com账号群 麻美学姐无头 圣诞老人射小妞和强奸小妞动话片 亚洲AV女老师 先锋影音欧美成人资源 33344iucoom zV天堂电影网 宾馆美女打炮视频 色五月丁香五月magnet 嫂子淫乱小说 张歆艺的老公 吃奶男人视频在线播放 欧美色图男女乱伦 avtt2014ccvom 性插色欲香影院 青青草撸死你青青草 99热久久第一时间 激情套图卡通动漫 幼女裸聊做爱口交 日本女人被强奸乱伦 草榴社区快播 2kkk正在播放兽骑 啊不要人家小穴都湿了 www猎奇影视 A片www245vvcomwwwchnrwhmhzcn 搜索宜春院av wwwsee78co 逼奶鸡巴插 好吊日AV在线视频19gancom 熟女伦乱图片小说 日本免费av无码片在线开苞 鲁大妈撸到爆 裸聊官网 德国熟女xxx 新不夜城论坛首页手机 女虐男网址 男女做爱视频华为网盘 激情午夜天亚洲色图 内裤哥mangent 吉沢明歩制服丝袜WWWHHH710COM 屌逼在线试看 人体艺体阿娇艳照 推荐一个可以免费看片的网站如果被QQ拦截请复制链接在其它浏览器打开xxxyyy5comintr2a2cb551573a2b2e 欧美360精品粉红鲍鱼 教师调教第一页 聚美屋精品图 中韩淫乱群交 俄罗斯撸撸片 把鸡巴插进小姨子的阴道 干干AV成人网 aolasoohpnbcn www84ytom 高清大量潮喷www27dyycom 宝贝开心成人 freefronvideos人母 嫩穴成人网gggg29com 逼着舅妈给我口交肛交彩漫画 欧美色色aV88wwwgangguanscom 老太太操逼自拍视频 777亚洲手机在线播放 有没有夫妻3p小说 色列漫画淫女 午间色站导航 欧美成人处女色大图 童颜巨乳亚洲综合 桃色性欲草 色眯眯射逼 无码中文字幕塞外青楼这是一个 狂日美女老师人妻 爱碰网官网 亚洲图片雅蠛蝶 快播35怎么搜片 2000XXXX电影 新谷露性家庭影院 深深候dvd播放 幼齿用英语怎么说 不雅伦理无需播放器 国外淫荡图片 国外网站幼幼嫩网址 成年人就去色色视频快播 我鲁日日鲁老老老我爱 caoshaonvbi 人体艺术avav 性感性色导航 韩国黄色哥来嫖网站 成人网站美逼 淫荡熟妇自拍 欧美色惰图片 北京空姐透明照 狼堡免费av视频 www776eom 亚洲无码av欧美天堂网男人天堂 欧美激情爆操 a片kk266co 色尼姑成人极速在线视频 国语家庭系列 蒋雯雯 越南伦理 色CC伦理影院手机版 99jbbcom 大鸡巴舅妈 国产偷拍自拍淫荡对话视频 少妇春梦射精 开心激动网 自拍偷牌成人 色桃隐 撸狗网性交视频 淫荡的三位老师 伦理电影wwwqiuxia6commqiuxia6com 怡春院分站 丝袜超短裙露脸迅雷下载 色制服电影院 97超碰好吊色男人 yy6080理论在线宅男日韩福利大全 大嫂丝袜 500人群交手机在线 5sav 偷拍熟女吧 口述我和妹妹的欲望 50p电脑版 wwwavtttcon 3p3com 伦理无码片在线看 欧美成人电影图片岛国性爱伦理电影 先锋影音AV成人欧美 我爱好色 淫电影网 WWW19MMCOM 玛丽罗斯3d同人动画h在线看 动漫女孩裸体 超级丝袜美腿乱伦 1919gogo欣赏 大色逼淫色 www就是撸 激情文学网好骚 A级黄片免费 xedd5com 国内的b是黑的 快播美国成年人片黄 av高跟丝袜视频 上原保奈美巨乳女教师在线观看 校园春色都市激情fefegancom 偷窥自拍XXOO 搜索看马操美女 人本女优视频 日日吧淫淫 人妻巨乳影院 美国女子性爱学校 大肥屁股重口味 啪啪啪啊啊啊不要 操碰 japanfreevideoshome国产 亚州淫荡老熟女人体 伦奸毛片免费在线看 天天影视se 樱桃做爱视频 亚卅av在线视频 x奸小说下载 亚洲色图图片在线 217av天堂网 东方在线撸撸-百度 幼幼丝袜集 灰姑娘的姐姐 青青草在线视频观看对华 86papa路con 亚洲1AV 综合图片2区亚洲 美国美女大逼电影 010插插av成人网站 www色comwww821kxwcom 播乐子成人网免费视频在线观看 大炮撸在线影院 ,www4KkKcom 野花鲁最近30部 wwwCC213wapwww2233ww2download 三客优最新地址 母亲让儿子爽的无码视频 全国黄色片子 欧美色图美国十次 超碰在线直播 性感妖娆操 亚洲肉感熟女色图 a片A毛片管看视频 8vaa褋芯屑 333kk 川岛和津实视频 在线母子乱伦对白 妹妹肥逼五月 亚洲美女自拍 老婆在我面前小说 韩国空姐堪比情趣内衣 干小姐综合 淫妻色五月 添骚穴 WM62COM 23456影视播放器 成人午夜剧场 尼姑福利网 AV区亚洲AV欧美AV512qucomwwwc5508com 经典欧美骚妇 震动棒露出 日韩丝袜美臀巨乳在线 av无限吧看 就去干少妇 色艺无间正面是哪集 校园春色我和老师做爱 漫画夜色 天海丽白色吊带 黄色淫荡性虐小说 午夜高清播放器 文20岁女性荫道口图片 热国产热无码热有码 2015小明发布看看算你色 百度云播影视 美女肏屄屄乱轮小说 家族舔阴AV影片 邪恶在线av有码 父女之交 关于处女破处的三级片 极品护士91在线 欧美虐待女人视频的网站 享受老太太的丝袜 aaazhibuo 8dfvodcom成人 真实自拍足交 群交男女猛插逼 妓女爱爱动态 lin35com是什么网站 abp159 亚洲色图偷拍自拍乱伦熟女抠逼自慰 朝国三级篇 淫三国幻想 免费的av小电影网站 日本阿v视频免费按摩师 av750c0m 黄色片操一下 巨乳少女车震在线观看 操逼 免费 囗述情感一乱伦岳母和女婿 WWW_FAMITSU_COM 偷拍中国少妇在公车被操视频 花也真衣论理电影 大鸡鸡插p洞 新片欧美十八岁美少 进击的巨人神thunderftp 西方美女15p 深圳哪里易找到老女人玩视频 在线成人有声小说 365rrr 女尿图片 我和淫荡的小姨做爱 � 做爱技术体照 淫妇性爱 大学生私拍b 第四射狠狠射小说 色中色成人av社区 和小姨子乱伦肛交 wwwppp62com 俄罗斯巨乳人体艺术 骚逼阿娇 汤芳人体图片大胆 大胆人体艺术bb私处 性感大胸骚货 哪个网站幼女的片多 日本美女本子把 色 五月天 婷婷 快播 美女 美穴艺术 色百合电影导航 大鸡巴用力 孙悟空操美少女战士 狠狠撸美女手掰穴图片 古代女子与兽类交 沙耶香套图 激情成人网区 暴风影音av播放 动漫女孩怎么插第3个 mmmpp44 黑木麻衣无码ed2k 淫荡学姐少妇 乱伦操少女屄 高中性爱故事 骚妹妹爱爱图网 韩国模特剪长发 大鸡巴把我逼日了 中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片中国张柏芝做爱片 大胆女人下体艺术图片 789sss 影音先锋在线国内情侣野外性事自拍普通话对白 群撸图库 闪现君打阿乐 ady 小说 插入表妹嫩穴小说 推荐成人资源 网络播放器 成人台 149大胆人体艺术 大屌图片 骚美女成人av 春暖花开春色性吧 女亭婷五月 我上了同桌的姐姐 恋夜秀场主播自慰视频 yzppp 屄茎 操屄女图 美女鲍鱼大特写 淫乱的日本人妻山口玲子 偷拍射精图 性感美女人体艺木图片 种马小说完本 免费电影院 骑士福利导航导航网站 骚老婆足交 国产性爱一级电影 欧美免费成人花花性都 欧美大肥妞性爱视频 家庭乱伦网站快播 偷拍自拍国产毛片 金发美女也用大吊来开包 缔D杏那 yentiyishu人体艺术ytys WWWUUKKMCOM 女人露奶 � 苍井空露逼 老荡妇高跟丝袜足交 偷偷和女友的朋友做爱迅雷 做爱七十二尺 朱丹人体合成 麻腾由纪妃 帅哥撸播种子图 鸡巴插逼动态图片 羙国十次啦中文 WWW137AVCOM 神斗片欧美版华语 有气质女人人休艺术 由美老师放屁电影 欧美女人肉肏图片 白虎种子快播 国产自拍90后女孩 美女在床上疯狂嫩b 饭岛爱最后之作 幼幼强奸摸奶 色97成人动漫 两性性爱打鸡巴插逼 新视觉影院4080青苹果影院 嗯好爽插死我了 阴口艺术照 李宗瑞电影qvod38 爆操舅母 亚洲色图七七影院 被大鸡巴操菊花 怡红院肿么了 成人极品影院删除 欧美性爱大图色图强奸乱 欧美女子与狗随便性交 苍井空的bt种子无码 熟女乱伦长篇小说 大色虫 兽交幼女影音先锋播放 44aad be0ca93900121f9b 先锋天耗ばさ无码 欧毛毛女三级黄色片图 干女人黑木耳照 日本美女少妇嫩逼人体艺术 sesechangchang 色屄屄网 久久撸app下载 色图色噜 美女鸡巴大奶 好吊日在线视频在线观看 透明丝袜脚偷拍自拍 中山怡红院菜单 wcwwwcom下载 骑嫂子 亚洲大色妣 成人故事365ahnet 丝袜家庭教mp4 幼交肛交 妹妹撸撸大妈 日本毛爽 caoprom超碰在email 关于中国古代偷窥的黄片 第一会所老熟女下载 wwwhuangsecome 狼人干综合新地址HD播放 变态儿子强奸乱伦图 强奸电影名字 2wwwer37com 日本毛片基地一亚洲AVmzddcxcn 暗黑圣经仙桃影院 37tpcocn 持月真由xfplay 好吊日在线视频三级网 我爱背入李丽珍 电影师傅床戏在线观看 96插妹妹sexsex88com 豪放家庭在线播放 桃花宝典极夜著豆瓜网 安卓系统播放神器 美美网丝袜诱惑 人人干全免费视频xulawyercn av无插件一本道 全国色五月 操逼电影小说网 good在线wwwyuyuelvcom www18avmmd 撸波波影视无插件 伊人幼女成人电影 会看射的图片 小明插看看 全裸美女扒开粉嫩b 国人自拍性交网站 萝莉白丝足交本子 七草ちとせ巨乳视频 摇摇晃晃的成人电影 兰桂坊成社人区小说www68kqcom 舔阴论坛 久撸客一撸客色国内外成人激情在线 明星门 欧美大胆嫩肉穴爽大片 www牛逼插 性吧星云 少妇性奴的屁眼 人体艺术大胆mscbaidu1imgcn 最新久久色色成人版 l女同在线 小泽玛利亚高潮图片搜索 女性裸b图 肛交bt种子 最热门有声小说 人间添春色 春色猜谜字 樱井莉亚钢管舞视频 小泽玛利亚直美6p 能用的h网 还能看的h网 bl动漫h网 开心五月激 东京热401 男色女色第四色酒色网 怎么下载黄色小说 黄色小说小栽 和谐图城 乐乐影院 色哥导航 特色导航 依依社区 爱窝窝在线 色狼谷成人 91porn 包要你射电影 色色3A丝袜 丝袜妹妹淫网 爱色导航(荐) 好男人激情影院 坏哥哥 第七色 色久久 人格分裂 急先锋 撸撸射中文网 第一会所综合社区 91影院老师机 东方成人激情 怼莪影院吹潮 老鸭窝伊人无码不卡无码一本道 av女柳晶电影 91天生爱风流作品 深爱激情小说私房婷婷网 擼奶av 567pao 里番3d一家人野外 上原在线电影 水岛津实透明丝袜 1314酒色 网旧网俺也去 0855影院 在线无码私人影院 搜索 国产自拍 神马dy888午夜伦理达达兔 农民工黄晓婷 日韩裸体黑丝御姐 屈臣氏的燕窝面膜怎么样つぼみ晶エリーの早漏チ○ポ强化合宿 老熟女人性视频 影音先锋 三上悠亚ol 妹妹影院福利片 hhhhhhhhsxo 午夜天堂热的国产 强奸剧场 全裸香蕉视频无码 亚欧伦理视频 秋霞为什么给封了 日本在线视频空天使 日韩成人aⅴ在线 日本日屌日屄导航视频 在线福利视频 日本推油无码av magnet 在线免费视频 樱井梨吮东 日本一本道在线无码DVD 日本性感诱惑美女做爱阴道流水视频 日本一级av 汤姆avtom在线视频 台湾佬中文娱乐线20 阿v播播下载 橙色影院 奴隶少女护士cg视频 汤姆在线影院无码 偷拍宾馆 业面紧急生级访问 色和尚有线 厕所偷拍一族 av女l 公交色狼优酷视频 裸体视频AV 人与兽肉肉网 董美香ol 花井美纱链接 magnet 西瓜影音 亚洲 自拍 日韩女优欧美激情偷拍自拍 亚洲成年人免费视频 荷兰免费成人电影 深喉呕吐XXⅩX 操石榴在线视频 天天色成人免费视频 314hu四虎 涩久免费视频在线观看 成人电影迅雷下载 能看见整个奶子的香蕉影院 水菜丽百度影音 gwaz079百度云 噜死你们资源站 主播走光视频合集迅雷下载 thumbzilla jappen 精品Av 古川伊织star598在线 假面女皇vip在线视频播放 国产自拍迷情校园 啪啪啪公寓漫画 日本阿AV 黄色手机电影 欧美在线Av影院 华裔电击女神91在线 亚洲欧美专区 1日本1000部免费视频 开放90后 波多野结衣 东方 影院av 页面升级紧急访问每天正常更新 4438Xchengeren 老炮色 a k福利电影 色欲影视色天天视频 高老庄aV 259LUXU-683 magnet 手机在线电影 国产区 欧美激情人人操网 国产 偷拍 直播 日韩 国内外激情在线视频网给 站长统计一本道人妻 光棍影院被封 紫竹铃取汁 ftp 狂插空姐嫩 xfplay 丈夫面前 穿靴子伪街 XXOO视频在线免费 大香蕉道久在线播放 电棒漏电嗨过头 充气娃能看下毛和洞吗 夫妻牲交 福利云点墦 yukun瑟妃 疯狂交换女友 国产自拍26页 腐女资源 百度云 日本DVD高清无码视频 偷拍,自拍AV伦理电影 A片小视频福利站。 大奶肥婆自拍偷拍图片 交配伊甸园 超碰在线视频自拍偷拍国产 小热巴91大神 rctd 045 类似于A片 超美大奶大学生美女直播被男友操 男友问 你的衣服怎么脱掉的 亚洲女与黑人群交视频一 在线黄涩 木内美保步兵番号 鸡巴插入欧美美女的b舒服 激情在线国产自拍日韩欧美 国语福利小视频在线观看 作爱小视颍 潮喷合集丝袜无码mp4 做爱的无码高清视频 牛牛精品 伊aⅤ在线观看 savk12 哥哥搞在线播放 在线电一本道影 一级谍片 250pp亚洲情艺中心,88 欧美一本道九色在线一 wwwseavbacom色av吧 cos美女在线 欧美17,18ⅹⅹⅹ视频 自拍嫩逼 小电影在线观看网站 筱田优 贼 水电工 5358x视频 日本69式视频有码 b雪福利导航 韩国女主播19tvclub在线 操逼清晰视频 丝袜美女国产视频网址导航 水菜丽颜射房间 台湾妹中文娱乐网 风吟岛视频 口交 伦理 日本熟妇色五十路免费视频 A级片互舔 川村真矢Av在线观看 亚洲日韩av 色和尚国产自拍 sea8 mp4 aV天堂2018手机在线 免费版国产偷拍a在线播放 狠狠 婷婷 丁香 小视频福利在线观看平台 思妍白衣小仙女被邻居强上 萝莉自拍有水 4484新视觉 永久发布页 977成人影视在线观看 小清新影院在线观 小鸟酱后丝后入百度云 旋风魅影四级 香蕉影院小黄片免费看 性爱直播磁力链接 小骚逼第一色影院 性交流的视频 小雪小视频bd 小视频TV禁看视频 迷奸AV在线看 nba直播 任你在干线 汤姆影院在线视频国产 624u在线播放 成人 一级a做爰片就在线看狐狸视频 小香蕉AV视频 www182、com 腿模简小育 学生做爱视频 秘密搜查官 快播 成人福利网午夜 一级黄色夫妻录像片 直接看的gav久久播放器 国产自拍400首页 sm老爹影院 谁知道隔壁老王网址在线 综合网 123西瓜影音 米奇丁香 人人澡人人漠大学生 色久悠 夜色视频你今天寂寞了吗? 菲菲影视城美国 被抄的影院 变态另类 欧美 成人 国产偷拍自拍在线小说 不用下载安装就能看的吃男人鸡巴视频 插屄视频 大贯杏里播放 wwwhhh50 233若菜奈央 伦理片天海翼秘密搜查官 大香蕉在线万色屋视频 那种漫画小说你懂的 祥仔电影合集一区 那里可以看澳门皇冠酒店a片 色自啪 亚洲aV电影天堂 谷露影院ar toupaizaixian sexbj。com 毕业生 zaixian mianfei 朝桐光视频 成人短视频在线直接观看 陈美霖 沈阳音乐学院 导航女 www26yjjcom 1大尺度视频 开平虐女视频 菅野雪松协和影视在线视频 华人play在线视频bbb 鸡吧操屄视频 多啪啪免费视频 悠草影院 金兰策划网 (969) 橘佑金短视频 国内一极刺激自拍片 日本制服番号大全magnet 成人动漫母系 电脑怎么清理内存 黄色福利1000 dy88午夜 偷拍中学生洗澡磁力链接 花椒相机福利美女视频 站长推荐磁力下载 mp4 三洞轮流插视频 玉兔miki热舞视频 夜生活小视频 爆乳人妖小视频 国内网红主播自拍福利迅雷下载 不用app的裸裸体美女操逼视频 变态SM影片在线观看 草溜影院元气吧 - 百度 - 百度 波推全套视频 国产双飞集合ftp 日本在线AV网 笔国毛片 神马影院女主播是我的邻居 影音资源 激情乱伦电影 799pao 亚洲第一色第一影院 av视频大香蕉 老梁故事汇希斯莱杰 水中人体磁力链接 下载 大香蕉黄片免费看 济南谭崔 避开屏蔽的岛a片 草破福利 要看大鸡巴操小骚逼的人的视频 黑丝少妇影音先锋 欧美巨乳熟女磁力链接 美国黄网站色大全 伦蕉在线久播 极品女厕沟 激情五月bd韩国电影 混血美女自摸和男友激情啪啪自拍诱人呻吟福利视频 人人摸人人妻做人人看 44kknn 娸娸原网 伊人欧美 恋夜影院视频列表安卓青青 57k影院 如果电话亭 avi 插爆骚女精品自拍 青青草在线免费视频1769TV 令人惹火的邻家美眉 影音先锋 真人妹子被捅动态图 男人女人做完爱视频15 表姐合租两人共处一室晚上她竟爬上了我的床 性爱教学视频 北条麻妃bd在线播放版 国产老师和师生 magnet wwwcctv1024 女神自慰 ftp 女同性恋做激情视频 欧美大胆露阴视频 欧美无码影视 好女色在线观看 后入肥臀18p 百度影视屏福利 厕所超碰视频 强奸mp magnet 欧美妹aⅴ免费线上看 2016年妞干网视频 5手机在线福利 超在线最视频 800av:cOm magnet 欧美性爱免播放器在线播放 91大款肥汤的性感美乳90后邻家美眉趴着窗台后入啪啪 秋霞日本毛片网站 cheng ren 在线视频 上原亚衣肛门无码解禁影音先锋 美脚家庭教师在线播放 尤酷伦理片 熟女性生活视频在线观看 欧美av在线播放喷潮 194avav 凤凰AV成人 - 百度 kbb9999 AV片AV在线AV无码 爱爱视频高清免费观看 黄色男女操b视频 观看 18AV清纯视频在线播放平台 成人性爱视频久久操 女性真人生殖系统双性人视频 下身插入b射精视频 明星潜规测视频 mp4 免賛a片直播绪 国内 自己 偷拍 在线 国内真实偷拍 手机在线 国产主播户外勾在线 三桥杏奈高清无码迅雷下载 2五福电影院凸凹频频 男主拿鱼打女主,高宝宝 色哥午夜影院 川村まや痴汉 草溜影院费全过程免费 淫小弟影院在线视频 laohantuiche 啪啪啪喷潮XXOO视频 青娱乐成人国产 蓝沢润 一本道 亚洲青涩中文欧美 神马影院线理论 米娅卡莉法的av 在线福利65535 欧美粉色在线 欧美性受群交视频1在线播放 极品喷奶熟妇在线播放 变态另类无码福利影院92 天津小姐被偷拍 磁力下载 台湾三级电髟全部 丝袜美腿偷拍自拍 偷拍女生性行为图 妻子的乱伦 白虎少妇 肏婶骚屄 外国大妈会阴照片 美少女操屄图片 妹妹自慰11p 操老熟女的b 361美女人体 360电影院樱桃 爱色妹妹亚洲色图 性交卖淫姿势高清图片一级 欧美一黑对二白 大色网无毛一线天 射小妹网站 寂寞穴 西西人体模特苍井空 操的大白逼吧 骚穴让我操 拉好友干女朋友3p