1 1
Molecular phylogenetics and morphology of Gambierdiscus yasumotoi from tropical eastern
2
Australia
3 4
Shauna Murray*1,2, Paolo Momigliano*3,4, , Kirsten Heimann3,5, David Blair3
5
*These authors contributed equally
6
1
7
NSW 2007, Australia
8
2
Sydney Institute of Marine Sciences, Chowder Bay Rd, Mosman NSW Australia
9
3
School of Marine and Tropical Biology, James Cook University, QLD 4810, Australia
10
4
Department of Biological Sciences, Macquarie University, North Ryde 2019 NSW,
11
Australia
12
5
13
4810, Australia
Plant Functional Biology and Climate Change Cluster, University of Technology, Sydney
Centre for Sustainable Tropical Fisheries and Aquaculture, James Cook University, QLD
14 15
Corresponding author:
16
Paolo Momigliano
17
email address:
[email protected]
18
phone: +61 2 9850 8331
19 20 21
2 22
Abstract
23
Here the occurrence of the species Gambierdiscus yasumotoi is reported for the first time
24
along a latitudinal gradient spanning more than 1550 km of the Australian Great Barrier Reef
25
(GBR), a region with endemic ciguatera fish poisoning. Gambierdiscus yasumotoi was found
26
at three tropical and sub tropical coral reef sites, Raine Island (northern GBR) , Nelly Bay
27
(central GBR) and Heron Island (southern GBR), indicating a wide-ranging distribution in
28
tropical and subtropical eastern Australia. Specimens from Australia broadly fitted the
29
original description of Gambierdiscus yasumotoi, but differed in some aspects, showing some
30
similarities to G. ruetzleri. Molecular phylogenetic analyses based on nuclear rRNA gene
31
sequences and morphological analyses showed specimens to be intermediate between the two
32
species Gambierdiscus yasumotoi and Gambierdiscus ruetzleri. The full intraspecific
33
diversity of these two species appears to be incompletely known, and these two species may
34
represent a species complex. Strains of this species from other sites around the world have
35
been found to produce an as yet unknown toxin, possibly an analogue of maitotoxin.
36
Keywords: Gambierdiscus, Ciguatera, Great Barrier Reef, Phylogenetics
37
3 38
1 Introduction
39
The genus Gambierdiscus Adachi and Fukuyo (Gonyaulacales, Dinoflagellata, Alveolata) is
40
one of more than 550 known genera of dinoflagellates. Species of this genus are the main
41
causative agents of Ciguatera Fish Poisoning (CFP), caused by the consumption of fish
42
contaminated with toxins produced by these dinoflagellates, which are biomagnified along
43
the food web. CFP is the most common non-bacterial seafood poisoning disease in the world
44
(Fleming et al., 1998) and mainly occurs in tropical countries. Cases have been reported in
45
coral reef areas in Australia (Gillespie et al., 1985; Lewis, 2006), the Pacific and Indian
46
oceans as well as the Caribbean (Tosteson, 2004). In Australia, there were more than 1,400
47
documented cases between 1965 and 2010, including two fatalities (Gillespie et al., 1985;
48
Hamilton et al., 2010).
49
Gambierdiscus are epibenthic on seagrass, macroalgae, sand and coral rubble, however, they
50
can also occur in the plankton (Nakahara et al., 1996; Parsons et al., 2011). Species of this
51
genus have the plate formula: Po, 3', 7", 6c,?s, 5’’', 1p, 2”” (Chinain et al., 1999; Faust, 1995;
52
Fraga et al., 2011; Holmes, 1998; Litaker et al., 2009) and are distinguished from one another
53
based on characters such as surface patterns, cell flattening (apical or dorso-ventral), size and
54
shapes of the plates and of the cells.
55
The type species, Gambierdiscus toxicus, was first identified as recently as 1979 (Adachi and
56
Fukuyo, 1979), and until 15 years ago, only 3 species were known in the genus. It was
57
recognised that the type description of Gambierdiscus toxicus probably incorporated several
58
species, and G. toxicus was redescribed more narrowly, based on a lectotype (Litaker et al.,
59
2009). Currently, twelve species of Gambierdiscus are recognised: Gambierdiscus australes,
60
G. belizeanus, G. caribaeus, G. carolinianus, G. carpenteri, G. pacificus, G. polynesiensis,
61
G. ruetzleri, G. toxicus, G scabrosus, G. yasumotoi, and G. excentricus (Chinain et al., 1999;
62
Faust, 1995; Fraga et al., 2011; Holmes, 1998; Litaker et al., 2009; Nishimura et al., 2014).
63
At least five additional unnamed clades that are clearly genetically distinguishable have been
64
found, indicating that further species may be distinguished in future (Litaker et al., 2009;
65
Richlen et al., 2008; Xu et al., 2014).
66
Published reports initially indicated the presence of only a single species, Gambierdiscus
67
toxicus in Australia, from tropical regions of Queensland (Holmes et al., 1990; Holmes et al.,
68
1991; Holmes et al., 1994), and more recently, from New South Wales (Murray, 2010).
69
Cultures of this taxon from sites in Queensland were made in the 1990s (Babinchak et al.,
4 70
1994; Holmes et al., 1990; Holmes et al., 1991; Holmes et al., 1994) and subsequently lost.
71
Since these initial reports, our understanding of the morphological and genetic diversity
72
present in the genus Gambierdiscus has increased significantly (Chinain et al., 1999; Fraga et
73
al., 2011; Litaker et al., 2009; Nishimura et al., 2014; Richlen et al., 2008; Xu et al., 2014),
74
and it is therefore necessary to revisit initial sites in order to assess the presence of these
75
species. Recently, a second species of Gambierdiscus, G. carpenteri, has been reported from
76
several sites in Australia for the first time (Kohli et al., 2014a; Kohli et al., 2014b).
77
The species Gambierdiscus yasumotoi was originally described from Singapore waters
78
(Holmes, 1998), and was found to produce an uncharacterised toxin that was lethal to mice,
79
producing symptoms similar to those of maitotoxin (Holmes, 1998). This species has since
80
been reported from additional sites around the world, including from the Pacific and
81
Caribbean waters. Strains isolated from New Zealand were found not to produce known
82
analogues of ciguatoxins (CTXs) or maitotoxin (MTX) 1, but were likely to produce an as yet
83
uncharacterised analogue of maitotoxin (Rhodes et al., 2014a).
84
In this study the presence of Gambierdiscus yasumotoi is shown for the first time along a
85
latitudinal gradient of more than 1550 km in waters of the northern, central and southern
86
Great Barrier Reef (GBR), Queensland, Australia. Gambierdiscus yasumotoi is described
87
from Australia based on scanning electron microscopy of field material, as well as scanning
88
electron microscopy and phylogenetic analysis of a culture of G. yasumotoi isolated from the
89
central GBR region. The occurrence of two other species of Gamberdiscus is also
90
documented: G. carpenteri and G. cf belizeanus.
91
2 Materials and Methods
92
2.1 Sample collection.
93
Between the 3rd-8th December 2003, eight samples were taken at low tide at sites on Raine
94
Island (11° 45’.501 S, 144° 02’.309 E), a remote uninhabited tropical coral cay in the
95
northern GBR, in the vicinity of the tower end of the beach (Fig 1), as part of annual
96
sampling coordinated by the Raine Island Corporation. On 15th June 2003, four samples were
97
collected from Heron Island (23°25’ S, 151° 55’ E), a tropical coral cay in the southern GBR
98
(Fig 1). Samples were taken from encrusting macroalgae in shallow (less than 30 cm deep)
99
tidal pools at low tide.
100
5 101
Macroalgal samples were obtained from central GBR waters on low tide in the area between
102
Townsville (Pallarenda at the mouth of Three Mile Creek, 19° 12' 34'' S, 146° 46' 36'' E) and
103
Magnetic Island (Nelly Bay, 1910’ S, 14650’ E; Great Barrier Reef Marine Park Authority
104
Permit G06/20234.1) in August 2004 and March 2008 (Fig 1). Macroalgal material was
105
washed in filtered (0.45 µm) seawater to flush epiphytic dinoflagellates off their surface and
106
the wash water was sequentially concentrated by filtration through 53 and 20 µm nylon mesh
107
filters.
108 109
2.2 Culture establishment.
110
Dinoflagellates were isolated from Pallarenda and Nelly Bay macroalgal wash waters at the
111
North Queensland Algal Identification/Culturing Facility (NQAIF) at James Cook University,
112
Townsville, Australia using the microcapillary‐capturing‐technique (Andersen and Kawachi,
113
2005) at 10X magnification on an inverted light microscope (Olympus CKK 41, Olympus
114
Australia Ltd, Mt Waverley VIC 3149). Cells were dispensed in autoclaved and filtered (0.45
115
μm) seawater. Cells were allowed to swim for ten minutes and were then recaptured. This
116
procedure was repeated ten times to ensure that nano‐ and pico‐plankton were no longer in
117
the vicinity of the cell to be isolated. Cultures of NQAIF116 (Pallarenda) and NQAIF210
118
(Nelly Bay) were established in L1 medium (Guillard and Hargraves, 1993) prepared in
119
natural seawater and maintained at 24°C, a 12:12 h photoperiod and light intensity of 45
120
μmol photons m‐2s‐1 in a Contherm cross‐flow phytoplankton growth chamber. Cultures were
121
subcultured in L1 medium every 4 weeks
122 123
2.3 Microscopy.
124
Samples were immediately examined using a dissecting light microscope. Individual cells
125
were isolated onto 5 µm millipore filtersor 1 cm round glass coverslips that had been treated
126
with poly-L-lysine (Marchant and Thomas, 1983), to make preparations to observe in the
127
scanning electron microscope. Cells were dehydrated in a series of increasing ethanol
128
concentrations (15, 30, 50, 70, 90, and 100%), followed by 50% and 100%
129
hexamethyldisilizane (HMDS). The filters were air dried before being sputter coated with
130
gold or gold/palladium. They were observed using a Philips 505 scanning electron
131
microscope (Philips, Eindhoven, the Netherlands) at 10–15kV and a JEOL JSM-5410LV
132
scanning electron microscope (JEOL, Sydney, Australia).
133
6 134
2.4 DNA extraction and sequencing.
135
DNA was extracted from cultures NQAIF210 and NQAIF116 using a modified Chelex®
136
protocol (Walsh et al., 1991) as outlined in Momigliano et al. (2013). The near-complete
137
nuclear 18S rRNA gene, and the D1-D3 and D8-D10 regions of the 28S rRNA gene were
138
amplified and sequenced using the primers listed in Table S1. PCR reactions were set up as
139
follows: 1µl of template DNA, 800 µM each DNTPs, 10 pmol of each primer, 1X Green
140
GoTaq® Flexi Buffer and 1 unit GoTaq® Flexi Polymerase (Promega, Inc.). PCR reactions
141
underwent an initial denaturation of 3 min at 94 °C, followed by 35 cycles of 30 s
142
denaturation (94 °C), 30 s annealing (see Table S1 for details), and 1 min extension (72 °C),
143
and a final extension of 5 min at 72 °C. PCR products were cleaned by isopropanol
144
precipitation and sequenced using a commercial service (Macrogen, Inc.). All sequences
145
were deposited in GenBank (accession numbers: KM272970-272974).
146
2.5 Sequence alignment and phylogenetic analysis.
147
Overlapping fragments of the D1-D3 and D8-D10 regions of the 28S and the near-complete
148
18S genes were assembled using the software Chromas Pro (Technelysium Pty Ltd.). The
149
sequences obtained from culture NQAIF210 and NQAIF116 were aligned with available
150
sequences of Gambierdiscus spp. (for a list of the strains used see Table S2) using Clustal W
151
(Thompson et al., 1994) and the alignments were visually refined using BioEdit (Hall, 1999).
152
Phylogenetic analysis was performed on the datasets including the D8-D10 region of the 28S
153
rRNA and the near complete 18S rRNA gene. Substitution models were selected for each
154
dataset in jModelTest2 (Darriba et al., 2012), using Bayesian Information Criterion as a
155
measure of the relative quality of the models. Phylogenetic reconstructions were carried out
156
using both a Maximum Likelihood (ML) and Bayesian Inference (BI) approach.
157
ML trees were produced in PhyML 3.1 (Guindon et al., 2010), using a TrN+G substitution
158
model with 4 gamma categories (alpha= 0.46) and a TIM2+G model with 4 gamma
159
categories (alpha=0.46) for the D8-D10 region of the 28S gene and the 18S dataset
160
respectively, tree improvement by using the best of nearest neighbour interchanges (NNI)
161
and subtree pruning and regrafting (NPR), starting with 5 random trees. Bayesian analysis
162
was performed in MrBayes 3.2.2 (Ronquist and Huelsenbeck, 2003). The substitution model
163
used for the D8-D10 region of the 28S gene was the same as in the ML analysis. However as
164
MrBayes does not implement the TIM2 substitution model, we used the next-best available
165
model, the HKY+G model. The number of Markov Chain Monte Carlo generations was set to
7 166
2 000 000, but set to automatically stop when the standard deviation of split frequencies fell
167
below 0.01. Trees were sampled every 100 generations, and convergence of the 2
168
independent runs was checked by analysing the stability of each parameter estimate in the
169
software Tracer v1.6.
170
In addition, for each locus p-distances for all available strains belonging to the G.
171
yasumotoi/ruetzleri clade were estimated in the R statistical environment using the ape
172
package (Paradis et al., 2004). The distance matrices were further analysed using Principal
173
Coordinates Analysis (PCoA) implemented in the package ade4 (Dray and Dufour, 2007).
174 175
3 Results
176
3.1 Phylogenetic analyses.
177
The final alignments of the D8-D10 region of the 28S rRNA gene based on 45 sequences
178
included 929 unambiguously aligned nucleotides and 348 variable sites, of which 286 were
179
parsimony informative. The 18S alignment included 36 sequences, 1746 unambiguously
180
aligned nucleotides and 696 variable sites of which 603 were parsimony informative. Both
181
phylogenies show strong support for all terminal branches (Figs 2 & 3), and largely support
182
previously published phylogenies of the genus Gambierdiscus. Small differences in the trees
183
presented here relative to previously published phylogenies are due to the fact that out-group
184
taxa were not used in phylogenetic reconstruction and disappear if the trees are rooted with
185
the most basal Gambierdiscus clade (G. yasumotoi clade).
186
NQAIF210 grouped with high support within the G. ruetzleri/G. yasumotoi clade - Clade I in
187
Nishimura et al. (2013) in both phylogenies (Figs 2 & 3 ). NQAIF210 could not be
188
unambiguously identified as either G. yasumotoi nor G. ruetzleri based on 18S sequences. In
189
the phylogeny reconstructed using the D8-D10 region of the 28S gene, NQAIF210 grouped
190
with high support with the New Zealand strain CAWD210, and both NQAIF210 and
191
CAWD210 grouped with G. ruetzleri with moderate bootstrap support (Fig 3). The culture
192
NQAIF116 consistently grouped with high support within the G. carpenteri clade (Figs 2 &
193
3). Representations of p-distance matrices in reduced space (Fig 4) showed that NQAIF210
194
and CAWD210 have identical sequences for the D1-D3 region, and have small p-distances
195
for the D8-D10 region but the analysis also failed to show any clear grouping of NQAIF210
196
with other previously identified strains in this clade, from which NQAIF210 is roughly
8 197
equidistant (Fig. 4) . However, it must be remembered that all sequences of G. yasumotoi
198
utilized in the phylogenies and p-distance analyses are clones of the same strain (Gyasu), and
199
therefore they cannot be used as a proxy of intra-specific variation. Furthermore, p-distances
200
within this clade, and between putative species within it, are very small when compared to
201
between-species distances elsewhere in the genus. P-distances between G. yasumotoi and G.
202
ruetzleri are as little as 0.01 (D1-D3 and D8-D10 regions of the 28S rRNA gene) and 0.004
203
(near complete 18S gene).
204
Interestingly, p-distances among clones of the same strain of G. yasumotoi are comparable to
205
those observed between clones of G. yasumotoi and strains of G. ruetzleri (Fig 4, particularly
206
4A, 4C and 4E).
207
3.2 Morphological analyses
208
The morphology of Gambierdiscus yasumotoi from field samples and from the cultured strain
209
is shown (Figs 5,6,7). Cells from field samples were laterally compressed in ventral view, and
210
oval to round in lateral view, with a depth of 57 µm (range 54- 59 µm), width 48 µm (range
211
46-50 µm) and length 56 µm (range 53-59 µm) (Figs 6,7). The depth to width ratio was 1.18,
212
while the length to width ratio was 1.16. Cultured cells were slightly smaller, more globular
213
and less obviously laterally compressed, with a depth of 49 µm (range 44- 54 µm), width 45
214
µm (range 40-49 µm) and length 51 µm (range 49-54 µm) (Fig 5). The depth to width ratio
215
was 1.08, while the length to width ratio was 1.13 (Table 1).
216
The thecal plates were smooth with many small round pores, approximately 0.3 µm diameter
217
(Figs 5b, 7a,b). The epicone had the plate formula of Gambierdiscus, with three apical plates
218
and seven precingular plates (Figs 5a, 6b,g,h,i). The apical pore plate was tear drop shaped,
219
while the pore was long and fish hook shaped at the cell apex, approximately 8-9 µm long,
220
with 33-45 pores surrounding (Figs 5b, 6g, 7a). The hypocone consisted of five postcingular
221
plates, two antapical plates and one posterior intercalary plate (Figs 5f, g, i, 6f).
222
The cingulum was deeply excavated and descending, displaced about 1-2 cingular widths
223
(Figs 5c, 6a, b, c, 7b, c, d). Prominent cingular lists were present in all cells,(Figs 5i, 6a, c,
224
7b, c, d). In cultured cells, the cingulum was positioned slightly higher than the mid point of
225
the cell in left lateral view, approximately in the mid point of the cell in dorsal view, and just
226
below the mid point of the cell in right lateral view (Figs 5c, d, e). In cells from field samples,
227
the cingulum was positioned higher, approximately 1/3 of the way down the cell from the
9 228
apex in left lateral view, and slightly above the mid point of the cell in right lateral view (Figs
229
6a, c, f).
230
The sulcus has been reported to consist of seven plates, the Sdp, Sda, Ssp, Ssa, t, Sma and
231
Smp (Litaker et al., 2009). Of these, the Ssp, Ssa, Sdp, t and Sda plates were clearly visible in
232
this study (Fig. 7b, c, d). The tiny, oval shaped Sma plate, which is adjacent to the Sda plate,
233
was also visible in several images (Fig. 7c, d). The Smp plate could not be observed in our
234
images. Prominent sulcal lists were present in all cells from the field (Figs 5i, 6a, c, 7b, c, d),
235
and some cultured cells (Fig. 5c), but absent in other cultured cells (Fig 5h). The sulcal lists
236
consisted of a long list on the left side of the Ssa and Ssp plates, and a shorter list on the right
237
side of the Sdp plate (Fig. 7b, c, d). The 2'''' plate was fork shaped and invaded the sulcus
238
(Figs 5g, h, i, 6a).
239
Substantial plate overlaps, particularly in the hypothecal cells, were found in some cells in the
240
cultured material (Fig 5g).
241
Two other species of Gambierdiscus were found in field samples: Gambierdiscus cf
242
belizeanus, and Gambierdiscus carpenteri. These species were only partially documented
243
from the sites, and the morphological information about them is shown in the Supplementary
244
Figs.
245
4 Discussion
246
Gambierdiscus yasumotoi has previously been reported from sub-tropical Japan (Nishimura
247
et al., 2013),the Persian Gulf and the Red Sea (Saburova et al., 2013), Singapore (Holmes,
248
1998), the Mexican Caribbean (Hernández-Becerril and Almazán Becerril, 2004), and even
249
the northern coast of New Zealand (Rhodes et al., 2014a), but not from Australian waters.
250
Here the occurrence of G. yasumotoi at three sites along a latitudinal gradient spanning more
251
than 1550 km of the Australian GBR is reported, a region with endemic CFP.
252
Records of G. yasumotoi in tropical and sub-tropical regions in the Indian, and Pacific
253
Oceans, both in the northern and southern hemispheres, suggest that this species has a
254
circum-global distribution in tropical and sub-tropical waters. Interesting is the occurrence of
255
G. yasumotoi in northern New Zealand, reported by Rhodes et al. (2014a), as sea surface
256
temperatures (SST) fall well below 15 º C in this region, the lowest thermal limit for any
257
species of Gambierdiscus reported to date (Kibler et al., 2012). However, as (Rhodes et al.,
258
2014a) noted, the specimens were collected in summer, when SST was 23 ºC. It is possible
10 259
that the occurrence of G. yasumotoi in New Zealand is seasonal, with cells being transported
260
from tropical and sub-tropical Australia during the warmer months via the East Australian
261
Current (EAC) and the Tasman Front. The EAC is known to affect phytoplankton
262
communities in the eastern Australian coast, transporting tropical species into sub-tropical
263
and temperate areas (Armbrecht et al., 2013).
264
4.1 Comparison between Gambierdiscus yasumotoi and G. ruetzleri
265
Interestingly, most of the strains of G. yasumotoi which have been recently isolated have
266
intermediate morphologies between G. yasumotoi and the more recently described G.
267
ruetzleri (Table 1). According to the original description of G. ruetzleri, the features that
268
distinguish this species from G. yasumotoi are: size (G. ruetzleri is significantly smaller and
269
has a cell width less than 42 µm); a smaller, narrower and differently shaped 2’’’’ plate; a
270
greater depth to width ratio, meaning they are proportionally narrower; a smaller length to
271
width ratio, and a smaller epitheca to hyptheca ratio, as the cingulum is closer to the centre of
272
the cell, rather than closer to the cell apex (Litaker et al., 2009).
273
There was some variability in these features in both individual cells from the culture, and
274
between cells found in field samples as compared to cultured cells (Table 1). Both the field
275
and cultured material had L:W and D:W ratios that were more similar to those in the original
276
description of Gambierdiscus yasumotoi than that of G. ruetzleri (Table 1). However, the
277
epitheca to hypotheca ratios appeared to differ between the field and cultured material. For
278
the cultured strain, this ratio seems to be smaller than that reported for G. yasumotoi, but the
279
difference is not nearly as pronounced as in the original description of G. ruetzleri (Litaker et
280
al., 2009). In previous studies, the D:W and L:W ratios of material from New Zealand were
281
most similar to those of G. ruetzleri. However in overall cell size, these specimens more
282
closely resembled G. yasumotoi (Rhodes et al. 2014a, Table 1). In the material from Kuwait
283
(Saburova et al., 2013), the D:W ratio and epitheca:hypotheca ratio appeared to be most
284
similar to those of Gambierdiscus ruetzleri, while the cell width and overall cell size were
285
most similar to those of G. yasumotoi.
286
In terms of the shape of the 2’’’’ plate, variability was found in the Australian samples. In
287
some cells this appeared to have an almost horizontal posterior side (Figs 5f, 6a), similar to
288
the shape as described for G. yasumotoi. Other cells had a 2’’’’ plate with a 45° angle on the
289
posterior side of this plate (Fig. 6c), a shape more similar to that of G. ruetzleri (Litaker et al.,
290
2009).
11 291
Several specimens of Gambierdiscus yasumotoi with aberrant plate pattern have been
292
observed among cultured cells (Saburova et al., 2013) and cells with substantial plate
293
overlaps were found in this study (Fig 5g), which would tend to distort the plate size or
294
shape. For this reason, interpreting species-level differences based on minor differences in
295
plate sizes or shapes could be misleading.
296 297
In our molecular analysis, NQAIF210 grouped with high support within the G.
298
yasumtoi/ruetzleri clade, but it was not clearly grouped with either based on 18S gene
299
sequences (Figs 3,4). The phylogeny constructed using the D8-D10 region of the 28S rRNA
300
gene suggests that NQAIF210 and the G. cf. yasumotoi isolated from New Zealand constitute
301
a highly supported sub-clade, and the two strains had identical sequences for the D1-D3
302
region of the same gene (Fig. 3). The genetic data support the morphological observations,
303
suggesting that these strains are indeed very similar. Interestingly, while both strains were
304
morphologically closer to the original description of G. yasumotoi, in the D8-D10 phylogeny
305
they grouped with high support with G. ruetzleri (Fig. 3). There are now reports of a number
306
of strains, including NQAIF210 (this study) , IR4G and Go3 (Nishimura et al., 2013) and
307
CAWD210 (Rhodes et al., 2014a), which always group within the G. yasumotoi /G. ruetzleri
308
clade, but cannot be clearly identified based on molecular and morphological data as either
309
species.
310 311
Our analyses of p-distances among strains (Fig 4) belonging to the Gambierdiscus yasumotoi
312
clade suggest that caution is needed in trying to delineate species boundaries within this
313
group. Firstly, the erection of G. ruetzleri as a new species was based on a comparison with a
314
single strain of G. yasumotoi, the strain used in the original description (Litaker et al., 2009).
315
The estimates of within-species p-distances for G. yasumotoi are based on multiple sequences
316
from this same strain, and therefore it could be considered a form of pseudo-replication to use
317
such differences to recognise G. yasumotoi and G. ruetzleri as distinct, well supported entities
318
in the phylogenetic analysis. Furthermore, our PCoA analyses of p-distance matrices for 18S
319
and 28S rRNA gene sequences seems to suggest that p-distances among clones of G.
320
yasumotoi are in some cases comparable to p-distances between strains of G. yasumotoi and
321
G. ruetzleri (Fig. 4). In addition, the genetic differences between these two species are less
322
than those determined for within-species divergence for other Gambierdiscus species
323
(Nishimura et al., 2013; Rhodes et al., 2014a). Morphological and genetic distances similar
324
to, or even larger than the ones reported between putative species within the G.
12 325
yasumoti/ruetzleri clade have been described within species of other dinoflagellate genera,
326
such as Coolia (Fraga et al., 2008; Momigliano et al., 2013; Rhodes et al., 2014b) and
327
Amphidinium carterae (Murray et al., 2012). Within the G. yasumotoi/ruetzleri clade, the
328
discovery of a number of strains which show intermediate morphologies and cannot be
329
unambiguously identified as either G. yasumotoi or G. ruetzleri (Nishimura et al., 2013;
330
Rhodes et al., 2014a) throws into doubt the status of these as distinct species. It is possible
331
that G. yasumotoi and G. ruetzleri are conspecific, and the minor genetic and morphological
332
differences reported are at the population level, or that they form an as yet unresolved species
333
complex.
334 335
4.2Toxicity
336
In the original description of Gambierdiscus yasumotoi, methanol extracts were found to be
337
toxic to mice with symptoms similar to those of maitotoxins (MTX) (Holmes, 1998). The
338
strain of Gambierdiscus yasumotoi from New Zealand was found by LC-MS analysis not to
339
produce known analogues of MTX or of the monitored analogues of CTXs at detectable
340
levels (Rhodes et al., 2014a). However, a putative MTX analogue which had the same mass
341
as MTX-3 (Holmes and Lewis, 1994; Lewis et al., 1994) was found (Rhodes et al., 2014a). A
342
study using the human erythrocyte lysis assay (HELA) showed significant toxicity in
343
Gambierdiscus ruetzleri (Holland et al., 2013). Experiments which used specific inhibitors
344
of the MTX pathway and purified MTX, Gambierdiscus whole cell extracts, and hydrophilic
345
cell extracts containing MTX, were consistent with MTX as the primary hemolytic
346
compound produced in those studies (Holland et al., 2013). MTX was generally considered
347
not to be related to CFP symptoms, which are thought to be solely caused by ciguatoxins
348
(CTX) (Lewis and Holmes, 1993). However, the propensity for MTX to accumulate in fish
349
flesh following exposure to Gambierdiscus spp has been re-examined experimentally, and its
350
role in fish toxicity remains incompletely known(Kohli et al., 2014c). A strain of G.
351
yasumotoi was recently found to increase mortality and the expression of stress related genes
352
in copepods (Lee et al., 2014)
353 354 355
Acknowledgements
13 356
We thank the Raine Island Corporation for funding the field trip to Raine Island (SM), the
357
Australian Biological Resources Study for additional funding (SM) and the University of
358
Sydney Australian Centre for Microscopy and Microanalysis for access to the scanning
359
electron microscopes (SM). We also thank the Marine and Tropical Sciences Research
360
Facility for funding field work and laboratory work in the central GBR (KH, DB, PM).
361 362
Author contributions
363
Conceived and designed the experiments: PM, SM, KH, DB. Performed the experiments:
364
PM, SM, KH. Analysed the data: PM, SM. Contributed reagents/materials/analytical tools:
365
KH, SM, DB. Wrote the paper: PM, SM. Provided feedback on the manuscript: KH, DB.
366 367 368
References
369
Adachi, R., Fukuyo, Y., 1979. The thecal structure of a marine dinoflagellate Gambierdiscus
370
toxicus gen. et sp. nov. collected in a ciguatera-endemic area. Bull. Jap. Soc. Sci. Fish. 45(1),
371
67-71.
372
Andersen, R.A., Kawachi, M., 2005. Traditional microalgae isolation techniques, In:
373
Anderson, R.A. (Ed.), Algal culturing techniques. Elsevier, Amsterdam, pp. 83-100.
374
Armbrecht, L.H., Roughan, M., Rossi, V., Schaeffer, A., Davies, P.L., Waite, A.M., Armand,
375
L.K., 2013. Phytoplankton composition under contrasting oceanographic conditions:
376
Upwelling and downwelling (Eastern Australia). Cont. Shelf Res.
377
Babinchak, J.A., Moeller, P.D.R., Van Dolah, F.M., Eyo, P.B., Ramsdell, J.S., 1994.
378
Production of ciguatoxins in cultured Gambierdiscus toxicus. Mem. Queensl. Mus. 34, 447-
379
453.
380
Chinain, M., Faust, M.A., Pauillac, S., 1999. Morphology and molecular analyses of three
381
toxic species of Gambierdiscus (Dinophyceae): G. pacificus, sp. nov., G. australes, sp. nov.,
382
and G. polynesiensis, sp. nov. J. Phycol. 35(6), 1282-1296.
14 383
Darriba, D., Taboada, G.L., Doallo, R., Posada, D., 2012. jModelTest 2: more models, new
384
heuristics and parallel computing. Nat. Methods 9(8), 772-772.
385
Dray, S., Dufour, A.-B., 2007. The ade4 package: implementing the duality diagram for
386
ecologist. J. Stat. Softw. 22(4), 1-20.
387
Faust, M.A., 1995. Observation of sand-dwelling toxic dinoflagellates (Dinophyceae) from
388
widely different sites, including two new species. J. Phycol. 31(6), 996-1003.
389
Fleming, L.E., Baden, D.G., Bean, J.A., Weisman, R., Blythe, D.G., 1998. Seafood toxin
390
diseases: issues in epidemiology and community outreach, In: Reguera, B., Blanco, J.,
391
Fernandez, M., Wyatt, T. (Eds.), Harmful Algae. Xunta de Galicia and Intergovernmental
392
Oceanographic Commission of UNESCO, pp. 245-248.
393
Fraga, S., Penna, A., Bianconi, I., Paz, B., Zapata, M., 2008. Coolia canariensis sp. nov.
394
(Dynophyceae), a new nontoxic epiphytic benthic dinoflagellate frm the Canary Islands. J.
395
Phycol. 44(4), 1060-1070.
396
Fraga, S., Rodríguez, F., Caillaud, A., Diogène, J., Raho, N., Zapata, M., 2011.
397
Gambierdiscus excentricus sp. nov.(Dinophyceae), a benthic toxic dinoflagellate from the
398
Canary Islands (NE Atlantic Ocean). Harmful Algae 11, 10-22.
399
Gillespie, N., Lewis, R., Pearn, J., Bourke, A., Holmes, M., Bourke, J., Shields, W., 1985.
400
Ciguatera in Australia. Occurrence, clinical features, pathophysiology and management. Med.
401
J. Aust. 145(11-12), 584-590.
402
Guillard, R., Hargraves, P., 1993. Stichochrysis immobilis is a diatom, not a chrysophyte.
403
Phycologia 32(3), 234-236.
404
Guindon, S., Dufayard, J.-F., Lefort, V., Anisimova, M., Hordijk, W., Gascuel, O., 2010.
405
New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the
406
performance of PhyML 3.0. Syst. Biol. 59(3), 307-321.
407
Hall, T., 1999. BioEdit: a user-friendly biological sequence alignment editor and analysis
408
program for Windows 95/98/NT. Nucl. Acid S. 41, 95-98.
409
Hamilton, B., Whittle, N., Shaw, G., Eaglesham, G., Moore, M.R., Lewis, R.J., 2010. Human
410
fatality associated with Pacific ciguatoxin contaminated fish. Toxicon 56(5), 668-673.
15 411
Hernández-Becerril, D., Almazán Becerril, A., 2004. Especies de dinoflagelados del género
412
Gambierdiscus (Dinophyceae) del Mar Caribe mexicano. Rev. Biol. Trop. 52, 77-87.
413
Holland, W.C., Litaker, R.W., Tomas, C.R., Kibler, S.R., Place, A.R., Davenport, E.D.,
414
Tester, P.A., 2013. Differences in the toxicity of six Gambierdiscus (Dinophyceae) species
415
measured using an in vitro human erythrocyte lysis assay. Toxicon 65, 15-33.
416
Holmes, M.J., 1998. Gambierdiscus yasumotoi sp. nov. (Dinophyceae), a toxic benthic
417
dinoflagellate from southeastern Asia. J. Phycol. 34(4), 661-668.
418
Holmes, M.J., Lewis, R.J., 1994. Purification and characterisation of large and small
419
maitotoxins from cultured Gambierdiscus toxicus. Nat. Toxins 2(2), 64-72.
420
Holmes, M.J., Lewis, R.J., Gillespie, N.C., 1990. Toxicity of Australian and French
421
Polynesian strains of Gambierdiscus Toxicus (Dinophyceae) grown in culture:
422
Characterization of a new type of maitotoxin. Toxicon 28(10), 1159-1172.
423
Holmes, M.J., Lewis, R.J., Poli, M.A., Gillespie, N.C., 1991. Strain dependent production of
424
ciguatoxin precursors (gambiertoxins) by Gambierdiscus toxicus (Dinophyceae) in culture.
425
Toxicon 29(6), 761-775.
426
Holmes, M.J., Lewis, R.J., Sellin, M., Street, R., 1994. The origin of ciguatera in Platypus
427
Bay. Mem. Queensl. Mus. 34, 505-512.
428
Kibler, S.R., Litaker, R.W., Holland, W.C., Vandersea, M.W., Tester, P.A., 2012. Growth of
429
eight Gambierdiscus (Dinophyceae) species: Effects of temperature, salinity and irradiance.
430
Harmful Algae 19, 1-14.
431
Kohli, G.S., Murray, S., Neilan, B.A., Rhodes, L., Harwood, D.T., Smith, K., Mayer, L.,
432
Capper, A., Hallegraeff, G., 2014a. Abundance of the potentially maitotoxic dinoflagellate
433
Gambierdiscus carpenteri in temperate waters of New South Wales, Australia. Harmful
434
Algae 39, 134-145.
435
Kohli, G.S., Neilan, B.A., Brown, M.V., Hoppenrath, M., Murray, S.A., 2014b. Cob gene
436
pyrosequencing enables characterization of benthic dinoflagellate diversity and
437
biogeography. Environ. Microbiol. 16(2), 467-485.
16 438
Kohli, G.S., Papiol, G.G., Rhodes, L.L., Harwood, D.T., Selwood, A., Jerrett, A., Murray,
439
S.A., Neilan, B.A., 2014c. A feeding study to probe the uptake of Maitotoxin by snapper
440
Pagrus auratus. Harmful Algae 37, 125-132.
441
Lee, K.-W., Kang, J.-H., Baek, S.H., Choi, Y.-U., Lee, D.-W., Park, H.-S., 2014. Toxicity of
442
the dinoflagellate Gambierdiscus sp. toward the marine copepod Tigriopus japonicus.
443
Harmful Algae 37, 62-67.
444
Lewis, R.J., 2006. Ciguatera: Australian perspectives on a global problem. Toxicon 48(7),
445
799-809.
446
Lewis, R.J., Holmes, M.J., 1993. Origin and transfer of toxins involved in ciguatera. Comp.
447
Biochem. Physiol., C: Toxicol. Pharmacol. 106(3), 615-628.
448
Lewis, R.J., Holmes, M.J., Alewood, P.F., Jones, A., 1994. Lonspray mass spectrometry of
449
ciguatoxin‐1, maitotoxin‐2 and‐3, and related marine polyether toxins. Nat. Toxins 2(2), 56-
450
63.
451
Litaker, R.W., Vandersea, M.W., Faust, M.A., Kibler, S.R., Chinain, M., Holmes, M.J.,
452
Holland, W.C., Tester, P.A., 2009. Taxonomy of Gambierdiscus including four new species,
453
Gambierdiscus caribaeus, Gambierdiscus carolinianus, Gambierdiscus carpenteri and
454
Gambierdiscus ruetzleri (Gonyaulacales, Dinophyceae). Phycologia 48(5), 344-390.
455
Marchant, H.J., Thomas, D.P., 1983. Polylysine as an adhesive for the attachment of
456
nanoplankton to substrates for electron microscopy. J. Microsc. 131(1), 127-129.
457
Momigliano, P., Sparrow, L., Blair, D., Heimann, K., 2013. The diversity of Coolia
458
spp.(Dinophyceae Ostreopsidaceae) in the central Great Barrier Reef region. PloS one 8(10),
459
e79278.
460
Murray, S., 2010. Benthic dinoflagellates, In: Hallegraeff, G.M., Bolch, C.J.S., Hil, l.D.R.A.,
461
Jameson, I., LeRoi, J.-M., McMinn, A., Murray, S., de Salas, M.F., Saunders, K. (Eds.),
462
Algae of Australia: Phytoplankton flora of temperate coastal waters. CSIRO Publishing,
463
Canberra.
464
Murray, S.A., Garby, T., Hoppenrath, M., Neilan, B.A., 2012. Genetic diversity,
465
morphological uniformity and polyketide production in dinoflagellates (Amphidinium,
466
Dinoflagellata). PloS one 7(6), e38253.
17 467
Nakahara, H., Sakami, T., Chinain, M., Ishida, Y., 1996. The role of macroalgae in
468
epiphytism of the toxic dinoflagellate Gambierdiscus toxicus (Dinophyceae). Phycol. Res.
469
44(2), 113-117.
470
Nishimura, T., Sato, S., Tawong, W., Sakanari, H., Uehara, K., Shah, M.M.R., Suda, S.,
471
Yasumoto, T., Taira, Y., Yamaguchi, H., 2013. Genetic diversity and distribution of the
472
ciguatera-causing dinoflagellate Gambierdiscus spp.(Dinophyceae) in coastal areas of Japan.
473
PloS one 8(4), e60882.
474
Nishimura, T., Sato, S., Tawong, W., Sakanari, H., Yamaguchi, H., Adachi, M., 2014.
475
Morphology of Gambierdiscus scabrosus sp. nov.(Gonyaulacales): a new epiphytic toxic
476
dinoflagellate from coastal areas of Japan. J. Phycol.
477
Paradis, E., Claude, J., Strimmer, K., 2004. APE: analyses of phylogenetics and evolution in
478
R language. Bioinformatics 20(2), 289-290.
479
Parsons, M.L., Settlemier, C.J., Ballauer, J.M., 2011. An examination of the epiphytic nature
480
of Gambierdiscus toxicus, a dinoflagellate involved in ciguatera fish poisoning. Harmful
481
Algae 10(6), 598-605.
482
Rhodes, L., Giménez Papiol, G., Smith, K., Harwood, T., 2014a. Gambierdiscus cf.
483
yasumotoi (Dinophyceae) isolated from New Zealand's sub-tropical northern coastal waters.
484
N. Z. J. Mar. Freshwat. Res.(ahead-of-print), 1-8.
485
Rhodes, L., Smith, K., Papiol, G.G., Adamson, J., Harwood, T., Munday, R., 2014b.
486
Epiphytic dinoflagellates in sub-tropical New Zealand, in particular the genus Coolia
487
Meunier. Harmful Algae 34, 36-41.
488
Richlen, M.L., Morton, S.L., Barber, P.H., Lobel, P.S., 2008. Phylogeography, morphological
489
variation and taxonomy of the toxic dinoflagellateGambierdiscus toxicus(Dinophyceae).
490
Harmful Algae 7(5), 614-629.
491
Ronquist, F., Huelsenbeck, J.P., 2003. MrBayes 3: Bayesian phylogenetic inference under
492
mixed models. Bioinformatics 19(12), 1572-1574.
493
Saburova, M., Polikarpov, I., Al-Yamani, F., 2013. New records of the genus Gambierdiscus
494
in marginal seas of the Indian Ocean. Mar. Biodivers. Rec. 6, e91.
18 495
Thompson, J.D., Higgins, D.G., Gibson, T.J., 1994. CLUSTAL W: improving the sensitivity
496
of progressive multiple sequence alignment through sequence weighting, position-specific
497
gap penalties and weight matrix choice. Nucleic Acids Res. 22(22), 4673-4680.
498
Tosteson, T., 2004. Caribbean ciguatera: a changing paradigm. Rev. Biol. Trop. 52, 109-113.
499
Walsh, P.S., Metzger, D.A., Higuchi, R., 1991. Chelex 100 as a medium for simple extraction
500
of DNA for PCR-based typing from forensic material. BioTechniques 10(4), 506-513.
501
Xu, Y., Richlen, M.L., Morton, S.L., Mak, Y.L., Chan, L.L., Tekiau, A., Anderson, D.M.,
502
2014. Distribution, abundance and diversity of Gambierdiscus spp. from a ciguatera-endemic
503
area in Marakei, Republic of Kiribati. Harmful Algae 34, 56-68.
504 505 506 507 508 509 510 511 512 513 514 515 516 517 518
19 519
Figures Legends
520 521
Figure 1. Map of the sampling sites in Great Barrier Reef (GBR) region, east coast,
522
Queensland, Australia. Inset box, showing the location of the sampling area within Australia.
523
The three regions in which sampling took place: Raine Island in the far northern GBR; Nelly
524
Bay and Pallarenda, Townsville, in the central GBR; Heron Island in the southern GBR.
20
525 526 527
Figure 2. Maximum Likelihood phylogeny of the genus Gambierdiscus based on near-
528
complete 18S gene sequences. Branch support values (ML/BI) represent bootstrap support
529
based on 1000 pseudo-replicate datasets and Bayesian clade credibility values respectively.
530
Branches with 100/100 support are indicated by an asterisk. Scale represents number of
531
changes per nucleotide. Strains isolated in this study are shown in bold.
21
532 533
Figure 3. Maximum Likelihood phylogeny of the genus Gambierdiscus based on the D8-D10
534
region of the 28S rRNA gene. Branch support values (ML/BI) represent bootstrap support
535
based on 1000 pseudo-replicate datasets and Bayesian clade credibility values respectively.
536
Branches with 100/100 support are indicated by an asterisk. Scale represents number of
537
changes per nucleotide. Strains isolated in this study are shown in bold.
538 539 540 541 542 543 544 545
22 546 547
548 549
Figure 4. Principal Coordinate Analyses (PCoA) based on p-distance matrices among
550
sequences within the G. yasumotoi/ G. ruetzleri clade. PCoA based on near complete 18S
551
gene sequences (A & B) and on the D8-D10 (C & D) and D1-D3 (E & F) regions of the 28S
552
rRNA gene. Clones from the same strain are represented by partially transparent symbols
553
with the same colour, allowing visualization of density of points via colour saturation .
554
Eigenvalues of the plotted axes are shown in black.
23
555 556
Figure 5. Scanning electron microscope images of Gambierdiscus yasumotoi from the
557
cultured strain, NQAIF210. Scale bars are shown. A. Apical view, showing the epithecal
558
plate pattern and shape of the epicone. B. The apical pore plate, showing the comma and
559
pores. C. Cell in left lateral view. D. Cell in dorsal view. E. Cell in dorsal view. F. Antapical-
560
ventral view, showing antapical plates and part of the sulcus. G. Antapical view, showing
561
substantial plate overlaps. H. Ventral view, showing the sulcus and a lack of sulcal lists. I.
562
Antapical/lateral view, showing the sulcal lists.
563 564
24
565 566
Figure 6. Scanning electron microscope images of Gambierdiscus yasumotoi from the field
567
material from Raine Island and Heron Island. Scale bars are 10 µm. A. Cell in ventral/lateral
568
view, showing sulcal lists. B. Cell in apical/ventral view, showing part of the epicone. C. Cell
569
in ventral view. D. Cell in right lateral view. E. Cell in left lateral view. F. Cell in left lateral
570
view. G. Cell in apical view, showing the epicone, including the apical pore plate and
571
comma. H. Cell in apical view, showing the epicone. I. Cell in apical view, showing the
572
epicone.
573 574 575
25
576 577
Figure 7. Scanning electron microscope images of Gambierdiscus yasumotoi from the field
578
material from Raine Island and Heron Island, showing the apical pore region, and sulcus.
579
Scale bars represent 10 µm (A) and (B).A. The apical pore plate. B. The sulcal plates,
580
showing the list, and ssa and ssp. C. Showing the sma and sda plates. D. Showing the sma
581
plate and sulcal list.
582
26
Tables Table 1. Cell sizes and comparison of other morphological features of G.yasumotoi, from our data and previous reports.
site
G.
G.
G. cf
G.
G.
G. yasumotoi
G. yasumotoi
G. ruetzleri
yasumotoi
yasumotoi
yasumotoi
yasumotoi
yasumotoi
(Holmes 1998),
(Litaker et al
(Litaker et al
(this study
(this study
(Rhodes et
(Saburova
(Saburova
type description
2009
2009)
field cells)
culture)
al in 2014a
et al 2013
et al 2013
Heron
Townsville northern
Kuwait
Jordan
Island,
New
Raine
Zealand
Singapore
Island Depth
Mean =
Mean = 49
Mean =
Mean=62.9 Mean= 62.6
Mean = 50
Mean = 56.8
Mean = 45.5
57.2
Range 44-
54.8
± 4.4
(n=17)
(n=14)
Range 41-55
Range 54-
54
± 5.7
(n=30)
Mean = 56
Mean= 51
Mean =
Mean=
Mean= 62.9
Mean = 53
Mean= 62.4
Mean = 51.6
Range 53-
Range 49-
59.8 ± 7.5
61.7 ± 6.2,
± 6.8, N=7
(n=17)
(n=14)
45-59.5
59
54
Mean=48
Mean=45
Mean =
Mean=
Mean=
Mean = 44
Mean = 51.7
Mean=37.5
Range 46-
Range 40-
42.5 ± 4.1
54.1 ± 5.1,
54.7± 1.5,
38-50
(n=14)
30.9-42.2
50
49
N=9
N=3
± 6.1, N =7
Range 43-61
59 Length
Width
N = 20
45-63
27
Epitheca
~0.5
~1.0
D:W ratio 1.18
1.08
L:W ratio
1.13
~1.0
0.5-0.6
~1.0
1.29
1.28
1.14
1.13
1.1
1.35
1.41
1.14
1.15
1.20
1.2
1.45
to hypotheca ratio
1.16
28