Variations in pelagic bacterial communities in the North Atlantic Ocean coincide with water bodies

Similar documents
Multivariate Analysis of Ecological Data

Multivariate Statistics 101. Ordination (PCA, NMDS, CA) Cluster Analysis (UPGMA, Ward s) Canonical Correspondence Analysis

Homework 5: Background Ocean Water Properties & Stratification

BIO 682 Multivariate Statistics Spring 2008

Lecture 2: Diversity, Distances, adonis. Lecture 2: Diversity, Distances, adonis. Alpha- Diversity. Alpha diversity definition(s)

Introduction to multivariate analysis Outline

1.3. Principal coordinate analysis. Pierre Legendre Département de sciences biologiques Université de Montréal

Short Cruise Report - RV Maria S. Merian, cruise MSM-58/2 -

Chad Burrus April 6, 2010

Algebra of Principal Component Analysis

Short Cruise Report R/V Maria S. Merian MSM-35T. Istanbul - Málaga Chief Scientist: Dr. Tomas Feseker Captain: Ralf Schmidt

Multivariate analysis of genetic data an introduction

Supplementary Materials for

Arctic Ocean Biology. from the surface to the deep sea

Supporting Information

4. Ordination in reduced space

Dimensionality Reduction Techniques (DRT)

Unconstrained Ordination

FIG S1: Rarefaction analysis of observed richness within Drosophila. All calculations were

INTRODUCTION TO MULTIVARIATE ANALYSIS OF ECOLOGICAL DATA

Lecture 5: Ecological distance metrics; Principal Coordinates Analysis. Univariate testing vs. community analysis

Physical Oceanography

An Introduction to Ordination Connie Clark

Supplement of Vegetation greenness and land carbon-flux anomalies associated with climate variations: a focus on the year 2015

The Ocean Floor THE VAST WORLD OCEAN

GeoComputation 2011 Session 4: Posters Discovering Different Regimes of Biodiversity Support Using Decision Tree Learning T. F. Stepinski 1, D. White

Lecture 5: Ecological distance metrics; Principal Coordinates Analysis. Univariate testing vs. community analysis

Physical Oceanography

Ordination & PCA. Ordination. Ordination

Characterizing and predicting cyanobacterial blooms in an 8-year

The California current is the eastern boundary current that lies to the west of

Indirect Gradient Analysis

An Introduction to R for the Geosciences: Ordination I

Chapter 1 Test on Geography Skills

Supplementary Information

Multivariate analysis

Habitat use patterns and edge effects across a seagrassunvegetated ecotone depend on species-specific behaviors and sampling methods

Multivariate Statistics (I) 2. Principal Component Analysis (PCA)

G E INTERACTION USING JMP: AN OVERVIEW

sphericity, 5-29, 5-32 residuals, 7-1 spread and level, 2-17 t test, 1-13 transformations, 2-15 violations, 1-19

Multivariate analysis of genetic data: an introduction

... Europe. Based on Bloom s Taxonomy. Environment Interactions Movement. Human & Location. Regions. Place

Outline Classes of diversity measures. Species Divergence and the Measurement of Microbial Diversity. How do we describe and compare diversity?

Lecture Topic Projects 1 Intro, schedule, and logistics 2 Applications of visual analytics, data types 3 Data sources and preparation Project 1 out 4

Statistics Toolbox 6. Apply statistical algorithms and probability models

Interannual changes in the zooplankton community structure on the southeastern Bering Sea shelf and Chukchi Sea during summers of

Canary. European shelf break. Canary current N Conjunction of Northern Equatorial Current and Countercurrent. from Cameroon to Congo

Jae-Bong Lee 1 and Bernard A. Megrey 2. International Symposium on Climate Change Effects on Fish and Fisheries

Other resources. Greengenes (bacterial) Silva (bacteria, archaeal and eukarya)

Level 3 Earth and Space Science, 2016

Correspondence Analysis & Related Methods

MAR 110 LECTURE #10 The Oceanic Conveyor Belt Oceanic Thermohaline Circulation

The North Atlantic Oscillation: Climatic Significance and Environmental Impact

C

Map Skills Unit. Note taking unit

Multivariate Data Analysis a survey of data reduction and data association techniques: Principal Components Analysis

RECENT STUDIES OF THE SEA OF OKHOTSK

Depth distributions of alkaline phosphatase and phosphonate utilization genes in the North Pacific Subtropical Gyre

