ALBERT

Description: Compositional variations of amphibole stratigraphically recovered from multiple eruptions at a given volcano have a great potential to archive long-term magmatic processes in its crustal plumbing system. Calcic amphibole is a ubiquitous yet chemically and texturally diverse mineral at Mount St. Helens (MSH), where it occurs in dacites and in co-magmatic enclaves throughout the Spirit Lake stage (last ~4000 years of eruptive history). It forms three populations with distinct geochemical trends in key major and trace elements, which are subdivided into a high-Al (11–14.5 wt% Al2O3), a medium-Al (10–12.5 wt% Al2O3), and a low-Al (7.5–10 wt% Al2O3) amphibole population. The oldest investigated tephra record (Smith Creek period, 3900–3300 years BP) yields a bimodal amphibole distribution in which lower-crustal, high-Al amphibole cores (crystallized dominantly from basaltic andesite to andesite melts) and upper-crustal, low-Al amphibole rims (crystallized from rhyolitic melt) document occasional recharge of a shallow silicic mush by a more mafic melt from a lower-crustal reservoir. The sudden appearance of medium-Al amphiboles enriched in incompatible trace elements in eruptive periods younger than 2900 years BP is associated with a change in reservoir conditions toward hotter and drier magmas, which indicates recharge of the shallow silicic reservoir by basaltic melt enriched in incompatible elements. Deep-crystallizing, high-Al amphibole, however, appears mostly unaffected by such incompatible-element-enriched basaltic recharge, suggesting that these basalts bypass the lower crustal reservoir. This could be the result of the eastward offset position of the lower crustal reservoir relative to the upper crustal storage zone underneath the MSH edifice. Amphibole has proven to be a sensitive geochemical archive for uncovering storage conditions of magmas at MSH. In agreement with geophysical observations, storage and differentiation have occurred in two main zones: an upper crustal and lower crustal reservoir (the lower one being chemically less evolved). The upper crustal silicic reservoir, offset to the west of the lower crustal reservoir, has captured compositionally unusual mafic recharge (drier, hotter, and enriched in incompatible trace elements in comparison to the typical parental magmas in the region), resulting in an increased chemical diversity of amphiboles and their carrier intermediate magmas, in the last ~3000 years of MSH’s volcanic record.

Description: Current global warming results in rising sea-water temperatures, and the loss of sea ice in arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon-sequencing data from sedimentary ancient DNA (sedaDNA) of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial–interglacial cycle, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by pico-sized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic pico-sized chlorophytes and Triparmaceae. In contrast, the glacial period is characterized by micro-sized phototrophic protists, including sea-ice associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, and resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift towards the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community’s potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea’s function as an effective carbon sink.

Description: Several species from various zooplankton taxa perform seasonal vertical migrations (SVM) of typically several hundred meters between the surface layer and overwintering depths, particularly in high-latitude regions. We use OPtimality-based PLAnkton (OPPLA) ecosystem model) to simulate SVM behavior in zooplankton in the Labrador Sea. Zooplankton in OPPLA is a generic functional group without life cycle, which facilitates analyzing SVM evolutionary stability and interactions between SVM and the plankton ecosystem. A sensitivity analysis of SVM-related parameters reveals that SVM can amplify the seasonal variations of phytoplankton and zooplankton and enhance the reduction of summer surface nutrient concentrations. SVM is often explained as a strategy to reduce exposure to visual predators during winter. We find that species doing SVM can persist and even dominate the summer-time zooplankton community, even in the presence of Stayers, which have the same traits as the migrators, but do not perform SVM. The advantage of SVM depends strongly on the timing of the seasonal migrations, particularly the day of ascent. The presence of higher (visual) predators tends to suppress the Stayers in our simulations, whereas the SVM strategy can persist in the presence of non-migrating species even without higher predators.

Description: Background and Aims Dysregulated mineral homeostasis is common in chronic kidney disease (CKD) and associated with bone demineralization and vascular calcification. The balance between bone formation and resorption, which reflect the bone calcium (Ca) balance (BCaB), cannot be determined without bone biopsy which is invasive and not easily repeatable. Recently, we have shown that stable (i.e. non-radioactive) Ca isotopes, 42Ca and 44Ca, can be measured in serum and their ratio (δ44/42Caserum) quantitatively determines net bone gain or loss of Ca. Thus, when bone formation exceeds bone resorption, the net BCaB is positive and δ44/42Caserum is high, and when bone resorption is the predominant process δ44/42Caserum is low compared to age-matched controls. In this study we compared δ44/42Caserum against δ44/42Cabone and arterial biopsy samples (δ44/42Caartery) and the sensitivity of δ44/42Ca in predicting changes in bone histology.

Description: The seismic receiver function (RF) technique is widely used as an economic method to image earth's deep interior in a large number of seismic experiments. P-wave receiver functions (RFs) constrain crustal thickness and average Vp/Vs in the crust by analysis of the Ps phase and multiples (reflected/converted waves) from the Moho. Regional studies often show significant differences between the Moho depth constrained by RF and by reflection/refraction methods. We compare the results from RF and controlled source seismology for the Baikal Rift Zone by calculating 1480 synthetic RFs for a seismic refraction/reflection velocity model and processing them with two common RF techniques [H–κ and Common Conversion Point (CCP) stacking]. We compare the resulting synthetic RF structure with the velocity model, a density model (derived from gravity and the velocity model), and with observed RFs. Our results demonstrate that the use of different frequency filters, the presence of complex phases from sediments and gradual changes in the properties of crustal layers can lead to erroneous interpretation of RFs and incorrect geological interpretations. We suggest that the interpretation of RFs should be combined with other geophysical methods, in particular in complex tectonic regions and that the long-wavelength Bouguer gravity anomaly signal may provide effective calibration for the determination of the correct Moho depth from RF results. We propose and validate a new automated, efficient method for this calibration.

Description: Eutrophication-driven harmful algal blooms (HABs) can have secondary effects on larval fishes that rely on estuaries as nurseries. However, few studies worldwide have quantified these effects despite the global rise in eutrophication. This study presents a novel approach using biochemical body condition analyses to evaluate the impact of HABs on the growth and body condition of the larvae of an estuarine resident fish. Recurrent phytoplankton blooms of Heterosigma akashiwo occur in the warm-temperate Sundays Estuary on the southeast coast of South Africa. The response in body condition and assemblage structure on larval estuarine roundherring (Gilchristella aestuaria) was measured in conjunction with bloom conditions, water quality and zooplanktonic prey and predators. Larvae and early juveniles were sampled during varying intensity levels, duration and frequency of hypereutrophic blooms. This study demonstrated that extensive HABs could significantly impact larval roundherring, G. aestuaria, by decreasing larval nutritional condition and limiting their growth, resulting in poor grow-out into the juvenile phase. Poor condition and growth may likely affect recruitment success to adult populations, and since G. aestuaria is an important forage fish and zooplanktivore, poor recruitment will hold consequences for estuarine food webs.

Description: The neritic-oceanic squid Illex argentinus supports one of the largest fisheries in the Southwest Atlantic. It is characterized by extensive migrations across the Patagonian Shelf and complex population structure comprising distinct seasonal spawning groups. To address uncertainty as to the demographic independence of these groups that may compromise sustainable management, a multidisciplinary approach was applied integrating statolith ageing with genome-wide single-nucleotide polymorphism (SNP) analysis. To obtain complete coverage of the spawning groups, sampling was carried out at multiple times during the 2020 fishing season and covered a large proportion of the species' range across the Patagonian Shelf. Statolith and microstructure analysis revealed three distinct seasonal spawning groups of winter-, spring-, and summer-hatched individuals. Subgroups were identified within each seasonal group, with statolith microstructure indicating differences in environmental conditions during ontogeny. Analysis of 〉10 000 SNPs reported no evidence of neutral or non-neutral genetic structure among the various groups. These findings indicate that I. argentinus across the Patagonian Shelf belong to one genetic population and a collaborative management strategy involving international stakeholders is required. The connectivity among spawning groups may represent a "bet-hedging" mechanism important for population resilience.

Description: The dynamics of marine systems at decadal scales are notoriously hard to predict—hence references to this timescale as the “grey zone” for ocean prediction. Nevertheless, decadal-scale prediction is a rapidly developing field with an increasing number of applications to help guide ocean stewardship and sustainable use of marine environments. Such predictions can provide industry and managers with information more suited to support planning and management over strategic timeframes, as compared to seasonal forecasts or long-term (century-scale) predictions. The most significant advances in capability for decadal-scale prediction over recent years have been for ocean physics and biogeochemistry, with some notable advances in ecological prediction skill. In this paper, we argue that the process of “lighting the grey zone” by providing improved predictions at decadal scales should also focus on including human dimensions in prediction systems to better meet the needs and priorities of end users. Our paper reviews information needs for decision-making at decadal scales and assesses current capabilities for meeting these needs. We identify key gaps in current capabilities, including the particular challenge of integrating human elements into decadal prediction systems. We then suggest approaches for overcoming these challenges and gaps, highlighting the important role of co-production of tools and scenarios, to build trust and ensure uptake with end users of decadal prediction systems. We also highlight opportunities for combining narratives and quantitative predictions to better incorporate the human dimension in future efforts to light the grey zone of decadal-scale prediction.

Description: Post-collisional volcanism contains important clues for understanding the processes that prevail in orogenic belts, including those in the mantle and the uplift and collapse of continents. Here we report new geochronological and geochemical data for a suite of post-collisional Miocene to Pleistocene volcanic rocks from northwest Iran. Four groups of volcanic rocks can be distinguished according to their geochemical and isotopic signatures, including: (1) Miocene depleted lavas with high Nd and Hf but low Pb and Sr isotopic ratios, (2) less depleted lavas with quite variable Pb isotopic composition, (3) lavas with non-radiogenic Nd and Hf isotopic values, but highly radiogenic Sr and Pb isotopic composition, and (4) Pleistocene adakitic rocks with depleted isotopic signatures. The isotopic data reveal that the Miocene rocks are derived from asthenospheric and highly heterogeneous sub-continental lithospheric mantle sources. Evidence suggests that the lithospheric mantle contains recycled upper continental material and is isotopically similar to the enriched mantle two (EMII) end-member. Analysis of Sr-Nd-Pb-Hf-O isotopes in both mineral and rock groundmass, in conjunction with energy-constrained assimilation and fractional crystallization (EC-AFC) numerical modeling, demonstrates that the incorporation of continental crust during magma fractionation via AFC had an insignificant impact on the isotopic composition of the Miocene lavas. Moreover, adakites are the youngest rocks and show a geochemical signature consistent with the partial melting of a young and mafic continental lower crust. Both seismological data and geochemical signatures on these Miocene to Pleistocene volcanic rocks indicate the initiation of asthenospheric upwelling and orogen uplift in the Arabia-Eurasia collision zone, which occurred after slab break-off, following the Neotethyan closure.

