Next Article in Journal
The Self-Efficacy of Preservice Physical Education Teachers in Disabilities Education in China
Previous Article in Journal
Information and Communication Technology Solutions for the Circular Economy
Previous Article in Special Issue
Variations in the Temporal and Spatial Distribution of Microalgae in Aquatic Environments Associated with an Artificial Weir
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Comparative Analysis of Sequence Polymorphism in Complete Organelle Genomes of the ‘Golden Tide’ Seaweed Sargassum horneri between Korean and Chinese Forms

1
Molecular Ecology and Evolution Laboratory, Department of Biological Science, Sangji University, Wonju 26339, Korea
2
Department of Biological Science, Sangji University, Wonju 26339, Korea
3
Oceanic Climate and Ecology Research Division, National Institute of Fisheries Science, Busan 46083, Korea
4
Estuarine and Coastal Ecology Laboratory, Department of Marine Life Sciences, Jeju National University, Jeju 63243, Korea
5
Department of Oceanography, Faculty of Earth Systems and Environmental Sciences, Chonnam National University, Gwangju 61186, Korea
6
Department of Marine Sciences, Incheon National University, Incheon 22012, Korea
*
Author to whom correspondence should be addressed.
Sustainability 2020, 12(18), 7280; https://0-doi-org.brum.beds.ac.uk/10.3390/su12187280
Submission received: 29 May 2020 / Revised: 28 August 2020 / Accepted: 30 August 2020 / Published: 4 September 2020
(This article belongs to the Special Issue Harmful Organisms and their Management for Sustainable Environment)

Abstract

:
Drifting and inundating brown seaweed Sargassum horneri biomass is called “golden tide”, as it resembles golden massive algal blooms like green tides. This phenomenon occurs globally and its serious ecological impacts on coastal ecosystems have recently begun to be paid attention to. In the present study, by sequencing whole organelle genomes of Korean indigenous S. horneri, we aimed to develop novel molecular markers that can be used for differentiating indigenous from nonindigenous individuals. To this end, we analyzed sequence polymorphisms in mitochondrial (mt) and chloroplast (cp) genomes of two Korean benthic samples in comparison to Chinese ones as a reference. We mapped mt genomes of 34,620~34,628 bp and cp genomes of 123,982~124,053 bp for the Korean samples. In comparative analyses, mtDNA cytochrome c oxidase subunit II (cox2) gene showed the highest number of single nucleotide polymorphisms (SNPs) between Korean and Chinese individuals. NADH dehydrogenase subunit 7 (Nad7)-proline tRNA (trnP) intergenic spacer (IGS) in the mt genome showed a 14 bp insertion or deletion (indel) mutation. For the cp genome, we found a total of 54 SNPs, but its overall evolution rate was approximately four-fold lower than the mt genome. Interestingly, analysis of Ka/Ks ratio in the cp genome revealed a signature of positive selection on several genes, although only negative selection prevalent in mt genome. The ‘candidate’ genetic markers that we found can be applied to discriminate between Korean indigenous and nonindigenous individuals. This study will assist in developing a molecular-based early detection method for effectively managing nonindigenous S. horneri in Korean waters.

1. Introduction

