(Table S1) Molecular analyses and morphological comparisons of deep-sea clams

Large vesicomyid clams are common inhabitants of sulphidic deep-sea habitats such as hydrothermal vents, hydrocarbon seeps and whale-falls. Yet, the species- and genus-level taxonomy of these diverse clams has been unstable due to insufficiencies in sampling and absence of detailed taxonomic studies that would consistently compare molecular and morphological characters. To clarify uncertainties about species-level assignments, we examined DNA sequences from mitochondrial cytochrome-c-oxidase subunit I (COI) in conjunction with morphological characters. New and published COI sequences were used to create a molecular database for 44 unique evolutionary lineages corresponding to species. Overall, the congruence between molecular and morphological characters was good. Several discrepancies due to synonymous species designations were recognized, and acceptable species names were rectified with published COI sequences in cases where morphological specimens were available. We identified seven species with trans-Pacific distributions, and two species with Indo-Pacific distributions. Presently, 27 species have only been documented from one region, which might reflect limited ranges, or insufficient geographical sampling. Vesicomyids exhibit the greatest species diversity along the northwest Pacific ridge systems and in the eastern Pacific, along the western America margin, where depth zonation typically results in segregation of closely related species. The broad distributions of several vesicomyid species suggest that their required chemosynthetic habitats might be more common than previously recognized and occur along most continental margins.

Geochemistry of ash and pumice of ODP Hole 132-810C

Air-fall volcanic ash and pumice were recovered from 22 intervals in upper Miocene-Pleistocene nannofossil oozes cored in Hole 810C on Shatsky Rise, northwest Pacific. Shatsky Rise is near the eastern limit of ash falls produced by explosive volcanism in arc systems in northern Japan and the Kuriles, more than 1600 km away. Electron probe analyses establish that the ash beds and pumice pebbles are andesitic to rhyolitic in composition, and belong to both tholeiitic and high-alumina lineages similar to tephra from Japanese volcanoes. High-speed winds in the polar-front and subtropical jets are evidently what propelled the ash for such a distance. The pumice arrived by flotation, driven from the same directions by winds, waves, and currents. It is not ice-rafted debris from the north. One thick pumice bed probably was deposited when a large pumice mat passed over Shatsky Rise. Far more abundant ash occurs in sediments cored at DSDP Sites 578 through 580, about 500 km west of Shatsky Rise. Most of the ash and pumice at Shatsky Rise can be correlated with specific ash beds at 1, 2, or all 3 of these sites by interpolating to precisely determined magnetic reversal sequences in the cores. Most of the correlations are to thick ash layers (5.7 +/- 3.0 cm) at one or more sites. These must represent extremely large eruptions that spread ash over very wide areas. Whereas several of the thicker correlative ashes fell from elongate east-trending plumes directed from central Japan, the majority of them - dating from about 2 Ma - came from the North Honshu and Kurile arc systems to the northwest. This direction probably was in response to both long-term and seasonal fluctuations in the location and velocity of the polar-front jet, and to more vigorous winter storm fronts originating over glaciated Siberia.

Geochemical data from pore water and solids from cores GS08-153-91GC12 and GS08-153-53-GC6, Arctic Mid-Ocean Ridge

