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: 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: 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: 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: 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: During the last ice age, the Northern Hemisphere experienced a series of abrupt millennial-scale climatic changes linked to variations in the strength of the Atlantic Meridional Overturning Circulation and sea-ice extent. However, our understanding of their impacts on decadal-scale climate variability in central Europe has been limited by the lack of high-resolution continental archives. Here, we present a near annual-resolution climate proxy record of central European temperature reconstructed from the Eifel maar lakes of Holzmaar and Auel in Germany, spanning the past 60,000 years. The lake sediments reveal a series of previously undocumented multidecadal climate cycles of around 20 to 150 years that persisted through the last glacial cycle. The periodicity of these cycles suggests that they are related to the Atlantic multidecadal climate oscillations found in the instrumental record and in other climate archives during the Holocene. Our record shows that multidecadal variability in central Europe was strong during all warm interstadials, but was substantially muted during all cold stadial periods. We suggest that this decrease in multidecadal variability was the result of the atmospheric circulation changes associated with the weakening of the Atlantic Meridional Overturning Circulation and the expansion of North Atlantic sea-ice cover during the coldest parts of the last ice age.

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: 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: 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: The Trans-European Suture Zone (TESZ) is the transition zone from the Precambrian East European Craton in the north and east to the younger Phanerozoic mobile belts to the south and west. It is the most prominent lithospheric tectonic feature of Europe. The term Trans-European Suture Zone was only adapted around year 2000 during the Pan-European EUROPROBE program of the European Science Foundation. Until then, parts of the zone were termed Teisseyre-Tornquist Zone, Sorgenfrei-Tornquist Zone, Trans-European Fault, and Tornquist Fan.

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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: 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: Infections arising from multidrug-resistant pathogenic bacteria are spreading rapidly throughout the world and threaten to become untreatable. The origins of resistance are numerous and complex, but one underlying factor is the capacity of bacteria to rapidly export drugs through the intrinsic activity of efflux pumps. In this Review, we describe recent advances that have increased our understanding of the structures and molecular mechanisms of multidrug efflux pumps in bacteria. Clinical and laboratory data indicate that efflux pumps function not only in the drug extrusion process but also in virulence and the adaptive responses that contribute to antimicrobial resistance during infection. The emerging picture of the structure, function and regulation of efflux pumps suggests opportunities for countering their activities.

Description: Microbial communities often exhibit incredible taxonomic diversity, raising questions regarding the mechanisms enabling species coexistence and the role of this diversity in community functioning. On the one hand, many coexisting but taxonomically distinct microorganisms can encode the same energy-yielding metabolic functions, and this functional redundancy contrasts with the expectation that species should occupy distinct metabolic niches. On the other hand, the identity of taxa encoding each function can vary substantially across space or time with little effect on the function, and this taxonomic variability is frequently thought to result from ecological drift between equivalent organisms. Here, we synthesize the powerful paradigm emerging from these two patterns, connecting the roles of function, functional redundancy and taxonomy in microbial systems. We conclude that both patterns are unlikely to be the result of ecological drift, but are inevitable emergent properties of open microbial systems resulting mainly from biotic interactions and environmental and spatial processes.

Description: Naturally produced polybrominated diphenyl ethers (PBDEs) pervade the marine environment and structurally resemble toxic man-made brominated flame retardants. PBDEs bioaccumulate in marine animals and are likely transferred to the human food chain. However, the biogenic basis for PBDE production in one of their most prolific sources, marine sponges of the order Dysideidae, remains unidentified. Here, we report the discovery of PBDE biosynthetic gene clusters within sponge-microbiome-associated cyanobacterial endosymbionts through the use of an unbiased metagenome-mining approach. Using expression of PBDE biosynthetic genes in heterologous cyanobacterial hosts, we correlate the structural diversity of naturally produced PBDEs to modifications within PBDE biosynthetic gene clusters in multiple sponge holobionts. Our results establish the genetic and molecular foundation for the production of PBDEs in one of the most abundant natural sources of these molecules, further setting the stage for a metagenomic-based inventory of other PBDE sources in the marine environment.

Description: Recognition that evolution operates on the same timescale as ecological processes has motivated growing interest in eco-evolutionary dynamics. Nonetheless, generating sufficient data to test predictions about eco-evolutionary dynamics has proved challenging, particularly in natural contexts. Here we argue that genomic data can be integrated into the study of eco-evolutionary dynamics in ways that deepen our understanding of the interplay between ecology and evolution. Specifically, we outline five major questions in the study of eco-evolutionary dynamics for which genomic data may provide answers. Although genomic data alone will not be sufficient to resolve these challenges, integrating genomic data can provide a more mechanistic understanding of the causes of phenotypic change, help elucidate the mechanisms driving eco-evolutionary dynamics, and lead to more accurate evolutionary predictions of eco-evolutionary dynamics in nature.

Description: Although nearly all 2 °C scenarios use negative CO2 emission technologies, only relatively small investments are being made in them, and concerns are being raised regarding their large-scale use. If no explicit policy decisions are taken soon, however, their use will simply be forced on us to meet the Paris climate targets.

Description: Marine microscopic plastic (microplastic) debris is a modern societal issue, illustrating the challenge of balancing the convenience of plastic in daily life with the prospect of causing ecological harm by careless disposal. Here we develop the concept of microplastic as a complex, dynamic mixture of polymers and additives, to which organic material and contaminants can successively bind to form an ‘ecocorona’, increasing the density and surface charge of particles and changing their bioavailability and toxicity. Chronic exposure to microplastic is rarely lethal, but can adversely affect individual animals, reducing feeding and depleting energy stores, with knock-on effects for fecundity and growth. We explore the extent to which ecological processes could be impacted, including altered behaviours, bioturbation and impacts on carbon flux to the deep ocean. We discuss how microplastic compares with other anthropogenic pollutants in terms of ecological risk, and consider the role of science and society in tackling this global issue in the future.

Description: Autophagy is a mechanism by which cellular material is delivered to lysosomes for degradation, leading to the basal turnover of cell components and providing energy and macromolecular precursors. Autophagy has opposing, context-dependent roles in cancer, and interventions to both stimulate and inhibit autophagy have been proposed as cancer therapies. This has led to the therapeutic targeting of autophagy in cancer to be sometimes viewed as controversial. In this Review, we suggest a way forwards for the effective targeting of autophagy by understanding the context-dependent roles of autophagy and by capitalizing on modern approaches to clinical trial design.

Description: Phytoplankton photosynthesis is a critical flux in the carbon cycle, accounting for approximately 40% of the carbon dioxide fixed globally on an annual basis and fuelling the productivity of aquatic food webs. However, rapid evolutionary responses of phytoplankton to warming remain largely unexplored, particularly outside the laboratory, where multiple selection pressures can modify adaptation to environmental change. Here, we use a decade-long experiment in outdoor mesocosms to investigate mechanisms of adaptation to warming (+4 °C above ambient temperature) in the green alga Chlamydomonas reinhardtii, in naturally assembled communities. Isolates from warmed mesocosms had higher optimal growth temperatures than their counterparts from ambient treatments. Consequently, warm-adapted isolates were stronger competitors at elevated temperature and experienced a decline in competitive fitness in ambient conditions, indicating adaptation to local thermal regimes. Higher competitive fitness in the warmed isolates was linked to greater photosynthetic capacity and reduced susceptibility to photoinhibition. These findings suggest that adaptive responses to warming in phytoplankton could help to mitigate projected declines in aquatic net primary production by increasing rates of cellular net photosynthesis.

Description: Nitrogen fixation — the reduction of dinitrogen (N2) gas to biologically available nitrogen (N) — is an important source of N for terrestrial and aquatic ecosystems. In terrestrial environments, N2-fixing symbioses involve multicellular plants, but in the marine environment these symbioses occur with unicellular planktonic algae. An unusual symbiosis between an uncultivated unicellular cyanobacterium (UCYN-A) and a haptophyte picoplankton alga was recently discovered in oligotrophic oceans. UCYN-A has a highly reduced genome, and exchanges fixed N for fixed carbon with its host. This symbiosis bears some resemblance to symbioses found in freshwater ecosystems. UCYN-A shares many core genes with the 'spheroid bodies' of Epithemia turgida and the endosymbionts of the amoeba Paulinella chromatophora. UCYN-A is widely distributed, and has diversified into a number of sublineages that could be ecotypes. Many questions remain regarding the physical and genetic mechanisms of the association, but UCYN-A is an intriguing model for contemplating the evolution of N2-fixing organelles.

Description: Ocean acidification (OA) is increasing due to anthropogenic CO2 emissions and poses a threat to marine species and communities worldwide. To better project the effects of acidification on organisms’ health and persistence, an understanding is needed of the 1) mechanisms underlying developmental and physiological tolerance and 2) potential populations have for rapid evolutionary adaptation. This is especially challenging in nonmodel species where targeted assays of metabolism and stress physiology may not be available or economical for large-scale assessments of genetic constraints. We used mRNA sequencing and a quantitative genetics breeding design to study mechanisms underlying genetic variability and tolerance to decreased seawater pH (-0.4 pH units) in larvae of the sea urchin Strongylocentrotus droebachiensis. We used a gene ontology-based approach to integrate expression profiles into indirect measures of cellular and biochemical traits underlying variation in larval performance (i.e., growth rates). Molecular responses to OA were complex, involving changes to several functions such as growth rates, cell division, metabolism, and immune activities. Surprisingly, the magnitude of pH effects on molecular traits tended to be small relative to variation attributable to segregating functional genetic variation in this species. We discuss how the application of transcriptomics and quantitative genetics approaches across diverse species can enrich our understanding of the biological impacts of climate change.

Description: To the Editor: We congratulate Halsey et al. for their important discussion of the random nature of P values in an earlier issue of this journal1. P values from identical experiments can differ greatly in a way that is surprising to many2, 3, 4, 5, 6. The failure to appreciate this…

Description: To understand how ocean acidification (OA) influences sediment microbial communities, naturally CO2-rich sites are increasingly being used as OA analogues. However, the characterization of these naturally CO2-rich sites is often limited to OA-related variables, neglecting additional environmental variables that may confound OA effects. Here, we used an extensive array of sediment and bottom water parameters to evaluate pH effects on sediment microbial communities at hydrothermal CO2seeps in Papua New Guinea. The geochemical composition of the sediment pore water showed variations in the hydrothermal signature at seep sites with comparable pH, allowing the identification of sites that may better represent future OA scenarios. At these sites, we detected a 60% shift in the microbial community composition compared with reference sites, mostly related to increases inChloroflexisequences. pH was among the factors significantly, yet not mainly, explaining changes in microbial community composition. pH variation may therefore often not be the primary cause of microbial changes when sampling is done along complex environmental gradients. Thus, we recommend an ecosystem approach when assessing OA effects on sediment microbial communities under natural conditions. This will enable a more reliable quantification of OA effects via a reduction of potential confounding effects

Description: Connectivity of pelagic, early life stages via transport by ocean currents may affect survival chances of offspring, recruitment success, and mixing of stocks across management units. Based on drift model studies, transport patterns of particles representing exogenously feeding cod larvae in the transition area between North Sea and Baltic were investigated to (i) determine long-term trends and variability in advective transport of larvae from spawning grounds to juvenile nursery areas, (ii) estimate the degree of exchange between different management areas, and (iii) compare the results with spatial distributions of juvenile cod. The transport of particles showed considerable intra- and interannual variability, but also some general patterns of retention within and dispersion to different management areas. Good spatial overlap of particle end positions, representing potential juvenile settlement areas, with observed distributions of juveniles in bottom trawl surveys suggests that the drift simulations provide reasonable estimates of early life stage connectivity between cod populations in the investigated areas. High exchange rates of particles between management areas of up to ca. 70% suggest that cod populations in the investigated areas are demographically correlated. Results are discussed in relation to their relevance for stock structure, fish stock assessment, and management.

Description: High-throughput techniques based on restriction site-associated DNA sequencing (RADseq) are enabling the low-cost discovery and genotyping of thousands of genetic markers for any species, including non-model organisms, which is revolutionizing ecological, evolutionary and conservation genetics. Technical differences among these methods lead to important considerations for all steps of genomics studies, from the specific scientific questions that can be addressed, and the costs of library preparation and sequencing, to the types of bias and error inherent in the resulting data. In this Review, we provide a comprehensive discussion of RADseq methods to aid researchers in choosing among the many different approaches and avoiding erroneous scientific conclusions from RADseq data, a problem that has plagued other genetic marker types in the past.

Description: While we know much about the evolutionary patterns of endosymbiotic organelle origins, we know less about how the actual process unfolded within each system. This is partly due to the massive changes endosymbiosis appears to trigger, and partly because most organelles evolved in the distant past. The dinotoms are dinoflagellates with diatom endosymbionts, and they represent a relatively recent but nevertheless obligate endosymbiotic association. We have carried out deep sequencing of both the host and endosymbiont transcriptomes from two dinotoms, Durinskia baltica and Glenodinium foliaceum, to examine how the nucleocytosolic compartments have functionally integrated. This analysis showed little or no functional reduction in either the endosymbiont or host, and no evidence for genetic integration. Rather, host and endosymbiont seem to be bound to each other via metabolites, such as photosynthate exported from the endosymbiont to the host as indicated by the presence of plastidic phosphate translocators in the host transcriptome. The host is able to synthesize starch, using plant-specific starch synthases, as a way to store imported photosynthate.

Description: Ubiquitous SAR11 Alphaproteobacteria numerically dominate marine planktonic communities. Because they are excruciatingly difficult to cultivate, there is comparatively little known about their physiology and metabolic responses to long- and short-term environmental changes. As surface oceans take up anthropogenic, atmospheric CO2, the consequential process of ocean acidification could affect the global biogeochemical significance of SAR11. Shipping accidents or inadvertent release of chemicals from industrial plants can have strong short-term local effects on oceanic SAR11. This study investigated the effect of 2.5-fold acidification of seawater on the metabolism of SAR11 and other heterotrophic bacterioplankton along a natural temperature gradient crossing the North Atlantic Ocean, Norwegian and Greenland Seas. Uptake rates of the amino acid leucine by SAR11 cells as well as other bacterioplankton remained similar to controls despite an instant ∼50% increase in leucine bioavailability upon acidification. This high physiological resilience to acidification even without acclimation, suggests that open ocean dominant bacterioplankton are able to cope even with sudden and therefore more likely with long-term acidification effects.

Description: Members of the archaeal phylum Bathyarchaeota are widespread and abundant in the energy-deficient marine subsurface sediments. However, their life strategies have remained largely elusive. Here, we provide genetic evidence that some lineages of Bathyarchaeota are acetogens, being capable of homoacetogenesis, a metabolism so far restricted to the domain Bacteria. Metabolic reconstruction based on genomic bins assembled from the metagenome of deep-sea subsurface sediments shows that the metabolism of some lineages of Bathyarchaeota is similar to that of bona fide bacterial homoacetogens, by having pathways for acetogenesis and for the fermentative utilization of a variety of organic substrates. Heterologous expression and activity assay of the acetate kinase gene ack from Bathyarchaeota, demonstrate further the capability of these Bathyarchaeota to grow as acetogens. The presence and expression of bathyarchaeotal genes indicative of active acetogenesis was also confirmed in Peru Margin subsurface sediments where Bathyarchaeota are abundant. The analyses reveal that this ubiquitous and abundant subsurface archaeal group has adopted a versatile life strategy to make a living under energy-limiting conditions. These findings further expand the metabolic potential of Archaea and argue for a revision of the role of Archaea in the carbon cycle of marine sediments.

Description: Remarkable technological advances have revealed ever more properties and behaviours of individual microorganisms, but the novel data generated by these techniques have not yet been fully exploited. In this Opinion article, we explain how individual-based models (IBMs) can be constructed based on the findings of such techniques and how they help to explore competitive and cooperative microbial interactions. Furthermore, we describe how IBMs have provided insights into self-organized spatial patterns from biofilms to the oceans of the world, phage–CRISPR dynamics and other emergent phenomena. Finally, we discuss how combining individual-based observations with IBMs can advance our understanding at both the individual and population levels, leading to the new approach of microbial individual-based ecology (μIBE).

Description: Bacteria use quorum sensing to orchestrate gene expression programmes that underlie collective behaviours. Quorum sensing relies on the production, release, detection and group-level response to extracellular signalling molecules, which are called autoinducers. Recent work has discovered new autoinducers in Gram-negative bacteria, shown how these molecules are recognized by cognate receptors, revealed new regulatory components that are embedded in canonical signalling circuits and identified novel regulatory network designs. In this Review we examine how, together, these features of quorum sensing signal–response systems combine to control collective behaviours in Gram-negative bacteria and we discuss the implications for host–microbial associations and antibacterial therapy.

Description: Mesoscale ocean eddies are ubiquitous coherent rotating structures of water with radial scales on the order of 100 kilometers. Eddies play a key role in the transport and mixing of momentum and tracers across the World Ocean. We present a global daily mesoscale ocean eddy dataset that contains ~45 million mesoscale features and 3.3 million eddy trajectories that persist at least two days as identified in the AVISO dataset over a period of 1993–2014. This dataset, along with the open-source eddy identification software, extract eddies with any parameters (minimum size, lifetime, etc.), to study global eddy properties and dynamics, and to empirically estimate the impact eddies have on mass or heat transport. Furthermore, our open-source software may be used to identify mesoscale features in model simulations and compare them to observed features. Finally, this dataset can be used to study the interaction between mesoscale ocean eddies and other components of the Earth System.

Description: Microorganisms produce a wealth of structurally diverse specialized metabolites with a remarkable range of biological activities and a wide variety of applications in medicine and agriculture, such as the treatment of infectious diseases and cancer, and the prevention of crop damage. Genomics has revealed that many microorganisms have far greater potential to produce specialized metabolites than was thought from classic bioactivity screens; however, realizing this potential has been hampered by the fact that many specialized metabolite biosynthetic gene clusters (BGCs) are not expressed in laboratory cultures. In this Review, we discuss the strategies that have been developed in bacteria and fungi to identify and induce the expression of such silent BGCs, and we briefly summarize methods for the isolation and structural characterization of their metabolic products.

Description: Cultivable Actinobacteria are the largest source of microbially derived bioactive molecules. The high demand for novel antibiotics highlights the need for exploring novel sources of these bacteria. Microbial symbioses with sessile macro-organisms, known to contain bioactive compounds likely of bacterial origin, represent an interesting and underexplored source of Actinobacteria. We studied the diversity and potential for bioactive-metabolite production of Actinobacteria associated with two marine lichens ( Lichina confinis and L. pygmaea ; from intertidal and subtidal zones) and one littoral lichen ( Roccella fuciformis ; from supratidal zone) from the Brittany coast (France), as well as the terrestrial lichen Collema auriforme (from a riparian zone, Austria). A total of 247 bacterial strains were isolated using two selective media. Isolates were identified and clustered into 101 OTUs (98% identity) including 51 actinobacterial OTUs. The actinobacterial families observed were: Brevibacteriaceae, Cellulomonadaceae, Gordoniaceae, Micrococcaceae, Mycobacteriaceae, Nocardioidaceae, Promicromonosporaceae, Pseudonocardiaceae, Sanguibacteraceae and Streptomycetaceae. Interestingly, the diversity was most influenced by the selective media rather than lichen species or the level of lichen thallus association. The potential for bioactive-metabolite biosynthesis of the isolates was confirmed by screening genes coding for polyketide synthases types I and II. These results show that littoral lichens are a source of diverse potentially bioactive Actinobacteria.

Description: Barnacles of the genus Galkinius occupy a large spectrum of host corals, making it one of the least host-specific genera within the Pyrgomatidae. Molecular analyses show that within the genus Galkinius there are highly supported clades, suggesting that the genus Galkinius is a complex of evolutionarily significant units (ESUs). The morphology of the opercular valves has been used as the basis for the separation of species of Galkinius. In this study, morphological variability was found both between specimens within ESUs extracted from different host species and between specimens extracted from the same colony. Identifications based on the opercular valves cannot therefore be assigned to different species despite being genetically distinguishable. It is proposed that in many cases the differences between valve morphology of different species of Galkinius are the outcome of ontogeny. Allometric growth of the valves has resulted in differences in the proportions of the parts of the valve.

Description: Large amounts of the greenhouse gas methane are released from the seabed to the water column1, where it may be consumed by aerobic methanotrophic bacteria2. The size and activity of methanotrophic communities, which determine the amount of methane consumed in the water column, are thought to be mainly controlled by nutrient and redox dynamics3–7. Here, we report repeated measurements of methanotrophic activity and community size at methane seeps west of Svalbard, and relate them to physical water mass properties and modelled ocean currents. We show that cold bottom water, which contained a large number of aerobic methanotrophs, was displaced by warmer water with a considerably smaller methanotrophic community within days. Ocean current simulations using a global ocean/sea-ice model suggest that this water mass exchange is consistent with short-term variations in the meandering West Spitsbergen Current. We conclude that the shift from an offshore to a nearshore position of the current can rapidly and severely reduce methanotrophic activity in the water column. Strong fluctuating currents are common at many methane seep systems globally, and we suggest that they affect methane oxidation in the water column at other sites, too.

