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: 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: Geological archives record multiple reversals of Earth’s magnetic poles, but the global impacts of these events, if any, remain unclear. Uncertain radiocarbon calibration has limited investigation of the potential effects of the last major magnetic inversion, known as the Laschamps Excursion [41 to 42 thousand years ago (ka)]. We use ancient New Zealand kauri trees (Agathis australis) to develop a detailed record of atmospheric radiocarbon levels across the Laschamps Excursion. We precisely characterize the geomagnetic reversal and perform global chemistry-climate modeling and detailed radiocarbon dating of paleoenvironmental records to investigate impacts. We find that geomagnetic field minima ~42 ka, in combination with Grand Solar Minima, caused substantial changes in atmospheric ozone concentration and circulation, driving synchronous global climate shifts that caused major environmental changes, extinction events, and transformations in the archaeological record.

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: 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: The last extended time period when climate may have been warmer than today was during the Last Interglacial (LIG; ca. 129 to 120 thousand years ago). However, a global view of LIG precipitation is lacking. Here, seven new LIG climate models are compared to the first global database of proxies for LIG precipitation. In this way, models are assessed in their ability to capture important hydroclimatic processes during a different climate. The models can reproduce the proxy-based positive precipitation anomalies from the preindustrial period over much of the boreal continents. Over the Southern Hemisphere, proxy-model agreement is partial. In models, LIG boreal monsoons have 42% wider area than in the preindustrial and produce 55% more precipitation and 50% more extreme precipitation. Austral monsoons are weaker. The mechanisms behind these changes are consistent with stronger summer radiative forcing over boreal high latitudes and with the associated higher temperatures during the LIG.

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: Variations in Earth’s orbit pace the glacial-interglacial cycles of the Quaternary, but the mechanisms that transform regional and seasonal variations in solar insolation into glacial-interglacial cycles are still elusive. Here, we present transient simulations of coevolution of climate, ice sheets, and carbon cycle over the past 3 million years. We show that a gradual lowering of atmospheric CO2 and regolith removal are essential to reproduce the evolution of climate variability over the Quaternary. The long-term CO2 decrease leads to the initiation of Northern Hemisphere glaciation and an increase in the amplitude of glacial-interglacial variations, while the combined effect of CO2 decline and regolith removal controls the timing of the transition from a 41,000- to 100,000-year world. Our results suggest that the current CO2 concentration is unprecedented over the past 3 million years and that global temperature never exceeded the preindustrial value by more than 2°C during the Quaternary.

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: To provide an observational basis for the Intergovernmental Panel on Climate Change projections of a slowing Atlantic meridional overturning circulation (MOC) in the 21st century, the Overturning in the Subpolar North Atlantic Program (OSNAP) observing system was launched in the summer of 2014. The first 21-month record reveals a highly variable overturning circulation responsible for the majority of the heat and freshwater transport across the OSNAP line. In a departure from the prevailing view that changes in deep water formation in the Labrador Sea dominate MOC variability, these results suggest that the conversion of warm, salty, shallow Atlantic waters into colder, fresher, deep waters that move southward in the Irminger and Iceland basins is largely responsible for overturning and its variability in the subpolar basin.

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: We quantify the oceanic sink for anthropogenic carbon dioxide (CO 2 ) over the period 1994 to 2007 by using observations from the global repeat hydrography program and contrasting them to observations from the 1990s. Using a linear regression–based method, we find a global increase in the anthropogenic CO 2 inventory of 34 ± 4 petagrams of carbon (Pg C) between 1994 and 2007. This is equivalent to an average uptake rate of 2.6 ± 0.3 Pg C year −1 and represents 31 ± 4% of the global anthropogenic CO 2 emissions over this period. Although this global ocean sink estimate is consistent with the expectation of the ocean uptake having increased in proportion to the rise in atmospheric CO 2 , substantial regional differences in storage rate are found, likely owing to climate variability–driven changes in ocean circulation.

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: As plastic waste pollutes the oceans and fish stocks decline, unseen below the surface another problem grows: deoxygenation. Breitburg et al. review the evidence for the downward trajectory of oxygen levels in increasing areas of the open ocean and coastal waters. Rising nutrient loads coupled with climate change—each resulting from human activities—are changing ocean biogeochemistry and increasing oxygen consumption. This results in destabilization of sediments and fundamental shifts in the availability of key nutrients. In the short term, some compensatory effects may result in improvements in local fisheries, such as in cases where stocks are squeezed between the surface and elevated oxygen minimum zones. In the longer term, these conditions are unsustainable and may result in ecosystem collapses, which ultimately will cause societal and economic harm.

