Open Access

DNA variation in the phenotypically-diverse brown alga Saccharina japonica

  • Evgeniy S Balakirev1, 2Email author,
  • Tatiana N Krupnova3 and
  • Francisco J Ayala1
BMC Plant Biology201212:108

DOI: 10.1186/1471-2229-12-108

Received: 25 October 2011

Accepted: 21 June 2012

Published: 11 July 2012



Saccharina japonica (Areschoug) Lane, Mayes, Druehl et Saunders is an economically important and highly morphologically variable brown alga inhabiting the northwest Pacific marine waters. On the basis of nuclear (ITS), plastid (rbcLS) and mitochondrial (COI) DNA sequence data, we have analyzed the genetic composition of typical Saccharina japonica (TYP) and its two common morphological varieties, known as the “longipes” (LON) and “shallow-water” (SHA) forms seeking to clarify their taxonomical status and to evaluate the possibility of cryptic species within S. japonica.


The data show that the TYP and LON forms are very similar genetically in spite of drastic differences in morphology, life history traits, and ecological preferences. Both, however, are genetically quite different from the SHA form. The two Saccharina lineages are distinguished by 109 fixed single nucleotide differences as well as by seven fixed length polymorphisms (based on a 4,286 bp concatenated dataset that includes three gene regions). The GenBank database reveals a close affinity of the TYP and LON forms to S. japonica and the SHA form to S. cichorioides. The three gene markers used in the present work have different sensitivity for the algal species identification. COI gene was the most discriminant gene marker. However, we have detected instances of interspecific COI recombination reflecting putative historical hybridization events between distantly related algal lineages. The recombinant sequences show highly contrasted level of divergence in the 5’- and 3’- regions of the gene, leading to significantly different tree topologies depending on the gene segment (5’- or 3’-) used for tree reconstruction. Consequently, the 5’-COI “barcoding” region (~ 650 bp) can be misleading for identification purposes, at least in the case of algal species that might have experienced historical hybridization events.


Taking into account the potential roles of phenotypic plasticity in evolution, we conclude that the TYP and LON forms represent examples of algae phenotypic diversification that enables successful adaptation to contrasting shallow- and deep-water marine environments, while the SHA form is very similar to S. cichorioides and should be considered a different species. Practical applications for algal management and conservation are briefly considered.


Brown algae Kelp Saccharina DNA polymorphism Interspecific COI recombination Phenotypic plasticity Saccharina japonica morphological forms


Genetically-based approaches (especially using multi-locus data) have resolved much phylogenetic uncertainty in large-scale algal trees of life (e.g., [13]). Less obvious progress, however, has been accomplished concerning the species level of algal taxonomy, when there is an abundance of morphological variants with uncertain taxonomic status. For instance, the genus Laminaria (that has been recently separated into two different genera, Laminaria and Saccharina; [1]) is rich with morphological types, varieties, ecotypes, and forms that have puzzled taxonomists since the very beginning of research on this subject [414]. Whether the distinct morphologies observed in many algae represent different species or only intraspecific forms is one of the main unresolved questions.

In this paper we focus on the most common and commercially important laminarialean species, Saccharina japonica[4, 5, 15] and its morphological forms, inhabiting contrasting depths of the northwest Pacific region. The species is morphologically and ecologically highly variable and some intraspecific forms have been described [1013]. Typical S. japonica (TYP) inhabits the littoral zone at the preferred depths of 5 – 11 m with wide distribution over the full species area. Its full distribution pattern, however, remains uncertain owing to the difficulty of reliable species identification. It grows along the sea shore on rocky substrate [5, 6, 1517]. (See Additional file 1 for a more detailed species description.) At present the populations of the TYP form are depressed due to overharvesting and possibly global climate changes [18, 19].

The deep-water (or “longipes”) S. japonica form (LON) inhabits the sublittoral zone at preferred depths of 14 – 25 m. This endemic form has a restricted distribution in the Sea of Japan and the Sea of Okhotsk. It grows in compact settlements at a significant distance (300 – 1000 m) from the seashore [6, 16, 17, 20, 21]. The LON form thallus is 2 – 3 times longer and 3 – 4 times heavier than the thallus of the TYP form. Differences in morphology and preferred ecology are detectable already among algae during the first year of growth. Originally the LON form was described as a deep-water variant of S. japonica[4, 6]; a new nomenclatural designation, S. japonica f. longipes, has been suggested recently [13]. The TYP and LON forms have significant differences in morphology, anatomy, reproductive biology, ecology, and other important features (Additional file 1). The opportunity for interbreeding between the LON and TYP forms is considered highly unlikely because the formation, release, and appearance of the mature gametes are separated in space and time. The zoospore spreading time of S. japonica (and the close species, S. angustata) is short (maximum 24 hours) leading to limited free dispersal gamete range (on average no more than 3.5 meters [19, 22]). Also the distribution areas of the LON and TYP forms are separated geographically (not less than three miles from each other) and their reproductive periods do not overlap (Additional file 1). The morphological differences between the TYP and LON forms of S. japonica are stable and easily recognizable and they are systematically reported from the northern Primorye coast region [20, 21]. The drastic differences between the TYP and LON forms have motivated changing the taxonomical status of the LON form, so that it would be considered a separate species [17]. However, transplant experiments [23, 24] show that zoospores of the LON form released in the habitat of TYP form produced algae morphologically indistinguishable from the TYP form, which did not support the biological species status of the LON form, but rather suggested high phenotypic plasticity for the species.

The shallow-water form of Saccharina japonica (SHA) inhabits the supralittoral zone (0.1 – 0.5 m depth) and it is widely distributed in the Primorye coast region, the Sea of Japan, including Peter the Great Bay. It is commonly accepted that the SHA form represents an intraspecific ecotype of the TYP S. japonica form. The zoospores of the SHA form are used to enhance the strongly overharvested resources of the TYP form [19]. Nevertheless, the SHA and TYP forms have clearly distinct morphological features (Additional file 1). There are also significant differences between the forms in heat and salinity tolerance [25]. The SHA form is more heat and desalination resistant and survives better than the TYP form under hyposalinity osmotic stress and increasing temperature conditions. It was suggested [26] that the SHA form may represent a subspecies of Laminaria (Saccharina) angustata adapted to the supralittoral zone.

We analyze the genetic composition of S. japonica and its common morphological forms, LON and SHA, seeking to clarify the taxonomic status of the forms and to evaluate the possibility of cryptic species within S. japonica, using DNA nucleotide sequences of three genes from mitochondrial, plastid, and nuclear genomes (4,286 bp total in a concatenated dataset). We also test for homogeneous phylogenetic signal within and among the marker loci (ITS, rbcLS, and COI). We have found that the TYP and LON morphological forms are genetically very similar; consequently they need not be considered as different biological species. However, taking into account the potential evolutionary roles of phenotypic plasticity, and its relevance to adaptive evolution and speciation in many organisms (e.g., [2729]) including algae (e.g., [30, 31]) we propose that the TYP and LON forms may represent an important example of algal diversification that enables successful adaptation to contrasting shallow- and deep-water marine environments. The SHA form has clear genetic differences from the TYP and LON forms and it might represent a morphological variant of S. cichorioides expanded to the supralittoral environment. The practical implications for mariculture are that all three forms should be considered separate evolutionary lineages and that restoration programs should include genetic data to avoid confounding.


Algae samples

The specimens of S. japonica morphological forms were collected from the Primorye coastal region, Sea of Japan. The TYP (five specimens) and SHA (four specimens) forms were collected near the Cape Dal'niy at depths of 6.0 m and 0.5 m, respectively. The LON form (three specimens) was collected near the Cape Zolotoi at depths of 15.0 m. We also analyzed first year individuals belonging to the S. japonica typical form (TYPF, four specimens), collected near the Cape Dal'niy; unidentified Saccharina species (CHE, four specimens; depth of 0.5 m), TYP form growing in close proximity with unidentified Saccharina species (TYPA, four specimens), and a sample of the TYP form without middle line (TYPW, a single specimen), all collected in the Bay Chernoruch'e near the Cape Khitrovo. The algae were identified by eye in the field.

DNA amplification and sequences

Total genomic DNA was extracted using the DNeasy Plant MiniKit protocol (Qiagen, Hilden, Germany). The procedures for DNA amplification and sequencing have been described previously [3234]. Plastid sequences for the partial rbcLS operon (RubisCO; 1,822 bp total) included 1,333 bp of the flanking rbcL gene, the complete spacer (282 bp), and 207 bp of the flanking rbcS gene. For the rbcLS operon, primers and PCR conditions followed [1]. Mitochondrial sequences of a 1,788-bp COI fragment were amplified using newly developed primers (annealing temperature 52°C): 5’-cttatcaaaaggtgcatctatgg-3’ (SacCOIF; forward) and 5’-acactctaccgctgagttacaag-3’ (SacCOIR; reverse) designed based on full mitochondrial genomes of 11 brown algal species [11, 3537]. Our sequences include the full coding region (1,602 bp, 534 codons) of the COI gene as well as 108 bp and 78 bp of the 5'- and 3'-flanking regions, respectively. The nuclear 676-bp fragment of the 18 S–ITS1–5.8 S–ITS2 rDNA region (ITS) was amplified using primers (annealing temperature 53°C): 5’-aggtccgaacgaaagtggta-3’ (SacITSF; forward) and 5’-acaaggtttccgtaggtgaac-3’ (SacITSR; reverse). The amplified fragment includes partial 18 S ribosomal RNA gene (17 bp), complete ITS1 (244 bp), complete 5.8 S ribosomal RNA gene (159 bp), and partial ITS2 (256 bp). Using all three gene regions we constructed a 4,286-bp concatenated alignment (with alignment gaps) of 25 algae samples. At least two independent PCR amplifications were sequenced for each polymorphic site in all 25 algae samples to correct for possible PCR or sequencing errors. (See Additional file 2 for the PCR details). The COI, ITS, and rbcLS sequences have been deposited in GenBank under accession numbers JN873222-JN873246, JN873247-JN873271, and JN873272-JN873296, respectively.

DNA sequence analysis

The sequences were assembled using the program SeqMan (Lasergene, DNASTAR, Inc.). Multiple alignment was carried out manually and using the program CLUSTAL W [38]. The computer programs DnaSP, version 5 [39] and PROSEQ, version 2.9 [40] were used for most intraspecific analyses. Saccharina coriacea, S.angustata, Laminaria digitata, and Agarum clathratum (see Additional file 3) were used as outgroup taxa; they were selected based on previous molecular evidence of close relationship to S. japonica[1, 2, 11, 4143] and screening of nucleotide sequences available in GenBank.

Model-based phylogeny reconstructions were performed with the rbcLS, COI, ITS, and concatenate sequence alignments using the neighbor-joining (NJ) and maximum-likelihood (ML) methods in MEGA, version 5 [44]. MEGA5 or jModelTest [45] were used to find the best model of substitution under the maximum likelihood criterion. Models with the lowest Bayesian information criterion [46] scores were considered most appropriate to describe evolution for each gene; the Akaike Information Criterion [47] was also applied. Maximum likelihood bootstrap analyses [48] consisted of 1000 replicates.

Partitioned analyses were performed with GARLI, version 2.0 [49] and MrBayes, version 3.2 [50] that allowed the overall rate to be different across partitions (rbcLS, COI, and ITS). Substitution model parameters were estimated separately for each gene. The most appropriate models were Tamura 3-parameter plus gamma for rbcLS, Hasegawa-Kishino-Yano plus I for COI, and Tamura 3-parameter for ITS.

The complete dataset was also analyzed by Bayesian inference using MrBayes, version 3.2 [50] under a general-time-reversible model plus gamma plus I. Proportion of invariable sites was uniformly distributed on the interval (0.00, 1.00). Gamma distribution was approximated using 4 categories. Analyses were performed as two independent runs, each with four incrementally heated Metropolis-coupled Monte-Carlo Markov Chains running for one million generations. Output trees and data were sampled every 500 generations. Likelihood values reached a plateau within approximately 1,000 generations. A total of 4002 trees in two files were read and 3002 of them were sampled. The numbers of unique site patterns were 100, 175, and 107 for rbcLS, COI, and ITS, respectively. It showed that COI as the most efficient gene marker consisting with the distance-based methods.

The average standard deviation of split frequencies at the end of the run was 0.0068 indicating stationary conditions. The log likelihood values increased from below −35292.434 to around −9652.102 in the first five thousand generations and then to around −9547.021 after one million generations. The likelihood of best state for “cold” chain of run 1 was −9523.80 and the likelihood of best state for “cold” chain of run 2 was −9524.14. At the end of run there were not any trends in a plot of generation versus the log probability of the observed data (the log likelihood values). A convergence diagnostic, the Potential Scale Reduction Factor (PSRF) [51] was between 1.00 and 1.003 indicating a good sample from the posterior probability distribution.

The topologies obtained with neighbor-joining and maximum-likelihood methods as well as with Bayesian inference were very similar. The close congruence could be explained by the fact that 1) the dataset was relatively straightforward and included just two main compared groups of sequences (TYP and SHA) and four outgroup species (S. coriacea, S. angustata, Laminaria digitata, and Agarum clathratum); 2) the extent of sequence divergence between the TYP and SHA forms for all gene regions was less than 5.0% (0.9% for rbcLS, 1.1% for ITS, and 4.8% for COI); and 3) the length of alignment was long enough (4.2 kb totally). As has been shown by many authors, the relative efficiencies of the ML and NJ methods in obtaining the correct tree topology were very close under these conditions [5254]. The only difference between the two methods in relation to our dataset concerned slightly different bootstrap values. To be conservative we show the lowest bootstrap values obtained with the ML method. When different methods produced similar or identical topologies the simplest, NJ method, was preferable.


Nucleotide diversity and divergence

We sequenced three gene regions (total 4,286-bp alignment) in 25 algae samples representing common morphological forms of S. japonica. All three gene regions (rbcLS, COI, and ITS) with comparable length were available in GenBank for some species related to S. japonica, which we have included in the analysis. The lengths of the GenBank concatenated sequences are (in parentheses): S. japonica (3,182 bp), S. diabolica (2,894 bp), S. longipedalis (2,896 bp), S. ochotensis (2,689 bp), S. religiosa (2,896 bp), S. angustata (4,256 bp), S. coriacea (2,711 bp), S. latissima (3,748 bp), and L. digitata (4,245 bp) (see Additional file 3 for GenBank accession numbers and full species names).

Figure 1 shows 114 nucleotide substitutions detected in three gene regions (18 sites in rbcLS, 88 sites in COI, and 8 sites in ITS) among the 25 algae samples. In addition, there were seven length polymorphisms in the rbcLS intergenic spacer and ITS regions (Additional file 4); no length polymorphisms were found in the COI gene. All algae samples sharply split up in two groups. The first group includes all samples of the TYP (including TYPF, TYPA, and TYPW) and LON forms; the second group includes the SHA form and one unknown species from the Chernoruch'e bay (CHE). The genetic structure is completely congruent for all three gene regions studied (Figure 1). The two groups (denoted as “TYP” and “SHA” lineages) differ by 109 fixed single nucleotide differences, as well as by seven fixed length polymorphisms. The lineages are differentially associated with indels. The TYP lineage is completely associated with three inferred deletions (13-, 4-, and 8-bp, within the rbcLS intergenic spacer and ITS region; 1, 2, and 3, Figure 1) and two inferred insertions (1 and 4) within the ITS region. The SHA lineage is completely associated with two inferred 1-bp insertions within the ITS region (2 and 3). The difference between the TYP and SHA lineages based on the full concatenate (including the three gene regions, rbcLS, COI, and ITS) is highly significant (Fst = 0.9949, P < 0.00001); the total sequence divergence (Dxy) between the lineages is 0.0261 (ignoring indels).
Figure 1

Nucleotide substitution sites in the rbc LS, COI , and ITS gene regions of the Saccharina japonica morphological forms. The numbers above the top sequence represent the position of segregating sites and the start of a deletion or insertion (see also Additional file 4). Nucleotides are numbered from the beginning of our sequence. Dots indicate the same nucleotide as the reference sequence. The hyphens represent deleted nucleotides. denotes a deletion; † denotes the absence of a deletion; denotes an insertion; ‡ denotes the absence of an insertion. The sample abbreviations are in the section "Methods". COR: S. coriacea. The underlined nucleotides in bold face show fixed differences between COR and the other sequences.

The three gene regions differed in levels of divergence between the TYP and SHA lineages (Figure 2). Total nucleotide divergence between TYP and SHA lineages is low for rbcLS (Dxy = 0.0099 ± 0.0025) and ITS (Dxy = 0.0106 ± 0.0037), but significantly higher (more than four times) for COI (Dxy = 0.0481 ± 0.0051). The same tendency is found for comparisons between all other laminarialean species pairs (data not shown). Thus the COI gene represents the most effective and sensitive gene marker of the three studied here for the comparison of close Saccharina species, which is in accordance with the data obtained for red algae and kelps (e.g., [5557]).
Figure 2

Sliding window plot of nucleotide divergence for three gene regions ( rbc LS, COI , and ITS ) between the TYP and SHA forms of Saccharina japonica. Window sizes are 100 nucleotides with one-nucleotide increments. A schematic representation of the regions is displayed at bottom. The thick black line marks the COI region investigated by McDevit and Saunders [43].

