Evidence of cryptic and pseudocryptic speciation in the Paracalanus parvus species complex (Crustacea, Copepoda, Calanoida), supplementary material

Introduction Many marine planktonic crustaceans such as copepods have been considered as widespread organisms. However, the growing evidence for cryptic and pseudo-cryptic speciation has emphasized the need of re-evaluating the status of copepod species complexes in molecular and morphological studies to get a clearer picture about pelagic marine species as evolutionary units and their distributions. This study analyses the molecular diversity of the ecologically important Paracalanus parvus species complex. Its seven currently recognized species are abundant and also often dominant in marine coastal regions worldwide from temperate to tropical oceans. Results COI and Cytochrome b sequences of 160 specimens of the Paracalanus parvus complex from all oceans were obtained. Furthermore, 42 COI sequences from GenBank were added for the genetic analyses. Thirteen distinct molecular operational taxonomic units (MOTU) and two single sequences were revealed with cladistic analyses (Maximum Likelihood, Bayesian Inference), of which seven were identical with results from species delimitation methods (barcode gaps, ABDG, GMYC, Rosenberg's P(AB)). In total, 10 to 12 putative species were detected and could be placed in three categories: (1) temperate geographically isolated, (2) warm-temperate to tropical wider spread and (3) circumglobal warm-water species. Conclusions The present study provides evidence of cryptic or pseudocryptic speciation in the Paracalanus parvus complex. One major insight is that the species Paracalanus parvus s.s. is not panmictic, but may be restricted in its distribution to the northeastern Atlantic.

A laboratory-scale experiment and a numerical simulation of unusual spiral plumes in a High-Prandtl-number fluid, links supplementary material

We experimentally and numerically investigated the generation of plumes from a local heat source (LHS) and studied the interaction of these plumes with cellular convective motion (CCM) in a rectangular cavity filled with silicon oil at a Prandtl number (Pr) of approximately two thousand. The LHS is generated using a 0.2-W green laser beam. A roll-type CCM is generated by vertically heating one side of the cavity. The CCM may lead to the formation of an unusual spiral convective plume that resembles a vertical Archimedes spiral. A similar plume is obtained in a direct numerical simulation. We discuss the physical mechanism for the formation of a spiral plume and the application of the results to mantle convection problems. We also estimate the Reynolds (Re) and Rayleigh (Ra) numbers and apply self-similarity theory to convection in the Earth's mantle. Spiral plumes can be used to interpret mantle tomography results over the last decade.

Lake Baikal amphipods under climate change: Thermal constraints and ecological consequences, links to supplementary material

Lake Baikal, the world's most voluminous freshwater lake, has experienced unprecedented warming during the last decades. A uniquely diverse amphipod fauna inhabits the littoral zone and can serve as a model system to identify the role of thermal tolerance under climate change. This study aimed to identify sublethal thermal constraints in two of the most abundant endemic Baikal amphipods, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, and Gammarus lacustris, a ubiquitous gammarid of the Holarctic. As the latter is only found in some shallow isolated bays of the lake, we further addressed the question whether rising temperatures could promote the widespread invasion of this non-endemic species into the littoral zone. Animals were exposed to gradual temperature increases (4 week, 0.8 °C/d; 24 h, 1 °C/h) starting from the reported annual mean temperature of the Baikal littoral (6 °C). Within the framework of oxygen- and capacity-limited thermal tolerance (OCLTT), we used a nonlinear regression approach to determine the points at which the changing temperature-dependence of relevant physiological processes indicates the onset of limitation. Limitations in ventilation representing the first limits of thermal tolerance (pejus (= "getting worse") temperatures (Tp)) were recorded at 10.6 (95% confidence interval; 9.5, 11.7), 19.1 (17.9, 20.2), and 21.1 (19.8, 22.4) °C in E. verrucosus, E. cyaneus, and G. lacustris, respectively. Field observations revealed that E. verrucosus retreated from the upper littoral to deeper and cooler waters once its Tp was surpassed, identifying Tp as the ecological thermal boundary. Constraints in oxygen consumption at higher than critical temperatures (Tc) led to an exponential increase in mortality in all species. Exposure to short-term warming resulted in higher threshold values, consistent with a time dependence of thermal tolerance. In conclusion, species-specific limits to oxygen supply capacity are likely key in the onset of constraining (beyond pejus) and then life-threatening (beyond critical) conditions. Ecological consequences of these limits are mediated through behavioral plasticity in E. verrucosus. However, similar upper thermal limits in E. cyaneus (endemic, Baikal) and G. lacustris (ubiquitous, Holarctic) indicate that the potential invader G. lacustris would not necessarily benefit from rising temperatures. Secondary effects of increasing temperatures remain to be investigated.

Wave-induced changes in seaweed toughness entail plastic modifications in snail traits maintaining consumption efficacy, link to supplementary material

1. Environmental stress can influence species traits and performance considerably. Using a seaweed-snail system from NW (Nova Scotia) and NE (Helgoland) Atlantic rocky shores, we examined how physical stress (wave exposure) modulates traits in the seaweed Fucus vesiculosus and indirectly in its main consumer, the periwinkle Littorina obtusata. 2. In both regions, algal tissue toughness increased with wave exposure. Reciprocal-transplant experiments showed that tissue toughness adjusts plastically to the prevailing level of wave exposure. 3. Choice experiments tested the feeding preference of snails from sheltered, exposed, and very exposed habitats for algae from such wave exposures. Snails from exposed and very exposed habitats consumed algal tissues at similar rates irrespective of the exposure of origin of the algae. However, snails from sheltered habitats consumed less algal tissues from very exposed habitats than tissues from sheltered and exposed habitats. Choice assays using reconstituted algal food (triturated during preparation) identified high thallus toughness as the explanation for the low preference of snails from sheltered habitats for algae from very exposed habitats. 4. Ultrastructural analyses of radulae indicated that rachidian teeth were longest and the number of cusps in lateral teeth (grazing-relevant traits) was highest in snails from very exposed habitats, suggesting that radulae are best suited to rupture tough algal tissues in such snails. 5. No-choice feeding experiments revealed that these radular traits are also phenotypically plastic, as they adjust to the toughness of the algal food. 6. Synthesis. This study indicates that the observed plasticity in the feeding ability of snails is mediated by wave exposure through phenotypic plasticity in the tissue toughness of algae. Thus, plasticity in consumers and their resource species may reduce the potential effects of physical stress on their interaction.

On the Southern Ocean CO2 uptake and the role of the biological carbon pump in the 21st century, links to supplementary material

We use a suite of eight ocean biogeochemical/ecological general circulation models from the MAREMIP and CMIP5 archives to explore the relative roles of changes in winds (positive trend of Southern Annular Mode, SAM) and in warming- and freshening-driven trends of upper ocean stratification in altering export production and CO2 uptake in the Southern Ocean at the end of the 21st century. The investigated models simulate a broad range of responses to climate change, with no agreement ona dominance of either the SAM or the warming signal south of 44° S. In the southernmost zone, i.e., south of 58° S, they concur on an increase of biological export production, while between 44 and 58° S the models lack consensus on the sign of change in export. Yet, in both regions, the models show an enhanced CO2 uptake during spring and summer. This is due to a larger CO 2 (aq) drawdown by the same amount of summer export production at a higher Revelle factor at the end of the 21st century. This strongly increases the importance of the biological carbon pump in the entire Southern Ocean. In the temperate zone, between 30 and 44° S all models show a predominance of the warming signal and a nutrient-driven reduction of export production. As a consequence, the share of the regions south of 44° S to the total uptake of the Southern Ocean south of 30° S is projected to increase at the end of the 21st century from 47 to 66% with a commensurable decrease to the north. Despite this major reorganization of the meridional distribution of the major regions of uptake, the total uptake increases largely in line with the rising atmospheric CO2. Simulations with the MITgcm-REcoM2 model show that this is mostly driven by the strong increase of atmospheric CO2, with the climate-driven changes of natural CO2 exchange offsetting that trend only to a limited degree (~10%) and with negligible impact of climate effects on anthropogenic CO2 uptake when integrated over a full annual cycle south of 30° S.

Treeline dynamics in Siberia under changing climates as inferred from an individual-based model for Larix, supplementary material

Siberian boreal forests are expected to expand northwards in the course of global warming. However, processes of the treeline ecotone transition, as well astiming and related climate feedbacks are still not understood. Here, we present 'Larix Vegetation Simulator' LAVESI, an individual-based spatially-explicit model that can simulate Larix gmelinii (RUPR.) RUPR. stand dynamics in an attempt to improve our understanding about past and future treeline movements under changing climates. The relevant processes (growth, seed production and dispersal, establishment and mortality) are incorporated and adjusted to observation data mainly gained from the literature. Results of a local sensitivity analysis support the robustness of the model's parameterization by giving relatively small sensitivity values. We tested the model by simulating tree stands under modern climate across the whole Taymyr Peninsula, north-central Siberia (c. 64-80° N; 92-119° E). We find tree densities similar to observed forests in the northern to mid-treeline areas, but densities are overestimated in the southern parts of the simulated region. Finally, from a temperature-forcing experiment, we detect that the responses of tree stands lag the hypothetical warming by several decades, until the end of 21st century. With our simulation experiments we demonstrate that the newly-developed model captures the dynamics of the Siberian latitudinal treeline.

Impacts of seawater acidification on mantle gene expression patterns of the Baltic Sea blue mussel: implications for shell formation and energy metabolism, link to supplementary material

Marine organisms have to cope with increasing CO2 partial pressures and decreasing pH in the oceans. We elucidated the impacts of an 8-week acclimation period to four seawater pCO2 treatments (39, 113, 243 and 405 Pa/385, 1,120, 2,400 and 4,000 µatm) on mantle gene expression patterns in the blue mussel Mytilus edulis from the Baltic Sea. Based on the M. edulis mantle tissue transcriptome, the expression of several genes involved in metabolism, calcification and stress responses was assessed in the outer (marginal and pallial zone) and the inner mantle tissues (central zone) using quantitative real-time PCR. The expression of genes involved in energy and protein metabolism (F-ATPase, hexokinase and elongation factor alpha) was strongly affected by acclimation to moderately elevated CO2 partial pressures. Expression of a chitinase, potentially important for the calcification process, was strongly depressed (maximum ninefold), correlating with a linear decrease in shell growth observed in the experimental animals. Interestingly, shell matrix protein candidate genes were less affected by CO2 in both tissues. A compensatory process toward enhanced shell protection is indicated by a massive increase in the expression of tyrosinase, a gene involved in periostracum formation (maximum 220-fold). Using correlation matrices and a force-directed layout network graph, we were able to uncover possible underlying regulatory networks and the connections between different pathways, thereby providing a molecular basis of observed changes in animal physiology in response to ocean acidification.