Setting Priorities for Eelgrass Conservation and Restoration. Robert Buchsbaum Massachusetts Audubon Society

Global Distribution of Northern Fur Seals (2013)

Activity #2 - Major Ocean Surface Currents

Bird Species richness per 110x110 km grid square (so, strictly speaking, alpha diversity) -most species live there!

2/15/2012. Earth System Science II EES 717 Spring 2012

ROBERTO BATTITI, MAURO BRUNATO. The LION Way: Machine Learning plus Intelligent Optimization. LIONlab, University of Trento, Italy, Apr 2015

Edwin A. Hernández-Delgado*

Marine Ecoregions. Marine Ecoregions. Slide 1. Robert G. Bailey. USDA Forest Service Rocky Mountain Research Station

Ocean Circulation- PART- I: In Class. Must be done inclass, and turned in before you leave for credit.

DIMENSION REDUCTION AND CLUSTER ANALYSIS

4/4/2018. Stepwise model fitting. CCA with first three variables only Call: cca(formula = community ~ env1 + env2 + env3, data = envdata)

Principal Components Analysis. Sargur Srihari University at Buffalo

Using Topological Data Analysis to find discrimination between microbial states in human microbiome data

Seasonal cycle of phytoplankton community composition off Newport, Oregon, in 2009

ROLES OF THE OCEAN MESOSCALE IN THE LATERAL SUPPLY OF MASS, HEAT, CARBON AND NUTRIENTS TO THE NORTHERN HEMISPHERE SUBTROPICAL GYRE

The fieldwork during the Polarstern cruise ANT XVI/2 as a contribution to the study of bottom water formation and sea ice transport in the Weddell Sea

4/2/2018. Canonical Analyses Analysis aimed at identifying the relationship between two multivariate datasets. Cannonical Correlation.

CHAPTER 7 Ocean Circulation Pearson Education, Inc.

Lecture 1. Equations of motion - Newton s second law in three dimensions. Pressure gradient + force force

SIO 210 Final examination Answer Key for all questions except Daisyworld. Wednesday, December 10, PM Name:

Fall equinox means 12 hours of sunlight for everyone, everywhere

Actual bathymetry (with vertical exaggeration) Geometry of the ocean 1/17/2018. Patterns and observations? Patterns and observations?

Ocean Circulation. In partnership with Dr. Zafer Top

Experimental Design and Data Analysis for Biologists

Proximity data visualization with h-plots

Gopalkrishna Veni. Project 4 (Active Shape Models)

Revision: Chapter 1-6. Applied Multivariate Statistics Spring 2012

Topic # 12 Natural Climate Processes

Multivariate Analysis of Ecological Data using CANOCO

Quantitative Genomics and Genetics BTRY 4830/6830; PBSB

Dissimilarity and transformations. Pierre Legendre Département de sciences biologiques Université de Montréal

Upper Ocean Circulation

Topic # 12 How Climate Works

Carlo Vittorio Cannistraci. Minimum Curvilinear Embedding unveils nonlinear patterns in 16S metagenomic data