The estimates of nucleotide diversity for all sequences of S. japonica and for the morphological forms or haplotype lineages separately are low. In the pooled sample, the total nucleotide diversity is similar for the TYP (π = 0.00011) and SHA (π = 0.00016) lineages (based on the three gene regions studied); the results are in accordance with literature data on intraspecific variability in S. japonica[58].

Figure 3 displays a maximum likelihood phylogram of the three gene regions separately and of the concatenate sequences obtained for the S. japonica forms in the present study, along with other Saccharina and Laminaria sequences obtained from GenBank. S. ochotensis, S. religiosa, S. longipedalis, and S. diabolica are conspesific with S. japonica (data not shown), as has been revealed previously [10, 11]. All trees showed that the sequences from the TYP and LON forms form a single clade, with no evidence of discrete species heterogeneity (Fst = 0.0201; P = 0.6600; Kst = − 0.0099; P = 0.7480); total sequence divergence (Dxy) is 0.0001. These data suggest that the TYP and LON morphological forms of S. japonica are not distinct biological species. In contrast, the sequences from the TYP and SHA lineages are distinct. The difference between the lineages is highly significant (Fst = 0.9949; P <0.0001; Kst = 0.8484; P <0.0001). The distance between them is 2.61% (2.77% including indels), which is in a range of the divergence between species observed within the genus Saccharina[1, 4143] suggesting that the TYP and SHA lineages would qualify as different species based on a phylogenetic species concept.
Figure 3

Phylogenetic trees of the rbc LS, COI, ITS , and concatenate sequences (4,286 bp total) of Saccharina japonica morphological forms. The topology of trees obtained with Maximum likelihood and Neighbor-joining methods were congruent. Numbers at the nodes are bootstrap percent support values based on 1,000 replications in Maximum likelihood analysis. Saccharina latissima was excluded from this analysis due to recombination in the COI gene (see text). Other comments as in Figure 1.

Species identity

Comparison of sequences obtained during this study with those available in the GenBank database suggests that the TYP lineage represents S. japonica, whereas the SHA lineage is closely related to a group of Saccharina species that include S. coriacea, S. cichorioides, S. sachalinensis, and S. yendoana (the last three species are not represented in Figure 3 because the COI sequences are not available for them). The distance between the SHA lineage and this group of Saccharina species (represented by S. coriacea in Figures 1 and 3) is low and not significant (Fst = 0.9242; P = 0.3300; Kst = 0.3890; P = 0.3300); the rbcLS and ITS regions are identical but the COI gene shows seven single nucleotide differences (Figure 1). Total sequence divergence is 0.38% (for the full 4,286-bp concatenate alignment), which is below the species divergence level observed in algae [1, 41] and other eukaryotes [59]. Yotsukura and colleagues [8, 9, 11, 12] have revealed high genetic similarity between S. coriacea, S. cichorioides, S. sachalinensis, and S. yendoana and suggested their conspecificity, in accordance with earlier results of Petrov [6], who treated L. coriacea and L. sachalinensis as intraspecific forms of L. cichorioides (the status of L. yendoana was not discussed by Petrov [6].) Following the priority rule of the International Code of Botanical Nomenclature, Selivanova and colleagues [13] introduced a new nomenclature combination: S. cichorioides with four intraspecific forms. The conspecificity of S. coriacea, S. cichorioides, S. sachalinensis, and S. yendoana was suggested on the basis of low genetic divergence detected with ITS-1, ITS-2, RubisCO spacer, and rbcL sequences [8, 9, 1113, 60]. In the present work we have detected low divergence between the SHA lineage and the group of Saccharina species listed above using all three gene markers, including the COI gene (Dxy = 0.0050), which is considerably more informative than the rbcLS and ITS markers (see above). Consequently, we conclude that the SHA form belongs to the S. cichorioides group (all “species” representing the group were previously synonymized based on genetic data as pointed out) and it may represent a morphological ecotype adapted to supralittoral environments. S. cichorioides is more heat resistant than S. japonica and its main populations are more southward in the Sea of Japan. However, this morphologically and ecologically versatile species (as it follows from our data) is capable of occupying different environments and successfully compete with local species in the northern Primorye coast region under hyposalinity osmotic stress and increased water temperature (see Background).

The unknown Saccharina collected in the Bay Chernoruch'e (CHE) are genetically identical to the SHA form, even though they have substantial differences in external morphologies (T. N. Krupnova, unpublished observations) but occupy similar depths (0.5 m). It has been suggested that laminarialean algae from the supralittoral zone represent S. angustata[26]. However, this is not the case (see Figure 3): the SHA form is close to the S. cichorioides group (represented by S. coriacea in Figure 3) and clearly different from S. angustata. The distance between the SHA form and S. angustata is 2.63% (which is consistent with interspecific divergence level), indicating species taxonomical range for these two evolutionary lineages. Thus, S. angustata was not found among 25 algal samples collected from distant areas of the Northern Primorye region. This species is common in the Southern Kuril Islands and Japan islands [4, 26], but it rarely occurs (if at all) along the continental coast of the Japan Sea.

Non-uniform pattern of divergence in the COIgene