Drifting and inundating brown macroalgae Sargassum biomass followed by its local blooming is called “golden tide”, as it resembles golden massive shoaling algal blooms like green tides. This phenomenon takes place globally and annually, and thus the urgency of its damaging impacts on the local economy and coastal ecosystems has recently begun to be paid attention to [1]. Its periodical outbreaks since 2013 have been reported in many parts of the world including Mexico, the Caribbean, West Africa, and Northeast Asia (northwestern Pacific), including Japan, China and Korea [1,2]. In the Atlantic Ocean and the Caribbean Sea, two pelagic Sargassum species, such as Sargassum natans and Sargassum fluitans, are the most predominant, continuously disperse with ocean currents, and occasionally accumulate massively on beaches, causing serious ecological and economic damages since the 1980s [1,2]. By comparison, in the northwestern Pacific, a benthic species Sargassum horneri is known to be only species contributing to the golden tides [3,4]. Since the first report of golden tides in the northwestern Pacific in early 1960s [5], there have been several golden tide formations observed in the East China Sea with no actual damage to coastal environments until early 2010s [6,7]. Sargassum horneri biomass drifting in the East China Sea and the Yellow Sea has recently inflowed into the southern coast of Korea with ocean currents. Small scale floating S. horneri patches were observed sporadically before 2015, but devastating effects of the golden tides (e.g., disturbing fisheries and tourism) from the East China Sea or the Yellow Sea have been ascertained since 2015, particularly on the South Sea and Jeju Island in Korea (Figure 1) [4,8,9].
The genus Sargassum C. Agardh (Sargassaceae, Phaeophyceae) is comprised of 537 species worldwide [10] and they typically occur from subtidal to intertidal environments. A majority of these species undertake a benthic life with holdfasts attached to the substrate. In a natural state, S. horneri and other Sargassum species provide valuable ecosystem functions and services for a diverse array of other coexisting marine species, serving as nurseries and spawning grounds and refuges. Sargassum horneri is dioecious and sexually reproduce with an obligate annual life cycle [11]. A recent demographic study of the two Korean populations of S. horneri showed that thalli recruited in autumn, grew in winter and senesced in summer [9]. When individuals of S. horneri detached from the substrate by the drag force of waves and/or currents, they can float and drift on the surface of the ocean by gas vesicles [12]. High mortality rate driven by S. horneri thalli detached from the substrate was found in winter and also in early spring [9]. Once detached, drifting individuals could damage coastal ecosystems if they form massive patches through rapid growth and propagation [13]. Nevertheless, drifting massive seaweed rafts in the open ocean sometimes offer refuge habitats for pelagic fishes and invertebrate fauna [14].
For possible management and control efforts of golden tides in Korean waters, early and accurate detection of nonindigenous individuals will be crucial for the rapid removal, in order to prevent potential hybridization incidences between these nonindigenous and Korean indigenous S. horneri individuals. Developing a molecular-based early detection technique would therefore be one of the options for their possible control. In our previous study [4], we found that there were significant differences in the genetic structure between floating (nonindigenous) and Korean benthic (indigenous) populations based on seven polymorphic microsatellite loci. Nevertheless, analyses of mitochondrial (mt) DNA cytochrome c oxidase subunit III (cox3) gene could not differentiate indigenous from nonindigenous individuals, due to the fact that this marker showed insufficient polymorphism. By considering geographic distribution of cox3 haplotypes determined in previous studies [4,15] along with the direction of major oceanic currents (e.g., the Kuroshio Current) surrounding this region, we were only able to speculate that S. horneri biomass floating to the Korean coast was likely to originate from Zhoushan, Zhejiang province of China [4]. Since the geographic origin(s) of floating S. horneri could not be firmly determined with mtDNA-based phylogeographic analysis, we, here, set out exploring genetic variation at whole organelle-genomic levels by comparing mt and chloroplast (cp) genomes between Korean and Chinese S. horneri. These analyses will help to differentiate between indigenous and nonindigenous individuals and further determine the origins of nonindigenous individuals more precisely.
The mt and cp genomes of the genus Sargassum have provided important insights into understanding their inter- and intraspecific phylogenetic relationships and evolutionary histories [16,17]. Given uniparental inheritance and accelerated evolutionary rates, mtDNA may provide a higher resolution in elucidating phylogenetic relationships among closely related taxa in brown algae [18]. The cpDNA would also help to advance our understanding of phylogenetic relationships and evolutionary processes that have not been fully resolved, particularly for terrestrial plants and seaweeds [19]. Sequence polymorphisms in mt and cpDNA genomes have been successfully applied to identifying different strains and tracking evolutionary origins of golden tide seaweed biomass in the Yellow Sea [20].
In this study, we analyzed sequence polymorphisms in the whole organelle genomes between two Korean indigenous samples, including one from Jindo Island in the South Sea and the other from Sokcho in the East Sea and nonindigenous Chinese sample (which was collected from Nanji Islands, Zhejiang Province) [21,22]. The aims of the present study were two-fold. First, we constructed and mapped the whole mt and cp genomes for the two Korean indigenous S. horneri individuals. Second, by comparing these organelle genomes with the previously published Chinese ones, we sought to develop novel genetic markers, such as indels (insertion or deletion mutations) and/or single nucleotide polymorphisms (SNPs) that can distinguish Korean indigenous from nonindigenous populations. This study will inform control and management efforts of golden tides and help to develop a “molecular-based early detection method” for nonindigenous S. horneri individuals in Korean waters, which is particularly essential for the rapid removal.

2. Materials and Methods

2.1. Sample Collection and DNA Extraction

Sargassum horneri individuals were sampled by scuba-diving from the rocky bed in approximately 3 m water depth at Jindo Island (JD; 34°22′35″N, 126°18′27″E) in the South Sea in January 2018 and at Sokcho (SC; 38°12′49.47″N, 128°36′8.48″E) in the East Sea in October 2018 (Figure 2). The sampled specimens were transported to the laboratory in ice box with frozen ice packs (approximately 4–10 °C) that blocked out light until DNA extraction. Using fresh leaf tissue, genomic DNA was extracted using a DNeasy Plant Mini Kit (Qiagen, Germany) according to the manufacturer’s directions after ground to fine power in liquid nitrogen. The DNA quality and quantity were checked by electrophoresis on 1 % agarose gel and also by a Qubit® 2.0 Fluorometer (Invitrogen, USA).

2.2. Genome Sequencing, Assembly, and Mapping