Description: Stratification of the deep Southern Ocean during the Last Glacial Maximum is thought to have facilitated carbon storage and subsequent release during the deglaciation as stratification broke down, contributing to atmospheric CO2 rise. Here, we present neodymium isotope evidence from deep to abyssal waters in the South Pacific that confirms stratification of the deepwater column during the Last Glacial Maximum. The results indicate a glacial northward expansion of Ross Sea Bottom Water and a Southern Hemisphere climate trigger for the deglacial breakup of deep stratification. It highlights the important role of abyssal waters in sustaining a deep glacial carbon reservoir and Southern Hemisphere climate change as a prerequisite for the destabilization of the water column and hence the deglacial release of sequestered CO2 through upwelling.

Description: Marine protected areas (MPAs) are increasingly used as a primary tool to conserve biodiversity. This is particularly relevant in heavily exploited fisheries hot spots such as Europe, where MPAs now cover 29% of territorial waters, with unknown effects on fishing pressure and conservation outcomes. We investigated industrial trawl fishing and sensitive indicator species in and around 727 MPAs designated by the European Union. We found that 59% of MPAs are commercially trawled, and average trawling intensity across MPAs is at least 1.4-fold higher as compared with nonprotected areas. Abundance of sensitive species (sharks, rays, and skates) decreased by 69% in heavily trawled areas. The widespread industrial exploitation of MPAs undermines global biodiversity conservation targets, elevating recent concerns about growing human pressures on protected areas worldwide.

Description: Hoffmann et al. (Reports, 23 February 2018, p. 912) report the discovery of parietal art older than 64,800 years and attributed to Neanderthals, at least 25 millennia before the oldest parietal art ever found. Instead, critical evaluation of their geochronological data seems to provide stronger support for an age of 47,000 years, which is much more consistent with the archaeological background in hand.

Description: Plate-boundary fault rupture during the 2004 Sumatra-Andaman subduction earthquake extended closer to the trench than expected, increasing earthquake and tsunami size. International Ocean Discovery Program Expedition 362 sampled incoming sediments offshore northern Sumatra, revealing recent release of fresh water within the deep sediments. Thermal modeling links this freshening to amorphous silica dehydration driven by rapid burial-induced temperature increases in the past 9 million years. Complete dehydration of silicates is expected before plate subduction, contrasting with prevailing models for subduction seismogenesis calling for fluid production during subduction. Shallow slip offshore Sumatra appears driven by diagenetic strengthening of deeply buried fault-forming sediments, contrasting with weakening proposed for the shallow Tohoku-Oki 2011 rupture, but our results are applicable to other thickly sedimented subduction zones including those with limited earthquake records.

Description: Arctic sea-ice loss is a leading indicator of climate change and can be attributed, in large part, to atmospheric forcing. Here, we show that recent ice reductions, weakening of the halocline, and shoaling of intermediate-depth Atlantic Water layer in the eastern Eurasian Basin have increased winter ventilation in the ocean interior, making this region structurally similar to that of the western Eurasian Basin. The associated enhanced release of oceanic heat has reduced winter sea-ice formation at a rate now comparable to losses from atmospheric thermodynamic forcing, thus explaining the recent reduction in sea-ice cover in the eastern Eurasian Basin. This encroaching “atlantification” of the Eurasian Basin represents an essential step toward a new Arctic climate state, with a substantially greater role for Atlantic inflows.

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 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: 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: 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: 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: 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: 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: Changes in the formation of dense water in the Arctic Ocean and Nordic Seas [the “Arctic Mediterranean” (AM)] probably contributed to the altered climate of the last glacial period. We examined past changes in AM circulation by reconstructing radiocarbon ventilation ages of the deep Nordic Seas over the past 30,000 years. Our results show that the glacial deep AM was extremely poorly ventilated (ventilation ages of up to 10,000 years). Subsequent episodic overflow of aged water into the mid-depth North Atlantic occurred during deglaciation. Proxy data also suggest that the deep glacial AM was ~2° to 3°C warmer than modern temperatures; deglacial mixing of the deep AM with the upper ocean thus potentially contributed to the melting of sea ice, icebergs, and terminal ice-sheet margins.

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: The profound influence of marine plankton on the global carbon cycle has been recognized for decades, particularly for photosynthetic microbes that form the base of ocean food chains. However, a comprehensive model of the carbon cycle is challenged by unicellular eukaryotes (protists) having evolved complex behavioral strategies and organismal interactions that extend far beyond photosynthetic lifestyles. As is also true for multicellular eukaryotes, these strategies and their associated physiological changes are difficult to deduce from genome sequences or gene repertoires - a problem compounded by numerous unknown function proteins. Here, we explore protistan trophic modes in marine food webs and broader biogeochemical influences. We also evaluate approaches that could resolve their activities, link them to biotic and abiotic factors, and integrate them into an ecosystems biology framework. © 2015, American Association for the Advancement of Science. All rights reserved.

Description: Spliceosomal introns are a hallmark of eukaryotic genes that are hypothesized to play important roles in genome evolution but have poorly understood origins. Although most introns lack sequence homology to each other, new families of spliceosomal introns that are repeated hundreds of times in individual genomes have recently been discovered in a few organisms. The prevalence and conservation of these introner elements (IEs) or introner-like elements in other taxa, as well as their evolutionary relationships to regular spliceosomal introns, are still unknown. Here, we systematically investigate introns in the widespread marine green alga Micromonas and report new families of IEs, numerous intron presence-absence polymorphisms, and potential intron insertion hot-spots. The new families enabled identification of conserved IE secondary structure features and establishment of a novel general model for repetitive intron proliferation across genomes. Despite shared secondary structure, the IE families from each Micromonas lineage bear no obvious sequence similarity to those in the other lineages, suggesting that their appearance is intimately linked with the process of speciation. Two of the new IE families come from an Arctic culture (Micromonas Clade E2) isolated from a polar region where abundance of this alga is increasing due to climate induced changes. The same two families were detected in metagenomic data from Antarctica-a system where Micromonas has never before been reported. Strikingly high identity between the Arctic isolate and Antarctic coding sequences that flank the IEs suggests connectivity between populations in the two polar systems that we postulate occurs through deep-sea currents. Recovery of Clade E2 sequences in North Atlantic Deep Waters beneath the Gulf Stream supports this hypothesis. Our research illuminates the dynamic relationships between an unusual class of repetitive introns, genome evolution, speciation, and global distribution of this sentinel marine alga. © 2015 The Author.

Description: The GeoChip 4.2 gene array was employed to interrogate the microbial functional gene repertoire of sponges and seawater collected from the Red Sea and the Mediterranean. Complementary amplicon sequencing confirmed the microbial community composition characteristic of high microbial abundance (HMA) and low microbial abundance (LMA) sponges. By use of GeoChip, altogether 20 273 probes encoding for 627 functional genes and representing 16 gene categories were identified. Minimum curvilinear embedding analyses revealed a clear separation between the samples. The HMA/LMA dichotomy was stronger than any possible geographic pattern, which is shown here for the first time on the level of functional genes. However, upon inspection of individual genes, very few specific differences were discernible. Differences were related to microbial ammonia oxidation, ammonification, and archaeal autotrophic carbon fixation (higher gene abundance in sponges over seawater) as well as denitrification and radiation-stress-related genes (lower gene abundance in sponges over seawater). Except for few documented specific differences the functional gene repertoire between the different sources appeared largely similar. This study expands previous reports in that functional gene convergence is not only reported between HMA and LMA sponges but also between sponges and seawater.