Microbial communities and their associated metabolic activity in marine sediments have a profound impact on global biogeochemical cycles. Their composition and structure are attributed to geochemical and physical factors, but finding direct correlations has remained a challenge. Here we show a significant statistical relationship between variation in geochemical composition and prokaryotic community structure within deep-sea sediments. We obtained comprehensive geochemical data from two gravity cores near the hydrothermal vent field Loki's Castle at the Arctic Mid-Ocean Ridge, in the Norwegian-Greenland Sea. Geochemical properties in the rift valley sediments exhibited strong centimeter-scale stratigraphic variability. Microbial populations were profiled by pyrosequencing from 15 sediment horizons (59,364 16S rRNA gene tags), quantitatively assessed by qPCR, and phylogenetically analyzed. Although the same taxa were generally present in all samples, their relative abundances varied substantially among horizons and fluctuated between Bacteria- and Archaea-dominated communities. By independently summarizing covariance structures of the relative abundance data and geochemical data, using principal components analysis, we found a significant correlation between changes in geochemical composition and changes in community structure. Differences in organic carbon and mineralogy shaped the relative abundance of microbial taxa. We used correlations to build hypotheses about energy metabolisms, particularly of the Deep Sea Archaeal Group, specific Deltaproteobacteria, and sediment lineages of potentially anaerobic Marine Group I Archaea. We demonstrate that total prokaryotic community structure can be directly correlated to geochemistry within these sediments, thus enhancing our understanding of biogeochemical cycling and our ability to predict metabolisms of uncultured microbes in deep-sea sediments.

(Table 1) Number of alleles per locus and heterozygosity of Arctic charr (Salvelinus alpinus) populations in Iceland

The ecological theory of adaptive radiation predicts that the evolution of phenotypic diversity within species is generated by divergent natural selection arising from different environments and competition between species. Genetic connectivity among populations is likely also to have an important role in both the origin and maintenance of adaptive genetic diversity. Our goal was to evaluate the potential roles of genetic connectivity and natural selection in the maintenance of adaptive phenotypic differences among morphs of Arctic charr, Salvelinus alpinus, in Iceland. At a large spatial scale, we tested the predictive power of geographic structure and phenotypic variation for patterns of neutral genetic variation among populations throughout Iceland. At a smaller scale, we evaluated the genetic differentiation between two morphs in Lake Thingvallavatn relative to historically explicit, coalescent-based null models of the evolutionary history of these lineages. At the large spatial scale, populations are highly differentiated, but weakly structured, both geographically and with respect to patterns of phenotypic variation. At the intralacustrine scale, we observe modest genetic differentiation between two morphs, but this level of differentiation is nonetheless consistent with strong reproductive isolation throughout the Holocene. Rather than a result of the homogenizing effect of gene flow in a system at migration-drift equilibrium, the modest level of genetic differentiation could equally be a result of slow neutral divergence by drift in large populations. We conclude that contemporary and recent patterns of restricted gene flow have been highly conducive to the evolution and maintenance of adaptive genetic variation in Icelandic Arctic charr.

The effect of repeated exposure to parasite genotypes of the eyefluke Diplostomum pseudospathaceum on infection rate and immunization in three-spined sticklebacks

Adaptive immunity in vertebrates can confer increased resistance against invading pathogens upon re-infection. But how specific parasite genotypes affect the transition from innate to adaptive immunity is poorly understood. Here, we investigated the effects of homologous and heterologous exposures of genetically distinct parasite lineages of the eye fluke Diplostomum pseudospathaceum on gene expression patterns of adaptive immunity in sticklebacks (Gasterosteus aculeatus). We showed that observable differences were largely attributable to final exposures and that there is no transcription pattern characteristic for a general response to repeated infections with D. pseudospathaceum. Final exposure did not unify expression patterns of heterologous pre-exposed fish. Interestingly, heterologous final exposures showed similarities between different treatment groups subjected to homologous pre-exposure. The observed pattern was supported by parasite infection rates and suggests that host immunization was optimized towards an adaptive immune response that favored effectiveness against parasite diversity over specificity.

Ichthyolith (fish tooth and denticle) counts from DSDP Site 91-596 and ODP Site 145-886

