The number of species of non-native and invasive marine algae is growing, with concomitant public concern about native ecosystems and coastlines.
Gelidium vagum
, recently introduced from northeast Asia to Europe and North America, commonly occurs from the intertidal to subtidal zones in Korea, China, and Japan. To investigate the level of genetic diversity of native populations, we analyzed mitochondrial
cox
1 and
cob
from 108 specimens of
G. vagum
from Korea, China, eastern Russia, including from the Netherlands and USA. The haplotype network of individual and
cox
1 +
cob
datasets revealed no genetic structure in local populations, suggesting genetic flow between Korean populations. Our results corroborate a typical pattern of genetic diversity for introduced species, with low levels in introduced populations and high levels in native populations. All haplotypes were shared between the Netherlands and USA, but not between Korea and the Netherlands / USA except
cox
1. Additional sampling will identify donor populations in native northeast Asian waters. This is the first report of the utility of the mitochondrial coding
cob
sequences in red algae.
INTRODUCTION
The number of species of non-native and invasive marine algae is growing, with concomitant public concern about native ecosystems and coastlines. One of the interesting biological questions is how these species are extending their natural ranges despite low genetic diversity (e.g.,
Jousson et al. 2000
,
Provan et al. 2005
,
Voisin et al. 2005
,
Kim et al. 2010
). However, the genetic structure of native populations has been less intensively studied than that of the non-native ones, and may provide clues to the mechanism of invasion.
The red alga
Gelidium vagum
Okamura commonly occurs in Korea, Japan, China and far eastern Russia (
Kang 1966
,
Zhang and Xia 1988
,
Yoshida 1998
,
Kozhenkova 2009
). The species was first described on the basis of intertidal specimens from western coast of Hokkaido to Iwaki Province (
Okamura 1934
). The type locality, although not noted in the protologue (
Okamura 1934
), was later designated as Onahama Harbor in Fukushima Prefecture, Japan (
Yoshida 1998
). In Korea,
G. vagum
was first reported at Kiseong on the east coast, Odongdo in south, and Huksando in west (
Kang 1966
); its distribution has since been extended to Jeju and Ocheongdo (
Sohn and Kang 1978
,
Lee 1988
,
1994
,
Lee and Kim 1995
).
Renfrew et al. (1989)
, relying morphological characters, reported the introduction of
Gelidium vagum
to British Columbia, Canada. Further introductions of this species in California and Washington, USA, and the Netherlands have been confirmed by comparing
rbc
L sequences (
Hughey et al. 1996
,
Mineur et al. 2012
). The Asian origin of
G. vagum
has been demonstrated by phylogenetic analyses of
rbc
L, and
cox
1 sequences (
Boo et al. 2014
).
Gelidium vagum
serves as an excellent model system for the study of heterosis in red algae because fronds can readily be grown under laboratory conditions and the sexual reproductive cycle can be completed in a few weeks under optimum conditions (
van der Meer and Patwary 1991
,
Patwary et al. 1993
,
Patwary and van der Meer 1994
). Although its mitochondrial genome was fully analyzed (
Yang et al. 2014
), little information on the genetic diversity of
G. vagum
in its native range is available.
Mitochondrial-encoded
cox
1 is a fast-evolving gene that has recently proved valuable for barcoding red algal species (
Saunders and Moore 2013
) and for revealing the biogeographic structure of populations (
Kim et al. 2010
,
2012
). Cytochrome
b
(
cob
), one of the 12 markers undergoing tests to determine its utility for phylogenetic studies by the Research Group for the Red Algal Tree of life (
Saunders and Moore 2013
), has potential for assessing genetic variation within species. The aim of this study was to contribute to an understanding of the phylogeography of this invasive alga by investigating haplotype and nucleotide diversities of both
cox
1 and
cob
genes in
Gelidium vagum
from Korea.
MATERIALS AND METHODS
Specimens in this study are deposited in the Herbarium of Chungnam National University, Daejeon, Korea. A total of 108 individuals was collected in Korea, China, far eastern Russia, the Netherlands, and USA (
Table 1
). Field-collected specimens were desiccated in silica gel. Epiphytes on specimens were removed under a dissecting microscope before DNA extraction.
Collection information, voucher code, GenBank accession number, and haplotypes for specimens ofGelidium vagumused in this study
Collection information, voucher code, GenBank accession number, and haplotypes for specimens of Gelidium vagum used in this study
DNA extraction, polymerase chain reaction amplification, and sequencing procedures followed
Boo et al. (2013)
. Specific primers for the amplification and sequencing were COXI43F and COXI1549R for
cox
1 (
Geraldino et al. 2006
) and CB44F (ATT TAA TWG ATT ATC CWA CDC C) and CB1006R (CAG GCA TWC CRC CWA TYC AKC) for
cob
(
Saunders and Moore 2013
). All sequences of the forward and reverse strands were determined for all taxa, and the electropherograms were edited using the program Chromas v.1.45 (Conor McCarthy, Griffith University, Southport, Queensland, Australia; http://www.technelysium.com.au/chromas.html) and rechecked manually for consistency. Nucleotide sequences were aligned using Se-Al v.2.0a11 (
Rambaut 2002
).
The minimum spanning network to analyze relationships and geographical distributions of each haplotype was constructed using Arlequin v.3.5.1.2 (
Excoffier and Lischer 2010
). Haplotype diversities (
H
d) and nucleotide diversities (π) were measured using DnaSP v.5.10.01 (
Librado and Rozas 2009
).
RESULTS
Our collections from Korea matched the morphological features of
Gelidium vagum
from Russia (
Fig. 1A
&
B
), the Netherlands and USA. Thalli are erect, up to 10 cm in height, pinnately to irregularly branched, with broad but very thin main axes and branches. The cortex is composed of round cells, and the medulla consists of a network of elongate cells with a few rhizoidal filaments in the outer medulla (figure not shown). Tetrasporangial stichidia, formed on terminal or lateral branches (
Fig. 1A
), are compressed and stipitate, with rounded tips. Tetrasporangia are in regular, acropetally developed rows, and are decussately or cruciately divided. Gametophytes are monoecious, forming cystocarps and spermatangia on the same individuals (
Fig. 1B
).
Morphology of Gelidium vagum. (A) Tetrasporangial thallus collected in Yangyang on the east coast on July 5, 2008. (B) Gametophyte bearing cystocarps and spermatangia in the same branchlets collected in Andreevka, Russia on July 3, 2013. Scale bars represent: A & B, 1 cm.
Thalli grew on rock or on the gastropod,
Acmaea pallida
at Oeyeondo Island (
Fig. 2
), in intertidal and subtidal zones, and sometimes were abundant on shaded rocks. Tetrasporangial and gametangial thalli were found from June to August in Korea and eastern Russia.
Gelidium vagum (approximately 4 cm in height) attached on different shells. (A) Cystocarpic thallus attached on Acmaea pallida at Oeyeondo Island on the west coast of Korea on July 28, 2013. (B) Thallus attached on Acmaea pallida at Uihangri, Taean on July 22, 2013.
Sequence comparisons of the 1,337 base pair (bp)
cox
1 fragment from 108 individuals revealed 16 polymorphic sites (1.2%), defining 17 haplotypes (
Table 2
). Ten haplotypes (59%) were ‘private’ (unique to a single region). Haplotype and nucleotide diversities of
cox
1 in the Mallipo population (n = 10) on the west coast were high (
H
d, 0.844; π, 0.173), while both indices were relatively low in Gyeokpo (
H
d, 0.673; π, 0.060). In the Russian population (n = 12), both diversities were zero (
Table 3
). Haplotype network of
cox
1 revealed no genetic structure in Korean populations of
G. vagum
, suggesting genetic flow among local populations of the species (Appendix 1). Haplotype 15 occurred in Yerseke, the Netherlands, Tomales Bay and Puget Sound, USA, as well as several places in Korea (
Table 1
).
Genetic information ofcox1,cob, andcox1 +cobdatasets
Genetic information of cox1, cob, and cox1 + cob datasets
Diversity measures of cox1, cob, and cox1 + cob datasets forGelidium vagumin Korea and Russia
n, number of individuals; h, number of haplotype; S, number of polymorphic sites; Hd, haplotype diversity; π, nucleotide diversity.
Sequence comparisons of the 761 bp
cob
fragment from 108 individuals revealed 11 polymorphic sites (1.5%), defining 11 haplotypes (
Table 2
). Haplotype and nucleotide diversities of
cox
1 coast were high in Namhaedo (
H
d, 0.711; π, 0.140), whereas both indices were low in Oeyeondo (
H
d, 0.133; π, 0.018) and in Gyeokpo (
H
d, 0.182; π, 0.024). Russian population (n = 12) was zero for both diversities (
Table 3
). Although the number of
cob
haplotypes was smaller than that of
cox
1 haplotypes, the network of
cob
haplotypes (Appendix 2) is similar to that of
cox
1, showing no genetic structure in Korean populations. Haplotype CB10 occurred in the Netherlands and USA (
Table 1
).
The combined
cox
1 +
cob
dataset (2,098 bp) revealed 27 polymorphic sites (1.3%), defining 28 haplotypes. The number of haplotypes was high on the west coast (32%, 21 haplotypes / 65 individuals) and relatively low on the south (26%, 4 / 15). Haplotype and nucleotide diversities were high in Mallipo (
H
d, 0.889; π, 0.127), while both indices were relatively low in Oeyeondo (
H
d, 0.781; π, 0.078) and Gyeokpo (
H
d, 0.764; π, 0.047). Both diversities in the Russian population (n = 12) were zero (
Table 3
).
Haplotype network of
cox
1 +
cob
(
Fig. 3
) was similar to that of individual networks, showing no genetic structure in Korean populations of
G. vagum
. Haplotype C1, placed in the center of a star-like structure, commonly occurred on the west and south coasts of Korea (
Fig. 4
). All specimens from Russia belonged to haplotype C1. Haplotype C27 occurred in the Netherlands and USA, but not in northeast Asian waters.
Haplotype network of the combined cox1 + cob dataset for Gelidium vagum.
Map of the combined cox1 + cob haplotypes for Gelidium vagum. Parentheses indicate the number of specimens with same haplotype.
DISCUSSION
Gelidium vagum
, a monoecious species (
Renfrew et al. 1989
,
van der Meer and Patwary 1991
), may be capable of selfing, resulting in reduced genetic variation within populations. However, Korean populations revealed a high number of haplotypes in
cox
1 (16%, 17 / 108 individuals),
cob
(10%, 11 / 108 individuals), and the combined
cox
1 +
cob
(26%, 28 / 108 individuals). Similarly, genetic diversities of
cox
1 (
H
d, 0.809 ± 0.026; π, 0.109 ± 0.009) for
G. vagum
were relatively high. Our results are consistent with previous reports (
Patwary et al. 1993
) that
G. vagum
isolates have sufficient genetic diversity for conducting a heterosis experiment using random amplified polymorphic DNA methods. In contrast, genetic diversities were relatively low (
H
d, 0.711) in
G. elegans
(
Kim et al. 2012
) and (
H
d, 0.757) in
Gracilaria vermiculophylla
(Ohmi) Papenfuss (
Kim et al. 2010
). The latter two species are dioecious, species, a trait which may account for increasing heterozygosity and variability. However, because red algae may often reproduce asexually, the relationship between genetic diversity and monoecy or dioecy is not clear within natural populations.
Recent studies on the genetic diversity of introduced marine algae have revealed invasive haplotype(s) that are shared between native and introduced populations (
Jousson et al. 2000
,
Provan et al. 2005
,
Voisin et al. 2005
,
Uwai et al. 2006
,
Kim et al. 2010
). Our
cox
1 dataset suggests that the introductions of
Gelidium vagum
to Europe and north America originated from populations in Korea and / or Japan, as shown in previous morphological and
rbc
L studies (
Renfrew et al. 1989
,
Hughey et al. 1996
,
Mineur et al. 2012
). In addition, the
cox
1 dataset revealed that the haplotype CX15, shared by populations in Korea, USA, and the Netherlands, is an invasive haplotype. However, because marine algal introductions might have occurred more than once (
Voisin et al. 2005
,
Hughey et al. 2009
), further sampling to search for other invasive haplotypes is necessary to determine whether
G. vagum
was introduced via single or multiple events.
The number of private haplotypes and both haplotype and nucleotide diversities were higher on the west coast of Korea than on the other coasts. Haplotype C27 was found in the Netherlands and USA, where haplotype and genetic diversities were zero. Our results corroborate the generalization that the genetic diversity of invasive species is low in introduced populations relative to that in native populations. Because the invasive haplotypes identified here, for example, CB10 in
cob
and C27 in
cob
+
cox
1, were not found in Korea, we cannot identify a specific donor region in the combined data without further study.
The vector for introduction of
Gelidium vagum
to Canada, USA, and the Netherlands may be the import of oysters, as shown in many cases of invasive macroalgae elsewhere (e.g.,
Voisin et al. 2005
,
Miller et al. 2011
), because the Hornby Island-Denman Island area, Canada and Oosterschelde estuary, the Netherlands, are hotspots for oyster aquaculture (
Renfrew et al. 1989
,
Mineur et al. 2012
). Other mollusks may prove to be vectors because
G. vagum
is often found on
Acmaea pallida
, as seen in
Fig. 2
.
Acmaea pallida
is commonly found on the hulls of ships, and may survive through long voyages.
This is the first report to analyze the mitochondrial
cob
gene as a molecular marker for red algae. The genetic distance values of
cob
(1.5%) were similar to those of
cox
1 (1.2%) (
Table 2
) and both demonstrated similar population structures for
G. vagum
(Appendices 1 & 2). In addition, our use of
cob
and
cox
1 genes reveals a very high number of haplotypes, compared to that of each individual gene, and thus increases the resolution of our haplotype network for
G. vagum
. The
cob
gene is a new marker for population studies of red algae, and its use can be extended to studies of other biogeographically interesting species of benthic algae. The analysis of more than one gene provides a more accurate picture of genetic diversity in marine algae and key information for conservation.
Acknowledgements
We thank Robert Waaland and Herre Stegenga for providing samples from Seattle and Netherlands and we appreciate the help of Tatiana Krupnova and Oksana Sergeeva during our collections in Russia. Il Hoe Kim assisted with the identification of gastropods and provided helpful discussion, for which we are grateful. This work was supported by Basic Science Grant (2013-0699) of Korean Research Foundation and Research Fund of Chungnam National University, Daejeon, Korea.
Boo G. H.
,
Kim K. M.
,
Nelson W. A.
,
Riosmena-Rodríguez R.
,
Yoon K. J.
,
Boo S. M.
2014
Taxonomy and distribution of selected species of the agarophyte genus Gelidium (Gelidiales, Rhodophyta)
J. Appl. Phycol.
http://dx.doi.org/10.1007/s10811-013-0111-7
Boo G. H.
,
Park J. K.
,
Boo S. M.
2013
Gelidiophycus (Rhodophyta: Gelidiales): a new genus of marine algae from East Asia
Taxon
http://dx.doi.org/10.12705/626.7
62
1105 -
1116
Excoffier L.
,
Lischer H. E. L.
2010
Arlequin suite ver. 3.5: a new series of program to perform population genetics analyses under Linux and Windows
Mol. Ecol. Resour.
http://dx.doi.org/10.1111/j.1755-0998.2010.02847.x
10
564 -
567
Geraldino P. J. L.
,
Yang E. C.
,
Boo S. M.
2006
Morphology and molecular phylogeny of Hypnea flexicaulis (Gigartinales, Rhodophyta) from Korea
Algae
http://dx.doi.org/10.4490/ALGAE.2006.21.4.417
21
417 -
423
Hughey J. R.
,
Kieldsen C. K.
,
Silva P. C.
,
Moe R. L.
,
DeCew T. C.
1996
Noteworthy collections: California
Madrono
43
432 -
Hughey J. R.
,
Miller K. A.
,
Lyman A.
2009
Noteworthy collections: Dasya sessilis, Grateloupia lanceolata, Grateloupia turuturu, and Neosiphonia harveyi
Madroño
http://dx.doi.org/10.3120/0024-9637-56.4.293
56
293 -
295
Jousson O.
,
Pawlowski J.
,
Zaninetti L.
,
Zechman F. W.
,
Dini F.
,
Di Guiseppe G.
,
Woodfield R.
,
Millar A.
,
Meinesz A.
2000
Invasive alga reaches California
Nature
http://dx.doi.org/10.1038/35041623
408
157 -
158
Kang J. W.
1966
On the geographical distribution of marine algae in Korea
Bull. Pusan Fish. Coll.
7
1 -
125
Kim K. M.
,
Hoarau G. G.
,
Boo S. M.
2012
Genetic structure and distribution of Gelidium elegans (Gelidiales, Rhodophyta) in Korea based on mitochondrial cox1 sequence data
Aquat. Bot.
http://dx.doi.org/10.1016/j.aquabot.2011.12.005
98
27 -
33
Kim S. Y.
,
Weinberger F.
,
Boo S. M.
2010
Genetic data hint at a common donor region for invasive Atlantic and Pacific populations of Gracilaria vermiculophylla (Gracilariales, Rhodophyta)
J. Phycol
http://dx.doi.org/10.1111/j.1529-8817.2010.00905.x
46
1346 -
1349
Kozhenkova S. I.
2009
Retrospective analysis of the marine flora of Vostok Bay, Sea of Japan
Russ. J. Mar. Biol.
http://dx.doi.org/10.1134/S1063074009040014
35
263 -
278
Lee H. B.
,
Abbott I. A.
1994
Taxonomy of Economic Seaweeds with Reference to Some Pacific and Caribbean Species
California Sea Grant College System
La Jolla, PA
Some species of Gelidium (Gelidiales, Rhodophyta) from Korea
67 -
79
Lee H. B.
,
Kim J. I.
,
Abbott I. A.
1995
Taxonomy of Economic Seaweeds with Reference to Some Pacific Species
California Sea Grant College System
La Jolla, PA
Notes on Gelidiales species from Korea
161 -
174
Lee Y. P.
1988
Taxonomic studies on the Gelidiaceae (Rhodophyta) in Cheju Island. I. Some members of Gelidium
Korean J. Plant Taxon
18
95 -
113
Librado P.
,
Rozas J.
2009
DnaSP v5: a software for comprehensive analysis of DNA polymorphism data
Bioinformatics
http://dx.doi.org/10.1093/bioinformatics/btp187
25
1451 -
1452
Miller K. A.
,
Aguilar-Rosas L. E.
,
Pedroche F. F.
2011
A review of non-native seaweeds from California, USA and Baja California, Mexico
Hidrobiológica
21
365 -
379
Mineur F.
,
Le Roux A.
,
Stegenga H.
,
Verlaque M.
,
Maggs C. A.
2012
Four new exotic red seaweeds on European shores
Biol. Invasions
http://dx.doi.org/10.1007/s10530-012-0186-0
14
1635 -
1641
Okamura K.
1934
On Gelidium and Pterocladia of Japan
J. Imp. Fish. Inst.
29
47 -
67
Patwary M. U.
,
MacKay R. M.
,
van der Meer J. P.
1993
Revealing genetic markers in Gelidium vagum (Rhodophyta) through the random amplified polymorphic DNA (RAPD) technique
J. Phycol.
http://dx.doi.org/10.1111/j.0022-3646.1993.00216.x
29
216 -
222
Patwary M. U.
,
van der Meer J. P.
1994
Application of RAPD markers in an examination of heterosis in Gelidium vagum (Rhodophyta)
J. Phycol.
http://dx.doi.org/10.1111/j.0022-3646.1994.00091.x
30
91 -
97
Provan J.
,
Murphy S.
,
Maggs C. A.
2005
Tracking the invasive history of the green alga Codium fragile ssp. tomentosoides
Mol. Ecol.
14
189 -
194
Rambaut A.
2002
Se-Al sequence alignment editor, version 2.0a11
Institute of Evolutionary Biology
Edinburgh, UK
Available from: http://tree.bio.ed.ac.uk/software/
Renfrew D. E.
,
Gabrielson P. W.
,
Scagel R. F.
1989
The marine algae of British Columbia, northern Washington, and southeast Alaska: division Rhodophyta (red algae), class Rhodophyceae, order Gelidiales
Can. J. Bot.
http://dx.doi.org/10.1139/b89-406
67
3295 -
3314
Saunders G. W.
,
Moore T. E.
2013
Refinements for the amplification and sequencing of red algal DNA barcode and RedToL phylogenetic markers: a summary of current primers, profiles and strategies
Algae
http://dx.doi.org/10.4490/algae.2013.28.1.031
28
31 -
43
Sohn C. H.
,
Kang J. W.
1978
The classification of the family Gelidiaceae (Rhodophyta) in Korea
Publ. Inst. Mar. Sci. Nat. Fish. Univ. Busan
11
29 -
40
Uwai S.
,
Nelson W.
,
Neill K.
,
Wang W. D.
,
Aguilar-Rosas L. E.
,
Boo S. M.
,
Kitayama T.
,
Kawai H.
2006
Genetic diversity in Undaria pinnatifida (Laminariales, Phaeophyceae) deduced from mitochondrial genes: origins and succession of introduced populations
Phycologia
http://dx.doi.org/10.2216/05-66.1
45
687 -
695
van der Meer J. P.
,
Patwary M. U.
1991
Genetic alleviation of the self-fertilization complication when hybridizing monoecious Gelidium vagum
Hydrobiologia
http://dx.doi.org/10.1007/BF00028373
221
167 -
179
Voisin M.
,
Engel C. R.
,
Viard F.
2005
Differential shuffling of native genetic diversity across introduced regions in a brown alga: aquaculture vs. maritime traffic effects
Proc. Natl. Acad. Sci. U. S. A.
http://dx.doi.org/10.1073/pnas.0501754102
102
5432 -
5437
Yang E. C.
,
Kim K. M.
,
Boo G. H.
,
Lee J. H.
,
Boo S. M.
,
Yoon H. S.
2014
Complete mitochondrial genome of the agarophyte red alga Gelidium vagum (Gelidiales)
Mitochondrial DNA
http://dx.doi.org/10.3109/19401736.2013.800499
Yoshida T.
1998
Marine algae of Japan
Uchida Rokakuho Publishing Co., Ltd.
Tokyo
1222 pp -
Zhang J. F.
,
Xia E.
,
Abbott I. A.
1988
Taxonomy of Economic Seaweeds with Reference to Some Pacific and Caribbean Species
California Sea Grant College System
La Jolla, CA
Chinese species of Gelidium Lamouroux and other Gelidiales (Rhodophyta) with key, list and distribution of common species
109 -
113