Description: Ferredoxins are electron shuttles harbouring iron–sulfur clusters that connect multiple oxido-reductive pathways in organisms displaying different lifestyles. Some prokaryotes and algae express an isofunctional electron carrier, flavodoxin, which contains flavin mononucleotide as cofactor. Both proteins evolved in the anaerobic environment preceding the appearance of oxygenic photosynthesis. The advent of an oxygen-rich atmosphere proved detrimental to ferredoxin owing to iron limitation and oxidative damage to the iron–sulfur cluster, and many microorganisms induced flavodoxin expression to replace ferredoxin under stress conditions. Paradoxically, ferredoxin was maintained throughout the tree of life, whereas flavodoxin is absent from plants and animals. Of note is that flavodoxin expression in transgenic plants results in increased tolerance to multiple stresses and iron deficit, through mechanisms similar to those operating in microorganisms. Then, the question remains open as to why a trait that still confers plants such obvious adaptive benefits was not retained. We compare herein the properties of ferredoxin and flavodoxin, and their contrasting modes of expression in response to different environmental stimuli. Phylogenetic analyses suggest that the flavodoxin gene was already absent in the algal lineages immediately preceding land plants. Geographical distribution of phototrophs shows a bias against flavodoxin-containing organisms in iron-rich coastal/freshwater habitats. Based on these observations, we propose that plants evolved from freshwater macroalgae that already lacked flavodoxin because they thrived in an iron-rich habitat with no need to back up ferredoxin functions and therefore no selective pressure to keep the flavodoxin gene. Conversely, ferredoxin retention in the plant lineage is probably related to its higher efficiency as an electron carrier, compared with flavodoxin. Several lines of evidence supporting these contentions are presented and discussed.

Description: Epibiotic biofilms have the potential to control major aspects of the biology and ecology of their hosts. Their composition and function may thus be essential for the health of the host. We tested the influence of salinity on the composition of epibacterial communities associated with the brown macroalga Fucus vesiculosus. Algal individuals were incubated at three salinities (5, 19, and 25) for 14days and nonliving reference substrata (stones) were included in the experiment. Subsequently, the composition of their surface-associated bacterial communities was analyzed by 454 pyrosequencing of 16S rRNA gene sequences. Redundancy analysis revealed that the composition of epiphytic and epilithic communities significantly differed and were both affected by salinity. We found that 5% of 2494 epiphytic operational taxonomic units at 97% sequence similarity were responsible for the observed shifts. Epibacterial -diversity was significantly lower at salinity 5 but did not differ between substrata. Our results indicate that salinity is an important factor in structuring alga-associated epibacterial communities with respect to composition and/or diversity. Whether direct or indirect mechanisms (via altered biotic interactions) may have been responsible for the observed shifts is discussed.

Description: In his In Depth News story “Warming may not swamp islands” (1 August, p. 496), C. Pala argues that “coral reefs supporting sandy atoll islands will grow and rise in tandem with the sea,” based largely on studies that showed stable Pacific-island area over recent decades (1–4). He suggests that recent land losses are driven mostly by bad choices and that islanders are being affected “for the same reason as millions of people on the continents: because they live too close to shore.” We disagree with these conclusions.

Description: Decadal trends in the properties of seawater adjacent to Antarctica are poorly known, and the mechanisms responsible for such changes are uncertain. Antarctic ice sheet mass loss is largely driven by ice shelf basal melt, which is influenced by ocean-ice interactions and has been correlated with Antarctic Continental Shelf Bottom Water (ASBW) temperature. We document the spatial distribution of long-term large-scale trends in temperature, salinity, and core depth over the Antarctic continental shelf and slope. Warming at the seabed in the Bellingshausen and Amundsen seas is linked to increased heat content and to a shoaling of the mid-depth temperature maximum over the continental slope, allowing warmer, saltier water greater access to the shelf in recent years. Regions of ASBW warming are those exhibiting increased ice shelf melt.

