Mitochondrial cox1 and cob sequence diversities in Gelidium vagum (Gelidiales, Rhodophyta) in Korea
Mitochondrial cox1 and cob sequence diversities in Gelidium vagum (Gelidiales, Rhodophyta) in Korea
ALGAE. 2014. Mar, 29(1): 15-25
Copyright © 2014, The Korean Society of Phycology
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License ( which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Received : February 20, 2014
  • Accepted : March 11, 2014
  • Published : March 15, 2014
Export by style
Cited by
About the Authors
Kyung Ju, Yoon
Department of Biology, Chungnam National University, Daejeon 305-764, Korea
Kyeong Mi, Kim
Department of Biology, Chungnam National University, Daejeon 305-764, Korea
Ga Hun, Boo
Department of Biology, Chungnam National University, Daejeon 305-764, Korea
Kathy Ann, Miller
University Herbarium, 1001 Valley Life Sciences Building 2465, University of California, Berkeley, CA 94720, USA
Sung Min, Boo
Department of Biology, Chungnam National University, Daejeon 305-764, Korea

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.
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.
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
PPT Slide
Lager Image
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; 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. ( Excoffier and Lischer 2010 ). Haplotype diversities ( H d) and nucleotide diversities (π) were measured using DnaSP v.5.10.01 ( Librado and Rozas 2009 ).
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 ).
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
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
PPT Slide
Lager Image
Genetic information of cox1, cob, and cox1 + cob datasets
Diversity measures of cox1, cob, and cox1 + cob datasets forGelidium vagumin Korea and Russia
PPT Slide
Lager Image
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.
PPT Slide
Lager Image
Haplotype network of the combined cox1 + cob dataset for Gelidium vagum.
PPT Slide
Lager Image
Map of the combined cox1 + cob haplotypes for Gelidium vagum. Parentheses indicate the number of specimens with same haplotype.
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.
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.
Boo G. H. , Park J. K. , Boo S. M. 2013 Gelidiophycus (Rhodophyta: Gelidiales): a new genus of marine algae from East Asia Taxon 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. 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 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 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 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. 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 46 1346 - 1349
Kozhenkova S. I. 2009 Retrospective analysis of the marine flora of Vostok Bay, Sea of Japan Russ. J. Mar. Biol. 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 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 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. 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. 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:
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. 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 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 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 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. 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
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