Phenotypic evolution in a venerid bivalve species lineage from the late Middle Miocene Central Paratethys Sea: a multi-approach morphometric analysis

A morphometric analysis was performed for the late Middle Miocene bivalve species lineage of Polititapes tricuspis (Eichwald, 1829) (Veneridae: Tapetini). Specimens from various localities grouped into two stratigraphically successive biozones, i.e. the upper Ervilia Zone and the Sarmatimactra Zone, were investigated using a multi-method approach. A Generalized Procrustes Analysis was computed for fifteen landmarks, covering characteristics of the hinge, muscle scars, and pallial line. The shell outline was separately quantified by applying the Fast Fourier Transform, which redraws the outline by fitting in a combination of trigonometric curves. Shell size was calculated as centroid size from the landmark configuration. Shell thickness, as not covered by either analysis, was additionally measured at the centroid. The analyses showed significant phenotypic differentiation between specimens from the two biozones. The bivalves become distinctly larger and thicker over geological time and develop circular shells with stronger cardinal teeth and a deeper pallial sinus. Data on the paleoenvironmental changes in the late Middle Miocene Central Paratethys Sea suggest the phenotypic shifts to be functional adaptations. The typical habitats for Polititapes changed to extensive, very shallow shores exposed to high wave action and tidal activity. Caused by the growing need for higher mechanical stability, the bivalves produced larger and thicker shells with stronger cardinal teeth. The latter are additionally shifted towards the hinge center to compensate for the lacking lateral teeth and improve stability. The deepening pallial sinus is related to a deeper burrowing habit, which is considered to impede being washed out in the new high-energy settings.

Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.