While the history of taxonomic diversification in open ocean lineages of ray-finned fish and elasmobranchs is increasingly known, the evolution of their roles within the open ocean ecosystem remains poorly understood. To assess the relative importance of these groups through time, we measured the accumulation rate of microfossil fish teeth and elasmobranch dermal denticles (ichthyoliths) in deep sea sediment cores from the North and South Pacific gyres over the past 85 million years. We find three distinct and stable open ocean ecosystem structures, each defined by the relative and absolute abundance of elasmobranch and ray-finned fish remains. The Cretaceous Ocean (pre-66 Ma), was characterized by abundant elasmobranch denticles, but low abundances of fish teeth. The Paleogene Ocean (66-20 Ma), initiated by the Cretaceous/Paleogene Mass Extinction, had nearly 4 times the abundance of fish teeth compared to elasmobranch denticles. This Paleogene Ocean structure remained stable during the Eocene greenhouse (50 Ma) and the Eocene-Oligocene glaciation (34 Ma), despite large changes in overall accumulation of both groups during those intervals, suggesting that climate change is not a primary driver of ecosystem structure. Dermal denticles virtually disappeared from open ocean ichthyolith assemblages about 20 Ma, while fish tooth accumulation increased dramatically in variability, marking the beginning of the Modern Ocean. Together, these results suggest that open ocean fish community structure is stable on long timescales, independent of total production and climate change. The timing of the abrupt transitions between these states suggests that the transitions may be due to interactions with other, non-preserved pelagic consumer groups.

Location, intervals and thickness of Cretaceous sections in Indian Ocean DSDP and ODP sites (Table 2)

The Indian Ocean covers approximately 73.5 * 10**6 km**3 from 25°N to 67°S and from 20° to 120°E. Several legs of the Deep Sea Drilling Project (DSDP) and the Ocean Drilling Program (ODP) have operated in its waters, many penetrating the Cretaceous. Most of the scientific drill sites are DSDP related and thus pre-dated modern biostratigraphic conventions. Foraminifers and calcareous nannoplankton were by far the dominant fossil groups studied in the earlier work, supplemented occasionally by studies of other fossil groups, The results of the Ocean Drilling Project phase are yet too young to be fully integrated but have been based on a broader range of techniques and fossil groups. During most of the Cretaceous, the proto-Indian Ocean basin lay in middle to high latitudes. Thus, it is unrealistic to expect successful routine application of low-latitude zonations. No planktonic foraminifer zonal scheme has been developed for the Indian Ocean basin for several reasons. There are no sections with complete or even significant partial sections to allow development of such a zonation. Carbonate compensation depth (CCD) effects have been marked in most sections, and significant intervals are devoid of planktonic foraminifers. The Indian Ocean now covers a great latitudinal range from tropics to polar regions and, at first glance, no scheme can be expected to be applicable over that entire range. In the Cretaceous the area was much smaller, though expanding progressively, and the paleolatitude range was quite small. Calcareous nannoplankton have proved valuable in dating Indian Ocean Cretaceous sediments and have, perhaps in contrast with the foraminifers, been consistently a more reliable means of applying zonal schemes developed elsewhere. For the Albian-Aptian, zonations based on well-known benthic foraminifer lineages (Scheibnerova, 1974) have been useful when nothing else was available or effective. Palynology has been used little, but where used, has proved excellent. It has the added value of providing valuable information on nearby terrestrial vegetation as the fossils were resistant to dissolution. Normally, when different fossil groups have been applied to a section, the results have been compatible or compatible to an acceptable degree. There are a few instances where incompatibility is noteworthy, and Site 263 is a classic example, as even two calcareous nannoplankton studies show irreconcilable differences here. All groups gave different results, but one benthic foraminifer analysis agreed with one calcareous nannoplankton study.

(Table 2) Collection sites and number of specimens of Usnea subgenus Neuropogon species