Description: Climate change is altering oceanic conditions in a complex manner, and the concurrent amendment of multiple properties will modify environmental stress for primary producers. So far, global modelling studies have focused largely on how alteration of individual properties will affect marine life. Here, we use global modelling simulations in conjunction with rotated factor analysis to express model projections in terms of regional trends in concomitant changes to biologically influential multi-stressors. Factor analysis demonstrates that regionally distinct patterns of complex oceanic change are evident globally. Preliminary regional assessments using published evidence of phytoplankton responses to complex change reveal a wide range of future responses to interactive multi-stressors with 〈20–300% shifts in phytoplankton physiological rates, and many unexplored potential interactions. In a future ocean, provinces will encounter different permutations of change that will probably alter the dominance of key phytoplankton groups and modify regional productivity, ecosystem structure and biogeochemistry. Consideration of regionally distinct multi-stressor patterns can help guide laboratory and field studies as well as the interpretation of interactive multi-stressors in global models.

Description: Reflection on the rapidly growing field of ocean acidification research highlights priorities for future research on the changing ocean.

Description: Thin oceanic crust is formed by decompression melting of the upper mantle at mid-ocean ridges, but the origin of the thick and buoyant continental crust is enigmatic. Juvenile continental crust may form from magmas erupted above intraoceanic subduction zones, where oceanic lithosphere subducts beneath other oceanic lithosphere. However, it is unclear why the subduction of dominantly basaltic oceanic crust would result in the formation of andesitic continental crust at the surface. Here we use geochemical and geophysical data to reconstruct the evolution of the Central American land bridge, which formed above an intra-oceanic subduction system over the past 70Myr. We find that the geochemical signature of erupted lavas evolved from basaltic to andesitic about 10Myr ago - coincident with the onset of subduction of more oceanic crust that originally formed above the Galápagos mantle plume. We also find that seismic P-waves travel through the crust at velocities intermediate between those typically observed for oceanic and continental crust. We develop a continentality index to quantitatively correlate geochemical composition with the average P-wave velocity of arc crust globally. We conclude that although the formation and evolution of continents may involve many processes, melting enriched oceanic crust within a subduction zone - a process probably more common in the Archaean - can produce juvenile continental crust.

Description: In this pilot study, we describe a high-pressure incubation system allowing multiple subsampling of a pressurized culture without decompression. The system was tested using one piezophilic (Photobacterium profundum), one piezotolerant (Colwellia maris) bacterial strain and a decompressed sample from the Mediterranean deep sea (3044 m) determining bacterial community composition, protein production (BPP) and cell multiplication rates (BCM) up to 27 MPa. The results showed elevation of BPP at high pressure was by a factor of 1.5 ± 1.4 and 3.9 ± 2.3 for P. profundum and C. maris, respectively, compared to ambient-pressure treatments and by a factor of 6.9 ± 3.8 fold in the field samples. In P. profundum and C. maris, BCM at high pressure was elevated (3.1 ± 1.5 and 2.9 ± 1.7 fold, respectively) compared to the ambient-pressure treatments. After 3 days of incubation at 27 MPa, the natural bacterial deep-sea community was dominated by one phylum of the genus Exiguobacterium, indicating the rapid selection of piezotolerant bacteria. In future studies, our novel incubation system could be part of an isopiestic pressure chain, allowing more accurate measurement of bacterial activity rates which is important both for modeling and for predicting the efficiency of the oceanic carbon pump.

Description: The stomachs of most vertebrates operate at an acidic pH of 2 generated by the gastric H+/K+-ATPase located in parietal cells. The acidic pH in stomachs of vertebrates is believed to aid digestion and to protect against environmental pathogens. Little attention has been placed on whether acidic gastric pH regulation is a vertebrate character or a deuterostome ancestral trait. Here, we report alkaline conditions up to pH 10.5 in the larval digestive systems of ambulacraria (echinoderm + hemichordate), the closest relative of the chordate. Microelectrode measurements in combination with specific inhibitors for acid-base transporters and ion pumps demonstrated that the gastric alkalization machinery in sea urchin larvae is mainly based on direct H+ secretion from the stomach lumen and involves a conserved set of ion pumps and transporters. Hemichordate larvae additionally utilized HCO 3- transport pathways to generate even more alkaline digestive conditions. Molecular analyses in combination with acidification experiments supported these findings and identified genes coding for ion pumps energizing gastric alkalization. Given that insect larval guts were also reported to be alkaline, our discovery raises the hypothesis that the bilaterian ancestor utilized alkaline digestive system while the vertebrate lineage has evolved a strategy to strongly acidify their stomachs.

Description: Symbiotic relationships between phytoplankton and N-2-fixing microorganisms play a crucial role in marine ecosystems. The abundant and widespread unicellular cyanobacteria group A (UCYN-A) has recently been found to live symbiotically with a haptophyte. Here, we investigated the effect of nitrogen (N), phosphorus (P), iron (Fe) and Saharan dust additions on nitrogen (N-2) fixation and primary production by the UCYN-A-haptophyte association in the subtropical eastern North Atlantic Ocean using nifH expression analysis and stable isotope incubations combined with single-cell measurements. N-2 fixation by UCYN-A was stimulated by the addition of Fe and Saharan dust, although this was not reflected in the nifH expression. CO2 fixation by the haptophyte was stimulated by the addition of ammonium nitrate as well as Fe and Saharan dust. Intriguingly, the single-cell analysis using nanometer scale secondary ion mass spectrometry indicates that the increased CO2 fixation by the haptophyte in treatments without added fixed N is likely an indirect result of the positive effect of Fe and/or P on UCYN-A N-2 fixation and the transfer of N-2-derived N to the haptophyte. Our results reveal a direct linkage between the marine carbon and nitrogen cycles that is fuelled by the atmospheric deposition of dust. The comparison of single-cell rates suggests a tight coupling of nitrogen and carbon transfer that stays balanced even under changing nutrient regimes. However, it appears that the transfer of carbon from the haptophyte to UCYN-A requires a transfer of nitrogen from UCYN-A. This tight coupling indicates an obligate symbiosis of this globally important diazotrophic association.

Description: Quasi-decadal variability in solar irradiance has been suggested to exert a substantial effect on Earth’s regional climate. In the North Atlantic sector, the 11-year solar signal has been proposed to project onto a pattern resembling the North Atlantic Oscillation (NAO), with a lag of a few years due to ocean-atmosphere interactions. The solar/NAO relationship is, however, highly misrepresented in climate model simulations with realistic observed forcings. In addition, its detection is particularly complicated since NAO quasi-decadal fluctuations can be intrinsically generated by the coupled ocean-atmosphere system. Here we compare two multi-decadal ocean-atmosphere chemistry-climate simulations with and without solar forcing variability. While the experiment including solar variability simulates a 1–2-year lagged solar/NAO relationship, comparison of both experiments suggests that the 11-year solar cycle synchronizes quasi-decadal NAO variability intrinsic to the model. The synchronization is consistent with the downward propagation of the solar signal from the stratosphere to the surface.

Description: Archaea associated with marine sponges are active and influence the nitrogen metabolism of sponges. However, we know little about their occurrence, specificity, and persistence. We aimed to elucidate the relative importance of host specificity and biogeographic background in shaping the symbiotic archaeal communities. We investigated these communities in sympatric sponges from the Mediterranean (Ircinia fasciculata and Ircinia oros, sampled in summer and winter) and from the Caribbean (Ircinia strobilina and Mycale laxissima). PCR cloning and sequencing of archaeal 16S rRNA and amoA genes showed that the archaeal community composition and structure were different from that in seawater and varied among sponge species. We found that the communities were dominated by ammonia-oxidizing archaea closely related to Nitrosopumilus. The community in M. laxissima differed from that in Ircinia spp., including the sympatric sponge I. strobilina; yet, geographical clusters within Ircinia spp. were observed. Whereas archaeal phylotypes in Ircinia spp. were persistent and belong to 'sponge-enriched' clusters, archaea in M. laxissima were closely related with those from diverse habitats (i.e. seawater and sediments). For all four sponge species, the expression of the archaeal amoA gene was confirmed. Our results indicate that host-specific processes, such as host ecological strategy and evolutionary history, control the sponge-archaeal communities.

Description: Dissolved oxygen in sea water affects marine habitats and biogeochemical cycles. Oceanic zones with oxygen deficits represent 7% of the volume and 8% of the area of the oceans, and are thought to be expanding. One of the most pronounced lies in the region off Peru, where mesoscale activity in the form of fronts and eddies is strong. Here, we study the dynamics of the Peruvian oxygen minimum zone in a Lagrangian framework, using a coupled physical-biogeochemical numerical model and finite-size Lyapunov exponent fields, to evaluate the role of mesoscale activity. We find that, at depths between 380 and 600 m, mesoscale structures have two distinct roles. First, their mean positions and paths delimit and maintain the oxygen minimum zone boundaries. Second, their high-frequency fluctuations inject oxygen across the oxygen minimum zone boundaries and eddy fluxes are one order of magnitude higher than mean oxygen fluxes. We conclude that these eddy fluxes contribute to the ventilation of the Peruvian oxygen minimum zone.

Description: The symbiotic association of corals and unicellular algae of the genus Symbiodinium in the southern Persian/Arabian Gulf (PAG) display an exceptional heat tolerance, enduring summer peak temperatures of up to 36 °C. As yet, it is not clear whether this resilience is related to the presence of specific symbiont types that are exclusively found in this region. Therefore, we used molecular markers to identify the symbiotic algae of three Porites species along 〉1000 km of coastline in the PAG and the Gulf of Oman and found that a recently described species, Symbiodinium thermophilum, is integral to coral survival in the southern PAG, the world’s hottest sea. Despite the geographic isolation of the PAG, we discovered that representatives of the S. thermophilum group can also be found in the adjacent Gulf of Oman providing a potential source of thermotolerant symbionts that might facilitate the adaptation of Indian Ocean populations to the higher water temperatures expected for the future. However, corals from the PAG associated with S. thermophilum show strong local adaptation not only to high temperatures but also to the exceptionally high salinity of their habitat. We show that their superior heat tolerance can be lost when these corals are exposed to reduced salinity levels common for oceanic environments elsewhere. Consequently, the salinity prevailing in most reefs outside the PAG might represent a distribution barrier for extreme temperature-tolerant coral/Symbiodinium associations from the PAG.

Description: Global warming was reported to cause growth reductions in tropical shallow water corals in both, cooler and warmer, regions of the coral species range. This suggests regional adaptation with less heat-tolerant populations in cooler and more thermo-tolerant populations in warmer regions. Here, we investigated seasonal changes in the in situ metabolic performance of the widely distributed hermatypic coral Pocillopora verrucosa along 12 degrees latitudes featuring a steep temperature gradient between the northern (28.5 degrees N, 21-27 degrees C) and southern (16.5 degrees N, 28-33 degrees C) reaches of the Red Sea. Surprisingly, we found little indication for regional adaptation, but strong indications for high phenotypic plasticity: Calcification rates in two seasons (winter, summer) were found to be highest at 28-29 degrees C throughout all populations independent of their geographic location. Mucus release increased with temperature and nutrient supply, both being highest in the south. Genetic characterization of the coral host revealed low inter-regional variation and differences in the Symbiodinium clade composition only at the most northern and most southern region. This suggests variable acclimatization potential to ocean warming of coral populations across the Red Sea: high acclimatization potential in northern populations, but limited ability to cope with ocean warming in southern populations already existing at the upper thermal margin for corals

Description: The ability of microorganisms to withstand long periods with extremely low energy input has gained increasing scientific attention in recent years. Starvation experiments in the laboratory have shown that a phylogenetically wide range of microorganisms evolve fitness-enhancing genetic traits within weeks of incubation under low-energy stress. Studies on natural environments that are cut off from new energy supplies over geologic time scales, such as deeply buried sediments, suggest that similar adaptations might mediate survival under energy limitation in the environment. Yet, the extent to which laboratory-based evidence of starvation survival in pure or mixed cultures can be extrapolated to sustained microbial ecosystems in nature remains unclear. In this review, we discuss past investigations on microbial energy requirements and adaptations to energy limitation, identify gaps in our current knowledge, and outline possible future foci of research on life under extreme energy limitation.

Description: Over the past 60 years, both average daily precipitation intensity and extreme precipitation have increased in many regions1, 2, 3. Part of these changes, or even individual events4, 5, have been attributed to anthropogenic warming6, 7. Over the Black Sea and Mediterranean region, the potential for extreme summertime convective precipitation has grown8 alongside substantial sea surface temperature increase. A particularly devastating convective event experienced in that region was the July 2012 precipitation extreme near the Black Sea town of Krymsk9. Here we study the effect of sea surface temperature (SST) increase on convective extremes within the region, taking the Krymsk event as a showcase example. We carry out ensemble sensitivity simulations with a convection-permitting atmospheric model and show the crucial role of SST increase in the extremeness of the event. The enhancement of lower tropospheric instability due to the current warmer Black Sea allows deep convection to be triggered, increasing simulated precipitation by more than 300% relative to simulations with SSTs characteristic of the early 1980s. A highly nonlinear precipitation response to incremental SST increase suggests that the Black Sea has exceeded a regional threshold for the intensification of convective extremes. The physical mechanism we identify indicates that Black Sea and Mediterranean coastal regions may face abrupt amplifications of convective precipitation under continued SST increase, and illustrates the limitations of thermodynamical bounds for estimating the temperature scaling of convective extremes.

Description: To date, the parameters that determine the rupture area of great subduction zone earthquakes remain contentious. On 1 April 2014, the Mw 8.1 Iquique earthquake ruptured a portion of the well-recognized northern Chile seismic gap but left large highly coupled areas un-ruptured. Marine seismic reflection and swath bathymetric data indicate that structural variations in the subducting Nazca Plate control regional-scale plate-coupling variations, and the limited extent of the 2014 earthquake. Several under-thrusting seamounts correlate to the southward and up-dip arrest of seismic rupture during the 2014 Iquique earthquake, thus supporting a causal link. By fracturing of the overriding plate, the subducting seamounts are likely further responsible for reduced plate-coupling in the shallow subduction zone and in a lowly coupled region around 20.5°S. Our data support that structural variations in the lower plate influence coupling and seismic rupture offshore Northern Chile, whereas the structure of the upper plate plays a minor role.

Description: Oceanic dissolved organic carbon (DOC) is an important carbon pool, similar in magnitude to atmospheric CO2, but the fate of its oldest forms is not well understood. Hot hydrothermal circulation may facilitate the degradation of otherwise un-reactive dissolved organic matter, playing an important role in the long-term global carbon cycle. The oldest, most recalcitrant forms of DOC, which make up most of oceanic DOC, can be recovered by solid-phase extraction. Here we present measurements of solid-phase extractable DOC from samples collected between 2009 and 2013 at seven vent sites in the Atlantic, Pacific and Southern oceans, along with magnesium concentrations, a conservative tracer of water circulation through hydrothermal systems. We find that magnesium and solid-phase extractable DOC concentrations are correlated, suggesting that solid-phase extractable DOC is almost entirely lost from solution through mineralization or deposition during circulation through hydrothermal vents with fluid temperatures of 212-401°C. In laboratory experiments, where we heated samples to 380°C for four days, we found a similar removal efficiency. We conclude that thermal degradation alone can account for the loss of solid-phase extractable DOC in natural hydrothermal systems, and that its maximum lifetime is constrained by the timescale of hydrothermal cycling, at about 40 million years

Description: Cold-water corals (CWC) are widely distributed around the world forming extensive reefs at par with tropical coral reefs. They are hotspots of biodiversity and organic matter processing in the world’s deep oceans. Living in the dark they lack photosynthetic symbionts and are therefore considered to depend entirely on the limited flux of organic resources from the surface ocean. While symbiotic relations in tropical corals are known to be key to their survival in oligotrophic conditions, the full metabolic capacity of CWC has yet to be revealed. Here we report isotope tracer evidence for efficient nitrogen recycling, including nitrogen assimilation, regeneration, nitrification and denitrification. Moreover, we also discovered chemoautotrophy and nitrogen fixation in CWC and transfer of fixed nitrogen and inorganic carbon into bulk coral tissue and tissue compounds (fatty acids and amino acids). This unrecognized yet versatile metabolic machinery of CWC conserves precious limiting resources and provides access to new nitrogen and organic carbon resources that may be essential for CWC to survive in the resource-depleted dark ocean.

Description: The interoceanic transfer of seawater between the Indian Ocean and the Atlantic, ‘Agulhas leakage’, forms a choke point for the overturning circulation in the global ocean. Here, by combining output from a series of high-resolution ocean and climate models with in situ and satellite observations, we construct a time series of Agulhas leakage for the period 1870–2014. The time series demonstrates the impact of Southern Hemisphere westerlies on decadal timescales. Agulhas leakage shows a correlation with the Atlantic Multi-decadal Oscillation on multi-decadal timescales; the former leading by 15 years. This is relevant for climate in the North Atlantic

Description: Sea-level rise1 is one of the most pressing aspects of anthropogenic global warming with far-reaching consequences for coastal societies. However, sea-level rise did2, 3, 4, 5, 6, 7 and will strongly vary from coast to coast8, 9, 10. Here we investigate the long-term internal variability effects on centennial projections of dynamic sea level (DSL), the local departure from the globally averaged sea level. A large ensemble of global warming integrations has been conducted with a climate model, where each realization was forced by identical CO2 increase but started from different atmospheric and oceanic initial conditions. In large parts of the mid- and high latitudes, the ensemble spread of the projected centennial DSL trends is of the same order of magnitude as the globally averaged steric sea-level rise, suggesting that internal variability cannot be ignored when assessing twenty-first-century DSL trends. The ensemble spread is considerably reduced in the mid- to high latitudes when only the atmospheric initial conditions differ while keeping the oceanic initial state identical; indicating that centennial DSL projections are strongly dependent on ocean initial conditions.

Description: Restructure data-gathering and evaluation networks to address climate change, energy, food, health and water provision, say Yonglong Lu and colleagues.

Description: The permafrost methane emission problem is the focus of attention on different climate models. Here, we present a mathematical model for permafrost lake methane emission and its influence on the climate system. We model this process using the theory of non-linear phase transitions. Further, we find that a climate catastrophe possibility depends on a value of feedback connecting the methane concentration in the atmosphere and temperature, and on the tundra permafrost methane pool.We note that the permafrost lake model that we developed for the methane emission positive feedback loop problem is a conceptual climate model.

Description: Global mean surface warming has stalled since the end of the twentieth century1, 2, but the net radiation imbalance at the top of the atmosphere continues to suggest an increasingly warming planet. This apparent contradiction has been reconciled by an anomalous heat flux into the ocean3, 4, 5, 6, 7, 8, induced by a shift towards a La Niña-like state with cold sea surface temperatures in the eastern tropical Pacific over the past decade or so. A significant portion of the heat missing from the atmosphere is therefore expected to be stored in the Pacific Ocean. However, in situ hydrographic records indicate that Pacific Ocean heat content has been decreasing9. Here, we analyse observations along with simulations from a global ocean–sea ice model to track the pathway of heat. We find that the enhanced heat uptake by the Pacific Ocean has been compensated by an increased heat transport from the Pacific Ocean to the Indian Ocean, carried by the Indonesian throughflow. As a result, Indian Ocean heat content has increased abruptly, which accounts for more than 70% of the global ocean heat gain in the upper 700 m during the past decade. We conclude that the Indian Ocean has become increasingly important in modulating global climate variability.

Description: An index of water-circulation strength in the North Atlantic Ocean has been derived from sea-level measurements. This provides fresh evidence of the ocean's leading role in multidecadal climate variability. See Letter p.508

Description: Increasingly, spatial geochemical zonation, present as geographically distinct, subparallel trends, is observed along hotspot tracks, such as Hawaii and the Galapagos. The origin of this zonation is currently unclear. Recently zonation was found along the last B70 Myr of the Tristan-Gough hotspot track. Here we present new Sr–Nd–Pb–Hf isotope data from the older parts of this hotspot track (Walvis Ridge and Rio Grande Rise) and re-evaluate published data from the Etendeka and Parana flood basalts erupted at the initiation of the hotspot track. We show that only the enriched Gough, but not the less-enriched Tristan, component is present in the earlier (70–132 Ma) history of the hotspot. Here we present a model that can explain the temporal evolution and origin of plume zonation for both the Tristan-Gough and Hawaiian hotspots, two end member types of zoned plumes, through processes taking place in the plume sources at the base of the lower mantle.