Mitochondrial cox1 and cob sequence diversities in Gelidium vagum (Gelidiales, Rhodophyta) in Korea

  • cc icon
  • ABSTRACT

    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 cox1 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 cox1 + 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 cox1. 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.


  • KEYWORD

    cob , cox1 , Gelidium , genetic diversity , haplotype network , invasive 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 rbcL sequences (Hughey et al. 1996, Mineur et al. 2012). The Asian origin of G. vagum has been demonstrated by phylogenetic analyses of rbcL, and cox1 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 cox1 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 cox1 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.

    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 cox1 (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 (Hd) 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).

    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.

    Sequence comparisons of the 1,337 base pair (bp) cox1 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 cox1 in the Mallipo population (n = 10) on the west coast were high (Hd, 0.844; π, 0.173), while both indices were relatively low in Gyeokpo (Hd, 0.673; π, 0.060). In the Russian population (n = 12), both diversities were zero (Table 3). Haplotype network of cox1 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).

    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 cox1 coast were high in Namhaedo (Hd, 0.711; π, 0.140), whereas both indices were low in Oeyeondo (Hd, 0.133; π, 0.018) and in Gyeokpo (Hd, 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 cox1 haplotypes, the network of cob haplotypes (Appendix 2) is similar to that of cox1, showing no genetic structure in Korean populations. Haplotype CB10 occurred in the Netherlands and USA (Table 1).

    The combined cox1 + 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 (Hd, 0.889; π, 0.127), while both indices were relatively low in Oeyeondo (Hd, 0.781; π, 0.078) and Gyeokpo (Hd, 0.764; π, 0.047). Both diversities in the Russian population (n = 12) were zero (Table 3).

    Haplotype network of cox1 + 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.

    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 cox1 (16%, 17 / 108 individuals), cob (10%, 11 / 108 individuals), and the combined cox1 + cob (26%, 28 / 108 individuals). Similarly, genetic diversities of cox1 (Hd, 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 (Hd, 0.711) in G. elegans (Kim et al. 2012) and (Hd, 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 cox1 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 rbcL studies (Renfrew et al. 1989, Hughey et al. 1996, Mineur et al. 2012). In addition, the cox1 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 + cox1, 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 cox1 (1.2%) (Table 2) and both demonstrated similar population structures for G. vagum (Appendices 1 & 2). In addition, our use of cob and cox1 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.

  • 1. Boo G. H., Kim K. M., Nelson W. A., Riosmena-Rodriguez R., Yoon K. J., Boo S. M. 2014 Taxonomy and distribution of selected species of the agarophyte genus Gelidium (Gelidiales, Rhodophyta) [J. Appl. Phycol.] google
  • 2. Boo G. H., Park J. K., Boo S. M. 2013 Gelidiophycus (Rhodophyta: Gelidiales): a new genus of marine algae from East Asia [Taxon] Vol.62 P.1105-1116 google
  • 3. 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.] Vol.10 P.564-567 google
  • 4. Geraldino P. J. L., Yang E. C., Boo S. M. 2006 Morphology and molecular phylogeny of Hypnea flexicaulis (Gigartinales, Rhodophyta) from Korea [Algae] Vol.21 P.417-423 google
  • 5. Hughey J. R., Kieldsen C. K., Silva P. C., Moe R. L., DeCew T. C. 1996 Noteworthy collections: California [Madrono] Vol.43 P.432 google
  • 6. Hughey J. R., Miller K. A., Lyman A. 2009 Noteworthy collections: Dasya sessilis, Grateloupia lanceolata, Grateloupia turuturu, and Neosiphonia harveyi [Madrono] Vol.56 P.293-295 google
  • 7. 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] Vol.408 P.157-158 google
  • 8. Kang J. W. 1966 On the geographical distribution of marine algae in Korea [Bull. Pusan Fish. Coll.] Vol.7 P.1-125 google
  • 9. 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.] Vol.98 P.27-33 google
  • 10. 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] Vol.46 P.1346-1349 google
  • 11. Kozhenkova S. I. 2009 Retrospective analysis of the marine flora of Vostok Bay, Sea of Japan [Russ. J. Mar. Biol.] Vol.35 P.263-278 google
  • 12. Lee H. B., Abbott I. A. 1994 Some species of Gelidium (Gelidiales, Rhodophyta) from Korea P.67-79 google
  • 13. Lee H. B., Kim J. I., Abbott I. A. 1995 Notes on Gelidiales species from Korea P.161-174 google
  • 14. Lee Y. P. 1988 Taxonomic studies on the Gelidiaceae (Rhodophyta) in Cheju Island. I. Some members of Gelidium [Korean J. Plant Taxon] Vol.18 P.95-113 google
  • 15. Librado P., Rozas J. 2009 DnaSP v5: a software for comprehensive analysis of DNA polymorphism data [Bioinformatics] Vol.25 P.1451-1452 google
  • 16. Miller K. A., Aguilar-Rosas L. E., Pedroche F. F. 2011 A review of non-native seaweeds from California, USA and Baja California, Mexico [Hidrobiologica] Vol.21 P.365-379 google
  • 17. Mineur F., Le Roux A., Stegenga H., Verlaque M., Maggs C. A. 2012 Four new exotic red seaweeds on European shores [Biol. Invasions] Vol.14 P.1635-1641 google
  • 18. Okamura K. 1934 On Gelidium and Pterocladia of Japan [J. Imp. Fish. Inst.] Vol.29 P.47-67 google
  • 19. 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.] Vol.29 P.216-222 google
  • 20. Patwary M. U., van der Meer J. P. 1994 Application of RAPD markers in an examination of heterosis in Gelidium vagum (Rhodophyta) [J. Phycol.] Vol.30 P.91-97 google
  • 21. Provan J., Murphy S., Maggs C. A. 2005 Tracking the invasive history of the green alga Codium fragile ssp. tomentosoides [Mol. Ecol.] Vol.14 P.189-194 google
  • 22. Rambaut A. 2002 Se-Al sequence alignment editor, version 2.0a11 google
  • 23. 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.] Vol.67 P.3295-3314 google
  • 24. 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] Vol.28 P.31-43 google
  • 25. Sohn C. H., Kang J. W. 1978 The classification of the family Gelidiaceae (Rhodophyta) in Korea [Publ. Inst. Mar. Sci. Nat. Fish. Univ. Busan] Vol.11 P.29-40 google
  • 26. 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] Vol.45 P.687-695 google
  • 27. van der Meer J. P., Patwary M. U. 1991 Genetic alleviation of the self-fertilization complication when hybridizing monoecious Gelidium vagum [Hydrobiologia] Vol.221 P.167-179 google
  • 28. 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.] Vol.102 P.5432-5437 google
  • 29. 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] google
  • 30. Yoshida T. 1998 Marine algae of Japan P.1222 pp google
  • 31. Zhang J. F., Xia E., Abbott I. A. 1988 Chinese species of Gelidium Lamouroux and other Gelidiales (Rhodophyta) with key, list and distribution of common species P.109-113 google
  • [Table 1.] Collection information, voucher code, GenBank accession number, and haplotypes for specimens of Gelidium vagum used in this study
    Collection information, voucher code, GenBank accession number, and haplotypes for specimens of Gelidium vagum used in this study
  • [Fig. 1.] 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.
    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.
  • [Fig. 2.] 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.
    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.
  • [Table 2.] Genetic information of cox1, cob, and cox1 + cob datasets
    Genetic information of cox1, cob, and cox1 + cob datasets
  • [Table 3.] Diversity measures of cox1, cob, and cox1 + cob datasets for Gelidium vagum in Korea and Russia
    Diversity measures of cox1, cob, and cox1 + cob datasets for Gelidium vagum in Korea and Russia
  • [Fig. 3.] Haplotype network of the combined cox1 + cob dataset for Gelidium vagum.
    Haplotype network of the combined cox1 + cob dataset for Gelidium vagum.
  • [Fig. 4.] Map of the combined cox1 + cob haplotypes for Gelidium vagum. Parentheses indicate the number of specimens with same haplotype.
    Map of the combined cox1 + cob haplotypes for Gelidium vagum. Parentheses indicate the number of specimens with same haplotype.