Description: In spite of considerable insights into the microbial diversity of marine sponges, quantitative information on microbial abundances and community composition remains scarce. Here, we established qPCR assays for the specific quantification of four bacterial phyla of representative sponge symbionts as well as the kingdoms Eubacteria and Archaea. We could show that the 16S rRNA gene numbers of Archaea, Chloroflexi, and the candidate phylum Poribacteria were 4–6 orders of magnitude higher in high microbial abundance (HMA) than in low microbial abundance (LMA) sponges and that actinobacterial 16S rRNA gene numbers were 1–2 orders higher in HMA over LMA sponges, while those for Cyanobacteria were stable between HMA and LMA sponges. Fluorescence in situ hybridization of Aplysina aerophoba tissue sections confirmed the numerical dominance of Chloroflexi, which was followed by Poribacteria. Archaeal and actinobacterial cells were detected in much lower numbers. By use of fluorescence-activated cell sorting as a primer- and probe-independent approach, the dominance of Chloroflexi, Proteobacteria, and Poribacteria in A. aerophoba was confirmed. Our study provides new quantitative insights into the microbiology of sponges and contributes to a better understanding of the HMA/LMA dichotomy.

Description: Dinoflagellates harboring diatom endosymbionts (termed “dinotoms”) have undergone a process often referred to as “tertiary endosymbiosis”—the uptake of algae containing secondary plastids and integration of those plastids into the new host. In contrast to other tertiary plastids, and most secondary plastids, the endosymbiont of dinotoms is distinctly less reduced, retaining a number of cellular features, such as their nucleus and mitochondria and others, in addition to their plastid. This has resulted in redundancy between host and endosymbiont, at least between some mitochondrial and cytosolic metabolism, where this has been investigated. The question of plastidial redundancy is particularly interesting as the fate of the host dinoflagellate plastid is unclear. The host cytosol possesses an eyespot that has been postulated to be a remnant of the ancestral peridinin plastid, but this has not been tested, nor has its possible retention of plastid functions. To investigate this possibility, we searched for plastid-associated pathways and functions in transcriptomic data sets from three dinotom species. We show that the dinoflagellate host has indeed retained genes for plastid-associated pathways and that these genes encode targeting peptides similar to those of other dinoflagellate plastid-targeted proteins. Moreover, we also identified one gene encoding an essential component of the dinoflagellate plastid protein import machinery, altogether suggesting the presence of a functioning plastid import system in the host, and by extension a relict plastid. The presence of the same plastid-associated pathways in the endosymbiont also extends the known functional redundancy in dinotoms, further confirming the unusual state of plastid integration in this group of dinoflagellates.

Description: The macroalga Fucus vesiculosus carries a specific community of surface bacteria. To identify chemical compounds that possibly mediate abundance and community composition of algae-associated bacteria, we tested the effect of surface-available algal compounds on bacterial settlement and community composition under field conditions. Compounds on algal thalli were separated from the surface by extraction with organic solvents and investigated on growth inhibition and settlement of bacterial isolates. Based on in vitro data, partially purified extract fractions were then exposed to bacterial colonizers in vivo followed by bacterial enumeration and community analysis. The polar fraction of the algal surface extract revealed a significant profouling effect for Vibrionales, whereas the nonpolar fraction containing the xanthophyll pigment fucoxanthin and other unidentified nonpolar surface compounds revealed a significant 80% reduction of surface colonizing bacteria. The analysis of bacterial surface communities by 454 pyrosequencing demonstrated that the antifouling activity of nonpolar algal surface compounds was targeting the abundance of natural bacterial colonizers rather than the relative composition of bacterial members within the community. Moreover, the bacterial community composition on F.vesiculosus was markedly different from artificial control substrates and chemically manipulated experimental treatments, suggesting that other, nonextractable surface features and/or physical properties render algal-specific epiphytic bacterial communities.

Description: Following their transmission from the human to the mosquito with the bloodmeal, malaria parasites have to persevere in the mosquito midgut for approximately 1 d. During this period the parasites are highly vulnerable to factors of the mosquito midgut, including bacteria. We here aimed at determining the microbial diversity of gut bacteria of the Asian malaria vector Anophebs stephensi (Liston) during development and under different feeding regimes, including feeds on malaria parasite-infected blood. 16S rRNA and denaturing gradient gel electrophoresis analyses demonstrated an increasing reduction in the microbial diversity during mosquito development from egg to adult and identified the gram-negative bacterium Elizabethkingia meningoseptica King as the dominant species in the midgut of lab-reared male and female mosquitoes. E. meningoseptica is transmitted between generations and its predominance in the mosquito midgut was not altered by diet, when the gut microbiota was compared between sugar-fed and blood-fed female mosquitoes. Furthermore, feeds on blood infected with malaria parasites did not impact the presence of E. meningoseptica in the gut. Extracts from cultured E. meningoseptica were active against gram-positive and negative bacteria and yeast and against the blood and gametocyte transmission stages of the malaria parasite Plasmodium falciparum Welch. The antimicrobial and antiplasmodial activities of E. meningoseptica may account for its dominance in the midgut of the malaria vector.

Description: The concentration of CO2 in global surface ocean waters is increasing due to rising atmospheric CO2 emissions, resulting in lower pH and a lower saturation state of carbonate ions. Such changes in seawater chemistry are expected to impact calcification in calcifying marine organisms. However, other physiological processes related to calcification might also be affected, including enzyme activity. In a mesocosm experiment, macroalgal communities were exposed to three CO2 concentrations (380, 665, and 1486 µatm) to determine how the activity of two enzymes related to inorganic carbon uptake and nutrient assimilation in Corallina officinalis, an abundant calcifying rhodophyte, will be affected by elevated CO2 concentrations. The activity of external carbonic anhydrase, an important enzyme functioning in macroalgal carbon-concentrating mechanisms, was inversely related to CO2 concentration after long-term exposure (12 weeks). Nitrate reductase, the enzyme responsible for reduction of nitrate to nitrite, was stimulated by CO2 and was highest in algae grown at 665 µatm CO2. Nitrate and phosphate uptake rates were inversely related to CO2, while ammonium uptake was unaffected, and the percentage of inorganic carbon in the algal skeleton decreased with increasing CO2. The results indicate that the processes of inorganic carbon and nutrient uptake and assimilation are affected by elevated CO2 due to changes in enzyme activity, which change the energy balance and physiological status of C. officinalis, therefore affecting its competitive interactions with other macroalgae. The ecological implications of the physiological changes in C. officinalis in response to elevated CO2 are discussed.

Description: The sheer scale of the metagenomic and metatranscriptomic datasets that are now available warrants the development of automated protocols for organizing, annotating and comparing the samples in terms of their metabolic profiles. We describe a user-friendly java program FROMP (Fragment Recruitment on Metabolic Pathways) for mapping and visualizing enzyme annotations onto the Kyoto Encyclopedia of Genes and Genomes (KEGG) metabolic pathways or custom-made pathways and comparing the samples in terms of their Pathway Completeness Scores,their relative Activity Scores or enzyme enrichment odds ratios. This program along with our fully-configurable PERL based annotation organization pipeline Meta2Pro (METAbolic PROfiling of META-omic data) offers a quick and accurate standalone solution for metabolic profiling of environmental samples or cultures from different treatments. Apart from pictorial comparisons, FROMP can also generate score matrices for multiple meta-omics samples which can be used directly by other statistical programs

Description: We investigated the gelatinous carnivore zooplankton community in the coastal waters of northeast Taiwan during the period 20072010. The community assemblage was composed of 45 species, of which only 14 appeared recurrent in Taiwanese waters. Although there was no clear seasonality, higher richness and abundances occurred in spring and autumn. Examination of potential physical drivers of coastal biomass accumulation did not show any link with water mass transport; instead, peak events were associated with typhoon disturbances, suggesting a potential resource pulse effect.

Description: The mechanics of great subduction earthquakes are influenced by the frictional properties, structure, and composition of the plate-boundary fault. We present observations of the structure and composition of the shallow source fault of the 2011 Tohoku-Oki earthquake and tsunami from boreholes drilled by the Integrated Ocean Drilling Program Expedition 343 and 343T. Logging-while-drilling and core-sample observations show a single major plate-boundary fault accommodated the large slip of the Tohoku-Oki earthquake rupture, as well as nearly all the cumulative interplate motion at the drill site. The localization of deformation onto a limited thickness (less than 5 meters) of pelagic clay is the defining characteristic of the shallow earthquake fault, suggesting that the pelagic clay may be a regionally important control on tsunamigenic earthquakes.

Description: Large coseismic slip was thought to be unlikely to occur on the shallow portions of plate-boundary thrusts, but the 11 March 2011 Tohoku-Oki earthquake [moment magnitude (Mw) = 9.0] produced huge displacements of ~50 meters near the Japan Trench with a resultant devastating tsunami. To investigate the mechanisms of the very large fault movements, we conducted high-velocity (1.3 meters per second) friction experiments on samples retrieved from the plate-boundary thrust associated with the earthquake. The results show a small stress drop with very low peak and steady-state shear stress. The very low shear stress can be attributed to the abundance of weak clay (smectite) and thermal pressurization effects, which can facilitate fault slip. This behavior provides an explanation for the huge shallow slip that occurred during the earthquake.

Description: Sediment-covered basalt on the flanks of mid-ocean ridges constitutes most of Earth’s oceanic crust, but the composition and metabolic function of its microbial ecosystem are largely unknown. By drilling into 3.5-million-year-old subseafloor basalt, we demonstrated the presence of methane- and sulfur-cycling microbes on the eastern flank of the Juan de Fuca Ridge. Depth horizons with functional genes indicative of methane-cycling and sulfate-reducing microorganisms are enriched in solid-phase sulfur and total organic carbon, host d13C- and d34S-isotopic values with a biological imprint, and show clear signs of microbial activity when incubated in the laboratory. Downcore changes in carbon and sulfur cycling show discrete geochemical intervals with chemoautotrophic d13C signatures locally attenuated by heterotrophic metabolism.

Description: This study aimed to examine interactive effects between ocean acidification and temperature on the photosynthetic and growth performance of Neosiphonia harveyi. N. harveyi was cultivated at 10 and 17.5 °C at present (~380 µatm), expected future (~800 µatm), and high (~1500 µatm) pCO2. Chlorophyll a fluorescence, net photosynthesis, and growth were measured. The state of the carbon-concentrating mechanism (CCM) was examined by pH-drift experiments (with algae cultivated at 10 °C only) using ethoxyzolamide, an inhibitor of external and internal carbonic anhydrases (exCA and intCA, respectively). Furthermore, the inhibitory effect of acetazolamide (an inhibitor of exCA) and Tris (an inhibitor of the acidification of the diffusive boundary layer) on net photosynthesis was measured at both temperatures. Temperature affected photosynthesis (in terms of photosynthetic efficiency, light saturation point, and net photosynthesis) and growth at present pCO2, but these effects decreased with increasing pCO2. The relevance of the CCM decreased at 10 °C. A pCO2 effect on the CCM could only be shown if intCA and exCA were inhibited. The experiments demonstrate for the first time interactions between ocean acidification and temperature on the performance of a non-calcifying macroalga and show that the effects of low temperature on photosynthesis can be alleviated by increasing pCO2. The findings indicate that the carbon acquisition mediated by exCA and acidification of the diffusive boundary layer decrease at low temperatures but are not affected by the cultivation level of pCO2, whereas the activity of intCA is affected by pCO2. Ecologically, the findings suggest that ocean acidification might affect the biogeographical distribution of N. harveyi.

Description: The 2011 moment magnitude 9.0 Tohoku-Oki earthquake produced a maximum coseismic slip of more than 50 meters near the Japan trench, which could result in a completely reduced stress state in the region. We tested this hypothesis by determining the in situ stress state of the frontal prism from boreholes drilled by the Integrated Ocean Drilling Program approximately 1 year after the earthquake and by inferring the pre-earthquake stress state. On the basis of the horizontal stress orientations and magnitudes estimated from borehole breakouts and the increase in coseismic displacement during propagation of the rupture to the trench axis, in situ horizontal stress decreased during the earthquake. The stress change suggests an active slip of the frontal plate interface, which is consistent with coseismic fault weakening and a nearly total stress drop.