Usnea species of the Neuropogon group are amongst the most widespread and abundant macrolichens in Antarctic regions. Four principal species, U. antarctica, U. aurantiaco-atra, U. sphacelata and U. subantarctica, have been described on morphological grounds. However, identification to species level is often difficult and atypical morphologies frequently arise. Over 400 specimens were collected on the Antarctic Peninsula and Falkland Islands. Both morphological and molecular characters (ITS and RPB1) were used to compare samples to clarify taxonomic relationships. Morphological characteristics used included presence of apothecia, apothecial rays, soredia, papillae, fibrils, pigmentation and the diameter of the central axis as a proportion of branch diameter. Results revealed a very close relationship between U. antarctica and U. aurantiaco-atra, suggesting that they might constitute a species pair or be conspecific. Usnea sphacelata was comprised of at least two genetically distinct groups with no clear differences in morphology. One group included the first reported fertile specimen of this species. Usnea subantarctica was phylogenetically distinct from the other main Antarctic Usnea species, but clustered with U. trachycarpa. Genetic variation was evident within all species although there was no clear correlation between geographic origin and genetic relatedness. Phylogenetic analyses indicated that species circumscription in the Neuropogon group needs revision, with the principal species being non-monophyletic. None of the morphological characters, or groups of characters, used in this study proved to be completely unambiguous markers for a single species. However, axis thickness was supported as being informative for the identification of monophyletic lineages within the group.

Seawater carbonate chemistry and biological processes of sand dollars Dendraster excentricus during experiments, 2011

Reduction in global ocean pH due to the uptake of increased atmospheric CO2 is expected to negatively affect calcifying organisms, including the planktonic larval stages of many marine invertebrates. Planktonic larvae play crucial roles in the benthic-pelagic life cycle of marine organisms by connecting and sustaining existing populations and colonizing new habitats. Calcified larvae are typically denser than seawater and rely on swimming to navigate vertically structured water columns. Larval sand dollars Dendraster excentricus have calcified skeletal rods supporting their bodies, and propel themselves with ciliated bands looped around projections called arms. Ciliated bands are also used in food capture, and filtration rate is correlated with band length. As a result, swimming and feeding performance are highly sensitive to morphological changes. When reared at an elevated PCO2 level (1000 ppm), larval sand dollars developed significantly narrower bodies at four and six-arm stages. Morphological changes also varied between four observed maternal lineages, suggesting within-population variation in sensitivity to changes in PCO2 level. Despite these morphological changes, PCO2 concentration alone had no significant effect on swimming speeds. However, acidified larvae had significantly smaller larval stomachs and bodies, suggesting reduced feeding performance. Adjustments to larval morphologies in response to ocean acidification may prioritize swimming over feeding, implying that negative consequences of ocean acidification are carried over to later developmental stages.

Sampling information, length measurements, DNA alignments, models and phylogenetics trees for Oithona similis s.l.

Traditionally, many small-sized copepod species are considered to be widespread, bipolar or cosmopolitan. However, these large-scale distribution patterns need to be re-examined in view of increasing evidence of cryptic and pseudo-cryptic speciation in pelagic copepods. Here, we present a phylogeographic study of Oithona similis s.l. populations from the Arctic Ocean, the Southern Ocean and its northern boundaries, the North Atlantic and the Mediterrranean Sea. O. similis s.l. is considered as one of the most abundant species in temperate to polar oceans and acts as an important link in the trophic network between the microbial loop and higher trophic levels such as fish larvae. Two gene fragments were analysed: the mitochondrial cytochrome oxidase c subunit I (COI), and the nuclear ribosomal 28S genetic marker. Seven distinct, geographically delimitated, mitochondrial lineages could be identified, with divergences among the lineages ranging from 8 to 24 %, thus representing most likely cryptic or pseudocryptic species within O. similis s.l. Four lineages were identified within or close to the borders of the Southern Ocean, one lineage in the Arctic Ocean and two lineages in the temperate Northern hemisphere. Surprisingly the Arctic lineage was more closely related to lineages from the Southern hemisphere than to the other lineages from the Northern hemisphere, suggesting that geographic proximity is a rather poor predictor of how closely related the clades are on a genetic level. Molecular clock application revealed that the evolutionary history of O. similis s.l. is possibly closely associated with the reorganization of the ocean circulation in the mid Miocene and may be an example of allopatric speciation in the pelagic zone.