Full length of the COI coding region in laminarialean algae is around 1.6 kb (1,602 bp in Saccharina). The 5’ part of the gene (~ 650 bp in length; 5’-COI) has been chosen as a “barcoding” region [61] and recommended for algal species identification [4143, 61]. We have shown above that the most pronounced differences between Saccharina species are detected within the COI gene. However, the level of divergence is highly non-uniform within the COI gene in the comparisons, including S. latissima and other Saccharina species (Figure 4A – D). The average divergence for the 5’-COI is 0.0380 ± 0.0106 but 4.2 times more, 0.1579 ± 0.0034 for the 3’-COI (paired samples t-test P = 0.0001; Mann–Whitney test P = 0.0090; Table 1). This difference is also obvious for the comparisons with other laminarialean species for which long enough COI sequences are available in GenBank [2, 11]: Laminaria digitata, Alaria esculenta, Ecklonia radiata, Undaria pinnatifida, and Agarum clathratum. The average divergence for the 5’-COI between these species is 0.1011 ± 0.0051 but 1.5 times more, 0.1516 ± 0.0030, for the 3’-COI (paired samples t-test P = 0.0022; Mann–Whitney test P = 0.0090 (Table 1). The pattern is the same for all species pair comparisons: the 5’-COI (~ 740 bp in length) was significantly less diverged between S. latissima and other laminarialean algae than the 3’-COI (~ 862 bp in length) (Figure 4A – D; Table 1).
Figure 4

Sliding-window plots of divergence along the COI gene region between S. latissima and other laminarialean algae: S. japonica (A), the SHA form (B), S. coriacea (C), and S. angustata (D). Window sizes are 50 nucleotides with ten-nucleotide increments. The thick black lines at the bottom mark the COI regions investigated by McDevit and Saunders [43] and Silberfeld with colleagues [2].

Table 1

Pairwise distances (with SE) between Saccharina latissima and close species, based on different segments of COI





Species pair

Dxy ± SE

Dxy ± SE

Dxy ± SE

S. latissimaS. coriacea

0.0141 ± 0.0049

0.1488 ± 0.0135

0.0774 ± 0.0081

S. latissima – SHA

0.0106 ± 0.0042

0.1508 ± 0.0136

0.0765 ± 0.0080

S. latissima – TYP

0.0528 ± 0.0091

0.1647 ± 0.0138

0.1054 ± 0.0093

S. latissimaS. japonica

0.0528 ± 0.0091

0.1647 ± 0.0138

0.1054 ± 0.0093

S. latissimaS. angustata

0.0599 ± 0.0100

0.1607 ± 0.0138

0.1073 ± 0.0094

S. latissimaLaminaria digitata

0.0915 ± 0.0115

0.1488 ± 0.0143

0.1185 ± 0.0101

S. latissimaAgarum clathratum

0.0951 ± 0.0115

0.1607 ± 0.0146

0.1259 ± 0.0101

S. latissimaEcklonia radiata

0.1039 ± 0.0120

0.1548 ± 0.0140

0.1278 ± 0.0103

S. latissimaAlaria esculenta

0.0951 ± 0.0118

0.1508 ± 0.0137

0.1213 ± 0.0097

S. latissimaUndaria pinnatifida

0.1197 ± 0.0130

0.1429 ± 0.0142

0.1306 ± 0.0103

S. latissimaLeathesia difformis

0.1655 ± 0.0144

0.1290 ± 0.0134

0.1483 ± 0.0110

S. latissimaAsperococcus bullosus

0.1532 ± 0.0143

0.1349 ± 0.0138

0.1446 ± 0.0104

S. latissimaPunctaria latifolia

0.1602 ± 0.0143

0.1349 ± 0.0138

0.1483 ± 0.0110

S. latissimaElachista fucicola

0.1989 ± 0.0166

0.1627 ± 0.0142

0.1819 ± 0.0115

S. latissimaHincksia granulosa

0.1549 ± 0.0144

0.1290 ± 0.0131

0.1427 ± 0.0109

S. latissimaScytosiphon lomentaria

0.1637 ± 0.0148

0.1349 ± 0.0136

0.1502 ± 0.0110

S. latissimaPetalonia fascia

0.1532 ± 0.0145

0.1468 ± 0.0137

0.1502 ± 0.0109

S. latissimaColpomenia peregrina

0.1495 ± 0.0138

0.1468 ± 0.0144

0.1483 ± 0.0107

S. latissimaEctocarpus siliculosus

0.1373 ± 0.0139

0.1627 ± 0.0143

0.1493 ± 0.0109

S. latissimaPylaiella littoralis

0.1585 ± 0.0143

0.1607 ± 0.0144

0.1595 ± 0.0113

5’-COI: the 658-bp fragment covering the 5’-flanking region of the COI gene. The region starts 123 bp downstream of the COI start codon and ends 822 bp upstream of the COI stop codon. This fragment has been recommended for algae “barcoding” (species identification) [4143, 61]. 3’-COI: the 597-bp fragment covering the 3’-flanking region of the COI gene that starts 781 bp downstream of the COI start codon and ends 225 bp upstream of the COI stop codon. This fragment was not investigated by the above authors, but has been recently used (along with the 5’-COI region) for algae phylogenetic reconstruction [2]. Full COI: the 1378-bp fragment covering most of the COI gene. The fragment starts 123 bp downstream of the COI start codon and ends 225 bp upstream of the COI stop codon (the full COI region represents the largest sequence available for laminarialean algae in GenBank, excluding species for which full mtDNA sequences have been obtained). See Additional file 3 for the GenBank accession numbers.

The difference in the level of divergence between the S. latissima and other species for the 5’-COI and 3’-COI regions is highly pronounced (Figure 4; Table 1). The BLAST procedure (limited by the current GenBank submissions) reveals, as expected, very high identity between the 5’-COI region of S. latissima and other Saccharina species; however, the 3’-COI region of S. latissima has paradoxical similarity (but not complete identity) to a number of species belonging to the order Ectocarpales. Indeed, there is a noticeable decrease in the divergence for the 3’-COI region between S. latissima and the following Ectocarpales species: Leathesia difformis, Asperococcus bullosus, Punctaria latifolia, Elachista fucicola, Hincksia granulosa, Scytosiphon lomentaria, and Petalonia fascia. The average divergence for the 5’-COI between S. latissima and these species is 0.1642 ± 0.0061, but significantly less, 0.1389 ± 0.0046, for the 3’-COI (paired samples t-test P = 0.0007; Mann–Whitney test P = 0.0087; Table 1). Other Ectocarpales species for which long enough COI sequences are available (Pylaiella littoralis, Ectocarpus siliculosus, and Colpomenia peregrina) do not show any difference in the divergence level between 5’-COI and 3’-COI regions (Table 1).

Phylogenetically discordant signals within the COIgene

Thus, the difference in the level of divergence between the 5’-COI and 3’-COI regions is significant and exhibits an opposite pattern for comparisons between S. latissima vs. Laminariales species and S. latissima vs. Ectocarpales species: the 5’-COI region of S. latissima has obvious similarity to Saccharina species, whereas the 3’-COI regions has unexpected similarity with some Ectocarpales species. As a consequence, the position of S. latissima on 5’- and 3’-COI based trees are sharply different (Figure 5); on the 5’-COI based tree, S. latissima is within the same cluster with other Saccharina species but not distinguishable from S. coriacea (Figure 5A). On the 3’-COI tree S. latissima is significantly different from Laminariales algae and clusters with some species of the order Ectocarpales (Figure 5B). On the full length COI tree, S. latissima is within the order Laminariales (Figure 5C), but significantly different from other Saccharina species. The high similarity of the S. latissima 3’-COI region with Ectocarpales is surprising considering the large evolutionary distance between the two brown algae orders, which separated around 100 Ma [2].
Figure 5

Phylogenetic trees based on different fragments of the COI coding region: (A) 5’- COI , (B) 3’- COI , and (C) full COI region. Representative sequences of the orders Laminariales (L) and Ectocarpales (E) included in these trees are marked by vertical lines. Red algae COI sequences of Callithamnion pikeanum and Hypnea nidulans (GenBank accession numbers EU194965 and FJ694907, respectively) are used as outgroups. Note the changed position of S. latissima (in bold) depending on the COI region used for the tree. The S. latissima COI sequence denoted as “S” is from [2]. See Additional file 3 for GenBank accession numbers. Other comments as in Table 1 and Figure 3.

Recombination in the COIgene

Based on these observations we suspected that the unusual patterns of the COI gene divergence in comparisons of S. latissima with other species might reflect historical hybridization event(s) between representatives of distantly related Ectocarpales and Laminariales algae, which might have resulted in interspecific recombination of mitochondrial DNA. We therefore analyzed the COI alignments for evidence of recombination (and breakpoints) using various recombination detection methods implemented in the program RDP3 [62] (Table 2). The parental and recombinant sequences were determined using the VisRD method [63], modified version of PHYLPRO [64], and EEEP [65] also implemented in RDP3 (default settings).
Table 2

Recombination assessed by six different methods (implemented in RDP3; see [ [62]])



Recomb. species

AverageP - value



S. latissima

2.739 × 10-12

C. retorta

3.167 × 10-13



S. latissima

3.815 × 10-11

C. retorta

1.179× 10-15


[68, 69]

S. latissima

2.170 × 10-06

C. retorta

1.416 × 10-13



S. latissima

3.182 × 10-10

C. retorta

5.361× 10-08



S. latissima

1.459 × 10-08

C. retorta

4.560× 10-11



S. latissima

1.071 × 10-23


C. retorta

4.285× 10-25

All six methods detected a single recombination event in S. latissima within the 1,257 bp COI region (Figure 6A), with high statistical support (Table 2). The COI sequences from Petalonia fascia (Ectocarpales) and the SHA form (Laminariales) were major and minor parents, respectively. The recombination event involved the 3’-COI region, encompassing a 542 bp-long segment, between positions 1227 and 685.
Figure 6

Schematic representation of recombination events in the COI gene region from Saccharina latissima (A) and Cystophora retorta (B). The parental sequences are from Petalonia fascia and the SHA form (for S. latissima ); and Bifurcaria bifurcata and Phyllospora comosa (for C. retorta ).

Thus, the parental COI sequences of S. latissima come from different algae orders, Ectocarpales and Laminariales. The recombination would not have been detected if the alignment would be limited to the 5’-COI barcoding region. An interesting feature of algal evolution would not have been detected. Moreover, incorrect phylogenies would be obtained for species with historical hybridization events (see Figure 5).

To understand the uniqueness of this case of historical hybridization revealed in S. latissima, we analyzed additional, long enough COI sequences (~ 1200 bp; Additional file 3) of brown algae obtained from GenBank. We analyzed 57 sequences (GeneBank’s plus our own sequences) and detected one more highly significant instance of COI recombination between representatives of brown algae belonging to the families Sargassacea and Seirococcaceae. A recombinant COI sequence was revealed in Cystophora retorta (order Fucales, family Sargassacea). The COI sequences from Phyllospora comosa (order Fucales, family Seirococcaceae) and Bifurcaria bifurcata (order Fucales, family Sargassacea) were the major and minor parents, respectively. The recombination spanned the 5’-COI region (Figure 6B) with a beginning breakpoint at position 540 and ending breakpoint position 16. This recombination event involves representatives of two brown algae families, Sargassacea and Seirococcaceae, which diverged around 73 Ma [2]. Also, as in case of S. latissima, different and misleading phylogenetic patterns would come about depending on the COI region used in tree reconstruction (data not shown).

The recombination analysis has confirmed that the evolutionary histories of S. latissima and C. retorta have been influenced by gene flow between two highly diverged algae lineages, the orders Ectocarpales and Laminariales (for S. latissima) and between the families Sargassacea and Seirococcaceae (for C. retorta). The COI gene in S. latissima and C. retorta represents historical recombinant products between two distantly related ancestors. The level of COI divergence is 14.3% and 14.7% between the parental sequences of C. retorta and S. latissima, respectively, which seems very high for the occurrence of successful recombination (see, for instance, [73]). However, it has been shown [74] that recombination may occur among highly divergent maternal and paternal mtDNA sequences of the sea mussel Mytilus galloprovincialis that differ by > 20%.

There are multiple evidences of mtDNA recombination in plant, fungi, and animals, including human (reviews in [7577]). Mitochondrial heteroplasmy (the presence of more than one mitochondrial genome in an individual) that could facilitate recombination has been reported in the brown alga Fucus[78, 79]. Hybrid speciation is widespread in plants [8082] and mosaic genomes due to hybrid speciation are common in many marine microorganisms, plants and animals (review in [83]). Interspecific hybridization has been detected in the brown algal genus Fucus[8487] and the red algal genus Porphyra[88]. However, instances of interspecific mtDNA recombination were not, to our knowledge, previously known for algae; it is an infrequent phenomenon detected in fish and primates [89, 90].

The two examples of interspecific recombination come from published work by Silberfeld with colleagues [2]. It might be suggested that these recombinations are artificial, resulting from PCR errors. However, the architecture of the recombinant sequences is not a mosaic composition from different parts of the parental sequences; the recombinants are similar but not identical to any other sequence used in this study. These observations suggest that PCR errors are unlikely to account for these recombinant sequences. Nevertheless, more sampling will be necessary in order to settle this issue of putative hybrid speciation in algae.

Our results are relevant concerning the DNA barcoding for algae and possibly other organisms. The 5’-COI “barcode” region is not representative and might be even misleading (at least in case of S. latissima and C. retorta) in resolving taxonomic relationships between algal species. Although 5’-COI barcoding is practically convenient (because 600–700 bp can readily be sequenced in a single run); it may not be suitable.


We have found that the TYP and LON morphological forms of Saccharina japonica are genetically very similar and, therefore, they might not be thought of as distinct biological species, as suggested by Gusarova and Ivanova [17]. The SHA form, however, exhibits distinctive genetic differences from the TYP and LON forms (2.6% total DNA divergence) and it is closely related to the S. cichorioides species group. The SHA form might be thought of as a morphological variety of S. cichorioides that inhabits an ecological niche – the supralittoral zone – new for this species. Our results are of practical consequence. Zoospores of the SHA form are used to replenish populations of the TYP form suffering from overharvesting [19]. This practice is inappropriate because the TYP and SHA forms represent different species. Genetic analysis can, indeed, help to select appropriate stocks for mariculture.

Thus, the genetic data suggest the existence of two close Saccharina lineages inhabiting the Primorye coastal region, S. japonica (TYP + LON) and S. cichorioides (SHA). Significant genetic differences between the SHA form and the TYP plus LON forms are not surprising: these two lineages are also highly different in morphology and ecological preferences. Less expected is the nearly complete genetic identity of the TYP and LON forms, because these forms have drastic differences in morphology, life history traits, and ecology. Taking into account that interbreeding between the LON and TYP forms is highly unlikely (see Background), and also the fact that phenotypic plasticity may have potential roles for adaptive evolution and speciation, as it has been shown for many species, including algae (e.g., [2731]), we propose that the TYP and LON forms represent important resources for algae diversification that enable successful adaptation to contrasting shallow- and deep-water marine environments. The present results for Saccharina are in accordance with our previous results concerning two morphological forms of the sea urchin Strongylocentrotus intermedius, which are clearly different morphologically but very similar genetically [91].

Transplant site experiments are not likely to be informative concerning the taxonomical status of the Saccharina morphological forms since both, the LON form, which is genetically similar to the TYP form, and the SHA form, which is genetically quite different from the TYP form, exhibit drastic morphological transformations when transplanted to the TYP habitat area, so that they become morphologically indistinguishable from the TYP form [23, 24]. Breeding experiments would not likely be informative either to resolve the taxonomic status of the TYP and LON forms, because even distantly related algal species produce highly viable first generation hybrids [88, 9295]. Indeed, the biological species concept is not particularly helpful in kelps [92, 96]. Kraan and Guiry [92] have, for example, shown that interspecific sequence divergence in Alaria is smaller than intraspecific sequence divergence. That is, there is more variation within a species than between two species of the same genus, which casts doubts on the morphological and biological species concepts used in Alaria and on the usefulness of hybridization studies in assessing species-level differences [97100]. Conflict between speciation decisions based on morphological or molecular characteristics comes about due to the drastically different rates at which molecular and morphological changes accumulate [100]. The use of a morphological and biological species concept to separate Alaria or Saccharina species is not satisfactory and does not fully reflect phylogenetic relationships.

The data show that there is no consistent relationship between morphological and genetic variation in the algae Saccharina (present data) and Alaria[96], as well as in sea urchins [91]. The results of Druehl and Saunders [96] and Kraan and Guiry [92] demonstrate that fertility barriers, which in most cases indicate complete reproductive isolation, may arise without affecting genetic divergence in the particular genes under investigation [101]. Our results support a relatively minor contribution of genetic factors to the TYP and LON morphological forms differences observed in the three gene markers, suggesting an important phenotypic plasticity basis for the variability in morphology, life history, and ecology in S. japonica.

The relative importance of genes and phenotypic plasticity in the generation of adaptive variation is a far from resolved issue [102]. The possibility that plasticity could play a major role, however, has recently become more commonly considered than it was traditionally believed [103105]. A body of evidence suggests that plasticity may promote adaptive divergence in various systems, often followed by genetic changes in the direction of the plastic response [106] (see however, [107]). Furthermore, phenotypic plasticity enhances the survival and reproductive success of individuals by contributing to their ability to cope with environmental changes. In this way, it enables potential adaptation to new niches [108], and therefore can promote important biological processes such as adaptation, divergence, reproductive isolation, and evolutionary innovation [28, 109, 110]. Environmentally induced phenotypic variation has been argued to be a more likely source of evolutionary novelty than new mutations [103, 104]. It has, indeed, been asserted that phenotypic plasticity can accelerate adaptive evolution by shortening the time a population needs to discover a new genotypic network [111] even in the case of recurrent rapid changes in gene regulation [112].

A number of experimental studies have revealed the mechanisms of phenotypic plasticity using an integrative biology approach (review in [113]). The overall results indicate that phenotypic plasticity has an important role in adaptive evolution. This consideration, in turn, points out the importance of characterizing diverging phenotypes and identifying relevant evolutionary forces acting on those phenotypes. A great number of intraspecific morphological varieties are known for laminarialean algae (review in [14]). Further studies of algae morphological forms are important and necessary to further understand how biological diversity is generated and maintained in evolution.


  1. 1.

    We have investigated the genetic make-up in three S. japonica forms, TYP, LON, and SHA. We have found that in spite of their drastic differences in morphology, life history and ecology, the TYP and LON forms are genetically indistinguishable and clearly different from the SHA form, which has a close genetic relationship with S. cichorioides. Taking into account the potential evolutionary significance of phenotypic plasticity, the TYP and LON forms may represent an important resource for algal evolutionary diversification and successful adaptation to contrasting marine environments.

  2. 2.

    The three gene markers used in the present work have different sensitivity for algal species identification; the mitochondrial COI gene is the most efficient gene marker. However, we found that the designated 5’-COI “barcoding” region is not enough and may even be misleading for purposes of identification, at least in comparisons that include species with evolutionary histories involving gene flow between distantly related algae lineages.

  3. 3.

    Populations of the TYP and LON forms are depressed owing to overharvesting; restriction programs are in play to enhance these populations. But the correct foundation stocks must be used. Zoospores of the SHA form should not be used to seed the fields of the TYP and LON forms because they belong to different evolutionary lineages. Genetic determination should be used to select correct foundation stocks in mariculture and management programs.


Note added in proof

Since the work described in this paper was completed and submitted for publication, sequences of the COI gene (most sensitive gene marker; see text) were obtained for typical individuals of Saccharina cichorioides. We detected no differences between the SHA form and S. cichorioides, which supports our conclusion about the conspecificity of the SHA form and S. cichorioides. The COI sequences have been deposited in GenBank under accession numbers JQ792007-JQ792010.



We thank Elena Balakireva for encouragement and help and Vladimir A. Pavlyuchkov for help with algal collection. This study was supported by Bren Professor Funds at the University of California Irvine to F. J. A. and E. S. B.

Authors’ Affiliations

Department of Ecology and Evolutionary Biology, University of California
A. V. Zhirmunsky Institute of Marine Biology, Far Eastern Branch of the Russian Academy of Science
Pacific Research Fisheries Centre (TINRO-Centre)


  1. Lane CE, Mayes C, Druehl LD, Saunders GW: A multi-gene molecular investigation of the kelp (Laminariales, Phaeophyceae) supports substantial reorganisation. J Phycol. 2006, 42: 493-512. 10.1111/j.1529-8817.2006.00204.x.
  2. Silberfeld T, Leigh JW, Verbruggen H, Cruaud C, de Reviers B, Rousseau F: A multi-locus time-calibrated phylogeny of the brown algae (Heterokonta, Ochrophyta, Phaeophyceae): Investigating the evolutionary nature of the "brown algal crown radiation". Mol Phylogenet Evol. 2010, 56: 659-674. 10.1016/j.ympev.2010.04.020.PubMed
  3. Verbruggen H, Maggs CA, Saunders GW, Le Gall L, Yoon HS, De Clerck O: Data mining approach identifies research priorities and data requirements for resolving the red algal tree of life. BMC Evol Biol. 2010, 10: 16-10.1186/1471-2148-10-16.PubMedPubMed Central
  4. Miуabe K: Laminariaceae of Hokkaido. Rep Fish Dept Hokkaido-tyo. 1902, 3: 1-60.
  5. Zinova ES: Marine cabbage (Laminaria) and other commercially important algae. Izv TINRO. 1928, 1: 7-42.
  6. Petrov YuE: Systematics of some species of the genus Laminariae Lamour. from the Far East. Novit Syst Plant Non Vascularium. 1972, 9: 47-59.
  7. Tokida J, Nakamura Y, Druehl LD: Typification of species of Laminaria (Phaeophyta, Laminariales) described by Miyabe, and taxonomic notes on the genus in Japan. Phycologia. 1980, 19: 317-328. 10.2216/i0031-8884-19-4-317.1.
  8. Yotsukura N, Denboh T, Motomura T, Horiguchi T, Coleman AW, Ichimura T: Little divergence in ribosomal DNA internal transcribed spacer −1 and −2 sequences among non-digitate species of Laminaria (Phaeophyceae) from Hokkaido, Japan. Phycol Res. 1999, 47: 71-80. 10.1111/j.1440-1835.1999.tb00286.x.
  9. Yotsukura N, Kawai T, Kawashima S, Ebata H, Ichimura T: Nucleotide sequence diversity of the 5 S rDNA spacer in the simple blade kelp genera Laminaria, Cymathaere and Kjellmaniella (Laminariales, Phaeophyceae) from northern Japan. Phycol Res. 2006, 54: 269-279. 10.1111/j.1440-1835.2006.00434.x.
  10. Yotsukura N, Kawashima S, Kawai T, Abe T, Druehl LD: A systematic re-examination of four Laminaria species: L. japonica, L. religiosa, L. ochotensis, and L. diabolica. J Jpn Bot. 2008, 83: 165-176.
  11. Yotsukura N, Shimizu T, Katayama T, Druehl LD: Mitochondrial DNA sequence variation of four Saccharina species in Japan. J Appl Phycol. 2010, 22: 243-251. 10.1007/s10811-009-9452-7.
  12. Yotsukura N: Molecular phylogeny of advanced kelps (Laminariales, Phaeophyceae) growing in Japan. Nat Hist Res Spec. 2005, 8: 69-81.
  13. Selivanova ON, Zhigadlova GG, Hansen GI: Revision of the systematics of algae in the order Laminariales (Phaeophyta) from the Far-Eastern seas of Russia on the basis of molecular-phylogenetic data. Russ J Mar Biol. 2007, 33: 278-289. 10.1134/S1063074007050021.
  14. Bartsch I, Wiencke C, Bischof K, Buchholz CM, Buck BH, Eggert A, Feuerpfeil P, Hanelt D, Jacobsen S, Karez R, Karsten U, Molis M, Roleda MY, Schubert H, Schumann R, Valentin K, Weinberger F, Wiese J: The genus Laminaria sensu lato: recent insights and developments. Eur J Phycol. 2008, 43: 1-86. 10.1080/09670260701711376.
  15. Gail GI: Laminariaceous algae of the Far Eastern seas. Vestnik DV FAN SSSR. 1936, 19: 31-64.
  16. Sukhoveeva MV: The algae distribution of along the Primorye coast region. Izv TINRO. 1967, 61: 255-260.
  17. Gusarova IS, Ivanova NV: Intraspecific systematic of Laminaria japonica at the continental coast of the Japan Sea. Izv TINRO. 2006, 147: 157-168.
  18. Krupnova TN: Influence of oceanology-climatic factors on dynamics of Laminaria japonica Aresch. fields in northwest part of the Sea of Japan. Sea Hydrology and Hydrochemistry: 8(2):162–166. Hydrometeoizdat: Japan Sea. Leningrad; 2004.
  19. Krupnova TN: Instruction on Cultivation and Recovery of Laminaria Fields. Vladivostok: Russia: Pacific Research Fisheries Centre (TINRO-Centre); 2008.
  20. Paimeeva LG: Distribution and growth of Laminaria japonica Aresch f. longipes (Miyabe et Tokida) Petr. in northern Primorye. In Commercial Algae and Their Use. 1987, Moscow: Nauka, 26-33.
  21. Paimeeva LG, Gusarova IS: Laminaria japonicaAresch. f.longipes(Miyabe et Tokida) Ju. Petr. bed condition in Northern Primorye. Komarov Reading. 1993, 38: 20-36. Vladivostok: Dal'nauka
  22. Hasegawa YJ: An ecological study of Laminaria angustat Kjellman on the coast of Hidaka Prov., Hokkaido. Bull Hokkaido Reg Fish Lab. 1962, 24: 116-138.
  23. Krupnova TN, Pavlyuchkov VA: Developing biotechnology to increase commercial stocks of sea urchins by creating a favorable food supply. TINRO archive. 1999, 23285: 38-
  24. Krupnova TN, Pavlyuchkov VA: Biotechnology materials for the combined cultivation of kelp and sea urchins. TINRO archive. 2000, 23647: 54-
  25. Krupnova TN: Reorganization of Laminaria fields under the influence of natural and anthropogenic factors. TINRO archive. 2010, 26962: 56-
  26. Petrov YuE, Sukhoveeva MV: Laminaria angustata Kjellm. at the coasts of Primorskii region. Novit Syst Plant Non Vascularium. 1972, 9: 44-47.
  27. Price TD, Qvarnström A, Irwin DE: The role of phenotypic plasticity in driving genetic evolution. Proc R Soc Lond B. 2003, 270: 1433-1440. 10.1098/rspb.2003.2372.
  28. Pfennig DW, Wund MA, Snell-Rood EC, Cruickshank T, Schlichting CD, Moczek AP: Phenotypic plasticity’s impacts on diversification and speciation. Trends Ecol Evol. 2010, 25: 459-467. 10.1016/j.tree.2010.05.006.PubMed
  29. Thibert-Plante X, Hendry AP: The consequences of phenotypic plasticity for ecological speciation. J Evol Biol. 2011, 24: 326-342. 10.1111/j.1420-9101.2010.02169.x.PubMed
  30. Hays CG: Adaptive phenotypic differentiation across the intertidal gradient in the alga Silvetia compressa. Ecology. 2007, 88: 149-157. 10.1890/0012-9658(2007)88[149:APDATI]2.0.CO;2.PubMed
  31. Monro K, Poore AGB: Performance benefits of growth-form plasticity in a clonal red seaweed. Biol J Linn Soc. 2009, 97: 80-89. 10.1111/j.1095-8312.2008.01186.x.
  32. Balakirev ES, Balakirev EI, Rodriguez-Trelles F, Ayala FJ: Molecular evolution of two linked genes, Est-6 and Sod, in Drosophila melanogaster. Genetics. 1999, 153: 1357-1369.PubMedPubMed Central
  33. Balakirev ES, Balakirev EI, Ayala FJ: Molecular evolution of the Est-6 gene in Drosophila melanogaster: Contrasting patterns of DNA variability in adjacent functional regions. Gene. 2002, 288: 167-177. 10.1016/S0378-1119(02)00477-8.PubMed
  34. Balakirev ES, Chechetkin VR, Lobzin VV, Ayala FJ: DNA polymorphism in the β-esterase gene cluster of Drosophila melanogaster. Genetics. 2003, 164: 533-544.PubMedPubMed Central
  35. Oudot-Le Secq M-P, Fontaine J-M, Rousvoal S, Kloareg B, Loiseaux-de Goër S: The complete sequence of a brown algal mitochondrial genome, the Ectocarpale Pylaiella littoralis (L.) Kjellm. J Mol Evol. 2001, 53: 80-88.PubMed
  36. Oudot-Le Secq M-P, Kloareg B, Loiseaux-de Goër S: The mitochondrial genome of the brown alga Laminaria digitata: a comparative analysis. Eur J Phycol. 2002, 37: 163-172. 10.1017/S0967026202003542.
  37. Loiseaux-de Goër S, Stam WT, Olsen JL, Oudot-Le Secq M-P: Complete mitochondrial genomes of the three brown algae (Heterokonta: Phaeophyceae) Dictyota dichotoma, Fucus vesiculosus and Desmarestia viridis. Curr Genet. 2006, 49: 47-58. 10.1007/s00294-005-0031-4.PubMed
  38. Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 1994, 22: 4673-4680. 10.1093/nar/22.22.4673.PubMedPubMed Central
  39. Librado P, Rozas J: DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009, 25: 1451-1452. 10.1093/bioinformatics/btp187.PubMed
  40. Filatov DA: PROSEQ: a software for preparation and evolutionary analysis of DNA sequence data sets. Mol Ecol Notes. 2002, 2: 621-624. 10.1046/j.1471-8286.2002.00313.x.
  41. Lane CE, Lindstrom SC, Saunders GW: A molecular assessment of northern Pacific Alaria species (Laminariales, Phaeophyceae) with reference to the utility of DNA barcoding. Mol Phylogenet Evol. 2007, 44: 634-648. 10.1016/j.ympev.2007.03.016.PubMed
  42. McDevit DC, Saunders GW: On the utility of DNA barcoding for species differentiation among brown macroalgae (Phaeophyceae) including a novel extraction protocol. Phycol Res. 2009, 57: 131-141. 10.1111/j.1440-1835.2009.00530.x.
  43. McDevit DC, Saunders GW: A DNA barcode examination of the Laminariaceae (Phaeophyceae) in Canada reveals novel biogeographical and evolutionary insights. Phycologia. 2010, 49: 235-248. 10.2216/PH09-36.1.
  44. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: Molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMedPubMed Central
  45. Posada D: jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008, 25: 1253-1256. 10.1093/molbev/msn083.PubMed
  46. Schwarz GE: Estimating the dimension of a model. Ann Stat. 1978, 6: 461-464. 10.1214/aos/1176344136.
  47. Akaike H: A new look at the statistical model identification. IEEE Trans Automat Contr. 1974, 19: 716-723. 10.1109/TAC.1974.1100705.
  48. Felsenstein J: Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985, 39: 783-791. 10.2307/2408678.
  49. Zwickl DJ: Genetic algorithm approaches for the phylogenetic analysis of large biological sequence datasets under the maximum likelihood criterion. 2006, Austin: The University of Texas, Ph.D. thesis
  50. Ronquist F, Teslenko M, van der Mark P, Ayres D, Darling A, Höhna S, Larget B, Liu L, Suchard MA, Huelsenbeck JP: 2012. MrBayes 3.2: Efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012, 61: 1-4. 10.1093/sysbio/syr074.
  51. Gelman A, Rubin DB: Inference from iterative simulation using multiple sequences. Stat Sci. 1992, 7: 457-511. 10.1214/ss/1177011136.
  52. Saitou N, Imanishi T: Relative efficiencies of the Fitch-Margoliash, Maximum-Parsimony, Maximum-Likelihood, Minimum-Evolution, and Neighbor-joining methods of phylogenetic tree construction in obtaining the correct tree. Mol Biol Evol. 1989, 6: 514-525.
  53. Tateno Y, Takezaki N, Nei M: Relative efficiencies of the Maximum-Likelihood, Neighbor-joining, and Maximum-Parsimony methods when substitution rate varies with site. Mol Biol Evol. 1994, 11: 261-277.PubMed
  54. Russo CAM, Takezaki N, Nei M: Efficiencies of different genes and different tree-building methods in recovering a known vertebrate phylogeny. Mol Biol Evol. 1996, 13: 525-536. 10.1093/oxfordjournals.molbev.a025613.PubMed
  55. Robba L, Russell SJ, Barker GL, Brodie J: Assessing the use of the mitochondrial COX1 marker for use in DNA barcoding of red algae (Rhodophyta). Amer J Botany. 2006, 93: 1101-1108. 10.3732/ajb.93.8.1101.
  56. Fraser CI, Nikula R, Spencer HG, Waters JM: Kelp genes reveal effects of subantarctic sea ice during the Last Glacial Maximum. Proc Natl Acad USA. 2009, 106: 3249-3253. 10.1073/pnas.0810635106.
  57. Clarkston BE, Saunders GW: A comparison of two DNA barcode markers for species discrimination in the red algal family Kallymeniaceae (Gigartinales, Florideophyceae), with a description of Euthora timburtonii sp. nov. Botany. 2010, 88: 119-131. 10.1139/B09-101.
  58. Shan TF, Pang SJ, Zhang YR, Yakovleva IM, Skriptsova AV: An AFLP-based survey of genetic diversity and relationships of major farmed cultivars and geographically isolated wild populations of Saccharina japonica (Phaeophyta) along the northwest coasts of the Pacific. J Appl Phycol. 2010, 23: 35-45.
  59. Avise JC: Phylogeography: the History and Formation of Species. 2000, Cambridge, MA: Harvard University Press
  60. Selivanova ON, Yotsukura N, Kawashima S: Comparison of some Laminaria species from Pacific coasts of Russia and Japan. Phycologia. 2005, 44: 92-93.
  61. Saunders GW: Applying DNA barcoding to red macroalgae: a preliminary appraisal holds promise for future applications. Phil Trans R Soc B. 2005, 60: 1879-1888.
  62. Martin DP, Lemey P, Lott M, Moulton V, Posada D, Lefeuvre P: RDP3: a flexible and fast computer program for analyzing recombination. Bioinformatics. 2010, 26: 2462-2463. 10.1093/bioinformatics/btq467.PubMedPubMed Central
  63. Lemey P, Lott M, Martin DP, Moulton V: Identifying recombinants in human and primate immunodeficiency virus sequence alignments using quartet scanning. BMC Bioinforma. 2009, 10: 126-10.1186/1471-2105-10-126.
  64. Weiller GF: Phylogenetic profiles: a graphical method for detecting genetic recombinations in homologous sequences. Mol Biol Evol. 1998, 15: 326-335. 10.1093/oxfordjournals.molbev.a025929.PubMed
  65. Beiko RG, Hamilton N: Phylogenetic identification of lateral genetic transfer events. BMC Evol Biol. 2006, 6: 15-10.1186/1471-2148-6-15.PubMedPubMed Central
  66. Martin D, Rybicki E: RDP: detection of recombination amongst aligned sequences. Bioinformatics. 2000, 16: 562-563. 10.1093/bioinformatics/16.6.562.PubMed
  67. Padidam M, Sawyer S, Fauquet CM: Possible emergence of new geminiviruses by frequent recombination. Virology. 1999, 265: 218-225. 10.1006/viro.1999.0056.PubMed
  68. Martin DP, Posada D, Crandall KA, Williamson C: A modified bootscan algorithm for automated identification of recombinant sequences and recombination breakpoints. Aids Res Hum Retrovir. 2005, 21: 98-102. 10.1089/aid.2005.21.98.PubMed
  69. Bredell H, Martin DP, Van Harmelen J, Varsani A, Sheppard HW, Donovan R, Gray CM, Williamson C, Team HS: HIV type 1 subtype C gag and nef diversity in southern Africa. Aids Res Hum Retrovir. 2007, 23: 477-481. 10.1089/aid.2006.0232.PubMed
  70. Smith JM: Analyzing the mosaic structure of genes. J Mol Evol. 1992, 34: 126-129.PubMed
  71. Posada D, Crandall KA: Evaluation of methods for detecting recombination from DNA sequences: computer simulations. Proc Natl Acad Sci USA. 2001, 98: 13757-13762. 10.1073/pnas.241370698.PubMedPubMed Central
  72. Gibbs MJ, Armstrong JS, Gibbs AJ: Sister-Scanning: a Monte Carlo procedure for assessing signals in recombinant sequences. Bioinformatics. 2000, 16: 573-582. 10.1093/bioinformatics/16.7.573.PubMed
  73. Balakirev ES, Ayala FJ: Pseudogenes: are they “junk” or functional DNA?. Annu Rev Genet. 2003, 37: 123-151. 10.1146/annurev.genet.37.040103.103949.PubMed
  74. Ladoukakis ED, Theologidis I, Rodakis GC, Zouros E: Homologous recombination between highly diverged mitochondrial sequences: examples from maternally and paternally transmitted genomes. Mol Biol Evol. 2011, 28: 1847-1859. 10.1093/molbev/msr007.PubMed
  75. Ballard JW, Whitlock MC: The incomplete natural history of mitochondria. Mol Ecol. 2004, 13: 729-744. 10.1046/j.1365-294X.2003.02063.x.PubMed
  76. Barr CM, Neiman M, Taylor DR: Inheritance and recombination of mitochondrial genomes in plants, fungi and animals. New Phytol. 2005, 168: 39-50. 10.1111/j.1469-8137.2005.01492.x.PubMed
  77. Galtier N, Nabholz B, Glémin S, Hurst GDD: Mitochondrial DNA as a marker of molecular diversity: a reappraisal. Mol Ecol. 2009, 18: 4541-4550. 10.1111/j.1365-294X.2009.04380.x.PubMed
  78. Coyer JA, Hoarau G, Stam WT, Olsen JL: Geographically specific heteroplasmy of mitochondrial DNA in the seaweed, Fucus serratus (Heterokontophyta: Phaeophyceae, Fucales). Mol Ecol. 2004, 13: 1323-1326. 10.1111/j.1365-294X.2004.02128.x.PubMed
  79. Hoarau G, Coyer JA, Olsen JL: Paternal leakage of mitochondrial DNA in a Fucus (Phaeophyceae) hybrid zone. J Phycol. 2009, 45: 621-624. 10.1111/j.1529-8817.2009.00679.x.
  80. Reisberg LH: Hybrid origins of plant species. Annu Rev Ecol Syst. 1997, 28: 359-389. 10.1146/annurev.ecolsys.28.1.359.
  81. Chapman MA, Burke JM: Genetic divergence and hybrid speciation. Evolution. 2007, 61: 1773-1780. 10.1111/j.1558-5646.2007.00134.x.PubMed
  82. Mallet J: Hybrid speciation. Nature. 2007, 446: 279-283. 10.1038/nature05706.PubMed
  83. Arnold ML, Fogarty ND: Reticulate evolution and marine organisms: the final frontier?. Int J Mol Sci. 2009, 10: 3836-3860. 10.3390/ijms10093836.PubMedPubMed Central
  84. Coyer JA, Peters AF, Hoarau G, Stam WT, Olsen JL: Hybridization of the marine seaweeds, Fucus serratus and F. evanescens (Heterokontophyta: Phaeophyceae) in a 100-year old zone of secondary contact. Proc R Soc Lond B. 2002, 269: 1829-1834. 10.1098/rspb.2002.2093.
  85. Coyer JA, Hoarau G, Pearson G, Serrão E, Stam WT, Olsen JL: Convergent adaptation to a marginal habitatby homoploid hybrids and polyploid ecads in the seaweed genus Fucus. Biol Lett. 2006, 2: 405-408. 10.1098/rsbl.2006.0489.PubMedPubMed Central
  86. Coyer JA, Hoarau G, Stam WT, Olsen JL: Hybridization and introgression in a mixed population of the intertidal seaweeds Fucus evanescens and F. serratus. J Evol Biol. 2007, 20: 2322-2333. 10.1111/j.1420-9101.2007.01411.x.PubMed
  87. Wallace A, Klein AS, Mathieson AC: Determining the affinities of salt marsh fucoids using microsatellite markers: evidence of hybridization and introgression between two species of Fucus (Phaeophyta) in a Maine estuary. J Phycol. 2004, 40: 1013-1027. 10.1111/j.1529-8817.2004.04085.x.
  88. Niwa K, Sakamoto T: Allopolyploidy in natural and cultivated populations of Porphyra (Bangiales, Rhodophyta). J Phycol. 2010, 46: 1097-1105. 10.1111/j.1529-8817.2010.00897.x.
  89. Piganeau G, Gardner M, Eyre-Walker A: A broad survey of recombination in animal mitochondria. Mol Biol Evol. 2004, 21: 2319-2325. 10.1093/molbev/msh244.PubMed
  90. Ciborowski KL, Consuegra S, García de Leániz C, Beaumont MA, Wang J, Jordan WC: Rare and fleeting: an example of interspecific recombination in animal mitochondrial DNA. Biol Lett. 2007, 3: 554-557. 10.1098/rsbl.2007.0290.PubMedPubMed Central
  91. Balakirev ES, Pavlyuchkov VA, Ayala FJ: DNA variation and endosymbiotic associations in phenotypically-diverse sea urchin Strongylocentrotus intermedius. Proc Nat Acad Sci USA. 2008, 105: 16218-16223. 10.1073/pnas.0807860105.PubMedPubMed Central
  92. Kraan S, Guiry MD: Sexual hybridization experiments and phylogenetic relationships as inferred from RUBISCO spacer sequences in the genus Alaria (Phaeophyceae). J Phycol. 2000, 35: 190-198.
  93. Liptack MK, Druehl LD: Molecular evidence for an interfamilial laminarialean cross. Eur J Phycol. 2000, 35: 135-142. 10.1080/09670260010001735721.
  94. Druehl LD, Collins JD, Lane CE, Saunders GW: An evaluation of methods used to assess intergeneric hybridization in kelp using Pacific Laminariales (Phaeophyceae). J Phycol. 2005, 41: 250-262. 10.1111/j.1529-8817.2005.04143.x.
  95. Xu P, Yang L, Zhu J, Xu H, Lu Q: Analysis of hybridization strains of Porphyra based on rbcL gene sequences. J Appl Phycol. 2011, 23: 235-241. 10.1007/s10811-010-9537-3.
  96. Druehl LD, Saunders GW: Molecular explorations in kelp evolution. Prog Phycol Res. 1992, 8: 47-83.
  97. Harrison RG: Molecular changes at speciation. Annu Rev Ecol Syst. 1991, 22: 281-308. 10.1146/
  98. Guiry MD: Species concepts in marine algae. Prog Phycol Res. 1992, 8: 251-278.
  99. Manhart JR, McCourt RM: Molecular data and species concepts in algae. J Phycol. 1992, 28: 730-737. 10.1111/j.0022-3646.1992.00730.x.
  100. Medlin LK, Lange M, Barker GLA, Hayes PK: Can molecular techniques change our ideas about the species concept?. In Molecular Ecology of Aquatic Microbes. Edited by: Joint I. Berlin: Springer;1995: 133-152. [NATO AS1 series, vol G38.]
  101. Hillis DM: Molecular versus morphological approaches to systematics. Annu Rev Ecol Syst. 1987, 18: 23-42. 10.1146/
  102. Schwander T, Leimar O: Genes as leaders and followers in evolution. Trends Ecol Evol. 2011, 26: 143-151. 10.1016/j.tree.2010.12.010.PubMed
  103. Schlichting CD, Pigliucci M: Phenotypic Evolution: a Reaction Norm Perspective. 1998, Sunderland, MA: Sinauer Associates
  104. West-Eberhard MJ: Developmental Plasticity and Evolution. New York: Oxford University Press;2003.
  105. West-Eberhard MJ: Developmental plasticity and the origin of species differences. Proc Natl Acad Sci USA. 2005, 102: 6543-6549. 10.1073/pnas.0501844102.PubMedPubMed Central
  106. Crispo E: Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow. J Evol Biol. 2008, 21: 1460-1469. 10.1111/j.1420-9101.2008.01592.x.PubMed
  107. Martínez-Fernández M, Páez de la Cadena M, Rolán-Alvarez E: The role of phenotypic plasticity on the proteome differences between two sympatric marine snail ecotypes adapted to distinct micro-habitats. BMC Evol Biol. 2010, 10: 65-10.1186/1471-2148-10-65.PubMedPubMed Central
  108. Prada C, Schizas NV, Yoshioka PM: Phenotypic plasticity or speciation?. A case from a clonal marine organism. BMC Evol Biol. 2008, 8: 47-PubMed
  109. Beldade P, Mateus AR, Keller RA: Evolution and molecular mechanisms of adaptive developmental plasticity. Mol Ecol. 2010, 20: 1347-1363.
  110. Moczek AP, Sultan S, Foster S, Ledón-Rettig C, Dworkin I, Nijhout HF, Abouheif E, Pfennig DW: The role of developmental plasticity in evolutionary innovation. Proc R Soc B. 2011, 278: 2705-2713. 10.1098/rspb.2011.0971.PubMedPubMed Central
  111. Espinosa-Soto C, Martin OC, Wagner A: Phenotypic plasticity can facilitate adaptive evolution in gene regulatory circuits. BMC Evol Biol. 2011, 11: 5-10.1186/1471-2148-11-5.PubMedPubMed Central
  112. Scoville AG, Pfrender ME: Phenotypic plasticity facilitates recurrent rapid adaptation to introduced predators. Proc Natl Acad Sci USA. 2010, 107: 4260-4263. 10.1073/pnas.0912748107.PubMedPubMed Central
  113. Aubin-Horth N, Renn SC: Genomic reaction norms: using integrative biology to understand molecular mechanisms of phenotypic plasticity. Mol Ecol. 2009, 18: 3763-3780. 10.1111/j.1365-294X.2009.04313.x.PubMed
  114. Krupnova TN: Development of sporogenous tissue in Laminaria japonica under the influence of changing environment. Izv TINRO. 2002, 130: 474-482.


© Balakirev et al.; licensee BioMed Central Ltd. 2012

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://​creativecommons.​org/​licenses/​by/​2.​0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.