( ( ( ( ( ( ( (! ! 1!

EDAMI DIMENSION REDUCTION BY PRINCIPAL COMPONENT ANALYSIS

Introduction to ordination. Gary Bradfield Botany Dept.

Lab #3 Background Material Quantifying Point and Gradient Patterns

On Sampling Errors in Empirical Orthogonal Functions

-Principal components analysis is by far the oldest multivariate technique, dating back to the early 1900's; ecologists have used PCA since the

Near-Isometry by Relaxation: Supplement

2. In words, what is electrical current? 3. Try measuring the current at various points of the circuit using an ammeter.

Transcription:

The following supplement accompanies the article Variations in pelagic bacterial communities in the North Atlantic Ocean coincide with water bodies Richard L. Hahnke 1, Christina Probian 1, Bernhard M. Fuchs 2, Jens Harder 1, * 1 Department of Microbiology and 2 Department of Molecular Ecology, Max Planck Institute for Marine Microbiology, Celsiusstr. 1, 28359 Bremen, Germany *Corresponding author. Email: jharder@mpi-bremen.de Aquatic Microbial Ecology 71: 131 140 (2013) Supplement 1. The supplement provides a map of the sampling area, the density of the samples in different water depths, k-dominance curves of T-RFLP (terminal restriction fragment length polymorphism) patterns, graphs of nmds (non-metric multidimensional scaling) and PCA (principal component) analyses, and a documentation of TRFs (terminal restriction fragments) relative fluorescence intensity along the sampling transect Fig. S1. Sampling stations and water density of the RV 'Maria S. Merian cruise 03/1 VISION in 2006. (A) Water samples were obtained from the North Atlantic Ocean along the 30 W meridian from the productive cold Greenland current (Boreal Polar [BPLR], 66 39 N) across the cold north (Atlantic Arctic [ARCT]) and warm south (North Atlantic Drift [NADR]) of the North Atlantic Current to the oligotrophic central Atlantic Ocean (North Atlantic Subtropical Gyre [NAST], 34 24 N). (B) The density of the seawater (σ θ ) from depths between 20 m and 200 m indicated the pycnocline between 50 m and 100 m

Fig. S2. Identification of T-RFLP patterns with low phylogenetic information was possible with the k-dominance plot of T-RFLP patterns. The normalized relative fluorescence intensities (RFI) were visualized in rank versus cumulated abundance curves. Two exemplary sampling sites are shown (A,B), in which each line represents the cumulative RFI of forward (black) or reverse (gray) terminal restriction fragments of one T-RFLP analysis. The species rank at 100% cumulative abundance (RFI) represents the richness of TRFs of the sampling site. (A) The T- RFLP patterns with a high amount of false positive signals (indicated by a star) originated from a fixed fluorescence intensity threshold and overall low fluorescence intensity in this T-RFLP pattern. These T-RFLP patterns were identified as outliers and excluded from further analysis. (B) In contrast, the cumulated abundance curves of T-RFLP patterns of comparable good quality had only slightly different fluorescence intensities between T-RFLP patterns Fig. S3. The diversity richness differed between terminal restriction fragments derived from the forward (black) and from the reverse (gray) primer. Each black and gray line represents the cumulative relative fluorescence intensity of the terminal restriction fragments of 1 T-RFLP analysis. The species rank at 100% cumulative abundance (RFI) represents the richness of TRFs of the sampling site. TRFs of the reverse primer had a lower species rank compared to the TRFs of the forward primer at the same cumulative abundance. Thus, the richness of the reverse primer TRFs was less than that of the forward primer TRFs 2

Fig. S4. Nonmetric multidimensional scaling (nmds) was applied to visualize (A) the Bray-Curtis similarity based on relative fluorescence intensities and (B) the Sørensen index based on presence/absence of terminal restriction fragments (TRFs) in TRF patterns of each sampling site (β diversity). Both biplots had a low stress value of 0.12, indicating a meaningful 2-dimensional visualization. The presentation of the differences revealed a clustering of bacterial communities from the individual water bodies (square: BPLR, dot: ARCT, diamond: NADR, triangle: NAST) along a latitudinal gradient.. A hierarchical clustering defined groups of bacterial communities at 48% Bray-Curtis similarity and 55% Sørensen index (indicated by solid gray lines) Fig. S5. Principal component analysis was applied to find terminal restriction fragments (TRFs) that cause changes in the overall community structure. The eigenvectors and eigenvalues revealed a differentiation between the northern sampling sites, BPLR (dot) and ARCT (square), above and below 75 m depth, and the southern sampling sites in NADR (diamond) and NAST (triangle), along the first principal component (PC 1). The second principal component (PC 2) distinguished the northern sampling sites (BPLR, ARCT) above and below 75 m depth and 3

distinguished the southern sampling sites in NADR and NAST. TRF_189nt had the largest eigenvector (gray arrow) parallel to the first principal component, indicating a large contribution to the population in the south of the transect. Numbers in gray represent the depth of the sampling site. A circle represents sampling sites that fall into 1 hierarchical cluster of 50% similarity 4

Fig. S6. Relative fluorescence intensity (RFI) pattern of significant terminal restriction fragments (TRFs) affiliated with (A) the Flavobacteriaceae Polaribacter and (B) group NS4. Contour lines indicate the relative abundance of TRFs (% RFI). A star indicates the sampling site with the largest RFI and the inset shows the corresponding fluorescence intensity (FI) in the T-RFLP profile 5

Fig. S7. Relative fluorescence intensity (RFI) pattern of terminal restriction fragments (TRFs) with a regional distribution and no affiliation to a single taxon in the in silico (itrf) calculation. Contour lines indicate the relative abundance of TRFs (% RFI). A star indicates the sampling site with the largest RFI, and the inset shows the corresponding fluorescence intensity (FI) in the T-RFLP profile 6