Paired-end sequencing for the total organelle genomes (mt and cpDNA genomes) of two Korean S. horneri samples was performed on the MiSeq (Illumina Inc., USA) platform with a library insert size of 300-500 bp at LabGenomics (http://www.labgenomics.co.kr/) in Korea. We obtained 168,510 and 3,407,076 raw read pairs with a length of 301 bp for JD and SC samples, respectively. These raw paired-end read pairs were trimmed using Geneious Prime (Biomatters Ltd., New Zealand) by error probability limit set to 0.05. The trimmed reads were aligned using Geneious in reference with previously determined S. horneri mt (GenBank accession no. KJ938300.1) and cp (accession no. KP881334.1) genome sequences sampled from Nanji Islands, Zhejiang province in China (CH; 27°27′N, 121°04′E) [21,22]. The number of mapped reads to the reference genomes was found to be 4885 (JD)–12,832 (SC) and 11,949 (JD)–34,730 (SC) for mt and cp genomes, respectively. Only using the mapped reads, de novo assembly was then performed with the default setting (assembler: Geneious, sensitivity: medium sensitivity/fast) of the Geneious. The Geneious de novo assembly is based on an overlap assembler with a blast-like algorithm.
The protein coding genes (PCGs), transfer RNA genes (tRNAs), and ribosomal RNA genes (rRNAs) in the organelle genomes were predicted and annotated using Dual Organellar GenoMe Annotator (DOGMA) using default parameters and then manually edited through a comparison with the published genome database [23]. tRNAs were confirmed using tRNAscan-SE [24]. The physical maps of the circular mt and cp genomes were drawn using OGDRAW program [25]. Mt and cp genome sequences of the Korean S. horneri samples can be accessed via GenBank with accession numbers MT795186–MT795189, respectively.

2.3. Comparative Analyses of Nucleotide Diversity and Synonymous (Ks) and Non-Synonymous Substitution (Ka) Rates

To explore genomic regions that can distinguish nonindigenous floating and indigenous Korean benthic populations, three mt and three cp genomes (JD, SC and CH) were aligned using MAFFT v1.4.0 alignment tool by set to default parameters (algorithm: auto, scoring matrix: 200PAM/K = 2, gap open penalty = 1.53, offset value = 0.123) [26]. We then manually searched for SNP or indel sites between two Korean and one Chinese genomes for each of the mt and cpDNA.
The mt and cpDNA genomes of three samples (JD, SC and CH) were aligned in a pairwise manner using a Geneious Alignment tool. Sliding window analysis of the mt and cp genomes was then conducted to calculate nucleotide diversity (Pi) values using DnaSP 5 [27] with a window length of 600 bp and a step size of 200 bp. To determine genomic regions that might have undergone a positive selection, we estimated synonymous (Ks) and non-synonymous (Ka) substitution rates for a total of 35 and 137 PCGs of mt and cpDNA, respectively using Geneious Alignment. The aligned total length, Pi, and Ka/Ks ratio for each separate gene were calculated in DnaSP 5.

3. Results

3.1. Mitochondrial and Chloroplast Genome Features of Korean S. horneri

The mt genomes of Korean S. horneri (JD, SC) were successfully assembled to show a size of 34,620 bp and 34,628 bp respectively (Figure 3). It turned out that Korean indigenous S. horneri have a 14–22 bp longer in a size of the whole mt genome than Chinese one (34,606 bp) [22]. The overall nucleotide compositions of JD and SC mt genomes were: A—27.3% (SC: 27.2%), T—36.6% (36.3%), C—15.1% (15.3%), and G—21.0% (21.2%), with an overall GC content of 36.1% (36.5%), respectively (Table 1). The Korean S. horneri mt genomes contained three rRNAs, 25 tRNAs and 35 PCGs including two open reading frames (ORFs).
The cp genomes of Korean S. horneri exhibited a size of 123,982 bp (JD) and 124,053 bp (SC), respectively. JD’s cp genome was mapped as a canonical quadripartite structure with two 5436 bp inverted repeat regions (IRs), which divided the circular molecule into a small single-copy region (SSC: 39,790 bp) and a large single-copy region (LSC: 73,320 bp). The overall nucleotide composition was: A—34.5%, T—34.9%, C—14.9%, and G—15.7%, with a total GC content of 30.6% (Table 1). Similarly, SC’s cp genome was mapped as a canonical quadripartite structure with two 5,435 bp IRs, which separated the circular molecule into a SSC (39,887 bp) and a LSC (73,296 bp). The overall nucleotide composition was: A—34.6%, T—34.8%, C—14.9%, and G—15.7%, with a total GC content of 30.6% (Table 1). The Korean S. horneri cp genomes had six rRNAs, 28 tRNAs and 137 PCGs including two ORFs.

3.2. Comparative Analysis of Genome Sequences between Korean and Chinese Individuals

By comparing the sequenced mt and cp genomes of Korean samples with the published genome database for a Chinese sample, the number and order of genes were observed to be the same as the reference genomes, suggesting no structural variation such as insertions or rearrangements exists. We identified 14 bp indel (insertion/deletion) mutation in nad7 (NADH dehydrogenase subunit 7)- trnP (proline tRNA) intergenic spacer (IGS) (Table 2). Genetic variation between JD and CH genomes revealed the highest number of seven SNPs on cytochrome c oxidase subunit II (cox2) in a total length of 3,147 bp. For a comparison of SC with CH genomes, cox2 gene also showed the highest number of 140 SNPs in a total of 3,154 bp. Although the mt genome had a shorter length in size than the cp genome, it had considerably more variable sites between Korean and Chinese genomes.
For the cp genome, we found 90 bp deletion in trnW (tryptophan tRNA) gene only for JD. While the mt genomes between JD and CH had a total of 52 SNPs in 34,620 bp, the cp genomes had only a total of 54 SNPs in 123,982 bp (Table 2). Overall, the mt genome (0.15%) had an approximately four-fold greater polymorphism than the cp genome (0.04%).

3.3. Nucleotide Diversity and Ka/Ks Ratio

We found 38 Korean unique SNPs in the mt genome and 45 SNPs in the cp genome, based on our comparative analysis of JD, SC, and CH genomes (Table 2).
An average Pi value among the three mt genomes was 0.0190, ranging from 0.0011 to 0.0522 (Figure 4A). Rnl (large subunit ribosomal RNA)-trnk (maturase K) IGS (JD vs SC: 0.04) and cox2 gene (JD vs. SC: 0.05) showed the highest Pi values (Figure 4). Although the cp genomes showed relatively lower levels of Pi values, psbJ (photosystem II reaction center protein J)-psaI (photosystem I reaction center subunit VIII) IGS and ilvH (acetolactate synthase small subunit)-ycf34 (conserved hypothetical plastid protein) genes showed the highest Pi values (Figure 5).
The Ka/Ks ratios in the 20 PCGs of mt genomes were always lower than 1, while Ka/Ks ratios in some of the 94 PCGs of the cp genomes were > 1 (Supplementary Material Tables S1 and S2). The Ka/Ks ratios of cp genomes between JD and SC showed always < 1 (Table S2), whereas that between JD and CH showed > 1 in only three genes, including dnaB (replicative DNA helicase), orf467 (open reading frame), secA (protein translocase subunit), and that between SC and CH showed > 1 in two genes, including petF (ferredoxin) and psaI (photosystem I reaction center subunit VIII) (Table 3). These results suggest that the genes with Ka/Ks ratio > 1 might be under a positive selection, in response to different regional environments (e.g., Korea vs. China).

4. Discussion

We successfully constructed the whole mt and cp genomes of two Korean benthic individuals using a reference genome available. From our organelle genomes-based comparative analyses between the Korean and Chinese Sargassum horneri forms, we identified several candidate markers, such as an indel and SNPs that can potentially be applied to distinguish Korean indigenous from Chinese nonindigenous populations. Although the mtDNA cox3 gene is known to harbor relatively high polymorphism and has often been used as a useful genetic marker for population genetic or phylogeographic studies of Sargassum species e.g., [15,28,29], we could not identify any sequence polymorphism between nonindigenous floating and indigenous Korean benthic populations of S. horneri in our previous study [4]. In the present investigation, we therefore sought to explore sequence variation between Korean indigenous and Chinese (that were sampled from Nanji Islands, Zhejiang Province where floating populations in Korean waters are presumed to be originated from there) at the whole organelle genome levels [4]. Excessive levels of within-species variation in DNA sequences are sometimes large challenges to the delineation of ‘species’ [30], while we sought to find a variation for the ‘hotspot’ region in organelle genomes that can distinguish Korean indigenous and nonindigenous individuals.
Although the results of intraspecific comparisons of S. horneri in our study showed fairly conserved levels of polymorphism in both mtDNA and cpDNA sequences, the cpDNA genomes showed Pi values approximately four times lower than the mtDNA genomes, supporting the previous findings that cpDNA is more conserved than mtDNA in the genus Sargassum [20]. At the species level, a three-fold greater evolutionary rate was detected in mtDNA relative to cpDNA [20]. In a previous study comparing cpDNA genomes of brown algae, the Order Fucales including the genus Sargassum was the most conservative taxa, especially for photosynthetic-related genes [21]. Depending on the importance of gene functions in connection with a functional constraint, mutations in the cpDNA genome may readily be eliminated if they occur, but mutations in mtDNA genome may accumulate as they would likely be selectively neutral [31], otherwise its non-neutral evolution is also feasible. We planned to test those ‘candidate’ genetic markers (an indel and SNPs) on both the mt and cp genomes, with an attempt to genetically discriminate between Korean indigenous and Chinese nonindigenous samples. The variable sites that we found can be exploited to identify differences between benthic and floating populations (e.g., application of SNPs for only Korean or nonindigenous samples). Although only two Korean benthic samples were analyzed, we suggest that those markers can be applied for distinguishing between nonindigenous floating and indigenous benthic samples.
We found that the highest number of SNPs was present in mtDNA cox2 gene, which was shown to accompany the largest changes in gene lengths when compared to other Phaeophyceae [32]. The cytochrome c oxidase (CO) is an enzyme that activates an enzyme that transfers electrons directly to molecular oxygen. It produces large amounts of ATP through cellular respiration in mitochondria. Therefore, it is thought that this particular gene has the highest mutation rate, among which cox2, a second subunit of CO genes can be strongly influenced by chemically unstable environments by cellular respiration. Also, we found an indel mutation in nad7-trnP IGS on mtDNA. An IGS is noncoding DNA sequences spanned between coding regions and is therefore known to have faster mutation rates than any other sites. Therefore, indel sites of IGS have frequently been used to identify invasive or alien species (populations) or lineages [33].
Comparison of SNPs on organelle genomes across PCGs can often provide important information for determining whether particular gene(s) may be a target by action of natural selection [34]. Our observed Ka/Ks ratio < 1 in all the genes of mt genomes suggests that purifying selection acts on (slightly) deleterious mutations. However, cp genomes showed a Ka/Ks ratio > 1 only in dnaB, orf467 and secA gene between JD and CH, and in the petF and psaI genes between SC and CH genomes, although Korean comparisons always showed a signal of negative selection. These results suggest the possibility that beneficial mutations might have evolved in these particular genes with different functions as a result of positive natural selection due to different habitat environments between Korea and China. DnaB’s function is known to unwind double-stranded DNA from chromosomal replication fork. Detected ORFs are regions of unknown functions that are potentially translated [35]; further analysis is needed to clarify their functions but the elevated Ka/Ks in the JD-CH comparison of orf467 may suggest that their functions are important to fitness. SecA gene delivers proteins through the cell membrane through ATP hydrolysis. PetF gene transports electrons in a variety of metabolisms. PsaI gene is involved in photosynthesis. The occurrence of natural selection, such as positive selection identified in five genes, can be under rapid evolutionary radiation, and is thus likely to influence speciation [36]. Among them, the psaI gene revealed the highest Ka/Ks ratio with a value close to six, and it is known that genes related to photosystem I are mostly conservative [37]. However, beneficial non-synonymous substitutions remain in the psaI gene of the SC and CH cp genomes, which might have been beneficial for survival and reproduction in different regional environments, but further studies are needed to verify this interesting hypothesis. However, there is a possibility that the observed signature of positive selection might be false inferences resulting from simultaneous non-independent substitutions within short stretches of DNA, called multinucleotide mutations (MNMs), which are known to prevail in eukaryote genomes [38,39]. Therefore, it is necessary to analyze more plastid genomes to test for the presence of MNMs to ensure whether our findings of positive selection are valid.
We have identified various SNP regions by conducting comparative analysis of genomes, however, nonetheless given that mtDNA and cpDNA of seaweeds are sometimes relatively more conservative, we should consider looking for variation in the nuclear genome [40] using more advanced genome sequencing technologies, such as double digest restriction-site associated DNA-sequencing (ddRAD-Seq) [41]. A recent study suggested that developed SNP markers in the nuclear genome of Sargassum muticum with ddRAD-Seq technique provide much higher marker resolution than any other markers, such as mtDNA and microsatellites [42]. It could also be possible that physiological or ecological characteristics (e.g., stable isotopic signature, morphology), rather than DNA sequences per se can allow us to differentiate Korean indigenous from nonindigenous populations. For example, the notion that the longer the drifting period, the more changes accompany in the thickness of thallus and stem, or even in color, needs to be studied. This would yield insights into how long plants have been floating around the sea, and also how the developmental morphology differs between floating and benthic populations [32].
Although only two Korean samples of S. horneri were analyzed and there will still be further efforts for the development of novel genetic markers, the obtained genome information will serve as reference database and thus will help to more precisely identify the origin of the floating S. horneri biomass that flows into the Korean coast. In the future, our study will inform management efforts including the development of “a DNA-based early detection technique for nonindigenous S. horneri”, which helps to mitigate ecological and economic damages on the Korean coastal ecosystem.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/2071-1050/12/18/7280/s1, Table S1: The Ka/Ks ratio calculated by pairwise protein-coding gene of each three mitochondrial genomes; Table S2: The Ka/Ks ratio calculated by pairwise protein-coding gene of each three chloroplast genomes. Sample abbreviations as in Figure 1.

Author Contributions

Conceptualization, S.K., S.-H.Y., H.-J.O., S.R.P, and H.J.L.; data curation, K.-S.C.; formal analysis, S.Y.B., K.-S.C., and H.J.L.; funding acquisition, S.R.P., T.-H.K., J.K.K., and H.J.L.; investigation, S.Y.B., and S.K.; methodology, S.Y.B., and K.-S.C.; project administration, S.-H.Y., H.-J.O., and S.R.P.; resources, S.K.; supervision, H.J.L.; validation, H.J.L.; visualization, T.-H.K., and J.K.K.; writing—original draft, S.Y.B., and H.J.L.; writing—review and editing, K.-S.C., S.K., S.-H.Y., H.-J.O., S.R.P., T.-H.K., and J.K.K. All authors read and approved the manuscript.

Funding

This research was a part of the project entitled “Improvement of management strategies on marine disturbing and harmful organisms (No. 20190518)” funded by the Ministry of Oceans and Fisheries, Korea. This study was also supported by a grant from the National Institute of Fisheries Science (R2020062), Korea and by a Sangji University Research Fund, 2018.

Acknowledgments

We thank members of the Molecular Ecology and Evolution Laboratory of Sangji University and Estuarine and Coastal Ecology Laboratory of Jeju National University for helping to collect samples in the field.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Smetacek, V.; Zingone, A. Green and golden seaweed tides on the rise. Nature 2013, 504, 84–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Wang, M.; Hu, C.; Barnes, B.B.; Mitchum, G.; Lapointe, B.; Montoya, J.P. The great Atlantic Sargassum belt. Science 2019, 365, 83–87. [Google Scholar] [CrossRef] [PubMed]
  3. Zhang, J.; Ding, X.; Zhuang, M.; Wang, S.; Chen, L.; Shen, H.; He, P. An increase in new Sargassum (Phaeophyceae) blooms along the coast of the East China Sea and Yellow Sea. Phycologia 2019, 58, 374–381. [Google Scholar] [CrossRef]
  4. Byeon, S.Y.; Oh, H.-J.; Kim, S.; Yun, S.H.; Kang, J.H.; Park, S.R.; Lee, H.J. The origin and population genetic structure of the ‘golden tide’ seaweeds, Sargassum horneri, in Korean waters. Sci. Rep. 2019, 9, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Yoshida, T. Studies on the distribution and drift of the floating seaweeds. Bull. Tohoku Reg. Fish. Res. Lab. 1963, 23, 141–186. [Google Scholar]
  6. Komatsu, T.; Tatsukawa, K.; Filippi, J.B.; Sagawa, T.; Matsunaga, D.; Mikami, A.; Ishida, K.; Ajisaka, T.; Tanaka, K.; Aoki, M. Distribution of drifting seaweeds in eastern East China Sea. J. Mar. Syst. 2007, 67, 245–252. [Google Scholar] [CrossRef]
  7. Komatsu, T.; Mizuno, S.; Natheer, A.; Kantachumpoo, A.; Tanaka, K.; Morimoto, A.; Hsiao, S.-T.; Rothäusler, E.A.; Shishidou, H.; Aoki, M. Unusual distribution of floating seaweeds in the East China Sea in the early spring of 2012. J. Appl. Phycol. 2014, 26, 1169–1179. [Google Scholar] [CrossRef] [Green Version]
  8. Hwang, E.K.; Lee, S.J.; Ha, D.S.; Park, C.S. Sargassum golden tides in the Shinan-gun and Jeju Island, Korea. Korean J. Fish. Aquat. Sci. 2016, 49, 689–693. [Google Scholar]
  9. Choi, S.K.; Oh, H.-J.; Yun, S.-H.; Lee, H.J.; Lee, K.; Han, Y.S.; Kim, S.; Park, S.R. Population dynamics of the ‘golden tides’ seaweed, Sargassum horneri, on the southwestern coast of Korea: The extent and formation of golden tides. Sustainability 2020, 12, 2903. [Google Scholar] [CrossRef] [Green Version]
  10. Guiry, M.; Guiry, G. Sargassum horneri (Turner) C.Agardh 1820. Available online: http://www.algaebase.org (accessed on 1 September 2020).
  11. Yoshida, G.; Arima, S.; Terawaki, T. Growth and maturation of the ‘autumn-fruiting type’of Sargassum horneri (Fucales, Phaeophyta) and comparisons with the ‘spring-fruiting type’. Phycol. Res. 1998, 46, 183–189. [Google Scholar] [CrossRef]
  12. Dawes, C.J.; Mathieson, A.C. The Seaweeds of Florida, 1st ed.; University Press of Florida: Gainesville, FL, USA, 2008. [Google Scholar]
  13. Komatsu, T.; Ariyama, H.; Nakahara, H.; Sakamoto, W. Spatial and temporal distributions of water temperature in a Sargassum forest. J. Oceanogr. Soc. Jpn. 1982, 38, 63–72. [Google Scholar] [CrossRef]
  14. Thiel, M.; Gutow, L. The ecology of rafting in the marine environment. I. The floating substrata. Oceanogr. Mar. Biol. 2005, 42, 181–264. [Google Scholar]
  15. Hu, Z.M.; Uwai, S.; Yu, S.H.; Komatsu, T.; Ajisaka, T.; Duan, D.L. Phylogeographic heterogeneity of the brown macroalga Sargassum horneri (Fucaceae) in the northwestern Pacific in relation to late Pleistocene glaciation and tectonic configurations. Mol. Ecol. 2011, 20, 3894–3909. [Google Scholar] [CrossRef] [PubMed]
  16. Liu, F.; Jin, Z.; Wang, Y.; Bi, Y.; Melton, J.T. Plastid genome of Dictyopteris divaricata (Dictyotales, Phaeophyceae): Understanding the evolution of plastid genomes in brown algae. Mar. Biotechnol. 2017, 19, 627–637. [Google Scholar] [CrossRef] [PubMed]
  17. Liu, F.; Li, X.; Che, Z. Mitochondrial genome sequences uncover evolutionary relationships of two Sargassum subgenera, Bactrophycus and Sargassum. J. Appl. Phycol. 2017, 29, 3261–3270. [Google Scholar] [CrossRef]
  18. Ševčíková, T.; Klimeš, V.; Zbránková, V.; Strnad, H.; Hroudová, M.; Vlček, Č.; Eliáš, M. A comparative analysis of mitochondrial genomes in eustigmatophyte algae. Genome Biol. Evol. 2016, 8, 705–722. [Google Scholar] [CrossRef] [Green Version]
  19. Keeling, P.J. The endosymbiotic origin, diversification and fate of plastids. Philos. Trans. R. Soc. B 2010, 365, 729–748. [Google Scholar] [CrossRef] [Green Version]
  20. Liu, F.; Pan, J.; Zhang, Z.; Moejes, F.W. Organelle genomes of Sargassum confusum (Fucales, Phaeophyceae): mtDNA vs cpDNA. J. Appl. Phycol. 2018, 30, 2715–2722. [Google Scholar] [CrossRef]
  21. Liu, F.; Pang, S. Chloroplast genome of Sargassum horneri (Sargassaceae, Phaeophyceae): Comparative chloroplast genomics of brown algae. J. Appl. Phycol. 2016, 28, 1419–1426. [Google Scholar] [CrossRef]
  22. Liu, F.; Pang, S.; Li, X.; Li, J. Complete mitochondrial genome of the brown alga Sargassum horneri (Sargassaceae, Phaeophyceae): Genome organization and phylogenetic analyses. J. Appl. Phycol. 2015, 27, 469–478. [Google Scholar] [CrossRef]
  23. Wyman, S.K.; Jansen, R.K.; Boore, J.L. Automatic annotation of organellar genomes with DOGMA. Bioinformatics 2004, 20, 3252–3255. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Chan, P.P.; Lowe, T.M. tRNAscan-SE: Searching for tRNA genes in genomic sequences. In Gene Prediction; Humana Press, Springer: New York, NY, USA, 2019; Volume 1962, pp. 1–14. [Google Scholar]
  25. Lohse, M.; Drechsel, O.; Bock, R. OrganellarGenomeDRAW (OGDRAW): A tool for the easy generation of high-quality custom graphical maps of plastid and mitochondrial genomes. Curr. Genet. 2007, 52, 267–274. [Google Scholar] [CrossRef] [PubMed]
  26. Katoh, K.; Standley, D.M. MAFFT multiple sequence alignment software version 7: Improvements in performance and usability. Mol. Biol. Evol. 2013, 30, 772–780. [Google Scholar] [CrossRef] [Green Version]
  27. Librado, P.; Rozas, J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics 2009, 25, 1451–1452. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Kogame, K.; Uwai, S.; Shimada, S.; Masuda, M. A study of sexual and asexual populations of Scytosiphon lomentaria (Scytosiphonaceae, Phaeophyceae) in Hokkaido, northern Japan, using molecular markers. Eur. J. Phycol. 2005, 40, 313–322. [Google Scholar] [CrossRef] [Green Version]
  29. Uwai, S.; Kogame, K.; Yoshida, G.; Kawai, H.; Ajisaka, T. Geographical genetic structure and phylogeography of the Sargassum horneri/filicinum complex in Japan, based on the mitochondrial cox3 haplotype. Mar. Biol. 2009, 156, 901–911. [Google Scholar] [CrossRef]
  30. Doebeli, M.; Dieckmann, U. Speciation along environmental gradients. Nature 2003, 421, 259–264. [Google Scholar] [CrossRef] [Green Version]
  31. Ballard, J.W.O.; Kreitman, M. Is mitochondrial DNA a strictly neutral marker? Trens. Ecol. Evol. 1995, 10, 485–488. [Google Scholar] [CrossRef]
  32. Liu, F.; Liu, X.; Wang, Y.; Jin, Z.; Moejes, F.W.; Sun, S. Insights on the Sargassum horneri golden tides in the Yellow Sea inferred from morphological and molecular data. Limnol. Oceanogr. 2018, 63, 1762–1773. [Google Scholar] [CrossRef]
  33. Voisin, M.; Engel, C.R.; Viard, F. Differential shuffling of native genetic diversity across introduced regions in a brown alga: Aquaculture vs. maritime traffic effects. Proc. Natl. Acad. Sci. USA 2005, 102, 5432–5437. [Google Scholar] [CrossRef] [Green Version]
  34. Liu, J.; Zhang, Y.; Lei, X.; Zhang, Z. Natural selection of protein structural and functional properties: A single nucleotide polymorphism perspective. Genome. Biol. 2008, 9, R69. [Google Scholar] [CrossRef] [Green Version]
  35. Sieber, P.; Platzer, M.; Schuster, S. The definition of open reading frame revisited. Trends Genet. 2018, 34, 167–170. [Google Scholar] [CrossRef]
  36. Nevado, B.; Wong, E.L.; Osborne, O.G.; Filatov, D.A. Adaptive evolution is common in rapid evolutionary radiations. Curr. Biol. 2019, 29, 3081–3086. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Scheller, H.V.; Jensen, P.E.; Haldrup, A.; Lunde, C.; Knoetzel, J. Role of subunits in eukaryotic Photosystem I. Biochim. Biophys. Acta 2001, 1507, 41–60. [Google Scholar] [CrossRef] [Green Version]
  38. Venkat, A.; Hahn, M.W.; Thornton, J.W. Multinucleotide mutations cause false inferences of lineage-specific positive selection. Nat. Ecol. Evol. 2018, 2, 1280–1288. [Google Scholar] [CrossRef]
  39. Schrider, D.R.; Hourmozdi, J.N.; Hahn, M.W. Pervasive multinucleotide mutational events in eukaryotes. Curr. Biol. 2011, 21, 1051–1054. [Google Scholar] [CrossRef] [Green Version]
  40. Lynch, M.; Koskella, B.; Schaack, S. Mutation pressure and the evolution of organelle genomic architecture. Science 2006, 311, 1727–1730. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Baird, N.A.; Etter, P.D.; Atwood, T.S.; Currey, M.C.; Shiver, A.L.; Lewis, Z.A.; Selker, E.U.; Cresko, W.A.; Johnson, E.A. Rapid SNP discovery and genetic mapping using sequenced RAD markers. PLoS ONE 2008, 3, e3376. [Google Scholar] [CrossRef] [PubMed]
  42. Le Cam, S.; Daguin-Thiébaut, C.; Bouchemousse, S.; Engelen, A.H.; Mieszkowska, N.; Viard, F. A genome-wide investigation of the worldwide invader Sargassum muticum shows high success albeit (almost) no genetic diversity. Evol. Appl. 2019, 13, 500–514. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Golden tide seaweed Sargassum horneri biomass introduced to Jeju Island in South Korea during 2019. (A) S. horneri piled up along Sinyang Seopji Beach; (B) floating biomass in seawater from Udo Island; (C) floating biomass from Seogwipo; (D) floating S. horneri entangled on the rope of a ship from Jagunae Port.
Figure 1. Golden tide seaweed Sargassum horneri biomass introduced to Jeju Island in South Korea during 2019. (A) S. horneri piled up along Sinyang Seopji Beach; (B) floating biomass in seawater from Udo Island; (C) floating biomass from Seogwipo; (D) floating S. horneri entangled on the rope of a ship from Jagunae Port.
Sustainability 12 07280 g001
Figure 2. Map of sampling sites for comparative genomic analysis of S. horneri. Two Korean benthic samples (filled circles in dark blue) were collected at Jindo Island (JD) in the South Sea and at Sokcho (SC) in the East Sea of Korea, and a sample of reference genome (filled circles in dark red) [21,22] was collected at Nanji Island (CH), Zhejiang province in China. Areas in red denote Zhoushan region, Zhejiang province in China, which was assumed to be geographic/genetic origin(s) of floating populations in our previous study [4]. Areas in blue denote Jeju Island and Haenam in the South Sea, respectively, which were main coastal regions in Korea influenced by floating S. horneri populations.
Figure 2. Map of sampling sites for comparative genomic analysis of S. horneri. Two Korean benthic samples (filled circles in dark blue) were collected at Jindo Island (JD) in the South Sea and at Sokcho (SC) in the East Sea of Korea, and a sample of reference genome (filled circles in dark red) [21,22] was collected at Nanji Island (CH), Zhejiang province in China. Areas in red denote Zhoushan region, Zhejiang province in China, which was assumed to be geographic/genetic origin(s) of floating populations in our previous study [4]. Areas in blue denote Jeju Island and Haenam in the South Sea, respectively, which were main coastal regions in Korea influenced by floating S. horneri populations.
Sustainability 12 07280 g002
Figure 3. Map of the mitochondrial and chloroplast genomes of two Korean benthic individuals of Sargassum horneri. The genes drawn outside the circles are transcribed counterclockwise, whereas genes drawn inside the circles are transcribed clockwise. The dark-gray inner circles indicate GC content, and the light-gray inner circles represent AT content.
Figure 3. Map of the mitochondrial and chloroplast genomes of two Korean benthic individuals of Sargassum horneri. The genes drawn outside the circles are transcribed counterclockwise, whereas genes drawn inside the circles are transcribed clockwise. The dark-gray inner circles indicate GC content, and the light-gray inner circles represent AT content.
Sustainability 12 07280 g003
Figure 4. Sliding window analysis of nucleotide diversity (Pi) for S. horneri mitochondrial genomes. A window size of 600 bp and a step size of 200 bp were applied. (A) Levels of Pi among the three genomes and (BD) pairwise Pi values between the two genomes. Individual abbreviations as in Figure 1.
Figure 4. Sliding window analysis of nucleotide diversity (Pi) for S. horneri mitochondrial genomes. A window size of 600 bp and a step size of 200 bp were applied. (A) Levels of Pi among the three genomes and (BD) pairwise Pi values between the two genomes. Individual abbreviations as in Figure 1.
Sustainability 12 07280 g004
Figure 5. Sliding window analysis of nucleotide diversity (Pi) for S. horneri chloroplast genomes. A window size of 600 bp and a step size of 200 bp were applied. (A) Levels of Pi among the three genomes and (BD) pairwise Pi values between the two genomes. Individual abbreviations as in Figure 1.
Figure 5. Sliding window analysis of nucleotide diversity (Pi) for S. horneri chloroplast genomes. A window size of 600 bp and a step size of 200 bp were applied. (A) Levels of Pi among the three genomes and (BD) pairwise Pi values between the two genomes. Individual abbreviations as in Figure 1.
Sustainability 12 07280 g005
Table 1. Molecular features of mitochondrial (mt) and chloroplast (cp) genomes of two Korean benthic Sargassum horneri. Sample abbreviations as in Figure 1.
Table 1. Molecular features of mitochondrial (mt) and chloroplast (cp) genomes of two Korean benthic Sargassum horneri. Sample abbreviations as in Figure 1.
RegionJDSC
Molecular FeaturesMt GenomeCp GenomeMt GenomeCp Genome
GC%36.130.636.530.6
rRNA3636
tRNA25282528
Protein coding genes3513735137
Table 2. Comparison of sequence polymorphisms between two Korean benthic samples and Chinese reference genomes to identify candidate marker sites that can be used for distinguishing Korean indigenous and nonindigenous populations. Sample abbreviations as in Figure 1.
Table 2. Comparison of sequence polymorphisms between two Korean benthic samples and Chinese reference genomes to identify candidate marker sites that can be used for distinguishing Korean indigenous and nonindigenous populations. Sample abbreviations as in Figure 1.
Mitochondrial GenomeChloroplast Genome
PositionLocusCHJDSCPositionLocusCHJDSC
1192rnlTGG1880rnlATT
1295rnlTGG1881rnlTAA
3842rps8TCC1911rnlTCC
4024rps8CAA4201rnsGAA
4524rpl6AGG9208ycf47TCC
5846rps4GAA12,192psbVCTT
8384trnL-CC17,056dnaKAGG
10,071rpl14GAA20,737secAGAA
11,590rpl16, rps3GAA26,794chlBCTT
13,118rpl2CTT45,360acsFACC
16,567nad2ATT46,202rpl35GAA
17,023nad2ATT47,715ycf34-trnL IGST--
17,224nad2AGG67,316orf467ACC
20,832cobCTT67,665orf467TCC
22,306cox2AGG75,251petLTCC
23,256cox2AGG75,509ycf4GCC
24,936nad4CTT77,134psaIATT
25,965nad4-nad5 IGSCAA77,137psaIGAA
26,158nad5ATT77,142psaICAA
29,327nad11CTT77,148psaIATT
29,522nad3CTT77,151psaIATT
30,348atp8GAA77,153psaICAA
30,415trnSCTT77,156psaIGCC
31,478rpl31T--77,158psaITAA
34,505*nad7-trnP IGS
(indel region)
-AC79,770rnsCTT
34,506-AA82,090rnlATT
34,507-TT82,091rnlTAA
34,508-AA89,686ycf33-ycf39 IGST--
34,509-TT91,872sufBAGG
34,510-TT92,695sufBAGG
34,511-TT98,854psbD-trnN IGS
(indel region)
A--
34,512-AA98,855T--
34,513-GG98,856A--
34,514-CC98,857T--
34,515-TT98,858A--
34,516-GG98,859T--
34,517-AA10,0113dnaBGAA
34,518-AA10,0384dnaBAGG
10,3700petFCAA
10,3801petFTCC
10,3851petFGAA
11,0124tufATGG
11,4149rps13AGG
11,7303rps8AGG
12,3529rpl3TAA
Table 3. Results of synonymous (Ks) and non-synonymous (Ka) substitution rates for the protein coding genes (PCGs) of the chloroplast genomes between two Korean benthic and Chinese reference genomes. Although Ka/Ks ratios were always < 1 between Korean samples, two (petF, psaI) and three (dnaB, orf467, secA) genes showed Ka/Ks ratios >1 between Korean and Chinese comparisons. Sample abbreviations as in Figure 1.
Table 3. Results of synonymous (Ks) and non-synonymous (Ka) substitution rates for the protein coding genes (PCGs) of the chloroplast genomes between two Korean benthic and Chinese reference genomes. Although Ka/Ks ratios were always < 1 between Korean samples, two (petF, psaI) and three (dnaB, orf467, secA) genes showed Ka/Ks ratios >1 between Korean and Chinese comparisons. Sample abbreviations as in Figure 1.
Genomic LocusJD-CH Ka/KsSC-CH Ka/KsJD-SC Ka/Ks
dnaB1.73330.18130.2180
orf4673.88890.17400.1662
petF0.00001.19050.2925
psaI0.68085.94010.0000
secA3.60000.03720.0174

Share and Cite

MDPI and ACS Style

Byeon, S.Y.; Cheon, K.-S.; Kim, S.; Yun, S.-H.; Oh, H.-J.; Park, S.R.; Kim, T.-H.; Kim, J.K.; Lee, H.J. Comparative Analysis of Sequence Polymorphism in Complete Organelle Genomes of the ‘Golden Tide’ Seaweed Sargassum horneri between Korean and Chinese Forms. Sustainability 2020, 12, 7280. https://0-doi-org.brum.beds.ac.uk/10.3390/su12187280

AMA Style

Byeon SY, Cheon K-S, Kim S, Yun S-H, Oh H-J, Park SR, Kim T-H, Kim JK, Lee HJ. Comparative Analysis of Sequence Polymorphism in Complete Organelle Genomes of the ‘Golden Tide’ Seaweed Sargassum horneri between Korean and Chinese Forms. Sustainability. 2020; 12(18):7280. https://0-doi-org.brum.beds.ac.uk/10.3390/su12187280

Chicago/Turabian Style

Byeon, Seo Yeon, Kyeong-Sik Cheon, Sangil Kim, Suk-Hyun Yun, Hyun-Ju Oh, Sang Rul Park, Tae-Hoon Kim, Jang Kyun Kim, and Hyuk Je Lee. 2020. "Comparative Analysis of Sequence Polymorphism in Complete Organelle Genomes of the ‘Golden Tide’ Seaweed Sargassum horneri between Korean and Chinese Forms" Sustainability 12, no. 18: 7280. https://0-doi-org.brum.beds.ac.uk/10.3390/su12187280

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop