DNA Barcode Examination of Bryozoa(Class: Gymnolaemata) in Korean Seawater
DNA Barcode Examination of Bryozoa(Class: Gymnolaemata) in Korean Seawater
Animal Systematics, Evolution and Diversity. 2011. Jul, 27(2): 159-163
Copyright ©2011, The korean Society of Systematic Zoology
This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License( permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
  • Published : July 31, 2011
Export by style
Cited by
About the Authors
Hyun-Jung Lee
Division of EcoScience, Ewha Womans University, Seoul 120-750, Korea
Ye-Seul Kwan
Division of EcoScience, Ewha Womans University, Seoul 120-750, Korea
So-Ra Kong
Division of EcoScience, Ewha Womans University, Seoul 120-750, Korea
Bum Sik Min
Department of Life Science, The University of Suwon, Hwaseong 445-743, Korea
Ji Eun Seo
Department of Rehabilitation, Woosuk University, Wanju 565-701, Korea
Yong-Jin Won
Division of EcoScience, Ewha Womans University, Seoul 120-750, Korea
DNA barcoding of Bryozoa or “moss animals” has hardly advanced and lacks reference sequences for correct species identification. To date only a small number of cytochrome c oxidase subunit I (COI) sequences from 82 bryozoan species have been deposited in the National Center for Biotechnology Information (NCBI)GenBank and Barcode of Life Data Systems (BOLD). We here report COI data from 53 individual samples of 29 bryozoan species collected from Korean seawater. To our knowledge this is the single largest gathering of COI barcode data of bryozoans to date. The average genetic divergence was estimated as 23.3% among species of the same genus, 25% among genera of the same family, and 1.7% at intraspecific level with a few rare exceptions having a large difference, indicating a possibility of presence of cryptic species. Our data show that COI is a very appropriate marker for species identification of bryozoans, but does not provide enough phylogenetic information at higher taxonomic ranks. Greater effort involving larger taxon sampling for the barcode analyses is needed for bryozoan taxonomy.
Bryozoans or “moss animals” are distributed around the world and presently more than 5,000 species are known to exist.They are fouling animals, which are easily introduced into alien regions via ships and/or ballast waters (Seo, 2005). The attachment of these animals to either naturally or artificially submerged surfaces, such as ship hulls, enables their frequent relocation. Although bryozoans are taxonomically and ecologically important marine invertebrates, DNA barcoding of them are not yet progressed. Studies based on the standard maker of DNA barcoding of most animal species, mitochondrial cytochrome c oxidase subunit I (COI) sequences, are very rare in bryozoans and confined to a few species with,for instance, concerns of phylogeography (Gomez et al., 2007)and invasion (Mackie et al., 2006).
Approximately, 130 species are known to be distributed in Korea; however, morphological identification of these taxa is difficult due to the lack of distinguishing characteristics,particularly among closely related species. For this reason, we expect that there should be a number of specific synonymy.DNA barcoding, therefore, should contribute significantly to the proper identification of bryozoan species. To date DNA barcodes of only 82 bryozoan species have been deposited in the National Center for Biotechnology Information(NCBI) and the database of Barcode of Life Data Systems(BOLD). Thus, there is an urgent need for augmenting the worldwide database of COI barcodes.
We attempted to construct COI barcodes based upon 53 individual samples from 29 species of Bryozoa collected from Korean seawater (Table 1). Genomic DNA was isolated from a small piece of each species using DNeasy Blood and Tissue Kits (QIAGEN Inc., USA) following the manufacturer’s protocol. A partial DNA fragment of mitochondrial COI was first amplified by PCR using the primers LCO1490 and HCO2198 (Folmer et al., 1994); however, the PCR products were not distinguishable on agarose gel. Thus, we designed a nested second PCR reaction using a degenerative BR_bugul_CO1_F (5′-CTG GDA TAR TTG GWA GAG G-3′) and BR_bugul_CO1_R (5′-GTR TTW AAA TTW
Classification of Bryozoa (class Gymnolaemata) and GenBank accession numbers
PPT Slide
Lager Image
Classification of Bryozoa (class Gymnolaemata) and GenBank accession numbers
CGR TCK GTT A-3′) primer set specifically designed for the order Cheilostomatida, of which COI sequences were relatively abundant in GenBank. Fortunately, this new primer set was successful for cross-species amplification of COI from two species of the order Ctenostomata. The PCR products were purified with a LaboPassTM PCR purification kit (Cosmo Genetech Inc., Korea) and sequenced using an ABI 3730xl(Applied Biosystems Inc., USA). When we found mixed COI haplotypes due to a mixture of different individuals of conspecies,a typical problem of colony forming tissues in some specimens, we performed clone sequencing to phase out the ambiguous DNA sequence sites. Finally, pairwise sequence divergences and a neighbor-joining (NJ) tree were investigated with a Kimura-2-parameter (K2P) model using the computer software MEGA 4.0 (Kumar et al., 2004). GenBank accession numbers of the sequences that we deposited from this study and their taxonomical information are provided in Table 1 .
Our COI data was built from the largest taxon sampling to date, representing 29 species, 18 genera, 15 families, and 2 orders of Bryozoa. Assuming that the present COI data were obtained from the recognized species without error, the estimates of pairwise sequence divergence between species are as follows: 25% at a confamilial level, 23.3% at a congeneric taxonomic level, and 1.7% at a conspecific level ( Fig.1 ).
PPT Slide
Lager Image
Distribution of pairwise sequence divergence based on the Kimura-2-parameter for cytochrome c oxidase subunit I sequences according to taxonomic levels. The horizontal axis represents intervals of genetic distance in percentage and the vertical axis is the number of individuals associated with each distance interval.
Although Fig. 1 highlights increasing genetic divergence as taxonomic rank increases, the difference between generic and familial levels is still very small. Therefore, the main purpose of DNA barcoding for species identification has been fulfilled given the current taxon sampling with some exceptionally high divergent conspecific cases which strongly indicates the possibility of the presence of cryptic species in Scrupocellaria diadema, Watersipora subtorquata and Parasmittina contraria . Further examinations are needed for these cases. The NJ tree clearly shows high bootstrapping supports in most terminal nodes, which indicates that COI sequences have sufficient information for the identification of species given the present taxon sampling. Among all 29 species, only two, Biflustra perfragilis and Biflustra crenulata , could not be distinguished due to identical sequences. In contrast, we could not find such high bootstrap supports in most basal nodes of the tree, suggesting a lack of phylogenetic signal at higher taxonomic levels ( Fig.2 ). Furthermore, most families with multiple genera were not grouped as monophyletic clades, and many basal branches show polytomies again indicating a lack of phylogenetic signal.
Our study of bryozoans demonstrates the classic problems of species identification by molecular markers, which suffers from insufficient samplings, superficial knowledge of the range of molecular divergence, cryptic species, and phenotypic plasticity of recognized species (Vrijenhoek, 2009). In some cases, the COI marker did not provide resolution of species identification. As an extreme example, four sequences of Biflustra perfragilis and Biflustra crenulata turned out to be identical, while they are easily classified by distinctive morphology. Therefore, it will be interesting to prove whether the taxonomic key characteristics for the two species result from phenotypic plasticity or the COI marker is just not appropriate for differentiate such closely related species. On the other hand, Scrupocellaria diadema showed a relatively large difference among different individuals (n=4) (Fig. 2).This species is known to have relatively large number of complex key characteristics and show a particularly high degree of morphological variability (Seo, 2005). Therefore, it is possible that this species may be a species complex. Further molecular examinations based on multiple individuals and localities are clearly warranted for this species. In addition,extremely divergent individual sequences were observed from Watersipora subtorquata and Parasmittina contraria, respectively. Interestingly the genus, Watersipora, includes cosmopolitan species and has been the subject of the studies of alien species (Seo, 1999). Because this taxon has a rather
PPT Slide
Lager Image
Neighbor-joining (NJ) tree of Bryozoa based on cytochrome c oxidase subunit I sequences from 54 bryozoan samples. The root was drawn using the midroot method and the asterisks (*) indicate two reference sequences (Securiflustra securifrons FJ196090;and Flustra foliacea FJ196080) of the order Cheilostomatida obtained from GenBank for within group comparison. Numbers on branches show bootstrap values (¤50%) in NJ analysis with 500 replicates.
simple morphology, morphological species identification is difficult and, accordingly, multiple cryptic species could be involved. In contrast, Parasmittina contraria is known to highly variable species, but our COI data strongly suggest a possible presence of cryptic species within this taxon. This species is currently known to distribute only in the Jejudo Islands. The high divergence of Watersipora subtorquata and Parasmittina contraria samples corresponds to the right tail of the distribution of pairwise sequence divergence at the conspecific level ( Fig.1 ). Our DNA barcoding clearly separate the morphologically distinctive species including Calyptotheca wasinensis, Emballotheca pacifica, Celleporaria brunnea, Bugula umbelliformis, Cellaria punctata , and Bantariella bocki ( Fig.2 ). One of the advantages of COI sequences and their comparison could be added by the case of Fenestrulina malusii . Particularly this species causes a considerable difficulty in taxonomic classification due to its ambiguous morphological characters. Thus the resultant divergent COI sequence of it to other recognized bryozoan species can be used as one of many references to identify it in future( Fig.2 ).
It should be noted that one of the most fundamental problems encountered with DNA barcoding of bryozoans was the lack of a universal COI primer set and insufficient reference sequences. During our study, we experienced frequent PCR failures with the universal Folmer primers (Folmer et al., 1994) for many bryozoans (not shown in the present results).This problem is not surprising if we take into account the fact that taxonomically broad bryozoans may have an enormous degree of COI divergence. In fact, the range of divergence and polymorphism of COI in bryozoans is yet unknown. The previously known COI sequences are limited to a very small portion of bryozoans, which actually impedes the design of universal primers. Hopefully, as COI data of bryozoans grow, development of universal primers should be possible. Such group specific oligonucleotide sequences might be desirable to minimize contamination due to nonbryozoan amplification, known as “the peril of universal primers” (Vrijenhoek, 2009). Another experimental difficulty was the frequent observation of mixed sequence reads from well-amplified PCR products, which indicates two possibilities:that is, either a mixture of variable haplotypes originating from a tissue specimen consisted of numerous tiny colonial individuals or a contamination of quite different species.Although the characteristic colonial structure of bryozoans imposes technical difficulties for extracting clean COI sequences from raw specimens. If we have suitable COI data from taxa across broad bryozoans taxa, the problem of mixed sequence reads will be more easily solved by comparing reference sequences well annotated with taxonomic information. In this regard, our current data, except for some suspicious cases,will add better discriminating power to the COI database in GenBank and BOLD for comparative analyses.
In summary, our DNA barcode study involving 29 bryozoan species is the first attempt of its kind in Korea. We found a high degree of COI sequence divergence among most species.The amount of differences was sufficient for species identification in most cases except for two species. We therefore believe that COI can be served as a standard molecular maker for DNA barcoding of bryozoans.
We would like to thank anonymous reviewers for their constructive criticisms that improved the manuscript. This research was supported by the Korea Ministry of Environment as “The Eco-technopia 21 project (Project No. 052-071-066).”
Folmer O , Black M , Hoeh W , Lutz R , Vrijenhoek R 1994 DNA primers for amplification of mitochondrial cytochromecoxidase subunit I from diverse metazoan invertebrates Molecular Marine Biology and Biotechnology 3 294 - 299
Gómez A , Hughes RN , Wright PJ , Carvalho GR , Lunt DH 2007 Mitochondrial DNA phylogeography and mating compatibility reveal marked genetic structuring and speciation in the NE Atlantic bryozoanCelleporella hyalina Molecular Ecology 16 2173 - 2188
Kumar S , Tamura K , Nei M 2004 MEGA3: Integrated software for Molecular Evolutionary Genetics Analysis and sequence alignment Briefings in Bioinformatics 5 150 - 163
Mackie JA , Keough MJ , Christidis L 2006 Invasion patterns inferred from cytochrome oxidase I sequences in three bryozoansBugula neritina Watersipora subtorquataandWatersipora arcuata Marine Biology 149 285 - 295
Seo JE 1999 Taxonomic review of KoreanWatersipora(BryozoaGymnolaemata Cheilostomata) Korean Journal of Systematic Zoology KISTI1.1003/JNL.JAKO199911920829369 15 221 - 229
Seo JE 2005 Illustrated encyclopedia of fauna and flora of Korea. Vol. 40. Bryozoa Ministry of Education and Human Resources Development Seoul 1 - 596
Vrijenhoek RC 2009 Cryptic species phenotypic plasticityand complex life histories: assessing deep-sea faunal diversity with molecular markers Deep Sea Research Part II: Topical Studies in Oceanography 56 1713 - 1723