Description: Primary andesitic magmas could be an important component of arc magma genesis and might have played a key role in the advent of continents. Recent studies hypothesized that primary andesitic magmas occur in the oceanic arc, where the crust is thin. The Kermadec arc has the thinnest crust among all the studied oceanic arcs (〈15 km in thickness); however, there are no studies that corroborate the formation of primary andesitic magmas in the arc. The aim of this study is to develop a better understanding of primary andesites in oceanic arcs through the petrology of the Kermadec arc. Here, we present the petrology of volcanic rocks dredged from the Kibblewhite Volcano in the Kermadec arc during the R/V SONNE SO-255 expedition in 2017. Magma types range from andesite to rhyolite at the Kibblewhite Volcano, but basalts dominate at the neighboring cones. This study focuses on magnesian andesites from the northeastern flank of this volcano. The magnesian andesites are nearly aphyric and plagioclase free but contain microphenocrysts of olivine (Fo84–86) and clinopyroxene (Mg# = 81–87). Using olivine addition models, the primary magmas were estimated to contain 55–56 wt % SiO2 and 10–12 wt % MgO, similar to the high-Mg andesites observed in other convergent plate margins, indicating the generation of primary andesitic magma beneath the Kibblewhite Volcano. The trace element and isotopic characteristics of the magnesian andesites are typical of volcanic rocks from the Kermadec arc. This indicates that the subduction of a young plate or melting of a pyroxenitic source is not necessary to produce magnesian andesites. Instead, we propose that the magnesian andesites were produced by the direct melting of the uppermost mantle of the Kermadec arc. The thin crust of the Kermadec arc should yield low-pressure conditions in the uppermost mantle, allowing the sub-arc mantle to generate primary andesitic melts. This study supports the hypothesis that primary andesitic magmas generate in the arc where the crust is thin and provides a new insight into the magma genesis of the Kermadec arc.

Description: Mitochondrial genomes (mitogenomes) of flowering plants are composed of multiple chromosomes. Recombination within and between the mitochondrial chromosomes may generate diverse DNA molecules termed isoforms. The isoform copy number and composition can be dynamic within and among individual plants due to uneven replication and homologous recombination. Nonetheless, despite their functional importance, the level of mitogenome conservation within species remains understudied. Whether the ontogenetic variation translates to evolution of mitogenome composition over generations is currently unknown. Here we show that the mitogenome composition of the seagrass Zostera marina is conserved among worldwide populations that diverged ca. 350,000 years ago. Using long-read sequencing, we characterized the Z. marina mitochondrial genome and inferred the repertoire of recombination-induced configurations. To characterize the mitochondrial genome architecture worldwide and study its evolution, we examined the mitogenome in Z. marina meristematic region sampled in 16 populations from the Pacific and Atlantic oceans. Our results reveal a striking similarity in the isoform relative copy number, indicating a high conservation of the mitogenome composition among distantly related populations and within the plant germline, despite a notable variability during individual ontogenesis. Our study supplies a link between observations of dynamic mitogenomes at the level of plant individuals and long-term mitochondrial evolution.

Description: Marine algae are central to global carbon fixation and their productivity is dictated largely by resource availability. Reduced nutrient availability is predicted for vast oceanic regions as an outcome of climate change, however there is much to learn regarding response mechanisms of the tiny picoplankton that thrive in these environments, especially eukaryotic phytoplankton. Here, we investigate responses of the picoeukaryote Micromonas commoda, a green alga found throughout subtropical and tropical oceans. Under shifting phosphate (P) availability scenarios, transcriptomic analyses revealed altered expression of transfer RNA (tRNA) modification enzymes and biased codon usage of transcripts more abundant during P-limiting versus P-replete conditions, consistent with the role of tRNA modifications in regulating codon recognition. To associate the observed shift in expression of the tRNA modification enzyme complement with the tRNAs encoded by M. commoda, we also determined the tRNA repertoire of this alga revealing potential targets of the modification enzymes. Codon usage bias was particularly pronounced in transcripts encoding proteins with direct roles in managing P-limitation and photosystem-associated proteins that have ill-characterized putative functions in “light stress”. The observed codon usage bias corresponds to a proposed stress response mechanism in which the interplay between stress-induced changes in tRNA modifications and skewed codon usage in certain essential response genes drives preferential translation of the encoded proteins. Collectively, we expose a potential underlying mechanism for achieving growth under enhanced nutrient limitation, that extends beyond the catalog of up- or down-regulated protein-encoding genes, to the cell biological controls that underpin acclimation to changing environmental conditions.

Description: Omic BON is a thematic Biodiversity Observation Network under the Group on Earth Observations Biodiversity Observation Network (GEO BON), focused on coordinating the observation of biomolecules in organisms and the environment. Our founding partners include representatives from national, regional, and global observing systems; standards organizations; and data and sample management infrastructures. By coordinating observing strategies, methods, and data flows, Omic BON will facilitate the co-creation of a global omics meta-observatory to generate actionable knowledge. Here, we present key elements of Omic BON's founding charter and first activities.

Description: Young grapevines (Vitis vinifera) suffer and eventually can die from the crown gall disease caused by the plant pathogen Allorhizobium vitis (Rhizobiaceae). Virulent members of A. vitis harbor a tumor-inducing plasmid and induce formation of crown galls due to the oncogenes encoded on the transfer DNA. The expression of oncogenes in transformed host cells induces unregulated cell proliferation and metabolic and physiological changes. The crown gall produces opines uncommon to plants, which provide an important nutrient source for A. vitis harboring opine catabolism enzymes. Crown galls host a distinct bacterial community, and the mechanisms establishing a crown gall–specific bacterial community are currently unknown. Thus, we were interested in whether genes homologous to those of the tumor-inducing plasmid coexist in the genomes of the microbial species coexisting in crown galls. We isolated 8 bacterial strains from grapevine crown galls, sequenced their genomes, and tested their virulence and opine utilization ability in bioassays. In addition, the 8 genome sequences were compared with 34 published bacterial genomes, including closely related plant-associated bacteria not from crown galls. Homologous genes for virulence and opine anabolism were only present in the virulent Rhizobiaceae. In contrast, homologs of the opine catabolism genes were present in all strains including the nonvirulent members of the Rhizobiaceae and non-Rhizobiaceae. Gene neighborhood and sequence identity of the opine degradation cluster of virulent and nonvirulent strains together with the results of the opine utilization assay support the important role of opine utilization for cocolonization in crown galls, thereby shaping the crown gall community.

Description: Eutrophication in marine waters is traditionally assessed by checking if nutrients, algal biomass and oxygen are below/above a given threshold. However, increased biomass, nutrient concentrations and oxygen demand do not lead to undesirable environmental effects if the flow of carbon/energy from primary producers toward high trophic levels is consistently preserved. Consequently, traditional indicators might provide a misleading assessment of the eutrophication risk. To avoid this, we propose to evaluate eutrophication by using a new index based on plankton trophic fluxes instead of biogeochemical concentrations. A preliminary, model-based, assessment suggests that this approach might give a substantially different picture of the eutrophication status of our seas, with potential consequences on marine ecosystem management. Given the difficulties to measure trophic fluxes in the field, the use of numerical simulations is recommended although the uncertainty associated with biogeochemical models inevitably affects the reliability of the index. However, given the effort currently in place to develop refined numerical tools describing the marine environment (Ocean Digital Twins), a reliable, model-based, eutrophication index could be operational in the near future.

Description: Macroalgal (seaweed) genomic resources are generally lacking as compared to other eukaryotic taxa, and this is particularly true in the red algae (Rhodophyta). Understanding red algal genomes is critical to understanding eukaryotic evolution given that red algal genes are spread across eukaryotic lineages from secondary endosymbiosis and red algae diverged early in the Archaeplastids. The Gracilariales is a highly diverse and widely distributed order including species that can serve as ecosystem engineers in intertidal habitats and several notorious introduced species. The genus Gracilaria is cultivated worldwide, in part for its production of agar and other bioactive compounds with downstream pharmaceutical and industrial applications. This genus is also emerging as a model for algal evolutionary ecology. Here, we report new whole genome assemblies for two species (G. chilensis and G. gracilis), a draft genome assembly of G. caudata, and genome annotation of the previously published G. vermiculophylla genome. To facilitate accessibility and comparative analysis, we integrated these data in a newly created web-based portal dedicated to red algal genomics (https://rhodoexplorer.sb-roscoff.fr). These genomes will provide a resource for understanding algal biology and, more broadly, eukaryotic evolution.

Description: Biological invasions are a global challenge that has received insufficient attention. Recently available cost syntheses have provided policy- and decision makers with reliable and up-to-date information on the economic impacts of biological invasions, aiming to motivate effective management. The resultant InvaCost database is now publicly and freely accessible and enables rapid extraction of monetary cost information. This has facilitated knowledge sharing, developed a more integrated and multidisciplinary network of researchers, and forged multidisciplinary collaborations among diverse organizations and stakeholders. Over 50 scientific publications so far have used the database and have provided detailed assessments of invasion costs across geographic, taxonomic, and spatiotemporal scales. These studies have provided important information that can guide future policy and legislative decisions on the management of biological invasions while simultaneously attracting public and media attention. We provide an overview of the improved availability, reliability, standardization, and defragmentation of monetary costs; discuss how this has enhanced invasion science as a discipline; and outline directions for future development.

Description: The development and physiology of herring larvae were monitored for individuals reared in control and combined warming-acidification crossed with different food quality treatments. The experiment revealed that warming and acidification triggers a stress response at the molecular level and decrease herring larvae size-at-stage. Global change puts coastal systems under pressure, affecting the ecology and physiology of marine organisms. In particular, fish larvae are sensitive to environmental conditions, and their fitness is an important determinant of fish stock recruitment and fluctuations. To assess the combined effects of warming, acidification and change in food quality, herring larvae were reared in a control scenario (11 & DEG;C*pH 8.0) and a scenario predicted for 2100 (14 & DEG;C*pH 7.6) crossed with two feeding treatments (enriched in phosphorus and docosahexaenoic acid or not). The experiment lasted from hatching to the beginning of the post-flexion stage (i.e. all fins present) corresponding to 47 days post-hatch (dph) at 14 & DEG;C and 60 dph at 11 & DEG;C. Length and stage development were monitored throughout the experiment and the expression of genes involved in growth, metabolic pathways and stress responses were analysed for stage 3 larvae (flexion of the notochord). Although the growth rate was unaffected by acidification and temperature changes, the development was accelerated in the 2100 scenario, where larvae reached the last developmental stage at a smaller size (-8%). We observed no mortality related to treatments and no effect of food quality on the development of herring larvae. However, gene expression analyses revealed that heat shock transcripts expression was higher in the warmer and more acidic treatment. Our findings suggest that the predicted warming and acidification environment are stressful for herring larvae, inducing a decrease in size-at-stage at a precise period of ontogeny. This could either negatively affect survival and recruitment via the extension of the predation window or positively increase the survival by reducing the larval stage duration.

Description: Machine learning covers a large set of algorithms that can be trained to identify patterns in data. Thanks to the increase in the amount of data and computing power available, it has become pervasive across scientific disciplines. We first highlight why machine learning is needed in marine ecology. Then we provide a quick primer on machine learning techniques and vocabulary. We built a database of & SIM;1000 publications that implement such techniques to analyse marine ecology data. For various data types (images, optical spectra, acoustics, omics, geolocations, biogeochemical profiles, and satellite imagery), we present a historical perspective on applications that proved influential, can serve as templates for new work, or represent the diversity of approaches. Then, we illustrate how machine learning can be used to better understand ecological systems, by combining various sources of marine data. Through this coverage of the literature, we demonstrate an increase in the proportion of marine ecology studies that use machine learning, the pervasiveness of images as a data source, the dominance of machine learning for classification-type problems, and a shift towards deep learning for all data types. This overview is meant to guide researchers who wish to apply machine learning methods to their marine datasets.

Description: The Observing Air–Sea Interactions Strategy (OASIS) is a new United Nations Decade of Ocean Science for Sustainable Development programme working to develop a practical, integrated approach for observing air–sea interactions globally for improved Earth system (including ecosystem) forecasts, CO2 uptake assessments called for by the Paris Agreement, and invaluable surface ocean information for decision makers. Our “Theory of Change” relies upon leveraged multi-disciplinary activities, partnerships, and capacity strengthening. Recommendations from 〉40 OceanObs’19 community papers and a series of workshops have been consolidated into three interlinked Grand Ideas for creating #1: a globally distributed network of mobile air–sea observing platforms built around an expanded array of long-term time-series stations; #2: a satellite network, with high spatial and temporal resolution, optimized for measuring air–sea fluxes; and #3: improved representation of air–sea coupling in a hierarchy of Earth system models. OASIS activities are organized across five Theme Teams: (1) Observing Network Design & Model Improvement; (2) Partnership & Capacity Strengthening; (3) UN Decade OASIS Actions; (4) Best Practices & Interoperability Experiments; and (5) Findable–Accessible–Interoperable–Reusable (FAIR) models, data, and OASIS products. Stakeholders, including researchers, are actively recruited to participate in Theme Teams to help promote a predicted, safe, clean, healthy, resilient, and productive ocean.

Description: We analyse the consistency of the delay time data in the most recent version of the ISC-EHB bulletin published by the International Seismological Centre covering the years 1964–2018. Considering that the delays are influenced by the lateral heterogeneity in the Earth’s mantle, we construct a tomographic matrix. We use singular value decomposition of the tomographic matrix for 19 707 dense clusters of earthquakes to compute objective estimates of the standard error from data that project into the null space and should be zero if there were no errors. Using a robust initial estimate of the standard deviation of the clustered delay times, we remove a small fraction of outliers before calculating the ultimate errors. We found that the errors depend on the type of body wave, depth of the earthquake (crust or mantle) and the number of decimals with which the arrival time was reported. Using these parameters, we distinguish 45 different classes of delay times for 11 different types of body waves. The errors of each class so divided generally follow a distribution that is approximately normal with a mean that ranges from 0.32 s for PKPbc waves from mantle earthquakes, to 2.82 s for S waves from shallow earthquakes bottoming in the upper mantle. The widths of the distributions of the errors themselves are small enough to serve in formal statistical quantification of the quality of fit in tomographic experiments.

Description: Strong anisotropy of seismic velocity in the Earth’scrust poses serious challenges for seismic imaging. Where in situ seismic properties are not available the anisotropy can be determined from velocity analysis of surface and borehole seismic profiles. This is well established for dense, long-offset reflection seismic data. However, it is unknown how applicable this approach is for sparse seismic reflection data with low fold and short offsets in anisotropic metamorphic rocks. Here we show that anisotropy parameters can be determined from a sparse 3D data set at the COSC-1 borehole site in the Swedish Caledonides and that the results agree well with the seismic anisotropy parameters determined from seismic laboratory measurements on core samples. Applying these anisotropy parameters during 3D seismic imaging improves the seismic image of the high amplitude reflections especially in the vicinity of the lower part of the borehole. Strong reflections in the resulting seismic data show good correlation with the borehole-derived lithology. Our results aid the interpretation and extrapolation of the seismic stratigraphy of the Lower Seve Nappe in Jämtland and other parts in the Caledonides.

Description: The southern boundary of the Cayman Trough in the Caribbean is marked by the Swan Islands transform fault (SITF), which also represents the ocean-continent transition of the Honduras continental margin. This is one of the few places globally where a transform continental margin is currently active. The CAYSEIS experiment acquired an ∼165-km-long seismic refraction and gravity profile (P01) running across this transform margin, and along the ridge-axis of the Mid-Cayman Spreading Centre (MCSC) to the north. This profile reveals not only the crustal structure of an actively evolving transform continental margin, that juxtaposes Mesozoic-age continental crust to the south against zero-age ultraslow spread oceanic crust to the north, but also the nature of the crust and uppermost mantle beneath the ridge-transform intersection (RTI). The traveltimes of arrivals recorded by ocean-bottom seismographs (OBSs) deployed along-profile have been inverse and forward modelled, in combination with gravity modelling, to reveal an ∼25-km-thick continental crust that has been continuously thinned over a distance of ∼65 km to ∼10 km adjacent to the SITF, where it is juxtaposed against ∼3-4-km-thick oceanic crust. This thinning is primarily accommodated within the lower crust. Since Moho reflections are only sparsely observed, and, even then, only by a few OBSs located on the continental margin, the 7.5 km s-1 velocity contour is used as a proxy to locate the crust-mantle boundary along-profile. Along the MCSC, the crust-mantle boundary appears to be a transition zone, at least at the seismic wavelengths used for CAYSEIS data acquisition. Although the traveltime inversion only directly constrains the upper crust at the SITF, gravity modelling suggests that it is underlain by a higher density (〉3000 kg m-3) region spanning the width (∼15 km) of its bathymetric expression, that may reflect a broad region of metasomatism, mantle hydration or melt-depleted lithospheric mantle. At the MCSC ridge-axis to the north, the oceanic crust appears to be forming in zones, where each zone is defined by the volume of its magma supply. The ridge tip adjacent to the SITF is currently in a magma rich phase of accretion. However, there is no evidence for melt leakage into the transform zone. The width and crustal structure of the SITF suggests its motion is currently predominantly orthogonal to spreading. Comparison to CAYSEIS Profile P04, located to the west and running across-margin and through 10 Ma MCSC oceanic crust, suggests that, at about this time, motion along the SITF had a left-lateral transtensional component, that accounts for its apparently broad seabed appearance westwards.

Description: Volcano seismology is an essential tool for monitoring volcanic processes in the advent and during eruptions. A variety of seismic signals can be recorded at volcanoes, of which some are thought to be related to the migration of fluids which is of primary importance for the anticipation of imminent eruptions. We investigate the volcanic crises at Villarrica volcano in 2015 and report on a newly discovered very-long-period (VLP) signal that accompanies phases of periodic long period (LP) signal burst. Despite their low amplitude emergent character, we can locate the source region of the 1 Hz LP signals to the close vicinity of the volcano using a network-based correlation method. The source of the VLP signal with a period of about 30–100 s appears to locate in the vicinity of two stations a few kilometres from the summit. Both stations record very similar VLP waveforms that are correlated with the envelope of the LP bursts. A shallow magma reservoir was inferred by Contreras from surface deformation as the source of inflation following the eruption in 2015. Cyclic volume changes of 6 m3 in this reservoir at 3 km depth can explain the observed amplitudes of the vertical VLP signal. We propose that the LP signal is generated by the migration of gas or gas-rich magma that is periodically released from the inflating reservoir through a non-linear valve structure which modulates the flux, and thereby causes bursts of flow-related LP signals and pressure changes observed as VLP deformation. Our model predicts that the correlated occurrence of LP bursts and VLP surface motion depends on the intensity of the fluid flux. A weaker flux of fluids may not exceed the opening pressure of valve structure, and higher rates might maintain pressure above the closing pressure. In both cases, the VLP signal vanishes. Our observation provides constrains for models of fluid transport inside volcanoes. At Villarrica the VLP signal, and its relation to the LP activity, reveal additional information about fluxes in the magmatic reservoir that might aide forecasting of volcanic activity.

Description: The female flowers (“cones”) of the hop plant (Humulus L.) produce compounds that contribute to the flavor and other properties of beer. Hop leaves and cones produce many of the same compounds, which also confer agronomic traits such as insect and disease resistance. Targeted and untargeted ultraperformance liquid chromatography-quadrupole time-of-flight–mass spectrometry with Waters MSE technology (UPLC-QTof-MSE) metabolomics were used to compare leaf phytochemical compositions of greenhouse-grown southwestern American wild Humulus neomexicanus (A. Nelson and Cockerell) Rydb. against a group of commercial hop cultivars consisting of both pure European Humulus lupulus L. and European–North American hybrids. Principal component analysis showed a clear distinction in chemical profiles between the two groups. H. neomexicanus leaves had a significantly higher content of total α acids (p = 4.4 × 10–9), total bitter acids (p = 2.6 × 10–6), cohumulone (p = 1.0 × 10–13), humulone + adhumulone (p = 9.1 × 10–4), and the prenylflavonoids xanthohumol (p = 0.013) and desmethylxanthohumol (p = 0.029) as well as significantly higher densities of glandular trichomes (p = 1.3 × 10–6), the biosynthetic site of those compounds. Most flavonol glycosides measured were also significantly more abundant in H. neomexicanus (p = 1.5 × 10–22 to 0.0027), whereas phenolic acids were consistently, but generally nonsignificantly (p 〉 0.05), more abundant in the cultivars. The higher bitter acid, prenylflavonoid, and flavonol glycoside content of H. neomexicanus leaves may help to confer more favorable insect and disease-resistance properties.

Description: Iron is an indispensable metabolic cofactor in both pro- and eukaryotes, which engenders a natural competition for the metal between bacterial pathogens and their human or animal hosts. Bacteria secrete siderophores that extract Fe3+ from tissues, fluids, cells, and proteins; the ligand gated porins of the Gram-negative bacterial outer membrane actively acquire the resulting ferric siderophores, as well as other iron-containing molecules like heme. Conversely, eukaryotic hosts combat bacterial iron scavenging by sequestering Fe3+ in binding proteins and ferritin. The variety of iron uptake systems in Gram-negative bacterial pathogens illustrates a range of chemical and biochemical mechanisms that facilitate microbial pathogenesis. This document attempts to summarize and understand these processes, to guide discovery of immunological or chemical interventions that may thwart infectious disease.

Description: While thousands of environmental metagenomes have been mined for the presence of novel biosynthetic gene clusters, such computational predictions do not provide evidence of their in vivo biosynthetic functionality. Using fluorescent in situ enzyme assay targeting carrier proteins common to polyketide (PKS) and nonribosomal peptide synthetases (NRPS), we applied fluorescence-activated cell sorting to tunicate microbiome to enrich for microbes with active secondary metabolic capabilities. Single-cell genomics uncovered the genetic basis for a wide biosynthetic diversity in the enzyme-active cells and revealed a member of marine Oceanospirillales harboring a novel NRPS gene cluster with high similarity to phylogenetically distant marine and terrestrial bacteria. Interestingly, this synthase belongs to a larger class of siderophore biosynthetic gene clusters commonly associated with pestilence and disease. This demonstrates activity-guided single-cell genomics as a tool to guide novel biosynthetic discovery.

Description: Full waveform inversion (FWI) is a data-fitting technique capable of generating high-resolution velocity models with a resolution down to half the seismic wavelength. FWI is applied typically to densely sampled seismic data. In this study, we applied FWI to 3D wide-angle seismic data acquired using sparsely spaced ocean bottom seismometers (OBSs) from the Deep Galicia Margin west of Iberia. Our dataset samples the S-reflector, a low-angle detachment present in this area. Here we highlight differences between 2D, 2.5D and 3D-FWI performances using a real sparsely spaced dataset. We performed 3D FWI in the time domain and compared the results with 2D and 2.5D FWI results from a profile through the 3D model. When overlaid on multichannel seismic images, the 3D FWI results constrain better the complex faulting within the pre- and syn-rift sediments and crystalline crust compared to the 2D result. Furthermore, we estimate variable serpentinisation of the upper mantle below the S-reflector along the profile using 3D FWI, reaching a maximum of 45 per cent. Differences in the data residuals of the 2D, 2.5D and 3D inversions suggest that 2D inversion can be prone to overfitting when using a sparse dataset. To validate our results, we performed tests to recover the anomalies introduced by the inversions in the final models using synthetic datasets. Based on our comparison of the velocity models, we conclude that the use of 3D data can partially mitigate the problem of receiver sparsity in FWI.

Description: Diatoms often dominate phytoplankton in temperate, polar and upwelling regions. Decreases in silicate availability or silicon to nitrogen (Si:N) ratios may induce silicon limitation in diatoms and lower their proportion within phytoplankton communities. The effects of such changes on the nutritional quality of phytoplankton are not well understood. To examine how changing Si:N ratios affect plankton nutritional value, we applied a range of Si:N ratios on a natural plankton community and manipulated grazing pressure to assess top-down effects of copepod selective grazing. Diatom proportion in phytoplankton increased with increasing Si:N ratios and so did phytoplankton nutritional quality in terms of major fatty acid concentrations, such as polyunsaturated fatty acids, docosahexaenoic (DHA) and eicosapentaenoic (EPA) acids. However, stoichiometric quality (carbon to nitrogen and carbon to phosphorus ratios), DHA:EPA and omega 3:6 (omega 3:omega 6) ratios declined with increasing Si:N ratios, suggesting that proportions between essential compounds in copepod diet may be more favorable in lowered Si:N ratios. Copepods had a negative effect on DHA contents, DHA:EPA and omega 3:omega 6 ratios, indicating possible selective grazing on more nutritious plankton. Our findings show that declining silicate concentrations can affect stoichiometric and biochemical quality of phytoplankton, which copepods can also moderate by selective grazing.

Description: Temperature and dehydration stress are two major co-occurring environmental stressors threatening the physiology, biochemistry, and ecology of insects. As such, understanding adaptive responses to desiccation stress is critical for predicting climate change impacts, particularly its influence on insect invasions. Here, we assessed water balance and desiccation resistance of the invasive Tuta absoluta (Meyrick, 1917) (Lepidoptera: Gelechiidae), and infer how eco-physiology shapes its niche. We measured basal body water and lipid content, water loss rates (WLRs), and desiccation resistance in larvae (second to fourth instars) and adults. Body -water, -lipid, and WLRs significantly varied across life stages. Second instars recorded the lowest while fourth instars exhibited the highest body water and lipid content. Adult body water and lipid content were higher than second and third instars and lower than fourth instars while proportion of body water and lipid contents were highest in adults and second larval instars respectively. Water loss rates were significantly highest in fourth-instar larvae compared to other life stages, but differences among stages were less apparent at longer exposure durations (48 h). Desiccation resistance assays showed that second instars had greatest mortality while fourth-instar larvae and adults were the most desiccation tolerant. Our results show that T. absoluta fourth-instar larvae and adults are the most resilient developmental stages and potentially contribute most to the invasion success of the pest in arid environments. Incorporation of these species-specific eco-physiological traits in predictive models can help refine invasive species potential spread under changing climates.

Description: The fishery for Northern Atlantic cod (Gadus morhua) off Newfoundland and Labrador, Eastern Canada, presents the most spectacular case of an exploited stock crashed in a few decades by an industrial bottom trawl fishery under a seemingly sophisticated management regime after half a millennium of sustainable fishing. The fishery, which had generated annual catches of 100000 to 200000 tonnes from the beginning of the 16th century to the 1950s, peaked in 1968 at 810000 tonnes, followed by a devastating collapse and closure 24 years later. Since then, stock recovery may have been hindered by premature openings, with vessels targeting the remains of the cod population. Previous research paid little attention towards using multicentury time series to inform sustainable catches and recovery plans. Here, we show that a simple stock assessment model can be used to model the cod population trajectory for the entire period from 1508 to 2019 for which catch estimates are available. The model suggests that if fishing effort and mortality had been stabilized in the 1980s, precautionary annual yields of about 200000 tonnes could have been sustained. Our analysis demonstrates the value of incorporating prior knowledge to counteract shifting baseline effects on reference points and contemporary perceptions of historical stock status.

Description: Probing seismic anisotropy of the lithosphere provides valuable clues on the fabric of rocks. We present a 3-D probabilistic model of shear wave velocity and radial anisotropy of the crust and uppermost mantle of Europe, focusing on the mountain belts of the Alps and Apennines. The model is built from Love and Rayleigh dispersion curves in the period range 5–149 s. Data are extracted from seismic ambient noise recorded at 1521 broad-band stations, including the AlpArray network. The dispersion curves are first combined in a linearized least squares inversion to obtain 2-D maps of group velocity at each period. Love and Rayleigh maps are then jointly inverted at depth for shear wave velocity and radial anisotropy using a Bayesian Monte Carlo scheme that accounts for the trade-off between radial anisotropy and horizontal layering. The isotropic part of our model is consistent with previous studies. However, our anisotropy maps differ from previous large scale studies that suggested the presence of significant radial anisotropy everywhere in the European crust and shallow upper mantle. We observe instead that radial anisotropy is mostly localized beneath the Apennines while most of the remaining European crust and shallow upper mantle is isotropic. We attribute this difference to trade-offs between radial anisotropy and thin (hectometric) layering in previous studies based on least-squares inversions and long period data (〉30 s). In contrast, our approach involves a massive data set of short period measurements and a Bayesian inversion that accounts for thin layering. The positive radial anisotropy (VSH 〉 VSV) observed in the lower crust of the Apennines cannot result from thin layering. We rather attribute it to ductile horizontal flow in response to the recent and present-day extension in the region.

Description: The freshwater sponge Ephydatia muelleri and its Chlorella-like algal partner is an emerging model for studying animal: algal endosymbiosis. The sponge host is a tractable laboratory organism, and the symbiotic algae are easily cultured. We took advantage of these traits to interrogate questions about mechanisms that govern the establishment of durable intracellular partnerships between hosts and symbionts in facultative symbioses. We modified a classical experimental approach to discern the phagocytotic mechanisms that might be co-opted to permit persistent infections, and identified genes differentially expressed in sponges early in the establishment of endosymbiosis. We exposed algal-free E. muelleri to live native algal symbionts and potential food items (bacteria and native heat-killed algae), and performed RNA-Seq to compare patterns of gene expression among treatments. We found a relatively small but interesting suite of genes that are differentially expressed in the host exposed to live algal symbionts, and a larger number of genes triggered by host exposure to heat-killed algae. The upregulated genes in sponges exposed to live algal symbionts were mostly involved in endocytosis, ion transport, metabolic processes, vesicle-mediated transport, and oxidation–reduction. One of the host genes, an ATP-Binding Cassette transporter that is downregulated in response to live algal symbionts, was further evaluated for its possible role in the establishment of the symbiosis. We discuss the gene expression profiles associated with host responses to living algal cells in the context of conditions necessary for long-term residency within host cells by phototrophic symbionts as well as the genetic responses to sponge phagocytosis and immune-driven pathways.

Description: Marine sponges play a major ecological role in recycling resources on coral reef ecosystems. The cycling of resources may largely depend on the stability of the host-microbiome interactions and their susceptibility to altered environmental conditions. Given the current coral to algal phase shift on coral reefs, we investigated whether the sponge-associated bacterial communities of four sponge species, with either high or low microbial abundances (HMA and LMA), remain stable at two reefs sites with different coral to algae cover ratios. Additionally, we assessed the bacterial community composition of two of these sponge species before and after a reciprocal transplantation experiment between the sites. An overall stable bacterial community composition was maintained across the two sites in all sponge species, with a high degree of host-specificity. Furthermore, the core bacterial communities of the sponges remained stable also after a 21-day transplantation period, although a minor shift was observed in less abundant taxa (〈 1%). Our findings support the conclusion that host identity and HMA-LMA status are stronger traits in shaping bacterial community composition than habitat. Nevertheless, long-term microbial monitoring of sponges along with benthic biomass and water quality assessments are needed for identifying ecosystem tolerance ranges and tipping points in ongoing coral reef phase shifts.

Description: Seismic full waveform inversion (FWI) is a powerful method for estimating quantitative subsurface physical parameters from seismic data. As the full waveform inversion is a non-linear problem, the linearized approach updates model iteratively from an initial model, which can get trapped in local minima. In the presence of a high velocity contrast, such as at Moho, the reflection coefficient and recorded waveforms from wide-aperture seismic acquisition are extremely non-linear around critical angles. The problem at the Moho is further complicated by the interference of lower crustal (Pg) and upper mantle (Pn) turning ray arrivals with the critically reflected Moho arrivals (PmP). In order to determine velocity structure near Moho, a non-linear method should be used. We propose to solve this strong non-linear FWI problem at Moho using a trans-dimensional Markov chain Monte Carlo (MCMC) method, where the earth model between lower crust and upper mantle is idealy parameterized with a 1-D assumption using a variable number of velocity interfaces. Different from common MCMC methods that require determining the number of unknown as a fixed prior before inversion, trans-dimensional MCMC allows the flexibility for an automatic estimation of both the model complexity (e.g. the number of velocity interfaces) and the velocity-depth structure from the data. We first test the algorithm on synthetic data using four representative Moho models and then apply to an ocean bottom seismometer (OBS) data from the Mid-Atlantic Ocean. A 2-D finite-difference solution of an acoustic wave equation is used for data simulation at each iteration of MCMC search, for taking into account the lateral heterogeneities in the upper crust, which is constrained from travel time tomography and is kept unchanged during inversion; the 1-D model parameterization near Moho enables an efficient search of the trans-dimensional model space. Inversion results indicate that, with very little prior and the wide-aperture seismograms, the trans-dimensional FWI method is able to infer the posterior distribution of both the number of velocity interfaces and the velocity-depth model for a strong nonlinear problem, making the inversion a complete data-driven process. The distribution of interface matches the velocity discontinuities. We find that the Moho in the study area is a transition zone of 0.7 km, or a sharp boundary with velocities from around 7 km/s in the lower crust to 8 km/s of the upper mantle; both provide nearly identical waveform match for the field data. The ambiguity comes from the resolution limit of the band-limited seismic data and limited offset range for PmP arrivals.

Description: The maximum sustainable yield (MSY) concept is widely considered to be outdated and misleading. In response, fisheries scientists have developed models that often diverge radically from the first operational version of the concept. We show that the original MSY concept was deeply rooted in ecology and that going back to that version would be beneficial for fisheries, not least because the various substitutes have not served us well.

Description: Food webs are central entities mediating processes and external pressures in marine ecosystems. They are essential to understand and predict ecosystem dynamics and provision of ecosystem services. Paradoxically, utilization of food web knowledge in marine environmental conservation and resource management is limited. To better understand the use of knowledge and barriers to incorporation in management, we assess its application related to the management of eutrophication, chemical contamination, fish stocks, and non-indigenous species. We focus on the Baltic, a severely impacted, but also intensely studied and actively managed semi-enclosed sea. Our assessment shows food web processes playing a central role in all four areas, but application varies strongly, from formalized integration in management decisions, to support in selecting indicators and setting threshold values, to informal knowledge explaining ecosystem dynamics and management performance. Barriers for integration are complexity of involved ecological processes and that management frameworks are not designed to handle such information. We provide a categorization of the multi-faceted uses of food web knowledge and benefits of future incorporation in management, especially moving towards ecosystem-based approaches as guiding principle in present marine policies and directives. We close with perspectives on research needs to support this move considering global and regional change.

Description: To constrain seismic anisotropy under and around the Alps in Europe, we study SKS shear wave splitting from the region densely covered by the AlpArray seismic network. We apply a technique based on measuring the splitting intensity, constraining well both the fast orientation and the splitting delay. Four years of teleseismic earthquake data were processed, from 723 temporary and permanent broad-band stations of the AlpArray deployment including ocean-bottom seismometers, providing a spatial coverage that is unprecedented. The technique is applied automatically (without human intervention), and it thus provides a reproducible image of anisotropic structure in and around the Alpine region. As in earlier studies, we observe a coherent rotation of fast axes in the western part of the Alpine chain, and a region of homogeneous fast orientation in the Central Alps. The spatial variation of splitting delay times is particularly interesting though. On one hand, there is a clear positive correlation with Alpine topography, suggesting that part of the seismic anisotropy (deformation) is caused by the Alpine orogeny. On the other hand, anisotropic strength around the mountain chain shows a distinct contrast between the Western and Eastern Alps. This difference is best explained by the more active mantle flow around the Western Alps. The new observational constraints, especially the splitting delay, provide new information on Alpine geodynamics.

Description: In this article, the concepts and background of regional climate modeling of the future Baltic Sea are summarized and state-of-the-art projections, climate change impact studies, and challenges are discussed. The focus is on projected oceanographic changes in future climate. However, as these changes may have a significant impact on biogeochemical cycling, nutrient load scenario simulations in future climates are briefly discussed as well. The Baltic Sea is special compared to other coastal seas as it is a tideless, semi-enclosed sea with large freshwater and nutrient supply from a partly heavily populated catchment area and a long response time of about 30 years, and as it is, in the early 21st century, warming faster than any other coastal sea in the world. Hence, policymakers request the development of nutrient load abatement strategies in future climate. For this purpose, large ensembles of coupled climate–environmental scenario simulations based upon high-resolution circulation models were developed to estimate changes in water temperature, salinity, sea-ice cover, sea level, oxygen, nutrient, and phytoplankton concentrations, and water transparency, together with uncertainty ranges. Uncertainties in scenario simulations of the Baltic Sea are considerable. Sources of uncertainties are global and regional climate model biases, natural variability, and unknown greenhouse gas emission and nutrient load scenarios. Unknown early 21st-century and future bioavailable nutrient loads from land and atmosphere and the experimental setup of the dynamical downscaling technique are perhaps the largest sources of uncertainties for marine biogeochemistry projections. The high uncertainties might potentially be reducible through investments in new multi-model ensemble simulations that are built on better experimental setups, improved models, and more plausible nutrient loads. The development of community models for the Baltic Sea region with improved performance and common coordinated experiments of scenario simulations is recommended.

Description: Marine macroalgae, seaweeds, are exceptionally prolific producers of halogenated natural products. Biosynthesis of halogenated molecules in seaweeds is inextricably linked to reactive oxygen species (ROS) signaling as hydrogen peroxide serves as a substrate for haloperoxidase enzymes that participate in the construction these halogenated molecules. Here, using red macroalga Asparagopsis taxiformis, a prolific producer of the ozone depleting molecule bromoform, we provide the discovery and biochemical characterization of a ROS-producing NAD(P)H oxidase from seaweeds. This discovery was enabled by our sequencing of Asparagopsis genomes, in which we find the gene encoding the ROS-producing enzyme to be clustered with genes encoding bromoform-producing haloperoxidases. Biochemical reconstitution of haloperoxidase activities establishes that fatty acid biosynthesis can provide viable hydrocarbon substrates for bromoform production. The ROS production haloperoxidase enzymology that we describe here advances seaweed biology and biochemistry by providing the molecular basis for decades worth of physiological observations in ROS and halogenated natural product biosyntheses.

Description: Fusarium oxysporum is a destructive root-infecting plant pathogen that causes significant yield losses in many economically important crop species. Hence, a deeper understanding of pathogen infection strategies is needed. With liquid chromatography-tandem mass spectrometry and gas chromatography-time of flight mass spectrometry platforms, we analyzed the metabolic changes in a time-course experiment with Arabidopsis accessions either resistant (Col-0) or susceptible (Ler-0) to isolates of Fusarium oxysporum forma specialis matthioli infection. We showed a concurrent effect of Fusarium-derived polyols and the mycotoxin beauvericin in the suppression of the immune response of susceptible hosts. A significant increase in oxidized glutathione in the resistant host was probably associated with effective reactive oxygen species-mediated resistance responses. Through a combination of targeted and untargeted metabolomics, we demonstrated the concurrent action of several Arabidopsis defense systems as well as the concurrent action of several virulence systems in the fungal attack of susceptible Arabidopsis.

Description: Intense bottom-ice algal blooms, often dominated by diatoms, are an important source of food for grazers, organic matter for export during sea ice melt, and dissolved organic carbon. Sea-ice diatoms have a number of adaptations, including accumulation of compatible solutes, that allows them to inhabit this highly variable environment characterized by extremes in temperature, salinity, and light. In addition to protecting them from extreme conditions, these compounds present a labile, nutrient-rich source of organic matter, and include precursors to climate active compounds (e.g., dimethyl sulfide [DMS]), which are likely regulated with environmental change. Here, intracellular concentrations of 45 metabolites were quantified in three sea-ice diatom species and were compared to two temperate diatom species, with a focus on compatible solutes and free amino acid pools. There was a large diversity of metabolite concentrations between diatoms with no clear pattern identifiable for sea-ice species. Concentrations of some compatible solutes (isethionic acid, homarine) approached 1 M in the sea-ice diatoms, Fragilariopsis cylindrus and Navicula cf. perminuta, but not in the larger sea-ice diatom, Nitzschia lecointei or in the temperate diatom species. The differential use of compatible solutes in sea-ice diatoms suggests different adaptive strategies and highlights which small organic compounds may be important in polar biogeochemical cycles.

Description: Deep-towed geophysical surveys require precise knowledge of navigational parameters such as instrument position and orientation because navigational uncertainties reflect in the data and therefore in the inferred geophysical properties of the subseafloor. We address this issue for the case of electrical conductivity inferred from controlled source electromagnetic data. We show that the data error is laterally variable due to irregular motion during deep towing, but also due to lateral variations in conductivity, including those resulting from topography. To address this variability and quantify the data error prior to inversion, we propose a 2-D perturbation study. Our workflow enables stable and geologically reliable results for multicomponent and multifrequency inversions. An error estimation workflow is presented, which comprises the assessment of navigational uncertainties, perturbation of navigational parameters, and forward modelling of electric field amplitudes for a homogeneous and then a heterogeneous subseafloor conductivity model. Some navigational uncertainties are estimated from variations of direct measurements. Other navigational parameters required for inversion are derived from the measured quantities and their error is calculated by means of error propagation. Some navigational parameters show direct correlation with the measured electric fields. For example, the antenna dip correlates with the vertical electric field and the depth correlates with the horizontal electric field. For the perturbation study each standard deviation is added to the navigational parameters. Forward models are run for each perturbation. Amplitude deviations are summed in quadrature with the stacking error for a total, laterally varying, data error. The error estimation is repeated for a heterogeneous subseafloor model due to the large conductivity range (several orders of magnitude), which affects the forward model. The approach enables us to utilize data from several components (multiple electric fields, frequencies and receivers) in the inversion to constrain the final model and reduce ambiguity. The final model is geologically reasonable, in this case enabling the identification of conductive metal sulphide deposits on the seafloor.

Description: Many coastal oceans experience not only increased loads of nutrients but also changes in the stoichiometry of nutrient supply. Excess supply of nitrogen and stable or decreased supply of silicon lower silicon to nitrogen (Si:N) ratios, which may decrease diatom proportion in phytoplankton. To examine how Si:N ratios affect plankton community composition and food web structure, we performed a mesocosm experiment where we manipulated Si:N ratios and copepod abundance in a Baltic Sea plankton community. In high Si:N treatments, diatoms dominated. Some of them were likely spared from grazing unexpectedly resulting in higher diatom biomass under high copepod grazing. With declining Si:N ratios, dinoflagellates became more abundant under low and picoplankton under high copepod grazing. This altered plankton food web structure: under high Si:N ratios, edible diatoms were directly accessible food for copepods, while under low Si:N ratios, microzooplankton and phago-mixotrophs (mixoplankton) were a more important food source for mesograzers. The response of copepods to changes in the phytoplankton community was complex and copepod density-dependent. We suggest that declining Si:N ratios favor microzoo- and mixoplankton leading to increased complexity of planktonic food webs. Consequences on higher trophic levels will, however, likely be moderated by edibility, nutritional value or toxicity of dominant phytoplankton species.

Description: Quantification and attribution of the food web changes associated with the invasion of non-indigenous species in the marine realm often remain a challenge. One of the pelagic non-indigenous species of concern in the recent history of aquatic bioinvasions is the predatory cladoceran Cercopagis pengoi, which invaded the Baltic Sea in the early 1990s. While several studies have reported immediate declines in abundances of its potential prey, the long-term effects of C. pengoi on the food webs remain to be examined. Based on the long-term time series (1968–2018) in the Gulf of Riga (Baltic Sea), we found significant declines in abundance of the cladoceran Pleopis spp. and copepod Eurytemora affinis by 90 and 80%, respectively, are associated with the invasion of C. pengoi as well as significant alterations in seasonal abundance patterns of Pleopis spp., E. affinis and cladoceran Bosmina spp. The invasion of the non-indigenous predator has led to the changed prey abundance–temperature relationships. Special caution was taken in data preprocessing, to minimize the likelihood that observed changes in the zooplankton prey could be associated with factors other than the invasion of C. pengoi, such as temperature and storminess.

Description: The copepod Acartia tonsa is a key component of a wide range of marine ecosystems, linking energy transfer from phytoplankton to higher trophic levels, and has a central role in productivity and biogeochemistry. The interaction of end-of-century global warming and ocean acidification scenarios with testing moderate temperature effects on a seminatural copepod community is needed to understand future community functioning. Here, we deployed a mesocosm experimental set-up with a full factorial design using two temperatures (13°C and 19°C) crossed with a pCO2 gradient ranging from ambient (550 μatm) to 3000 μatm. We used the natural bacteria, phyto- and microzooplankton species composition and biomass of the Kiel Bight and tested the response of A. tonsa development, carbon growth, mortality, size and condition. The tested traits were differently affected by the interaction of temperature and acidification. Ocean acidification increased development, carbon growth, size and mortality under the warming scenario of 19°C. At 13°C mortality rates decreased, while carbon growth, size and condition increased with acidification. We conclude from our experimental approach that a single species shows a variety of responses depending on the focal functional trait. Trait-specific mesozooplankton responses need to be further investigated and compared between geographical regions, seasons and taxonomic groups.

Description: Dinoflagellates possess many cellular characteristics with unresolved evolutionary histories. These include nuclei with greatly expanded genomes and chromatin packaged using histone-like proteins and dinoflagellate-viral nucleoproteins instead of histones, highly reduced mitochondrial genomes with extensive RNA editing, a mix of photosynthetic and cryptic secondary plastids, and tertiary plastids. Resolving the evolutionary origin of these traits requires understanding their ancestral states and early intermediates. Several early-branching dinoflagellate lineages are good candidates for such reconstruction, however these cells tend to be delicate and environmentally sparse, complicating such analyses. Here, we employ transcriptome sequencing from manually isolated and microscopically documented cells to resolve the placement of two cells of one such genus, Abedinium, collected by remotely operated vehicle in deep waters off the coast of Monterey Bay, CA. One cell corresponds to the only described species, Abedinium dasypus, whereas the second cell is distinct and formally described as Abedinium folium, sp. nov. Abedinium has classically been assigned to the early-branching dinoflagellate subgroup Noctilucales, which is weakly supported by phylogenetic analyses of small subunit ribosomal RNA, the single characterized gene from any member of the order. However, an analysis based on 221 proteins from the transcriptome places Abedinium as a distinct lineage, separate from and basal to Noctilucales and the rest of the core dinoflagellates. The transcriptome also contains evidence of a cryptic plastid functioning in the biosynthesis of isoprenoids, iron–sulfur clusters, and heme, a mitochondrial genome with all three expected protein-coding genes (cob, cox1, and cox3), and the presence of some but not all dinoflagellate-specific chromatin packaging proteins.

Description: Background: The giant squid (Architeuthis dux; Steenstrup, 1857) is an enigmatic giant mollusc with a circumglobal distribution in the deep ocean, except in the high Arctic and Antarctic waters. The elusiveness of the species makes it difficult to study. Thus, having a genome assembled for this deep-sea-dwelling species will allow several pending evolutionary questions to be unlocked. Findings: We present a draft genome assembly that includes 200 Gb of Illumina reads, 4 Gb of Moleculo synthetic long reads, and 108 Gb of Chicago libraries, with a final size matching the estimated genome size of 2.7 Gb, and a scaffold N50 of 4.8 Mb. We also present an alternative assembly including 27 Gb raw reads generated using the Pacific Biosciences platform. In addition, we sequenced the proteome of the same individual and RNA from 3 different tissue types from 3 other species of squid (Onychoteuthis banksii, Dosidicus gigas, and Sthenoteuthis oualaniensis) to assist genome annotation. We annotated 33,406 protein-coding genes supported by evidence, and the genome completeness estimated by BUSCO reached 92%. Repetitive regions cover 49.17% of the genome. Conclusions: This annotated draft genome of A. dux provides a critical resource to investigate the unique traits of this species, including its gigantism and key adaptations to deep-sea environments.

Description: Parasites are arguably among the strongest drivers of natural selection, constraining hosts to evolve resistance and tolerance mechanisms. Although, the genetic basis of adaptation to parasite infection has been widely studied, little is known about how epigenetic changes contribute to parasite resistance and eventually, adaptation. Here, we investigated the role of host DNA methylation modifications to respond to parasite infections. In a controlled infection experiment, we used the three-spined stickleback fish, a model species for host-parasite studies, and their nematode parasite Camallanus lacustris. We showed that the levels of DNA methylation are higher in infected fish. Results furthermore suggest correlations between DNA methylation and shifts in key fitness and immune traits between infected and control fish, including respiratory burst and functional trans-generational traits such as the concentration of motile sperm. We revealed that genes associated with metabolic, developmental and regulatory processes (cell death and apoptosis) were differentially methylated between infected and control fish. Interestingly, genes such as the neuropeptide FF receptor 2 and the integrin alpha 1 as well as molecular pathways including the Th1 and Th2 cell differentiation were hypermethylated in infected fish, suggesting parasite-mediated repression mechanisms of immune responses. Altogether, we demonstrate that parasite infection contributes to genome-wide DNA methylation modifications. Our study brings novel insights into the evolution of vertebrate immunity and suggests that epigenetic mechanisms are complementary to genetic responses against parasite-mediated selection.

Description: The Law of the Sea as well as regional and national laws and agreements require exploited populations or stocks to be managed so that they can produce maximum sustainable yields. However, exploitation level and stock status are unknown for most stocks because the data required for full stock assessments are missing. This study presents a new method (AMSY) that estimates relative population size when no catch data are available using time-series of catch-per-unit-effort or other relative abundance indices as the main input. AMSY predictions for relative stock size were not significantly different from the “true” values when compared with simulated data. Also, they were not significantly different from relative stock size estimated by data-rich models in 88% of the comparisons within 140 real stocks. Application of AMSY to 38 data-poor stocks showed the suitability of the method and led to the first assessments for 23 species. Given the lack of catch data as input, AMSY estimates of exploitation come with wide margins of uncertainty which may not be suitable for management. However, AMSY seems to be well suited for estimating productivity as well as relative stock size and may, therefore, aid in the management of data-poor stocks.

Description: International Ocean Discovery Program (IODP) Expedition 351 “Izu–Bonin–Mariana (IBM) Arc Origins” drilled Site U1438, situated in the north-western region of the Philippine Sea. Here volcaniclastic sediments and the igneous basement of the proto-IBM volcanic arc were recovered. To gain a better understanding of the magmatic processes and evolution of the proto-IBM arc, we studied melt inclusions hosted in fresh igneous minerals and sampled from 30- to 40-Ma-old deposits, reflecting the maturation of arc volcanism following subduction initiation at 52 Ma. We performed a novel statistical analysis on the major element composition of 237 representative melt inclusions selected from a previously published dataset, covering the full age range between 30 and 40 Ma. In addition, we analysed volatiles (H2O, S, F and Cl) and P2O5 by Secondary Ion Mass Spectrometry (SIMS) for a subset of 47 melt inclusions selected from the dataset. Based on statistical analysis of the major element composition of melt inclusions and by considering their trace and volatile element compositions, we distinguished five main clusters of melt inclusions, which can be further separated into a total of eight subclusters. Among the eight subclusters, we identified three major magma types: (1) enriched medium-K magmas, which form a tholeiitic trend (30–38 Ma); (2) enriched medium-K magmas, which form a calc-alkaline trend (30–39 Ma); and (3) depleted low-K magmas, which form a calc-alkaline trend (35–40 Ma). We demonstrate that (1) the eruption of depleted low-K calc-alkaline magmas occurred prior to 40 Ma and ceased sharply at 35 Ma; (2) the eruption of depleted low-K calc-alkaline magmas, enriched medium-K calc-alkaline magmas and enriched medium-K tholeiitic magmas overlapped between 35 and 38 − 39 Ma; and (3) the eruption of enriched medium-K tholeiitic and enriched medium-K calc-alkaline magmas became predominant thereafter at the proto-IBM arc. Identification of three major magma types are distinct from the previous work, in which enriched medium-K calc-alkaline magmas and depleted low-K calc-alkaline magmas were not identified. This indicates the usefulness of our statistical analysis as a powerful tool to partition a mixture of multivariable geochemical datasets, such as the composition of melt inclusions in this case. Our data suggest that a depleted mantle source had been replaced by an enriched mantle source due to convection beneath the proto-IBM arc from >40 to 35 Ma. Finally, thermodynamic modelling indicates that the overall geochemical variation of melt inclusions assigned to each cluster can be broadly reproduced either by crystallisation differentiation assuming P = 50 MPa (∼2-km deep) and ∼2 wt % H2O (almost saturated H2O content at 50 MPa) or P = 300 MPa (∼15-km deep) and ∼6 wt % H2O (almost saturated H2O content at 300 MPa). Assuming oxygen fugacity (fO2) of log fO2 equal to + 1 relative to nickel-nickel oxide (NNO) buffer best reproduces the overall geochemical variation of melt inclusions, but assuming a more oxidising conditions (log fO2 = +1 to + 2 NNO) likely reproduces the geochemical variation of enriched medium-K and calc-alkaline melt inclusions (30 − 39 Ma).

Description: The present study quantified prey preferences by adult males and females of the mysid Mesopodopsis wooldridgei fed the calanoid copepods Pseudodiaptomus hessei and Paracartia longipatella at varying proportions. Both sexes of M. wooldridgei showed a lack of prey switching and a strong preference for the smaller, less active P. longipatella irrespective of density. Given a lack of low-density prey refuge, this finding may have important implications for the distribution of P. longipatella in estuaries along the eastern seaboard of South Africa. Results of the present study contribute to a growing body of literature that suggests that selective predation may play an important role in structuring plankton prey populations in shallow water ecosystems.

Description: The two toothed jaws of cichlid fishes provide textbook examples of convergent evolution. Tooth phenotypes such as enlarged molar-like teeth used to process hard-shelled molluscs have evolved numerous times independently during cichlid diversification. While the ecological benefit of molar-like teeth to crush prey is known, it is unclear whether the same molecular mechanisms underlie these convergent traits. To identify genes involved in the evolution and development of enlarged cichlid teeth, we performed RNA-seq on the serially homologous toothed oral and pharyngeal jaws as well as the fourth toothless gill arch of Astatoreochromis alluaudi. We identified 27 genes that are highly upregulated on both tooth-bearing jaws compared to the toothless gill arch. Most of these genes have never been reported to play a role in tooth formation. Two of these genes (unk, rpfA) are not found in other vertebrate genomes but are present in all cichlid genomes. They also cluster genomically with two other highly expressed tooth genes (odam, scpp5) that exhibit conserved expression during vertebrate odontogenesis. Unk and rpfA were confirmed via in situ hybridization to be expressed in developing teeth of Astatotilapia burtoni. We then examined expression of the cluster's four genes in six evolutionarily independent and phylogenetically disparate cichlid species pairs each with a large- and a small-toothed species. Odam and unk commonly and scpp5 and rpfA always showed higher expression in larger-toothed cichlid jaws. Convergent trophic adaptations across cichlid diversity are associated with the repeated developmental deployment of this genomic cluster containing conserved and novel cichlid-specific genes.

Description: Polynoidae contains ~900 species within 18 subfamilies, some of them restricted to the deep sea. Macellicephalinae is the most diverse among these deep-sea subfamilies. In the abyssal Equatorial Pacific Ocean, the biodiversity of benthic communities is at stake in the Clarion-Clipperton Fracture Zone (CCFZ) owing to increased industrial interest in polymetallic nodules. The records of polychaetes in this region are scarce. Data gathered during the JPI Oceans cruise SO239 made a significant contribution to fill this gap, with five different localities sampled between 4000 and 5000 m depth. Benthic samples collected using an epibenthic sledge or a remotely operated vehicle resulted in a large collection of polynoids. The aims of this study are as follows: (1) to describe new species of deep-sea polynoids using morphology and molecular data (COI, 16S and 18S); and (2) to evaluate the monophyly of Macellicephalinae. Based on molecular and morphological phylogenetic analyses, ten subfamilies are synonymized with Macellicephalinae in order to create a homogeneous clade determined by the absence of lateral antennae. Within this clade, the Anantennata clade was well supported, being determined by the absence of a median antenna. Furthermore, 17 new species and four new genera are described, highlighting the high diversity hidden in the deep. A taxonomic key for the 37 valid genera of the subfamily Macellicephalinae is provided.

Description: Metabolites give us a window into the chemistry of microbes and are split into two subclasses: primary and secondary. Primary metabolites are required for life whereas secondary metabolites have historically been classified as those appearing after exponential growth and are not necessarily needed for survival. Many microbial species are estimated to produce hundreds of metabolites and can be affected by differing nutrients. Using various analytical techniques, metabolites can be directly detected in order to elucidate their biological significance. Currently, a single experiment can produce anywhere from megabytes to terabytes of data. This big data has motivated scientists to develop informatics tools to help target specific metabolites or sets of metabolites. Broadly, it is imperative to identify clear biological questions before embarking on a study of metabolites (metabolomics). For instance, studying the effect of a transposon insertion on phenazine biosynthesis in Pseudomonas is a very different from asking what molecules are present in a specific banana-derived strain of Pseudomonas. This review is meant to serve as a primer for a ‘choose your own adventure’ approach for microbiologists with limited mass spectrometry expertise, with a strong focus on liquid chromatography mass spectrometry based workflows developed or optimized within the past five years.

Description: Three monoterpenoid indole alkaloids (MIAs), tabernabovines A–C (1–3), were isolated from Tabernaemontana bovina. They were elucidated by spectroscopic data and computational calculations. Unlike precursors of MIAs, strictosidine and alstrostine A, alkaloid 1 consists of tryptamine and secologanin in a 2:1 ratio. Alkaloid 2 is a cage compound, and 3 possesses a bridged ring. Tabernabovine A exhibited inhibitory activity against NO production with IC50 44.1 μM compared to l-NMMA with IC50 of 48.6 μM.

Description: Euphorkanlide A (1), a highly modified ingenane diterpenoid with a C24 appendage forming an additional hexahydroisobenzofuran-fused 19-membered macrocyclic bis-lactone ring system was isolated from the roots of Euphorbia kansuensis. Its structure was determined by extensive spectroscopic analysis and quantum-chemical calculations. Compound 1 showed significant cytotoxicities against a panel of cancer cell lines (IC50s 〈 5 μM). Mechanistic study revealed that 1 could induce the generation of ROS, leading to cell cycle arrest and cell apoptosis in drug-resistant cancer cell line HCT-15/5-FU.

Description: The vast amount of plastic waste emitted into the environment and the increasing concern of potential harm to wildlife has made microplastic and nanoplastic pollution a growing environmental concern. Plastic pollution has the potential to cause both physical and chemical harm to wildlife directly or via sorption, concentration, and transfer of other environmental contaminants to the wildlife that ingest plastic. Small particles of plastic pollution, termed microplastics (〉100 nm and 〈5 mm) or nanoplastics (〈100 nm), can form through fragmentation of larger pieces of plastic. These small particles are especially concerning because of their high specific surface area for sorption of contaminants as well as their potential to translocate in the bodies of organisms. These same small particles are challenging to separate and identify in environmental samples because their size makes handling and observation difficult. As a result, our understanding of the environmental prevalence of nanoplastics and microplastics is limited. Generally, the smaller the size of the plastic particle, the more difficult it is to separate from environmental samples. Currently employed passive density and size separation techniques to isolate plastics from environmental samples are not well suited to separate microplastics and nanoplastics. Passive flotation is hindered by the low buoyancy of small particles as well as the difficulty of handling small particles on the surface of flotation media. Here we suggest exploring alternative techniques borrowed from other fields of research to improve separation of the smallest plastic particles. These techniques include adapting active density separation (centrifugation) from cell biology and taking advantage of surface-interaction-based separations from analytical chemistry. Furthermore, plastic pollution is often challenging to quantify in complex matrices such as biological tissues and wastewater. Biological and wastewater samples are important matrices that represent key points in the fate and sources of plastic pollution, respectively. In both kinds of samples, protocols need to be optimized to increase throughput, reduce contamination potential, and avoid destruction of plastics during sample processing. To this end, we recommend adapting digestion protocols to match the expected composition of the nonplastic material as well as taking measures to reduce and account for contamination. Once separated, plastics in an environmental sample should ideally be characterized both visually and chemically. With existing techniques, microplastics and nanoplastics are difficult to characterize or even detect. Their low mass and small size provide limited signal for visual, vibrational spectroscopic, and mass spectrometric analyses. Each of these techniques involves trade-offs in throughput, spatial resolution, and sensitivity. To accurately identify and completely quantify microplastics and nanoplastics in environmental samples, multiple analytical techniques applied in tandem are likely to be required.

Description: Noninvasive diagnostic by imaging combined with a contrast agent (CA) is by now the most used technique to get insight into human bodies. X-ray and magnetic resonance imaging (MRI) are widely used technologies providing complementary results. Nowadays, it seems clear that bimodal CAs could be an emerging approach to increase the patient compliance, accessing different imaging modalities with a single CA injection. Owing to versatile designs, targeting properties, and high payload capacity, nanocarriers are considered as a viable solution to reach this goal. In this study, we investigated efficient superparamagnetic iron oxide nanoparticle (SPION)-loaded iodinated nano-emulsions (NEs) as dual modal injectable CAs for X-ray imaging and MRI. The strength of this new CA lies not only in its dual modal contrasting properties and biocompatibility, but also in the simplicity of the nanoparticulate assembling: iodinated oily core was synthesized by the triiodo-benzene group grafting on vitamin E (41.7% of iodine) via esterification, and SPIONs were produced by thermal decomposition during 2, 4, and 6 h to generate SPIONs with different morphologies and magnetic properties. SPIONs with most anisotropic shape and characterized by the highest r2/r1 ratio once encapsulated into iodinated NE were used for animal experimentation. The in vivo investigation showed an excellent contrast modification because of the presence of the selected NEs, for both imaging techniques explored, that is, MRI and X-ray imaging. This work provides the description and in vivo application of a simple and efficient nanoparticulate system capable of enhancing contrast for both preclinical imaging modalities, MRI, and computed tomography.

Description: Understanding mechanical interactions between hydrate and hosting sediments is critical for evaluating formation stability and associated environmental impacts of hydrate-bearing sediments during gas production. While core-scale studies of hydrate-bearing sediments are readily available and some explanations of observed results rely on pore-scale behavior of hydrate, actual pore-scale observations supporting the larger-scale phenomena are rarely available for hydrate-bearing sediments, especially with methane as guest molecules. The primary reasons for the scarcity include the challenge of developing tools for small-scale testing apparatus and pore-scale visualization capability. We present a testing assembly that combines pore-scale visualization and triaxial test capability of methane hydrate-bearing sediments. This testing assembly allows temperature regulation and independent control of four pressures: influent and effluent pore pressure, confining pressure, and axial pressure. Axial and lateral effective stresses can be applied independently to a 9.5 mm diameter and 19 mm long specimen while the pore pressure and temperature are controlled to maintain the stability of methane hydrate. The testing assembly also includes an X-ray transparent beryllium core holder so that 3D computed tomography scanning can be conducted during the triaxial loading. This testing assembly permits pore-scale exploration of hydrate-sediment interaction in addition to the traditional stress-strain relationship. Exemplary outcomes are presented to demonstrate applications of the testing assembly on geomechanical property estimations of methane-hydrate bearing sediments.

Description: Marine prey and predators will respond to future climate through physiological and behavioral adjustments. However, our understanding of how such direct effects may shift the outcome of predator–prey interactions is still limited. Here, we investigate the effects of ocean warming and acidification on foraging behavior and biomass of a common prey (shrimps, Palaemon spp.) tested in large mesocosms harboring natural resources and habitats. Acidification did not alter foraging behavior in prey. Under warming, however, prey showed riskier behavior by foraging more actively and for longer time periods, even in the presence of a live predator. No effects of longer-term exposure to climate stressors were detected on prey biomass. Our findings suggest that ocean warming may increase the availability of some prey to predators via a behavioral pathway (i.e., increased risk-taking by prey), likely by elevating metabolic demand of prey species.

Description: Climate change influences the Baltic Sea ecosystem via its effects on oceanography and biogeochemistry. Sea surface temperature has been projected to increase by 2 to 4 °C until 2100 due to global warming; the changes will be more significant in the northern areas and less so in the south. The warming up will also diminish the annual sea ice cover by 57% to 71%, and ice season will be one to three months shorter than in the early 21st century, depending on latitude. A significant decrease in sea surface salinity has been projected because of an increase in rainfall and decrease of saline inflows into the Baltic Sea. The increasing surface flow has, in turn, been projected to increase leaching of nutrients from the soil to the watershed and eventually into the Baltic Sea. Also, acidification of the seawater and sea-level rise have been predicted. Increasing seawater temperature speeds up metabolic processes and increases growth rates of many secondary producers. Species associated with sea ice, from salt brine microbes to seals, will suffer. Due to the specific salinity tolerances, species’ geographical ranges may shift by tens or hundreds of kilometres with decreasing salinity. A decrease in pH will slow down calcification of bivalve shells, and higher temperatures also alleviate establishment of non-indigenous species originating from more southern sea areas. Many uncertainties still remain in predicting the couplings between atmosphere, oceanography and ecosystem. Especially projections of many oceanographic parameters, such as wind speeds and directions, the mean salinity level, and density stratification, are still ambiguous. Also, the effects of simultaneous changes in multiple environmental factors on species with variable preferences to temperature, salinity, and nutrient conditions are difficult to project. There is, however, enough evidence to claim that due to increasing runoff of nutrients from land and warming up of water, primary production and sedimentation of organic matter will increase; this will probably enhance anoxia and release of phosphorus from sediments. Such changes may keep the Baltic Sea in an eutrophicated state for a long time, unless strong measures to decrease nutrient runoff from land are taken. Changes in the pelagic and benthic communities are anticipated. Benthic communities will change from marine to relatively more euryhaline communities and will suffer from hypoxic events. The projected temperature increase and salinity decline will contribute to maintain the pelagic ecosystem of the Central Baltic and the Gulf of Finland in a state dominated by cyanobacteria, flagellates, small-sized zooplankton and sprat, instead of diatoms, large marine copepods, herring, and cod. Effects vary from area to area, however. In particular the Bothnian Sea, where hypoxia is less common and rivers carry a lot of dissolved organic carbon, primary production will probably not increase as much as in the other basins. The coupled oceanography-biogeochemistry ecosystem models have greatly advanced our understanding of the effects of climate change on marine ecosystems. Also, studies on climate associated “regime shifts” and cascading effects from top predators to plankton have been fundamental for understanding of the response of the Baltic Sea ecosystem to anthropogenic and climatic stress. In the future, modeling efforts should be focusing on coupling of biogeochemical processes and lower trophic levels to the top predators. Also, fine resolution species distribution models should be developed and combined with 3-D modelling, to describe how the species and communities are responding to climate-induced changes in environmental variables.

Description: UPLC–TOF/MS profiling, followed by the recently reported differential off-line LC–NMR (DOLC–NMR) and quantitative 1H NMR spectroscopy (qHNMR), led to the differential qualitative analysis and accurate quantitation of l-tryptophan-induced metabolome alterations of Penicillium roqueforti, which is typically used in making blue-mold cheese. Among the 24 metabolites identified, two tetrapeptides, namely, d-Phe-l-Val-d-Val-l-Tyr and d-Phe-l-Val-d-Val-l-Phe, as well as cis-bis(methylthio)silvatin, are reported for the first time as metabolites of P. roqueforti. Antimicrobial activity tests showed strong effects of the catabolic l-tryptophan metabolites 3-hydroxyanthranilic acid, anthranilic acid, and 3-indolacetic acid against Saccharomyces cerevisiae, with IC50 values between 15.6 and 24.0 μg/mL, while roquefortine C and cis-bis(methylthio)silvatin inhibited the growth of Gram-negative Escherichia coli and Gram-positive Bacillus subtilis with IC50 values between 30.0 and 62.5 μg/mL.

Description: This review shows that the presence of seagrass microbial community is critical for the development of seagrasses; from seed germination, through to phytohormone production and enhanced nutrient availability, and defence against pathogens and saprophytes. The tight seagrass-bacterial relationship highlighted in this review supports the existence of a seagrass holobiont and adds to the growing evidence for the importance of marine eukaryotic microorganisms in sustaining vital ecosystems. Incorporating a micro-scale view on seagrass ecosystems substantially expands our understanding of ecosystem functioning and may have significant implications for future seagrass management and mitigation against human disturbance.

Description: Diurnal vertical migration (DVM) is a widespread phenomenon in the upper ocean, but it remains unclear to what degree it also involves passively transported micro- and meso-zooplankton. These organisms are difficult to monitor by in situ sensing and observations from discrete samples are often inconclusive. Prime examples of such ambiguity are planktonic foraminifera, where contradictory evidence for DVM continues to cast doubt on the stability of species vertical habitats, which introduces uncertainties in geochemical proxy interpretation. To provide a robust answer, we carried out highly replicated randomised sampling with 41 vertically resolved plankton net hauls taken within 26 hours in a confined area of 400 km2 in the tropical North Atlantic, where DVM in larger plankton occurs. Manual enumeration of planktonic foraminifera cell density consistently reveals the highest total cell concentrations in the surface mixed layer (top 50 m) and analysis of cell density in seven individual species representing different shell sizes, life strategies and presumed depth habitats reveals consistent vertical habitats not changing over the 26 hours sampling period. These observations robustly reject the existence of DVM in planktonic foraminifera in a setting where DVM occurs in other organisms.

Description: Sponges are one of the most dominant organisms in marine ecosystems. One reason for their success is their association with microorganisms that are besides the host itself responsible for the chemical defence. Sponge abundances have been increasing on coral reefs in the Western Indian Ocean (WIO) and are predicted to increase further with rising anthropogenic impacts on coral reefs. However, there is a paucity of information on chemical ecology of sponges from the WIO and their prokaryotic community composition. We used a combination of Illumina sequencing and a predictive metagenomic analysis to (i) assess the prokaryotic community composition of sponges from Zanzibar, (ii) predict the presence of KEGG metabolic pathways responsible for bioactive compound production and (iii) relate their presence to the degree of observed chemical defence in their respective sponge host. We found that sponges from Zanzibar host diverse prokaryotic communities that are host species-specific. Sponge-species and respective specimens that showed strong chemical defences in previous studies were also predicted to be highly enriched in various pathways responsible for secondary metabolite production. Hence, the combined sequencing and predictive metagenomic approach proved to be a useful indicator for the metabolic potential of sponge holobionts

Description: Microplastics (MPs) in aquatic organisms are raising increasing concerns regarding their potential damage to ecosystems. To date, Raman and Fourier transform infrared spectroscopy techniques have been widely used for detection of MPs in aquatic organisms, which requires complex protocols of tissue digestion and MP separation and are time- and reagentconsuming. This novel approach directly separates, identifies, and characterizes MPs from the hyperspectral image (HSI) of the intestinal tract content in combination with a support vector machine classification model, instead of using the real digestion/separation protocols. The procedures of HSI acquisition ( 1 min) and data analysis (5 min) can be completed within 6 min plus the sample preparation and drying time (30 min) where necessary. This method achieved a promising efficiency (recall 〉98.80%, precision 〉96.22%) for identifying five types of MPs (particles 〉0.2 mm). Moreover, the method was also demonstrated to be effective on field fish from three marine fish species, revealing satisfying detection accuracy (particles 〉0.2 mm) comparable to Raman analysis. The present technique omits the digestion protocol (reagent free), thereby significantly reducing reagent consumption, saving time, and providing a rapid and efficient method for MP analysis.

Description: Repeated and independent emergence of trait divergence that matches habitat differences is a sign of parallel evolution by natural selection. Yet, the molecular underpinnings that are targeted by adaptive evolution often remain elusive. We investigate this question by combining genome-wide analyses of copy number variants (CNVs), single nucleotide polymorphisms (SNPs), and gene expression across four pairs of lake and river populations of the three-spined stickleback (Gasterosteus aculeatus). We tested whether CNVs that span entire genes and SNPs occurring in putative cis-regulatory regions contribute to gene expression differences between sticklebacks from lake and river origins. We found 135 gene CNVs that showed a significant positive association between gene copy number and gene expression, suggesting that CNVs result in dosage effects that can fuel phenotypic variation and serve as substrates for habitat-specific selection. Copy number differentiation between lake and river sticklebacks also contributed to expression differences of two immune-related genes in immune tissues, cathepsin A and GIMAP7. In addition, we identified SNPs in cis-regulatory regions (eSNPs) associated with the expression of 1,865 genes, including one eSNP upstream of a carboxypeptidase gene where both the SNP alleles differentiated and the gene was differentially expressed between lake and river populations. Our study highlights two types of mutations as important sources of genetic variation involved in the evolution of gene expression and in potentially facilitating repeated adaptation to novel environments.

Description: Transparent exopolymer particles (TEP) are polysaccharide-rich microgels that are prevalent in the marine environment and have important roles in the aggregation of organic matter and carbon export from the euphotic zone. TEP are readily colonized by bacteria and utilized by specialized taxa, such as Alteromonadaceae. However, bacterial community composition specifically attached to natural TEP remains largely unknown. In this study, we isolated individual TEP from Plymouth Sound (UK) and performed DNA sequencing of the TEP-attached bacterial communities. We also sampled the cognate bulk seawater total bacterial communities for comparison. The bacterial communities associated with individual TEP showed distinct differences compared to the total bulk bacterioplankton communities, with Alteromonadaceae significantly more abundant on TEP. The TEP-associated Alteromonadaceae consisted of two operational taxonomic units that were closely related to Marinobacter and Glaciecola, both previously associated with biogenic aggregates and microgel-rich habitats. This study provides novel insight into marine bacterial–microgel interactions.

Description: Cichlid fishes provide textbook examples of explosive phenotypic diversification and sympatric speciation, thereby making them ideal systems for studying the molecular mechanisms underlying rapid lineage divergence. Despite the fact that gene regulation provides a critical link between diversification in gene function and speciation, many genomic regulatory mechanisms such as microRNAs (miRNAs) have received little attention in these rapidly diversifying groups. Therefore, we investigated the posttranscriptional regulatory role of miRNAs in the repeated sympatric divergence of Midas cichlids (Amphilophus spp.) from Nicaraguan crater lakes. Using miRNA and mRNA sequencing of embryos from five Midas species, we first identified miRNA binding sites in mRNAs and highlighted the presences of a surprising number of novel miRNAs in these adaptively radiating species. Then, through analyses of expression levels, we identified putative miRNA/gene target pairs with negatively correlated expression level that were consistent with the role of miRNA in downregulating mRNA. Furthermore, we determined that several miRNA/gene pairs show convergent expression patterns associated with the repeated benthic/limnetic sympatric species divergence implicating these miRNAs as potential molecular mechanisms underlying replicated sympatric divergence. Finally, as these candidate miRNA/gene pairs may play a central role in phenotypic diversification in these cichlids, we characterized the expression domains of selected miRNAs and their target genes via in situ hybridization, providing further evidence that miRNA regulation likely plays a role in the Midas cichlid adaptive radiation. These results provide support for the hypothesis that extremely quickly evolving miRNA regulation can contribute to rapid evolutionary divergence even in the presence of gene flow.

Description: Microplastics are ubiquitous across ecosystems, yet the exposure risk to humans is unresolved. Focusing on the American diet, we evaluated the number of microplastic particles in commonly consumed foods in relation to their recommended daily intake. The potential for microplastic inhalation and how the source of drinking water may affect microplastic consumption were also explored. Our analysis used 402 data points from 26 studies, which represents over 3600 processed samples. Evaluating approximately 15% of Americans’ caloric intake, we estimate that annual microplastics consumption ranges from 39000 to 52000 particles depending on age and sex. These estimates increase to 74000 and 121000 when inhalation is considered. Additionally, individuals who meet their recommended water intake through only bottled sources may be ingesting an additional 90000 microplastics annually, compared to 4000 microplastics for those who consume only tap water. These estimates are subject to large amounts of variation; however, given methodological and data limitations, these values are likely underestimates.

Description: The Arctic Limnocalanus macrurus is a prominent representative of large copepods which performs several essential functions in freshwater and marine ecosystems. Being a cold stenotherm species, its distribution is primarily confined to deeper water layers. Based on the long-term observations from one of the largest spatially confined natural populations of this species in the Baltic Sea, we detected profound long-term variability of L. macrurus during 1958–2016: high abundances before the 1980s, then nearly disappearance in the 1990s and recovery in the 2000s. The main environmental parameters explaining the interannual variability of L. macrurus in spring were herring spawning stock biomass in preceding year, winter severity, and bottom water temperature in preceding summer. The effect of winter severity and water temperature was also non-linear. The sliding window correlation analysis pointed to a non-stationary relationship between the abundance of L. macrurus and the key variables. Given the observed pronounced seasonality in the population structure of L. macrurus (young stages dominated in the beginning of the year and only adults were left in the population in summer and autumn) we identified the dynamics of key environmental variables to understand this species under different ecosystem configurations and different combinations of drivers of change.

Description: The Swan Islands Transform Fault (SITF) marks the southern boundary of the Cayman Trough and the ocean–continent transition of the North American–Caribbean Plate boundary offshore Honduras. The CAYSEIS experiment acquired a 180-km-long seismic refraction and gravity profile across this transform margin, ∼70 km to the west of the Mid-Cayman Spreading Centre (MCSC). This profile shows the crustal structure across a transform fault system that juxtaposes Mesozoic-age continental crust to the south against the ∼10-Myr-old ultraslow spread oceanic crust to the north. Ocean-bottom seismographs were deployed along-profile, and inverse and forward traveltime modelling, supported by gravity analysis, reveals ∼23-km-thick continental crust that has been thinned over a distance of ∼70 km to ∼10 km-thick at the SITF, juxtaposed against ∼4-km-thick oceanic crust. This thinning is primarily accommodated within the lower crust. Since Moho reflections are not widely observed, the 7.0 km s−1 velocity contour is used to define the Moho along-profile. The apparent lack of reflections to the north of the SITF suggests that the Moho is more likely a transition zone between crust and mantle. Where the profile traverses bathymetric highs in the off-axis oceanic crust, higher P-wave velocity is observed at shallow crustal depths. S-wave arrival modelling also reveals elevated velocities at shallow depths, except for crust adjacent to the SITF that would have occupied the inside corner high of the ridge-transform intersection when on axis. We use a Vp/Vs ratio of 1.9 to mark where lithologies of the lower crust and uppermost mantle may be exhumed, and also to locate the upper-to-lower crustal transition, identify relict oceanic core complexes and regions of magmatically formed crust. An elevated Vp/Vs ratio suggests not only that serpentinized peridotite may be exposed at the seafloor in places, but also that seawater has been able to flow deep into the crust and upper mantle over 20–30-km-wide regions which may explain the lack of a distinct Moho. The SITF has higher velocities at shallower depths than observed in the oceanic crust to the north and, at the seabed, it is a relatively wide feature. However, the velocity–depth model subseabed suggests a fault zone no wider than ∼5–10 km, that is mirrored by a narrow seabed depression ∼7500 m deep. Gravity modelling shows that the SITF is also underlain, at 〉2 km subseabed, by a ∼20-km-wide region of density 〉3000 kg m−3 that may reflect a broad region of metamorphism. The residual mantle Bouguer anomaly across the survey region, when compared with the bathymetry, suggests that the transform may also have a component of left-lateral trans-tensional displacement that accounts for its apparently broad seabed appearance, and that the focus of magma supply may currently be displaced to the north of the MCSC segment centre. Our results suggest that Swan Islands margin development caused thinning of the adjacent continental crust, and that the adjacent oceanic crust formed in a cool ridge setting, either as a result of reduced mantle upwelling and/or due to fracture enhanced fluid flow.

Description: Latest knowledge on the reactivity of charged nanoparticulate complexants toward aqueous metal ions is discussed in mechanistic detail. We present a rigorous generic description of electrostatic and chemical contributions to metal ion binding by nanoparticulate complexants, and their dependence on particle size, particle type (i.e., reactive sites distributed within the particle body or confined to the surface), ionic strength of the aqueous medium, and the nature of the metal ion. For the example case of soft environmental particles such as fulvic and humic acids, practical strategies are delineated for determining intraparticulate metal ion speciation, and for evaluating intrinsic chemical binding affinities and heterogeneity. The results are compared with those obtained by popular codes for equilibrium speciation modeling (namely NICA-Donnan and WHAM). Physicochemical analysis of the discrepancies generated by these codes reveals the a priori hypotheses adopted therein and the inappropriateness of some of their key parameters. The significance of the characteristic time scales governing the formation and dissociation rates of metal−nanoparticle complexes in defining the relaxation properties and the complete equilibration of the metal− nanoparticulate complex dispersion is described. The dynamic features of nanoparticulate complexes are also discussed in the context of predictions of the labilities and bioavailabilities of the metal species.

Description: Standard seismic acquisition and processing require appropriate source-receiver offsets. P-cable technology represents the opposite, namely, very short source-receiver offsets at the price of increased spatial and lateral resolution with a high-frequency source. To use this advantage, a processing flow excluding offset information is required. This aim can be achieved with a processing tuned to diffractions because point diffractions scatter the same information in offset and midpoint direction. Usually, diffractions are small amplitude events and a careful diffraction separation is required as a first step. We suggest the strategy to use a multiparameter stacking operator, e.g, common-reflection surface, and stack along the midpoint direction. The obtained kinematic wavefront attributes are used to calculate time-migration velocities. A diffractivity map serves as filter to refine the velocities. This strategy is applied to a 3D P-cable data set to obtain a time-migrated image.

Description: There is a recognized need for new methods with modest data requirements to provide preliminary estimates of stock status for data-limited stocks (e.g. Rudd and Thorson, 2018). Froese et al. (2018) provide such a method, which derives estimates of relative stock size from length frequency (LF) data of exploited stocks. They show that their length-based Bayesian biomass estimation method (LBB) can reproduce the “true” parameters used in simulated data and can approximate the relative stock size as estimated independently by more data-demanding methods in 34 real stocks. However, in a comment on LBB, Hordyk et al. (2019) claim (i) that the master equation of LBB is incomplete because it does not correct for the pile-up effect caused by aggregating length measurements into length classes or “bins”, (ii) that LBB is highly sensitive to equilibrium assumptions and wrongly uses maximum observed length (Lmax) for guidance in setting a prior for the estimation of asymptotic length (Linf), and (iii) that the default prior used by LBB for the ratio between natural mortality and somatic growth rate (M/K) of 1.5 (SD = 0.15) is inadequate for many exploited species. These comments are addressed below

Description: On 25 December 2016, a Mw 7.6 earthquake broke a portion of the Southern Chilean subduction zone south of Chiloé Island, located in the central part of the Mw 9.5 1960 Valdivia earthquake. This region is characterized by repeated earthquakes in 1960 and historical times with very sparse interseismic activity due to the subduction of a young (~15 Ma), and therefore hot, oceanic plate. We estimate the co-seismic slip distribution based on a kinematic finite fault source model, and through joint inversion of teleseismic body waves and strong motion data. The coseismic slip model yields a total seismic moment of 3.94×1020 Nm that occurred over ~30 s, with the rupture propagating mainly downdip, reaching a peak-slip of ~4.2 m. Regional moment tensor inversion of stronger aftershocks reveals thrust type faulting at depths of the plate interface. The fore- and aftershock seismicity is mostly related to the subduction interface with sparse seismicity in the overriding crust. The 2016 Chiloé event broke a region with increased locking and most likely broke an asperity of the 1960 earthquake. The updip limit of the main event, aftershocks, foreshocks and interseismic activity are spatially similar, located ~15 km offshore and parallel to Chiloé Islands west coast. The coseismic slip model of the 2016 Chiloé earthquake suggests a peak slip of 4.2 m that locally exceeds the 3.38 m slip deficit that has accumulated since 1960. Therefore, the 2016 Chiloé earthquake possibly released strain that has built up prior to the 1960 Valdivia earthquake.

Description: We used a molecular phylogenetic approach to investigate species delimitation and diversification in the northeastern Atlantic and Mediterranean musseldrills of the Ocinebrina aciculata complex, based on molecular data from topotypical material of many of the nominal taxa. The complex is shown to consist of at least five species: Ocinebrina aciculata (Lamarck, 1822) from the Atlantic and western Mediterranean; O. cf. corallina (Scacchi, 1836) from the central Mediterranean Sea; O. reinai Bonomolo & Crocetta, 2012 from the Tyrrhenian Sea; O. corallinoides Pallary, 1912 from the Gulf of Gabès; and O. aegeensis n. sp. currently known from the Aegean Sea only. The new species is differentiated from the other taxa by very subtle morphological diagnostic features, although it is clearly identified by genetic distance and apomorphic DNA-sequence characters. The identity of Murex corallinus Scacchi, 1836 (type species of Ocinebrina Jousseaume, 1880) could not be defined with certainty, pending genetic comparison of specimens of the â €? large Tyrrhenian morphotype' (corresponding to the neotype, but not assayed herein) with the assayed â €? small morphotype'.

Description: Marine sponges are early-branched metazoans known to harbor dense and diverse microbial communities. Yet the role of the so far uncultivable alphaproteobacterial lineages that populate these sessile invertebrates remains unclear. We applied a sequence composition-dependent binning approach to assemble one Rhodospirillaceae genome from the Spongia officinalis microbial metagenome and contrast its functional features with those of closely related sponge-associated and free-living genomes. Both symbiotic and free-living Rhodospirillaceae shared a suite of common features, possessing versatile carbon, nitrogen, sulfur and phosphorus metabolisms. Symbiotic genomes could be distinguished from their free-living counterparts by the lack of chemotaxis and motility traits, enrichment of genes required for the uptake and utilization of organic sulfur compounds—particularly taurine—, higher diversity and abundance of ABC transporters, and a distinct repertoire of genes involved in natural product biosynthesis, plasmid stability, cell detoxification and oxidative stress remediation. These sessile symbionts may more effectively contribute to host fitness via nutrient exchange, and also host detoxification and chemical defense. Considering the worldwide occurrence and high diversity of sponge-associated Rhodospirillaceae verified here using a tailored in silico approach, we suggest that these organisms are not only relevant to holobiont homeostasis but also to nutrient cycling in benthic ecosystems.

Description: Continental hyperextension during magma-poor rifting at the Deep Galicia Margin is characterised by a complex pattern of faulting, thin continental fault blocks, and the serpentinisation, with local exhumation, of mantle peridotites along the S-reflector, interpreted as a detachment surface. In order to understand fully the evolution of these features, it is important to image seismically the structure and to model the velocity structure to the greatest resolution possible. Travel-time tomography models have revealed the long-wavelength velocity structure of this hyperextended domain, but are often insufficient to match accurately the short-wavelength structure observed in reflection seismic imaging. Here we demonstrate the application of two-dimensional (2D) time-domain acoustic full-waveform inversion to deep water seismic data collected at the Deep Galicia Margin, in order to attain a high resolution velocity model of continental hyperextension. We have used several quality assurance procedures to assess the velocity model, including comparison of the observed and modelled waveforms, checkerboard tests, testing of parameter and inversion strategy, and comparison with the migrated reflection image. Our final model exhibits an increase in the resolution of subsurface velocities, with particular improvement observed in the westernmost continental fault blocks, with a clear rotation of the velocity field to match steeply dipping reflectors. Across the S-reflector there is a sharpening in the velocity contrast, with lower velocities beneath S indicative of preferential mantle serpentinisation. This study supports the hypothesis that normal faulting acts to hydrate the upper mantle peridotite, observed as a systematic decrease in seismic velocities, consistent with increased serpentinisation. Our results confirm the feasibility of applying the full-waveform inversion method to sparse, deep water crustal datasets.

Description: Overfishing and rapid environmental shifts pose severe challenges to the resilience and viability of marine fish populations. To develop and implement measures that enhance species’ adaptive potential to cope with those pressures while, at the same time, ensuring sustainable exploitation rates is part of the central goal of fisheries management. Here, we argue that a combination of biophysical modelling and population genomic assessments offer ideal management tools to define stocks, their physical connectivity and ultimately, their short-term adaptive potential. To date, biophysical modelling has often been confined to fisheries ecology whereas evolutionary hypotheses remain rarely considered. When identified, connectivity patterns are seldom explored to understand the evolution and distribution of adaptive genetic variation, a proxy for species’ evolutionary potential. Here, we describe a framework that expands on the conventional seascape genetics approach by using biophysical modelling and population genomics. The goals are to identify connectivity patterns and selective pressures, as well as putative adaptive variants directly responding to the selective pressures and, ultimately, link both to define testable hypotheses over species response to shifting ecological conditions and overexploitation.

Description: This study presents a new method (LBB) for the analysis of length frequency data from commercial catches. LBB works for species that grow throughout their lives, such as most commercially-important fish and invertebrates, and requires no input in addition to length frequency data. It estimates asymptotic length, length at first capture, relative natural mortality, and relative fishing mortality. Standard fisheries equations can then be used to approximate current exploited biomass relative to unexploited biomass. In addition, these parameters allow the estimation of length at first capture that would maximize catch and biomass for a given fishing effort, and estimation of a proxy for the relative biomass capable of producing maximum sustainable yields. Relative biomass estimates of LBB were not significantly different from the “true” values in simulated data and were similar to independent estimates from full stock assessments. LBB also presents a new indicator for assessing whether an observed size structure is indicative of a healthy stock. LBB results will obviously be misleading if the length frequency data do not represent the size composition of the exploited size range of the stock or if length frequencies resulting from the interplay of growth and mortality are masked by strong recruitment pulses.

Description: We present high-resolution resistivity imaging of gas hydrate pipe-like structures, as derived from marine controlled-source electromagnetic (CSEM) inversions that combine towed and ocean-bottom electric field receiver data, acquired from the Nyegga region, offshore Norway. Two-dimensional CSEM inversions applied to the towed receiver data detected four new prominent vertical resistive features that are likely gas hydrate structures, located in proximity to a major gas hydrate pipe-like structure, known as the CNE03 pockmark. The resistivity model resulting from the CSEM data inversion resolved the CNE03 hydrate structure in high resolution, as inferred by comparison to seismically constrained inversions. Our results indicate that shallow gas hydrate vertical features can be delineated effectively by inverting both ocean-bottom and towed receiver CSEM data simultaneously. The approach applied here can be utilised to map and monitor seafloor mineralisation, freshwater reservoirs, CO2 sequestration sites and near-surface geothermal systems.

Description: In this study, the complex frequency-shifted perfectly matched layer (CFS-PML) in stretching Cartesian coordinates, is successfully applied to three-dimensional (3D) frequency-domain marine controlled-source electromagnetic (CSEM) field modelling. The Dirichlet boundary, which is usually used within the traditional framework of EM modeling algorithms, assumes the electric or magnetic field values are zero at the boundaries. This requires the boundaries be sufficiently far away from the sources in the area of interest. To mitigate the boundary artifacts, a large modelling area may be necessary even though cell sizes are allowed to grow toward the boundaries due to the diffusion of the electromagnetic wave propagation. Compared with the conventional Dirichlet boundary, the PML boundary is preferred as the modelling area of interest could be restricted to the target region and only a few absorbing layers surrounding can effectively depress the artificial boundary effect without losing the numerical accuracy. Furthermore, for joint inversion of seismic and marine CSEM data, if we used the PML for CSEM field simulation instead of the conventional Dirichlet, the modeling area for these two different geophysical data collected from the same survey area could the same, which is convenient for joint inversion grid matching. We apply the CFS-PML boundary to 3D marine CSEM modelling by using the staggered finite-difference (SFD) discretization. Numerical test indicates that the modeling algorithm using the CFS-PML also shows good accuracy compared to the Dirichlet. Furthermore, the modeling algorithm using the CFS-PML shows advantages in computational time and memory saving than that using the Dirichlet boundary. For the 3D example in this study, the memory saving using the PML is nearly 42 % and the time saving is around 48% compared to using the Dirichlet.

Description: Pseudovibrio is a marine bacterial genus members of which are predominantly isolated from sessile marine animals, and particularly sponges. It has been hypothesised that Pseudovibrio spp. form mutualistic relationships with their hosts. Here, we studied Pseudovibrio phylogeny and genetic adaptations that may play a role in host colonization by comparative genomics of 31 Pseudovibrio strains, including 25 sponge isolates. All genomes were highly similar in terms of encoded core metabolic pathways, albeit with substantial differences in overall gene content. Based on gene composition, Pseudovibrio spp. clustered by geographic region, indicating geographic speciation. Furthermore, the fact that isolates from the Mediterranean Sea clustered by sponge species suggested host-specific adaptation or colonization. Genome analyses suggest that Pseudovibrio hongkongensis UST20140214-015BT is only distantly related to other Pseudovibrio spp., thereby challenging its status as typical Pseudovibrio member. All Pseudovibrio genomes were found to encode numerous proteins with SEL1 and tetratricopeptide repeats, which have been suggested to play a role in host colonization. For evasion of the host immune system, Pseudovibrio spp. may depend on type III, IV and VI secretion systems that can inject effector molecules into eukaryotic cells. Furthermore, Pseudovibrio genomes carry on average seven secondary metabolite biosynthesis clusters, reinforcing the role of Pseudovibrio spp. as potential producers of novel bioactive compounds. Tropodithietic acid, bacteriocin and terpene biosynthesis clusters were highly conserved within the genus, suggesting an essential role in survival e.g. through growth inhibition of bacterial competitors. Taken together, these results support the hypothesis that Pseudovibrio spp. have mutualistic relations with sponges.

Description: The BONUS symposium “Science delivery for sustainable use of the Baltic Sea living resources” held in Tallinn, Estonia, in October 2017 was an opportunity for the presentation and discussion of 107 papers that examined the state and dynamics of living resources of the Baltic Sea, and associated management challenges. The symposium included a half-day stakeholder panel discussion that addressed the main challenges related to sustainable management and matching research and policy/management needs. Based on the five symposium papers published in this Special Issue as well as the stakeholder panel discussion, it can be concluded that (i) new observations about the feeding ecology of clupeids supports a more complete understanding of trophic interactions in the pelagic realm and improved calibration of ecosystem models, (ii) to safequard sustainable and diverse fisheries resources, one should take into account the specific local characteristics of the fish community, (iii) to safeguard sustainable use of marine resources and mitigate cross-sectoral and transboundary conflicts, a risk-based approach should be adopted, and (iv) incorporation of scientific advice into management faces several obstacles including the reality that not all readily available knowledge is currently being incorporated into the decision-making process.

Description: Fisheries and marine ecosystem-based management requires a holistic understanding of the dynamics of fish communities and their responses to changes in environmental conditions. Environmental conditions can simultaneously shape the spatial distribution and the temporal dynamics of a population, which together can trigger changes in the functional structure of communities. Here, we developed a comprehensive framework based on complementary multivariate statistical methodologies to simultaneously investigate the effects of environmental conditions on the spatial, temporal and functional dynamics of species assemblages. The framework is tested using survey data collected during more than 4000 fisheries hauls over the Baltic Sea between 2001 and 2016. The approach revealed the Baltic fish community to be structured into three sub-assemblages along a strong and temporally stable salinity gradient decreasing from West to the East. Additionally, we highlight a mismatch between species and functional richness associated with a lower functional redundancy in the Baltic Proper compared with other sub-areas, suggesting an ecosystem more susceptible to external pressures. Based on a large dataset of community data analysed in an innovative and comprehensive way, we could disentangle the effects of environmental changes on the structure of biotic communities-key information for the management and conservation of ecosystems.

Description: High-resolution 3D (HR3D) seismic data are important for hydrocarbon exploration of shallow reservoirs, site characterization, and geohazard assessments. The goal of this contribution is to identify and quantify the parameters to increase the resolution of HR3D seismic data to meter scale. The main acquisition parameters controlling the resolution of the collected data are the spectrum of the seismic source, source-receiver offset range, and trace density. An evolution to one-meter-scale resolution of 3D seismic will rely on combining a reproducible seismic source with high frequencies up to at least 600 Hz, a high uniform trace density of more than 4 million traces per square kilometer, and an offset range shorter than approximately 200 m. The resulting 3D seismic data volume will reach meter-scale resolution for water and target depths of less than 600 m. The proposed HR3D system will be suitable for 3D and 4D characterization of seabed properties and shallow stratigraphy, the identification of geohazards and hydrocarbon leakage, and monitoring the environmental impact of offshore activities. The P-Cable 3D system is an excellent starting point for achieving one-meter-scale resolution due to its flexible and tight meter-scale shot and receiver spacing.

Description: The European and American eels spawn in the subtropical convergence zone (STCZ) in the Sargasso Sea, a dynamic and relatively productive area that is strongly influenced by front and eddy formations and subducted high-saline water masses. To understand how the physical and biological environments may affect the early life history of eels, we conducted a detailed bio-physical investigation of the water column at a site of high eel larvae abundance. Diel measurements and sampling in the upper 300 m revealed strong variations in hydrographic conditions and mean depths of different taxonomic groups; however, characteristics patterns of distribution were apparent. Most species showed diel vertical migrations, ascending about 20-30 m at night, whereas examples of night-time downward migration were also seen. European eel larvae were among the species showing more extensive diel vertical migration: their population mean depth changed from 160 m at day to 100 m at night where abundance peaked at 45 m depth. Distribution and migration of eel larvae corresponded to patterns observed for small hydrozoans, supporting a proposed predator-prey linkage. The study demonstrates the diverse and vertically strongly structured plankton community of STCZ where larvae of eel and other fish find a wide range of potential niches.

Description: An indoor mesocosm experiment was carried out to investigate the combined effects of ocean acidification and warming on the species composition and biogeochemical element cycling during a winter/spring bloom with a natural phytoplankton assemblage from the Kiel fjord, Germany. The experimental setup consisted of a "Control" (ambient temperature of similar to 4.8 degrees C and similar to 535 +/- 25 mu atm pCO(2)), a "High-CO2" treatment (ambient temperature and initially 1020 +/- 45 mu atm pCO(2)) and a "Greenhouse" treatment (similar to 8.5 degrees C and initially 990 +/- 60 mu atm pCO(2)). Nutrient replete conditions prevailed at the beginning of the experiment and light was provided at in situ levels upon reaching pCO(2) target levels. A diatom-dominated bloom developed in all treatments with Skeletonema costatum as the dominant species but with an increased abundance and biomass contribution of larger diatom species in the Greenhouse treatment. Conditions in the Greenhouse treatment accelerated bloom development with faster utilization of inorganic nutrients and an earlier peak in phytoplankton biomass compared to the Control and High CO2 but no difference in maximum concentration of particulate organic matter (POM) between treatments. Loss of POM in the Greenhouse treatment, however, was twice as high as in the Control and High CO2 treatment at the end of the experiment, most likely due to an increased proportion of larger diatom species in that treatment. We hypothesize that the combination of warming and acidification can induce shifts in diatom species composition with potential feedbacks on biogeochemical element cycling.

Description: Light and nutrients are essential resources for phytoplankton growth and considered to shape the size structure and other morphometric traits (surface:volume ratio, deviation from spherical shape) of phytoplankton communities. If morphometric traits influence the growth and resource use, shifts by one of the two factors should influence the capability to utilize the other factor. We performed a two-step experiment, where a natural phytoplankton community was first exposed to three different light levels (supposed to be limiting, saturating and slightly inhibiting for the majority of species) and grown until stationary phase. Then, the pre-conditioned communities were split into two nutrient treatments (control and saturating nutrient pulse) and again grown until stationary phase under the medium light intensity. During the experimental light phase, community mean cell-size increased with light, but surface:volume ratio and deviation from spherical shape decreased. Moreover, in response to the following nutrient pulse, the low light pre-conditioned communities showed the highest initial growth rates in response to the nutrient pulse. The high light pre-conditioned communities showed the highest conversion of the nutrient pulse into biomass during the stationary phase. These results demonstrate how the imprint of one environmental factor on trait distribution influences the ability to cope with another.

Description: Aurelia aurita (Linneaus, 1758) is a cosmopolitan scyphozoan, probably the most investigated jellyfish in temperate and highly productive coastal ecosystems. Despite a prominent top-down control in plankton food webs, a mechanistic understanding of A. aurita population dynamics and trophic interactions has been barely addressed. Here we develop a food web dynamic model to assess A. aurita role in the seasonal plankton dynamics of the Kiel Fjord, southwestern Baltic Sea. The model couples low trophic level dynamics, based on a classical Nutrient Phytoplankton Zooplankton Detritus (NPZD) model, to a stage-resolved copepod model (referencing Pseudocalanus sp.) and a jellyfish model (A. aurita ephyra and medusa) as consumers and predators, respectively. Simulations showed the relevance of high abundances of A. aurita, which appear related with warm winter temperatures, promoting a shift from a copepod-dominated food web to a ciliate and medusa dominated one. The model captured the intraspecific competition triggered by the medusae abundance and characterized by a negative relationship between population density and individual size/weight. Our results provide a mechanistic understanding of an emergent trait such as size shaping the food web functioning, driving predation rates and population dynamics of A. aurita, driving its sexual reproductive strategy at the end of the pelagic phase.

Description: Atlantic herring (Clupea harengus) is a benthic spawner, therefore its eggs are prone to encounter different water conditions during embryonic development, with bottom waters often depleted of oxygen and enriched in CO2. Some Atlantic herring spawning grounds are predicted to be highly affected by ongoing Ocean Acidification and Warming with water temperature increasing by up to +3°C and CO2 levels reaching ca. 1000 μatm (RCP 8.5). Although many studies investigated the effects of high levels of CO2 on the embryonic development of Atlantic herring, little is known about the combination of temperature and ecologically relevant levels of CO2. In this study, we investigated the effects of Ocean Acidification and Warming on embryonic metabolic and developmental performance such as mitochondrial function, respiration, hatching success (HS) and growth in Atlantic herring from the Oslo Fjord, one of the spawning grounds predicted to be greatly affected by climate change. Fertilized eggs were incubated under combinations of two PCO2 conditions (400 μatm and 1100 μatm) and three temperatures (6, 10 and 14°C), which correspond to current and end-of-the-century conditions. We analysed HS, oxygen consumption (MO2) and mitochondrial function of embryos as well as larval length at hatch. The capacity of the electron transport system (ETS) increased with temperature, reaching a plateau at 14°C, where the contribution of Complex I to the ETS declined in favour of Complex II. This relative shift was coupled with a dramatic increase in MO2 at 14°C. HS was high under ambient spawning conditions (6–10°C), but decreased at 14°C and hatched larvae at this temperature were smaller. Elevated PCO2 increased larval malformations, indicating sub-lethal effects. These results indicate that energetic limitations due to thermally affected mitochondria and higher energy demand for maintenance occur at the expense of embryonic development and growth.

Description: The ocean load in glacial isostatic adjustment (GIA) modelling is represented by the so-called sea level equation (SLE). The SLE describes the mass redistribution of water between ice sheets and oceans on a deforming Earth. Despite various teams independently investigating GIA, there has been no systematic intercomparison among the numerical solvers of the SLE through which the methods may be validated. The goal of this paper is to present a series of synthetic examples designed for testing and comparing the numerical implementations of the SLE in GIA modelling. The 10 numerical codes tested combine various temporal and spatial parametrizations. The time-domain or Laplace-domain discretizations are used to solve the SLE through time, while spherical harmonics, finite differences or finite elements parametrize the GIA-related field variables spatially. The surface ice-water load and solid Earth’s topography are represented spatially either on an equiangular grid, a Gauss–Legendre or an equiarea grid with icosahedron-shaped spherical pixels. Comparisons are made in a series of five benchmark examples with an increasing degree of complexity. Due to the complexity of the SLE, there is no analytical solution to it. The accuracy of the numerical implementations is therefore assessed by the differences of the individual solutions with respect to a reference solution. While the benchmark study does not result in GIA predictions for a realistic loading scenario, we establish a set of agreed-upon results that can be extended in the future by including more complex case studies, such as solutions with realistic loading scenarios, the rotational feedback in the linear-momentum equation, and by considering a 3-D viscosity structure of the Earth’s mantle. The test computations performed so far show very good agreement between the individual results and their ability to capture the main features of sea-surface variation and the surface vertical displacement. The differences found can often be attributed to the different approximations inherent in the various algorithms. This shows the accuracy that can be expected from different implementations of the SLE, which helps to assess differences noted in the literature between predictions for realistic loading cases.