Description: The frictional resistance on a fault during slip controls earthquake dynamics. Friction dissipates heat during an earthquake; therefore, the fault temperature after an earthquake provides insight into the level of friction. The Japan Trench Fast Drilling Project (Integrated Ocean Drilling Program Expedition 343 and 343T) installed a borehole temperature observatory 16 months after the March 2011 moment magnitude 9.0 Tohoku-Oki earthquake across the fault where slip was ~50 meters near the trench. After 9 months of operation, the complete sensor string was recovered. A 0.31°C temperature anomaly at the plate boundary fault corresponds to 27 megajoules per square meter of dissipated energy during the earthquake. The resulting apparent friction coefficient of 0.08 is considerably smaller than static values for most rocks.

Description: In the Arctic, under-ice primary production is limited to summer months and is restricted not only by ice thickness and snow cover but also by the stratification of the water column, which constrains nutrient supply for algal growth. Research Vessel Polarstern visited the ice-covered eastern-central basins between 82° to 89°N and 30° to 130°E in summer 2012, when Arctic sea ice declined to a record minimum. During this cruise, we observed a widespread deposition of ice algal biomass of on average 9 grams of carbon per square meter to the deep-sea floor of the central Arctic basins. Data from this cruise will contribute to assessing the effect of current climate change on Arctic productivity, biodiversity, and ecological function.

Description: The vast majority of microbes inhabiting the subseafloor remain uncultivated and their energy sources unknown. Thus, a focus of ocean drilling expeditions over the past decade has been to characterize the distribution of microbes associated with specific metabolic reactions. An important question has been whether microbes involved in key microbial processes, such as sulfate reduction and methanogenesis, differ fundamentally from their counterparts in surface environments. To this end, functional genes of anaerobic methane cycling (mcrA), sulfate reduction (dsrAB), acetogenesis (fhs), and dehalorespiration (rdhA) have been examined. A compilation of existing functional gene data suggests that subseafloor microbes involved in anaerobic methane cycling, sulfate reduction, acetogenesis, and dehalorespiration are not fundamentally different from their counterparts in the surface world. Moreover, quantifications of mcrA and dsrAB suggest that, unless the majority of subseafloor microbes involved in methane cycling and sulfate reduction are too genetically divergent to be detected with conventional methods, these processes only support a small fraction (〈 1%) of total microbial biomass in the deep biosphere. Ecological explanations for the observed trends, target processes and methods for future investigations, and strategies for tackling the unresolved issue of microbial contamination in samples obtained by ocean drilling are discussed.

Description: Trophodynamics of meso-zooplankton in the North Sea (NS) were assessed at a site in the southern NS, and at a shallow and a deep site in the central NS. Offshore and neritic species from different ecological niches, including Calanus spp., Temora spp. and Sagitta spp., were collected during seven cruises over 14 months from 2007 to 2008. Bulk stable isotope (SI) analysis, phospholipid-derived fatty acid (PLFA) compositions, and δ 13CPLFA data of meso-zooplankton and particulate organic matter (POM) were used to describe changes in zooplankton relative trophic positions (RTPs) and trophodynamics. The aim of the study was to test the hypothesis that the RTPs of zooplankton in the North Sea vary spatially and seasonally, in response to hydrographic variability, with the microbial food web playing an important role at times. Zooplankton RTPs tended to be higher during winter and lower during the phytoplankton bloom in spring. RTPs were highest for predators such as Sagitta sp. and Calanus helgolandicus and lowest for small copepods such as Pseudocalanus elongatus and zoea larvae (Brachyura). δ 15NPOM-based RTPs were only moderate surrogates for animals’ ecological niches, because of the plasticity in source materials from the herbivorous and the microbial loop food web. Common (16:0) and essential (eicosapentaenoic acid, EPA and docosahexaenoic acid, DHA) structural lipids showed relatively constant abundances. This could be explained by incorporation of PLFAs with δ 13C signatures which followed seasonal changes in bulk δ 13CPOM and PLFA δ 13CPOM signatures. This study highlighted the complementarity of three biogeochemical approaches for trophodynamic studies and substantiated conceptual views of size-based food web analysis, in which small individuals of large species may be functionally equivalent to large individuals of small species. Seasonal and spatial variability was also important in altering the relative importance of the herbivorous and microbial food webs.