Multiple origins and the population genetic structure of Rubus takesimensis (Rosaceae) on Ulleung Island: Implications for the genetic consequences of anagenetic speciation
aff001; Jae-Hong Pak
aff001; Masayuki Maki
aff002; Seung-Chul Kim
Authors place of work:
Department of Biology, Research Institute for Dok-do and Ulleung-do Island, Kyungpook National University, Daegu, Gyeongsangbuk-do, Korea
aff001; Botanical Gardens, Tohoku University, Sendai, Miyagi, Japan
aff002; Department of Biological Sciences, Sungkyunkwan University, Suwon, Gyeonggi-do, Korea
Published in the journal:
PLoS ONE 14(9)
To determine the origin and genetic consequences of anagenesis in Rubus takesimensis on Ulleung Island, Korea, we compared the genetic diversity and population structure of R. takesimensis with those of its continental progenitor R. crataegifolius. We broadly sampled a total of 315 accessions in 35 populations and sequenced five noncoding regions of chloroplast DNA. Rubus takesimensis emerged as nonmonophyletic and several geographically diverse continental populations were likely responsible for the origin of R. takesimensis; the majority of R. takesimensis accessions were sisters to the clade containing accessions of R. crataegifolius, primarily from the Korean peninsula, while rare accessions from three populations shared common ancestors with the ones from the southern part of the Korean peninsula, Jeju Island, and Japan. A few accessions from the Chusan population originated independently from the Korean peninsula. Of 129 haplotypes, 81 and 48 were found exclusively in R. crataegifolius and R. takesimensis, respectively. We found unusually high genetic diversity in two regions on Ulleung Island and no geographic population structure. For R. crataegifolius, two major haplotype groups were found; one for the northern mainland Korean peninsula, and the other for the southern Korean peninsula and the Japanese archipelago. Compared with populations of R. crataegifolius sampled from Japan, much higher haplotype diversity was found in populations from the Korean peninsula. The patterns of genetic consequences in R. takesimensis need to be verified for other endemic species based on chloroplast DNA and independent nuclear markers to synthesize emerging patterns of anagenetic speciation on Ulleung Island.
Biology and life sciences – Genetics – Heredity – Genetic mapping – Haplotypes – Population genetics – Phylogeography – Genetic polymorphism – Evolutionary biology – Population genetics – Phylogeography – Genetic polymorphism – Population biology – Population genetics – Phylogeography – Genetic polymorphism – Biogeography – Phylogeography – Ecology – Ecological metrics – Species diversity – Earth sciences – Geomorphology – Topography – Landforms – Islands – Geography – Biogeography – Phylogeography – Ecology and environmental sciences – Biogeography – Phylogeography – Ecology – Ecological metrics – Species diversity – People and places – Geographical locations – Asia – Japan – Korea
On oceanic islands, approximately 75% of the insular endemic plant species have originated presumably via cladogenesis, which is typically associated with the term adaptive radiation in insular setting [1,2]. This mode of speciation involves a founding population splitting into several phenotypically and ecologically distinct lineages that can be recognized at a species (or even generic) level; consequently, the initial founding gene pool is partitioned and subjected to intense selection in different habitats, resulting in a low level of genetic variation within each species [1,3]. To reveal and elucidate this mode of speciation, numerous independent research groups have focused on a diverse assemblage of endemic plants on several archipelagoes, such as Hawaii and the Galápagos in the Pacific [4–7] and the Canaries in the Atlantic [8–13]. Also, similar interests have been focused on relatively smaller oceanic archipelagos, such as the Juan Fernández Islands [14–17] and the Bonin Islands [18–19]. The remaining 25% of island endemic plant species, however, likely are derived via anagenetic speciation . This mode of speciation involves an initial founding population on an island, which gradually accumulates genetic variation by mutation and recombination and eventually, without any splitting events, transforms into morphologically and genetically distinct taxonomic forms [1,2]. Unlike the numerous examples of adaptive radiation, the mode of anagenetic speciation on oceanic islands have recently been investigated by a few research groups, primarily in the Juan Fernández Islands in the Pacific [17,20–22] and Ulleung Island in East Asia [23–26].
Ulleung Island, an oceanic volcanic island in the East Sea/Sea of Japan, is located about 137 km off the east coast of the Korean peninsula. The age of the island is approximately 1.8 million years old (myo) and has never been connected to the adjacent continents . A high level of endemism in oceanic archipelagos, such as the Hawaiian Islands (828 species) and the Canary Islands (429 species) is often associated with cladogenetic speciation (95% for Hawaii and 84% for the Canaries), but Ulleung Island is known for its exceptionally high level of anagenetic speciation (33 species; 88% of the total endemic species) . It has been hypothesized that low-elevation islands with low habitat heterogeneity are highly correlated with high levels of anagenetic speciation [1,26]. As of today, the outstanding studies focused on anagenetic speciation on Ulleung Island have been limited to two core eudicot families, i.e., Dystaenia takesimana (asteroid II, Apiales, Apiaceae)  and Acer okamotoanum and Acer takesimense (rosid II/Malvidae, Sapindales, Sapindaceae) [24,25]. In a recent review of anagenetic speciation in the endemic angiosperms of Ulleung Island, Takayam et al.  summarized that, unlike some previous results based on isozyme studies, AFLP (amplified fragment length polymorphisms) and SSR (simple sequence repeats) molecular markers revealed much higher levels of genetic variation in some species, highlighting the need for the appropriate use of molecular markers. They also pointed out that the anagenetically derived species on Ulleung Island is expected to show high levels of genetic variation similar to (or even higher than) those of the continental progenitor species. Furthermore, given low habitat heterogeneity, no geographic partitioning of the genetic variation owing to gene flow among subpopulations is expected over the landscape . Although these expected genetic consequences of anagenetic speciation on the Ulleung Islands are supported by the results of the few aforementioned studies, we know very little about whether similar patterns can be found in other diverse endemic angiosperm lineages, hindering a synthesis of any emerging patterns on Ulleung Island.
Of nearly 40 endemic species on Ulleung Island, the family Rosaceae (rosid I, Fabidae) includes the greatest number of endemic species (i.e., Cotoneaster wilsonii, Potentilla dickinsii var. glabrata, Prunus takesimensis, Rubus takesimensis, and Spiraea insularis). However, the lack of detailed population genetic studies based on broad samplings of both continental progenitors and anagenetically derived insular species hampers any generalizations regarding the genetic consequences of endemic species-rich angiosperm lineage. In this study, we determined the genetic consequences of anagenetic speciation using a species pair: an insular endemic Rubus takesimensis Nakai on Ulleung Island and a continental progenitor R. crataegifolius Bunge, based on broad geographical sampling and chloroplast DNA sequences. The genus Rubus, with worldwide distribution, is a very complex and challenging group owing to agamospermy, hybridization, and polyploidy [28,29], and regional variations in morphology [30–32]. Most species of Rubus in Korea belong to the subgenus Idaeobatus, and the phylogenetic relationships among them were preliminarily assessed [33,34]. Of approximately 21 species of Rubus recognized in Korea, Rubus takesimensis is the only endemic species on Ulleung Island, which is morphologically most closely related to R. crataegifolius found in northeastern Asia (Korea, Japan, China, and the Russian Far East) [34,35]. Rubus takesimensis, the insular derivative species, differs from the continental progenitor R. crataegifolius in its lack of prominent prickles on the stems (i.e., an example of a loss of a defense mechanism) and a general trend towards increasing leaf and flower size (i.e., an example of insular gigantism). Molecular phylogenetic studies established the progenitor-derivative species relationship, allowing us to use them as a model system to study anagenetic speciation on Ulleung Island [34,36,JY Yang unpublished data]. Our preliminary study based on limited insular derivative populations provided some possibilities of multiple introductions of R. takesimensis on Ulleung Island . This study, however, was based on a very limited sampling of R. takesimensis (i.e., 37 individuals) and one variable noncoding region of cpDNA (i.e., trnL-trnF; a total of 391 aligned sites). The initial sampling strategy of the previous study was to broadly sample the continental progenitor species, R. crataegifolius, to determine the source populations of R. takesimensis. Therefore, in this study, we made special efforts to broadly sample R. takesimensis on Ulleung Island (a total of 106 individuals from 12 populations) and also to sample other highly variable noncoding sequences of cpDNA (a total of 4,212 aligned sites). Specifically, we attempted to (1) determine the geographic origin of R. takesimensis, (2) compare levels of genetic diversity between continental progenitor and insular derivative species and the degree of genetic differentiation between them, and finally (3) examine population genetic structures within insular derivative species. The results are expected to allow us to synthesize any emerging patterns in the genetic consequences of anagenetic speciation on Ulleung Island.
Materials and methods
To determine the geographic origin of R. takesimensis and reveal the genetic structure among populations on Ulleung Island, we sampled the continental progenitor species R. crataegifolius throughout East Asia (i.e., the Korean peninsula and Jeju Island, China, the Japanese archipelago, and the Russian Far East) and the insular derivative species R. takesimensis throughout its entire range on Ulleung Island. There were no permit requirements to collect both species in east Asia. For R. crataegifolius, we sampled a total of 23 populations (209 individuals), including 14 populations (135 individuals) from Korea, 6 populations (51 individuals) from Japan, 1 population (8 individuals) from China, and 2 populations (15 individuals) from the Russian Far East (Table 1). The number of individuals per population ranged from 3 (ND, Chungcheongnam-do, Korea) to 13 (DH, Gangwon-do, Korea). In the case of R. takesimensis, we sampled a total of 12 populations (106 individuals) throughout Ulleung Island. The number of individuals per population ranged from 7 (THR) to 13 (JR). All individuals were sampled at least 10 m apart after considering the possibility of clonal reproduction in two species. Thus, we sampled a total of 35 populations (315 individuals) for our study species, and R. trifidus (two individuals) was also sampled as an outgroup for phylogenetic analysis (Table 1). All leaf samples were dried with silica gel until DNA extraction, and voucher specimens for representative individuals of each population were deposited in the SKK (the Ha Eun Herbarium, Sungkyunkwan University) and KNU (the Kyungpook National University Herbarium).
DNA isolation, PCR amplification, and sequencing
Total genomic DNA was extracted from dried leaves using a DNeasy plant mini kit (Qiagen, Valencia, California). After screening several noncoding regions of chloroplast DNA, we selected five highly variable intergenic spacers, i.e., atpH/atpI, psaI/accD, psbD/trnT, trnG/trnS, and trnL/trnF [37,38]. PCR amplification was performed in a final volume of 50 μL, with each primer at a concentration of 1.0 pmol, along with 0.2 mM dNTPs, 2.5 mM MgCl2, and 1.25 U of IncloneTM Taq DNA polymerase (InClone Biotech Co., Seoul, Korea). PCR conditions were an initial denaturation of 95°C for 2 min, 35 cycles of 95°C denaturation for 20 s, 56°C annealing for 40 s, 72°C extension for 1 min, and a final extension at 72°C for 5 min. PCR products were purified using the InCloneTM gel & PCR purification kit (InClone Biotech Co., Seoul, Korea). Sequencing reactions were performed for the purified PCR products using BigDye Terminator V3.1 Cycle Sequencing Kit (Applied Biosystems, Foster City, California) at the GenoTech Corp. (Daejeon, Korea).
Phylogenetic analysis, haplotype network, genetic diversity and differentiation, and demographic change
All sequences were assembled and edited using Sequencher 4.2.2. (Gene Codes, Ann Arbor, MI, USA). Sequences were aligned using Clustal W version 2 , and minor adjustments were made manually using MacClade version 4.08 . Based on the previous phylogenetic study of Rubus in Korea and phylogeographic analysis, R. trifidus was selected as an outgroup [33,34,36]. We performed maximum parsimony (MP) and maximum likelihood (ML) analysis using PAUP* 4.0b10  and IQ-TREE version 1.5.5 , respectively, to determine the source population of R. takesimensis on Ulleung Island. All gaps were treated as missing data or coded as simple binary characters . All characters were treated as unordered, and transformation of all characters was weighted equally. For ML analysis, the best-fit model “K3Pu+F+R2” was used to find the best scored ML tree using IQ-TREE program. For MP analysis, search options included heuristic searching with multistate taxa interpreted as uncertain, a simple stepwise addition of taxa, tree-bi-section reconnection (TBR) branch swapping, MulTrees options on, and zero branches collapsed. The bootstrap support (BS) value for each node was also calculated based on bootstrap analysis with 1,000 replicates using the same heuristic search options .
For DNA polymorphism, the number of haplotype, nucleotide diversity (π), gene diversity, and θπ were estimated by Arlequin version 3.5 . To determine haplotype relationships between R. crataegifolius and R. takesimensis, we used the haplotype network program TCS version 1.21 . Gaps were treated as missing data, and the connection limit excluding homoplastic changes was set to 95% in accordance with Hart and Sunday . An analysis of molecular variance (AMOVA) implemented by Arlequin version 3.5  was performed to estimate generic variation among and within populations. The AMOVA analyses were performed using three categories: (1) all samples from the two species; (2) the continental progenitor R. crataegifolius, using three groups from Korea, China and Russia, and Japan; and (3) the insular derivative, R. takesimensis. The significance of the variance components was calculated using 1,023 permutations. For demographic history, we used mismatch distribution analysis using Arlequin 3.5  and DnaSP V. 5 . The goodness-of-fit under a sudden expansion model was tested with the sum of squared deviations (SSD; ) and Harpending’s raggedness index (RI) . A neutrality test, Tajima’s D , and Fu and Li’s test [52,53] were also performed to determine any evidence of recent range expansion. A parametric bootstrap approach with 1,000 replicates was used to test the observed mismatch distribution fit to the sudden expansion model. Lastly, we conducted PERMUT analysis  to test for the existence of geographical structures within populations of R. takesimensis on Ulleung Island and compared two coefficients for gene differentiation, GST and NST, with 1,000 permutations. Nine populations with less than 3 individuals (i.e., ND, YW, HG, GWH, JJ-DANG, JJ-DRS, Shiko, BHD, and DMD) were excluded for the mismatch distribution, PERMUT analysis, and in the distribution of chloroplast haplotypes among populations.
Phylogenetic relationship between R. crataegifolius and R. takesimensis
The combined aligned sites of the five intergenic spacers were 4,212 base pairs (bp). The accD/psaI intergenic region was 706 bp long, including nine single nucleotide polymorphisms (SNPs, 1.2%) and four indels, while the atpH/atpI region was 576 bp (14 SNPs, 2.4%; one indel). The psbD/trnT region was 1,430 bp (30 SNPs, 2.0%; three indels), trnG/trnS was 560 bp (six SNPs, 1.07%; no indel), and trnL/trnF was 940 bp (five SNPs, 0.5%; three indels). Of 4,212 total aligned sites, 4,017 (95.37%) were constant, 57 (1.35%) were variable but parsimony-uninformative, and 138 (3.28%) were parsimony-informative characters. The MP analysis found 605 equally most parsimonious trees with a tree length (TL) of 296, a consistency index (CI) of 0.67, and a retention index (RI) of 0.96. Similar tree topology was found when the gaps were treated as either missing data or coded as binary characters, and thus, the strict consensus tree based on the treatment of gaps as missing data is shown (Fig 1). The ML tree (not shown) also showed essentially the same tree topology as the MP tree. Based on extensive sampling of R. takesimensis throughout the entire range on the island, our results were similar to what we showed earlier but provided new insights into the origin and evolution of R. takesimensis . Rubus takesimensis emerged as non-monophyletic, with 12 individuals from Chusan (CS), Seokpo (SP), and Dodong (DDHN, DDR16) populations closely related to R. crataegifolius individuals. The group containing the most R. takesimensis individuals throughout the entire range of Ulleung Island shared the most recent common ancestor with populations of R. crataegifolius sampled from the northern part of the Korean peninsula, China, and the Russian Far East (hereafter called group A). The other group (hereafter called group B) included R. crataegifolius, primarily sampled from the southern part of the Korean peninsula, Jeju Island, and Japan. Within this group, southern populations sampled from Jeju Island (JJ-GMO, SPpy, and YAO), Korea, were closely related to the ones sampled from Hiroshima (Hiro) and Yamagata (Yama) prefectures in Japan. In group A, some individuals sampled from China (YR1), were closely related to the ones from Russia (Rus1), while other individuals were closely related to populations sampled from Korea (SC, YS, and SA). Some individuals from Russia (Rus2) were more closely related to the ones from the mountains of Seolak (SA) and Hambaek (HB) in the Gwangwoan-do province of Korea. It was noticeable that some individuals of R. takesimensis were more closely related to R. crataegifolius than to other conspecific populations on the island; three of eight individuals sampled from Dodong (DDHN), two of ten individuals sampled from Jeodong (DDR16), and two of eight individuals from Seokpo (SP). Furthermore, five of eight individuals of R. takesimensis from the Chusan population (CS) were closely related to group A of R. crataegifolius; two individuals, 217 and 218, were closely related to populations sampled from Gwangwon-do province (DH) and Gwangwon-do/Gyeongsangbuk-do province (HB, SA, and BS), respectively.
Genealogical relationships among haplotypes
The TCS haplotype network (Fig 2) showed essentially the same relationships as shown in the phylogenetic tree (Fig 1). We found a total of 129 haplotypes between the continental progenitor R. crataegifolius and the insular-derived R. takesimensis. Of 129 haplotypes, 81 and 48 haplotypes were found exclusively in R. crataegifolius and R. takesimensis, respectively. Most haplotypes were found at low frequency, and no haplotype was shared between the two species pair. For the insular derivative R. takesimensis, most haplotypes, except for individuals sampled from Chusan (H3, H4, and H5), were connected by one mutational step from major haplotypes, suggesting recent origins from the ancestral haplotype. Haplotypes H1 (18 individuals), H24 (eight individuals), and H29 (ten individuals) showed high frequency in R. takesimensis on Ulleung Island.
In most cases, each haplotype was represented by one individual, but rarely by two to three individuals. Eight haplotypes, each one represented by a single individual, were found in the Chusan population (CS). Specifically, haplotypes H43 and H58 were more closely related to the haplotype found in R. crataegifolius sampled from the northern part of the Korean peninsula. Two haplotypes found in Dodong-ri (DDHN), H103 and H119, were also closely related to R. crataegifolius sampled from the southern part of the Korean peninsula. Of seven haplotypes found in Seokpo (SP), two haplotypes (H108 and H109) were closely connected to haplotypes of R. crataegifolius sampled from Jeju Island, Korea (H105 and H114). These cases represent several rare haplotypes on the island more closely related to those of R. crataegifolius, indicating their independent origins from different continental source areas.
In the case of the continental progenitor R. crataegifolius shown in the phylogenetic tree, we found two divergent clusters: group A included major haplotypes sampled from the northern part of the Korean peninsula, while group B included major haplotypes from the southern part of the Korea peninsula, Jeju Island, and Japan. Higher haplotype diversity was found in populations sampled from the Korean peninsula compared with the ones in Japan. The population sampled from Baekryeong Island (BR) in the Yellow Sea contained two divergent haplotype groups, one with H40 and H41, and the other with H97 and H98, and two somewhat divergent haplotypes (H40 and H41) were connected to the most common haplotype (H1) of R. takesimensis. The most common haplotype within this group was H42 (12 individuals), which was found in Gwangwon-do province (HB and DH), Gyeonggi-do province (CB), and Gyeongsangbuk-do province (BS) in Korea. Unlike haplotypes found in R. takesimensis, haplotypes found in this group A were connected by one to seven missing or inferred haplotypes, suggesting the divergent nature of the haplotype group. Of the five haplotypes found in one population from China (YR1), three (H60, H61, and H62) were connected to the haplotype H59 sampled from Russia (Rus1), while the remaining two (H46 and H84) were related to the southern part of the Korean peninsula (H47 from SC and H82 from BS). Most haplotypes found in Russia, except for two haplotypes, were connected to the ones from the northern part of the Korean peninsula (SA and BS). Two haplotypes (H43 and H58) of R. takesimensis sampled from Chusan (CS) on Ulleung Island were connected to the ones from Gangwon-do province (H44 from DH and H57 from HB, respectively).
The group B included haplotypes found primarily in the southern part of the Korean peninsula (Gyeongsangnam-do and Jeollanam-do province), Jeju Island, and Japan. The most common haplotype was H105 (23 individuals), which occurred in the southern part of Korea (YS and BU), Jeju Island, and Japan (Hiro and Miya). The haplotypes found in R. takesimensis were connected to these southern populations via several missing haplotypes. The haplotype found in Jeju Island (H105) showed regional sharing across the Halla Mountains between the northern and southern parts of the island. Two haplotypes (H108 and H109) of R. takesimensis found from Seokpo (SP) on Ulleung Island were connected to H105 and H114 found on the northern part of Jeju Island.
Diversity and distribution of chloroplast haplotypes among populations
The number of haplotypes found in each population of R. takesimensis on Ulleung Island ranged from four (DDO, NSJ, and SDR) to nine (JR) (Table 1 and Fig 3). Some populations sampled from the northern part of the island, such as the Nari Basin (AB; seven haplotypes), Chusan (CS; eight haplotypes), and Seokpo (SP; seven haplotypes) harbored diverse haplotypes. In addition, some haplotypes (e.g., H1, H24, and H29) occurred widely in different populations, whereas others (e.g., H2, H3, H4, H5, among others) were narrowly restricted to specific populations. A total of 81 haplotypes were found in R. crataegifolius, and two divergent groups were observed (Fig 4). Group A included populations primarily sampled from the mainland part of the Korean peninsula, Russia, and China. For the 14 populations sampled from Korea, the number of haplotypes per population ranged from three (JJ-BSG) to eight (HB, SA, and SC). Of the 14 populations, six populations (BR, CB, BS, SA, HB, and DH) contained primarily group A haplotypes, while the remaining three from the southern peninsula (BU, SC, and YS) and four populations from Jeju Island contained the group B haplotypes. For the three populations sampled from Russia and China, the number of haplotypes ranged from four (Rus2) to five (Rus1 and YR1), and they are closely related to the ones from the northern part of the Korean peninsula.
Contrary to group A, with its diverse haplotypes, group B included populations primarily from the southern part of the Korean peninsula and Japan. For the six populations sampled from Japan, all but two populations (Oita and TK) contained two haplotypes, showing highly reduced haplotype diversity. Two populations, Oita and TK, contained H123 (shared with the nearby Naga population) and H104 (specific to the TK population), respectively.
Genetic diversity and structure between R. crataegifolius and R. takesimensis
With regard to the partitioning of genetic variations between the two taxa (the continental progenitor R. crataegifolius and the insular derivative R. takesimensis), the majority of variations (73%) existed within populations, while approximately 27% of the variations existed among populations (Table 2). Within a continental progenitor species R. crataegifolius, slightly more than half of the total variations (i.e., 57%) occurred among populations, while the remaining 17% and 26.5% of the variations existed among countries and within populations, respectively. In the case of the insular derivative species R. takesimensis, total variations were partitioned more-or-less equally among populations (57%) and within populations (43%).
In the case of the genetic diversity between continental progenitors and insular derivative species, we found significantly higher diversity statistics in R. crataegifolius than in R. takesimensis: the number of haplotypes (81 vs. 48), nucleotide diversity (π) (0.015951 vs. 0.007045), and θπ (65.0013 vs. 27.67857) (Table 3). However, we cautiously interpreted these diversity statistics given the unequal sample size between two species: 209 individuals in R. crataegifolius vs. 106 individuals in R. takesimensis.
Demographic change and geographical structure
Mismatch distributions showed that R. crataegifolius had multimodal distributions with an SSD of 0.007 (P<0.05) and RI of 0.0074 (P<0.05) (Fig 5 and Table 4). Also, populations sampled from R. crataegifolius on the Korean peninsula only showed multimodal distributions. Mismatch distributions of R. takesimensis were also multimodal, with an SSD of 0.0149 (not significant) and an RI of 0.0355 (P<0.05). Therefore, it is likely that the population of R. takesimensis did not undergo recent population expansion or a bottleneck after introduction to Ulleung Island [53,54]. Demographic analysis showed negative values for Fu’s Fs and Tajima’s D test (Table 4), with significant values for R. crataegifolius only (Fs = -2.75529, P <0.05). The Harpending’s raggedness index (RI) derived from the mismatch distribution was significant (Table 4). In R. takesimensis on Ulleung Island, the total genetic diversity (hT = 0.885), based on cp DNA data across all populations, was much higher than the average within the population diversity (hS = 0.701). The GST value (0.208) of R. takesimensis on Ulleung Island was higher than its corresponding NST (0.139), suggesting no significant genetic structure in R. takesimensis on Ulleung Island.
Multiple origins of R. takesimensis and patterns of anagenetic speciation on Ulleung Island
This study represents the first attempt to utilize chloroplast DNA based on extensive sampling (315 individuals) and sequence data (aligned sequence of 4,212 bp) to characterize the genetic consequences of anagenetic speciation in one endemic species-rich group, Rosaceae, on Ulleung Island, Korea. Previous studies were limited to the Sapindaceae (i.e., Acer takesimensis and Acer okamotoanum) and the Apiaceae (Dystaenia takesimana), primarily utilizing nuclear markers (i.e., SSRs, AFLPs, and nrDNA ITS) with complimentary use of a chloroplast marker (trnL-F) [23–25,55]. It may not be possible to make direct comparisons between our study, which was based on maternally inherited cpDNA and thus focused on seed-mediated gene flow, and previous studies, which primarily utilized biparentally inherited nuclear DNA. However, the general trend and pattern found in this study could be compared and evaluated as a proxy of nuclear genetic diversity, given the unusually high level of chloroplast diversity in Rubus species.
First of all, we found that endemic species of R. takesimensis on Ulleung Island are not the result of a single colonization event. Rather, several geographically diverse continental source populations were likely involved in the origin of R. takesimensis. As shown in the phylogenetic tree (Fig 1) and the haplotype network (Fig 2), a major group of R. takesimensis accessions were closely related to the group containing R. crataegifolius accessions primarily sampled from the northern part of the Korean peninsula, while rare accessions from Dodong (DDR16 and DDHN) and Seokpo (SP) most likely originated from the southern part of the Korean peninsula, Jeju Island, and Japan. Furthermore, even within a population, especially Chusan (CS) in the north-facing part of the island, we found evidence for a few accessions originating separately from the Korean peninsula. These populations are located either nearby a harbor used for commercial fishing and passenger ferries or near typical bird nesting sites, thus increasing the likelihood of anthropogenic or natural independent colonization events from different source populations. It seems less likely though that human-mediated colonization is responsible for the origin of some accessions in those populations because the archaeological evidence indicates that the island has been inhabited since the 1st millennium BC, with the first confirmed historical reference for the year of 512 CE [56–58]. Considering the geographical proximity of Ulleung Island to possible source areas, i.e., 137 km from the Korean peninsula and approximately 240 km from Oki Island, Japan, multiple introductions may not be exceptional . Pfosser et al. , based on one noncoding cpDNA and AFLPs, initially raised the possibility of multiple origins for A. okamotoanum on Ulleung Island, but later suggested a single origin from the continental progenitor A. mono based on microsatellite or SSR markers . This discrepancy based on different types of molecular markers may be caused by incomplete lineage sorting and shared ancestral polymorphism and/or gene introgression between closely related species, as well as limited genetic information revealed by low evolutionary rate of the plastid genome [59–62]. It may be necessary to use highly variable SSRs or genome-wide SNPs to confirm our current results. Our earlier study  also raised the possibility of colonization of the most common haplotype in R. takesimensis back to the southern part of the Korean peninsula. However, current studies (JY Yang, unpublished data) involving extensive sampling from southern areas did not find any haplotype sharing between continental and insular species, suggesting Ulleung Island as a primary sink area for plant colonization despite its geographical proximity to the continent. Nevertheless, this study raises the possibility of multiple introductions for Ulleung Island endemics from nearby source areas, and additional investigations on diverse phylogenetic lineages with different seed dispersal modes will be helpful to determine the evolutionary history of each group.
Based on the few available studies of anagenetically derived species on Ulleung Island, some emerging patterns of the genetic consequences include a clear genetic distinction between continental progenitor and island-derived species and high levels of genetic variation within the island population with no geographic genetic partitioning [23,24,63]. Furthermore, a recent comprehensive review confirmed the same patterns, i.e., clear genetic differences between continental progenitors and insular derivative species, different levels of genetic variation depending on species, and no geographical partitioning of the genetic variation . Based on the chloroplast DNA noncoding sequences, we found a significant reduction in haplotype number and nucleotide diversity in R. takesimensis, the insular derivative, compared with its continental progenitor R. crataegifolius (Table 3). This finding, however, could be due to uneven sampling between the two species. Rubus takesimensis harbors unique and diverse chlorotypes despite its occurrence in a young (<1.8 million years old) and small isolated island (a total area of 73 km2). Although we used five highly variable noncoding regions of cpDNA in our study, this level of haplotype diversity (i.e., 48 chlorotypes) contrasts with that of other endemic species, such as Fagus multinervis (one haplotype based on the psbA-trnH noncoding region)  and Sedum takesimense (two haplotypes based on the trnL-F noncoding region) . Therefore, we suggest that R. takesimensis is genetically unique and contains comparable genetic diversity relative to its continental progenitor species.
Factors affecting patterns of genetic variation in R. takesimensis
As suggested by Takayama et al. , several factors, such as population size, breeding system, and species age can contribute to the genetic consequences of anagenetically derived species on Ulleung Island. The age of R. takesimensis on Ulleung Island was estimated to be 1.7 ± 3.02 million years old, suggesting its origin soon after the formation of the island (JY Yang, unpublished data). In addition, species of Rubus are known for bird-mediated seed dispersal, employ diverse modes of reproduction (sexual and asexual), and play as early forest succession frontiers. When a single individual of Rubus is transplanted into a bare and open area, it expands quickly as a population with asexual creeping stems and roots in less than five years and undergoes sexual reproduction with the production of flowers and fruits (JY Yang, personal observation). In addition, Rubus species can utilize both asexual (obligate apomixes to facultative apomixes) and sexual (obligate sexuality) reproductive modes as the population expands . As a consequence, these strategies allow fast formation and stabilization of a population via asexual reproduction, and maintain high genetic diversity from sexual reproduction [67–72]. High levels of genetic variation in vegetatively propagated species have been reported in the North American quaking aspen, Populus tremuloidess . Sarhanova et al.  pointed out that anthropogenic disturbances, such as road construction, deforestation, and landscape changes, lead to the fast evolution of bramble species in Europe. Rubus takesimensis can vigorously expand its habitat to open areas as frontier plants, accumulate mutations in its populations, and receive a continuous supply of new genetic materials from continental progenitors through migrating birds. A combination of these factors, i.e., large population size, dual mode of reproduction, bird-mediated seed dispersal, and the long presence on the island likely resulted in high genetic diversity and no geographical genetic structuring in R. takesimensis.
As for bird-mediated seed dispersal, variations in fruit and seed size have potential effects on seed dispersal. For example, it has been demonstrated that frugivory of many passerine birds is limited by their gape width [75,76]. Birds also prefer fruits with low seed to fruit ratios [75,77]. Moreover, large seeds are often regurgitated after a relatively short time, whereas small seeds pass through the gut and are defecated after longer periods . Therefore, small fruits and small seeds are more likely to be dispersed over long distances . Brambles are typical bird-dispersed plants with a wide spectrum of dispersers [80,81]. European blackberries are mainly dispersed by frugivorous passerines [76,80]. Also, Jordano  and Rejmanek  reported that migrant birds are the main seed dispersers of Spanish populations of R. ulmifolius and a southern Spanish population of blackberry. The Ulleung and Dok-do Islands are known as refuge islands for migratory birds . In the case of avian fauna on Ulleung Island, the Japanese wood pigeon (Columba janthina) is one notable bird species, and several other birds have also been reported; the black-tailed gull (Larus crassiroostris), the russet sparrow (Passer rutilans), the barn swallow (Hirundo rustica), and the oriental greenfinch (Carduelis sinica) . Other permanent residents on Ulleung Island include the brown-eared bulbul (Hypsipetes amaurotis), the great tit (Parus major), the pheasant (Phasianus colchicus), and the blue rock thrush (Monticola solitarius) . Although a total of 112 bird species in 38 families and 16 orders have been recorded , little is known regarding which species feed on R. takesimensis on Ulleung Island. Considering the habit of R. takesimensis (i.e., subshrubs), it is plausible that species belonging to the order Passeriformes could feed on aggregate druplets of R. takesimensis. The Passeriformes, suborders Passeri are small or medium-sized birds, mostly insectivorous or omnivorous, and known raspberry feeding birds include the catbird, the mocking bird, the robin, the tanager, the thrasher, the turkey, and the warbler . The Passeriformes, Carduelis sinica (oriental green finch), Passer rutilans (russet sparrow), and Zosterops japonicas (warbling white-eye) are found on Ulleung Island. Carduelis sinica is assumed to be a permanent resident on the island, while Zosterops japonicas occur in southern coastal areas of the Korean peninsula, Jeju Island, and Ulleung Island. Passer rutilans, a winter migratory bird that visits southern parts of the Korean peninsula for the summer, occurs much more broadly in eastern Asia, Sakhalin, Jeju Island, and Ulleung Island. Considering some R. takesimensis haplotypes (e.g., Chusan) are closely related to R. crataegifolius haplotypes from the northern part of the Korean peninsula, China, and Russian Far East and others (i.e., Dodong, Jeodong, and Seokpo) to the southern part of the Korean peninsula, Jeju Island, and Japan, it is conceivable that such migratory birds are responsible for introducing a diverse pool of R. crataegifolius haplotypes to Ulleung Island, eventually evolving into R. takesimensis.
Haplotype diversity in continental R. crataegifolius
For the patterns of genetic variation in the continental progenitor species R. crataegifolius, we found that populations sampled from the inland part of the Korean peninsula showed much higher haplotype diversity compared with those from Japan. In addition, haplotype composition among the populations from Japan was very simple, i.e., just one fixed haplotype, or up to two closely related haplotypes (except for H123). In Korea, R. crataegifolius frequently occurs in lowland open habitats and along mountainous trail edges all over the Korean peninsula, while it is rarely found in the mountainous areas of Japan. Considering the wide geographic natural distribution of R. crataegifolius in Japan, it was interesting to find very low genetic variation within populations and differentiation among populations compared with conspecific populations on the continent. Several factors (e.g., breeding system, demographic history, gene dispersal, among others) likely contribute to neutral genetic diversity within populations and genetic differentiation among populations. In particular, vegetative propagation and the selfing rate have profound outcomes for the genetic diversity of species [84,85]. In the genus Rubus, it has been demonstrated that sexual R. ideanus has much higher genetic variation than apomictic R. nessensis . Although the mode of reproduction in populations of R. crataegifolius is unknown, it is conceivable that populations in Korea are sexual reproducing diploid or facultative apomictic, while those in Japan are predominantly obligate apomictic. The breeding system of R. crataegifoloius and other possible factors for shaping genetic variation and differentiation requires further investigation.
Given the rarity of sharing certain haplotypes between and among populations, it was noticeable to find the haplotype H105, which occurs commonly in the southern part of the Korean peninsula (i.e., Gyeongsangnam-do, Jeollanam-do, and Jeju-do Provinces), is shared with populations (northern Miya and southern Hiro) throughout the Japanese archipelago. Three derived haplotypes from H105 by one mutational step, i.e., H104, H106, and H107, were shared among populations (H107) or unique (H104 and H106) to a specific population. Therefore, the sharing of haplotype H105 between southern Korea and Japan may suggest a recent migration of R. crataegifolius from the Korean peninsula to the Japanese archipelago via land bridges, which existed during the last glacial maximum (LGM; the period between 25,000 and 18,000 before present) and served as a dispersal corridor during that time [36,87–89]. Alternatively, it is also conceivable that the haplotype H105 of R. crataegifolius colonized Japan via long-distance dispersal most likely by birds during or even before the Quaternary climatic oscillation, and survived during the LGM. Given all other haplotypes are confined within each region (the Korean peninsula and the Japanese archipelago), the Korean/Tsushima Strait served as an effective barrier to seed-mediated gene flow [90–93].
We provided the first comprehensive results regarding the genetic consequences of the anagenetic speciation of R. takesimensis on Ulleung Island from the continental progenitor R. crataegifolius based on chloroplast DNA sequences. The insular derivate endemic species, R. takesimensis on Ulleung Island, originated multiple times from diverse geographical source populations. A distinct pool of chlorotypes was found between the progenitor and derivative species, suggesting significant genetic divergence between the two species. A highly diverse haplotype composition within the population and a lack of significant geographic structuring of R. takesimensis were found. Two major groups of chlorotypes in R. crataegifolius were found, with high diversity in the mainland continent populations than in the Japanese archipelago. Based on the results of anagenetic speciation of R. takesimensis, it seems necessary to test rigorously the monophyly of diverse insular endemic lineages on Ulleung Island based on extensive sampling strategy (i.e., covering diverse geographical source areas for continental progenitor species and anagenetically derived insular populations) and sequence data (i.e., including genome-wide SNPs and fastly evolving genes) and compare and contrast the emerging patterns in the genetic consequences of anagenetic speciation shaped by different types of molecular markers. In addition, the discovery of cryptic species on oceanic archipelagos in collaboration with others in the botanical disciplines is critical to the discovery and conservation of the subtle diversity in oceanic island flora .
1. Stuessy TF, Jakubowsky G, Gomez RS, Pfosser M, Schluter PM, Fer T, et al. Anagenetic evolution in island plants. J Biogeo. 2006;33: 1259−1265.
2. Stuessy TF. Evolution of specific and genetic diversity during ontogeny of island floras: the importance of understanding process for interpreting island biogeographic patterns. In: Ebach MC, Tangney RS, editors. Biogeography in a changing world. Boca Raton: CRC Press; 2007. pp. 117−133.
3. Stuessy TF, Takayama K, López-Sepúlveda P, Crawford DJ. Interpretation of patterns of genetic variation in endemic plant species of oceanic islands. Bot J Linn Soc. 2014;174: 276−288. doi: 10.1111/boj.12088 26074627
4. Eliasson U. Studies in Galapagos Plants. XIV. The genus Scalesia Arn. Opera Botanica. 1974;36: 1−117.
5. Schilling EE, Panero JL, Eliasson UH. Evidence from chloroplast DNA restriction site analysis on the relationships of Scalesia (Asteraceae: Heliantheae). Am J Bot. 1994;81: 248−254.
6. Wagner WL, Herbst DR, Sohmer SH. Manual of the flowering plants of Hawaii, vol. 1. 2nd ed. Honolulu: University of Hawaii and Bishop Museum Press; 1999.
7. Nielsen LR. Molecular differentiation within and among island populations of the endemic plant Scalesia affinis (Asteraceae) from the Galapagos Islands. Heredity. 2004;93: 434−442. doi: 10.1038/sj.hdy.6800520 15280895
8. Böhle UR, Hilger HH, Martin WF. Island colonization and evolution of the insular woody habit in Echium L. (Boraginaceae). Proc Natl Acad Sci USA. 1996;93: 11740−11745. doi: 10.1073/pnas.93.21.11740 8876207
9. Kim S-C, Crawford DJ, Francisco-Ortega J, Santos-Guerra A. A common origin for woody Sonchus and five related genera in the Macaronesian islands: molecular evidence for extensive radiation. Proc Natl Acad Sci USA. 1996;93: 7743−7748. doi: 10.1073/pnas.93.15.7743 8755546
10. Kim S-C, Crawford DJ, Jansen RK. Phylogenetic relationships among the genera of subtribe Sonchinae (Asteraceae): Evidence from ITS sequences. Syst Bot. 1996;221: 417−432.
11. Jorgensne TH, Olesen JM. Adaptive radiation of island plants: evidence from Aeonium (Crassulaceae) of the Canary Islands. Persepct Plant Evol. 2001;4: 29−42.
12. Mort ME, Soltis DE, Soltis PS, Francisco-Ortega J, Santos-Guerra A. Phylogenetic and evolution of the Macaronesian clade of Crassulaceae inferred from nuclear and chloroplast sequence data. Syst Bot. 2002;27: 271−288.
13. Kim S-C, McGowen MR, Lubinsky P, Barber JC, Mort ME, Santos-Guerra A. Timing and tempo of early and successive adaptive radiations in Macaronesia. PLoS ONE. 2008;3: 32139.
14. Sang T, Crawford DJ, Kim S-C, Stuessy TF. Radiation of the endemic genus Dendroseris (Asteraceae) in the Juan Fernandez Islands: evidence from sequences of the ITS regions of nuclear ribosomal DNA. Am J Bot. 1994;81: 1494−1501.
15. Sang T, Crawford DJ, Stuessy TF, Silva M. ITS sequences and the phylogeny of the genus Robinsonia (Asteraceae). Syst Bot. 1995;20: 55−64.
16. López-Sepúlveda P, Takayama K, Greimler J, Crawford DJ, Peñailillo P, Baeza M, et al. Speciation and biogeography of Erigeron (Asteraceae) in the Juan Fernández Archipelago, Chile, based on AFLPs and SSRs. Syst Bot. 2015;40: 888−899.
17. Takayama K., Lopez-sequlveda P, Greimler J, Crawford DJ, Penailillo P, Baeza M, et al. Relationships and genetic consequences of contrasting modes of speciation among endemic species of Robinsonia (Asteraceae, Senecioneae) of the Juan Fernandez Archipelago, Chile, based on AFLPS and SSRs. New Phytol. 2015;205: 415−428. doi: 10.1111/nph.13000 25209139
18. Ito M, Ono M. Allozyme diversity and the evolution of Crepidiastrum (Composiate) on the Bonin (Ogasawara) Islands. Bot Mag (Tokyo). 1990;103: 449−459.
19. Setoguchi H, Watanabe I. Intersecional gene flow between insular endemics of Ilex (Aquifoliaceae) on the Bonin Islands and the Ryukyu Islands. Am J Bot. 2000;87: 793−810. 10860910
20. López-Sepúlveda P, Takayama K, Greimler J, Peñailillo P, Crawford DJ, Baeza M, et al. Genetic variation (AFLPs and nuclear microsatellites) in two anagenetically derived endemic species of Myrceugenia (Myrtaceae) on the Juan Fernández Islands, Chile. Am J Bot. 2013;100: 722−734. doi: 10.3732/ajb.1200541 23510759
21. López-Sepúlveda P, Takayama K, Greimler J, Crawford DJ, Peñailillo P, Baeza M, et al. Progressive migration and anagenesis in Drimys confertifolia of the Juan Fernández Archipelago, Chile. J Plant Res. 2015;128: 73−90. doi: 10.1007/s10265-014-0666-7 25292282
22. López-Sepúlveda P, Takayama K, Crawford DJ, Greimler J, Peñailillo P, Baeza M, et al. Biogeography and genetic consequences of anagenetic speciation of Rhaphithamnus venustus (Verbenaceae) in the Jun Fernández Archipelago, Chile: insight from AFPL and SSR markers. Plant Spec Biol. 2017;32: 223−237.
23. Pfosser M, Jakubowsky G, Schlüter PM, Fer T, Kato H, Stuessy TF, et al. Evolution of Dystaenia takesimana (Apiaceae), endemic to Ullung Island, Korea. Plant Syst Evol. 2006;256: 159−170.
24. Takyama K, Sun B-Y, Stuessy TF. Genetic consequences of anagenetic speciation in Acer okamotoanum (Sapindaceae) on Ulleung Island, Korea. Ann Bot. 2012;109: 321−330. doi: 10.1093/aob/mcr280 22056411
25. Takyama K, Sun B-Y, Stuessy TF. Anagenetic speciation in Ullung Island, Korea: genetic diversity and structure in the island endemic species, Acer takesimense (Sapindaceae). J Plant Res. 2013;126: 323−333. doi: 10.1007/s10265-012-0529-z 23090156
26. Takayama K, Sun B-Y, Stuessy TF. Genetic consequences of anagenetic speciation in endemic angiosperms of Ullung Island. J Jpn Bot Suppl. 2016;91: 83−98.
27. Kim YK. Petrology of Ulreung volcanic island, Korea–part 1. Geology. J Japanese Assoc Mineral Petrol Econ Geol. 1985;80: 128−135.
28. Thompson MM. Chromosome numbers of Rubus species at the national clonal germplasm repository. HortScience. 1995;30: 1447−1452.
29. Alice LA, Campbell CS. Phylogeny of Rubus (Rosaceae) based on nuclear ribosomal DNA internal transcribed spacer region sequences. Am J Bot. 1999;86: 81−97. 21680348
30. Naruhashi N. Morphology of 34 Japanese Rubus species. Acta Hortic. 1980;112:177−181.
31. Suzuki W. Comparative ecology of Rubus Species (Rosaceae) II. Reproductive characteristics of three Rubus species, R. palmatus var. coptophyllus, R. microphyllus and R. crataegifolius. Pl Spec Biol. 1990;5: 263−275.
32. Weber HE. Former and modern taxonomic treatment of the apomictic Rubus complex. Folia Geobot Phytotax. 1996;31: 373−380.
33. Yang JY, Pak JH. Phylogeney of Korean Rubus (Rosaceae) based on ITS (nrDNA) and trnL/F intergenic region (cpDNA). J Plant Biol. 2006;49: 44−54.
34. Yang JY, Yoon HS, Pak JH. Phylogeny of Korean Rubus (Rosaceae) based on the second intron of the LEAFY gene. Can J Plant Sci. 2012;92: 461−472.
35. Nakai T. Notulae ad plants Japoniae et Koreae XVII. Bot Mag (Tokyo). 1918; 32: 104−105.
36. Lee W, Yang JY, Jung K-S, Pak J-H, Maki M, Kim S-C. Chloroplast DNA assessment of anagenetic speciation in Rubus takesimensis (Rosaceae) on Ulleung Island, Korea. J Plant Biol. 2017;60: 163−174.
37. Taberlet P, Gielly L, Pautou G, Bouvet J. Universal primers for amplification of three noncoding regions of chloroplast DNA. Plant Mol Biol. 1991;17: 1105−1109. doi: 10.1007/bf00037152 1932684
38. Shaw J, Lickey EB, Schilling EE, Small RL. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: The tortoise and the hare III. Am J Bot. 2007;94: 275−288. doi: 10.3732/ajb.94.3.275 21636401
39. Larkin MA, Blackshields G, Brown NP, Chenna R, McGettigan PA, McWilliam H, et al. "ClustalW and ClustalX version 2". Bioinformatics. 2007;23: 2947–2948. doi: 10.1093/bioinformatics/btm404 17846036
40. Maddison WP, Maddison DR. MacClade version 3.04 Analysis of phylogeny and character evolution. Massachusetts: Sinauer Associates; 1992.
41. Swofford DL. PAUP*: Phylogenetic analysis using parsimony (and other method). Version 4.0b10. Massachusetts: Sinauer Associates; 2002.
42. Nguyen L-T, Schmidt H-A, Haeseler Av, Minh BQ. IQ-TREE: A fast and effective stochastic algorithm for estimating maximum likelihood phylogenies. Mol Biol Evol. 2015;32: 268−274. doi: 10.1093/molbev/msu300 25371430
43. Simmon MP, Ochoterena H. Gaps as characters in sequence based phylogenetic analyses. Syst Biol. 2000;49: 369−381. 12118412
44. Felsenstein J. Confidence limits on phylogenies: an approach using the bootstrap. Evolution. 1985;39: 783−791. doi: 10.1111/j.1558-5646.1985.tb00420.x 28561359
45. Excoffier L, Lischer HEL. Arlequin suite ver 3.5: a new series of programs to perform population genetics analyses under Linux and Windows. Mol Ecol Resour. 2010;10: 564−567. doi: 10.1111/j.1755-0998.2010.02847.x 21565059
46. Clement M, Posada D, Crandall KA. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9: 1657−1659. 11050560
47. Hart MW, Sunday J. Things fall apart: biological species form unconnected parsimony networks. Biol Lett. 2007;3: 509−512. doi: 10.1098/rsbl.2007.0307 17650475
48. Librado P, Rozas J. DnaSP v5: A software for comprehensive analysis of DNA polymorphism data. Bioinformatics. 2009;25: 1451−1452. doi: 10.1093/bioinformatics/btp187 19346325
49. Schneider S, Excoffier L. Estimation of past demographic parameters from the distribution of pairwise differences when the mutation rates vary among sites: application to human mitochondrial DNA. Genetics 1999;152: 1079−1089. 10388826
50. Harpending HC. Signature of ancient population growth in a low-resolution mitochondrial DNA mismatch distribution. Hum Biol. 1994;66: 591−600. 8088750
51. Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics. 1989;123: 585−595. 2513255
52. Fu XY, Li WH. Statistical tests of neutrality of mutations. Genetics. 1993;133: 693−709. 8454210
53. Fu XY. Statistical tests of neutrality of mutations against population growth, hitchhiking and background selection. Genetics. 1997;147: 915−925. 9335623
54. Pons O, Petit R. Measuring and testing genetic differentiation with ordered versus unordered alleles. Genetics. 1996;144: 1237−1245. 8913764
55. Pfosser M, Guzy-Wróbelska J, Sun BY, Stuessy TF, Sugawara T, Fujii N. The origin of species of Acer (Sapindaceae) endemic to Ullung Island, Korea. Syst Bot. 2002;27: 351−367.
56. Kim P. (Lee BG; translator). The History of the Three Kingdoms (Samguk Sagi). Seoul: Withersbook; 2015.
57. Kim P. (Shultz EJ, Kang HHW; translator). The Silla Annals of the Samguk Sagi. Seoul: Academy of Korean Studies Press; 2012.
58. Han J-S. A study on the local government system in the middle period of Silla dynasty focused on ancient tomb in Ulleungdo. A Collection of Treatises of Korean Studies. 2014; 41: 131–167.
59. Barker NP, Howis S, Nordenstam B, Källersjö M, Eldenäs P, Griffioen C, et al. Nuclear and chloroplast DNA-based phylogenies of Chrysanthemoides Tourn. ex Medik. (Calenduleae; Asteraceae) reveal extensive incongruence and generic paraphyly, but support the recognition of infraspecific taxa in C. monilifera. S Afr J Bot. 2009;75: 560−572.
60. Manos PS, Zhou ZK, Cannon CH. Systematics of Fagaceae: phylogenetic tests of reproductive trait evolution. Int J. Plant Sci. 2001;162: 1361−1379.
61. Denk T, Grimm GW. The oaks of western Eurasia: traditional classifications and evidence from two nuclear markers. Taxon. 2010;59: 351−366.
62. Simeone MC, Piredda R, Papini A, Vessella F, Schirone B. Application of plastid and nuclear markers to DNA barcoding of Euro-Mediterranean oaks (Quercus, Fagaceae): problems, prospects and phylogenetic implications. Bot J Linn Soc.2013;172: 478−499.
63. Takayama K, López-Sepúlveda P, Greimler J, Crawford DJ, Peñailillo P, Baeza M, et al. Genetic consequences of cladogenetic vs. anagenetic speciation in endemic plants on oceanic islands. AoB Plants. 2015;7: plv102.
64. Oh S-H. Sea, wind, or bird: origin of Fagus multinervis (Fagaceae) inferred from chloroplast DNA sequences. Korean J Pl Taxon. 2015;45: 213−220.
65. Lee W, Pak JH. Intraspecific sequence variation of trnL/F intergenic region (cpDNA) in Sedum takesimense Nakai (Crassulaceae) and aspects of geographic distribution. Korean J Pl Taxon. 2010;40: 157−162.
66. Sochor M, Vasut RJ, Sharbel TF, Travnicek B. How just a few makes a lot: Speciation via reticulation and apomixes on example of European brambles (Rubus subgen. Rubus, Rosaceae). Mol Phylogenet Evol. 2015;899: 13−27.
67. Asker SE, Jerling L. Apomixis in Plants. London: CRC press; 1992.
68. Birky CW Jr. Heterozygosity, heteromorphy, and phylogenetic trees in asexual eukaryotes. Genetics. 1996;144: 427−237. 8878706
69. Horandle E, Paun O. Patterns and sources of genetic diversity in apomictic plants: implications for evolutionary potentials. In: Horandl E, Grossniklaus U, van Dijk PJ, Sharbel TF, editors. Apomixis: evolution, mechanisms and perspectives. Rugell: A.R.G. Gantner Verlag; 2007. pp. 169−174.
70. Majesky L, Vasut RJ, Kitner M, Schranz ME, Koch M, Kiefer C, et al. The pattern of genetic variability in apomictic clones of Taraxacum officinale indicates the alternation of asexual and sexual histories of apomicts. PLoS ONE. 2012;7: e41868. doi: 10.1371/journal.pone.0041868 22870257
71. Majesky L, Vasut RJ, Kitner M. Genotypic diversity of apomictic microspecies of the Taraxacum scanicum group (Taraxacum sect. Erythrosperma). Plant Syst Evol. 2015;301: 2105−2124.
72. Sochor M, Travnicek B. Melting pot of biodiversity: first insights into the evolutionary patterns of the Colchic bramble flora (Rubus subgenus Rubus, Rosaceae). Bot J Linn Soc. 2016;181: 610−620.
74. Sarhanova P, Sharbel TF, Sochor M, Vasut RJ, Dancak M, Travnicek B. Hybridization drives evolution of apomicts in Rubus subgenus Rubus: evidence from microsatellite markers. Ann Bot. 2017;120: 317−328. doi: 10.1093/aob/mcx033 28402390
75. Herrera CM. A study of avian frugivores, bird-dispersed plants, and their interaction in Mediterranean scrublands. Ecol Monogr. 1984;54: 1–23.
76. Snow B, Snow D. Birds and berries: A study of an ecological interaction. London: T&AD Poyser Ltd; 1988. pp. 60−63.
77. Murray KG, Winnett-Murray K, Cromie EA, Minor M, Meyers E. The influence of seed packaging and fruit color on feeding preferences of American robins. Vegetatio. 1993;107/108: 217–226.
78. Levey DJ. Methods of seed processing by birds and seed deposition patterns. In: Estrada A, Fleming TH, editors. Frugivores and seed dispersal, Dordrecht: Junk; 1986. pp. 147–158.
79. Kollmann J, Steinger T, Ray BA. Evidence of sexuality in European Rubus (Rosaceae) species based on AFLP and allozyme analysis. Am J Bot. 2000;87: 1592−1598. 11080109
80. Jordano P. Migrant birds are the main seed dispersers of blackberries in southern Spain. Oikos. 1982;38: 183–193.
81. Rejmanek M. Invasion of Rubus praecox (Rosaceae) is promoted by the native tree Aristotelia chilensis (Elaeocarpaceae) due to seed dispersal facilitation. Gayana Bot. 2015;72: 27–33.
82. Cha JS, Kim KA, Kim HJ, Park HC. The Avifauna in Ulleungdo. J Korean Nat. 2010;3: 151−157.
83. Yu J-P, Jin S-D, Kim W-B, Kang J-H, Kim I-K, Kang T-H, et al. Characteristics of birds community in Ulleung Island, Korea. J Asia Pac Biodivers. 2013;6: 175−187.
84. Hamrick JL, Godt MJW. Allozyme diversity in plant species. In: Brown AHD, Clegg MT, Kahler AL, Weir BS, editors. Plant population genetics, breeding and genetic resources. Massachusetts: Sinauer Associates; 1989. pp. 43−63.
85. Hamrick JL, Godt MJW. Effects of life history traits on genetic diversity in plant species. Philos Trans R Soc. 1996;351: 1291−1298.
86. Antonius K, Nybom H. DNA Fingerprinting reveals significant amounts of genetic variation in a wild raspberry Rubus idaeus population. Mol Ecol. 1994;3: 177−180.
87. Qi XS, Chen C, Comes HP, Sakaguchi S, Liu YH, Tanaka N, et et al. Molecular data and ecological niche modelling reveal a highly dynamic evolutionary history of the East Asian Tertiary relict Cercidiphylluum (Cercidiphyllaceae). New Phytol. 2012;196: 617−630. doi: 10.1111/j.1469-8137.2012.04242.x 22845876
88. Sakaguchi S, Qiu YX, Liu YH, Qi XS, Kim SH, Han J, et al. Climate oscillation during the Quaternary associated with landscape heterogeneity prompted allopatric lineage divergence of a temperate tree Kalopanax septemlobus (Araliaceae) in East Asia. Mol Ecol. 2012;21: 3823−3838. doi: 10.1111/j.1365-294X.2012.05652.x 22646502
89. Lee JH, Lee DH, Choi BH. Phylogeography and genetic diversity of East Asian Neolisttea sericea (Lauraceae) based on variations in chloroplast DNA sequences. J Plant Res. 2013;126: 193−202. doi: 10.1007/s10265-012-0519-1 22990429
90. Li EX, Yi S, Qiu YX, Guo JT, Comes HP, Fu CX. Phylogeography of two East Asian species in Croomia (Stemonaceae) inferred from chloroplast DNA and ISSR fingerprinting variation. Mol Phylogenet Evol. 2008;49: 702−714. doi: 10.1016/j.ympev.2008.09.012 18849001
91. Qiu YX, Sun Y, Zhang XP, Lee J, Fu CX, Comes HP. Molecular phylogeography of East Asian Kirengeshoma (Hydrangeaceae) in relation to Quaternary climate change and landbridge configurations. New Phytol. 2009;183: 480−495. doi: 10.1111/j.1469-8137.2009.02876.x 19496955
92. Qi XS, Yuan N, Comes HP, Sakaguchi S, Qiu YX. A strong ‘filter’ effect of the East China Sea land bridge for East Asia’s temperate plant species: inferences from molecular phylogeography and ecological niche modeling of Platycrater arguta (Hydrangeaceae). BMC Evol Biol. 2014;14: 41. doi: 10.1186/1471-2148-14-41 24593236
93. Wang JF, Gong X, Chiang YC, Kuroda C. Phylogenetic patterns and disjunct distribution in Ligularia hodgsonii Hook. (Asteraceae). J Biogeogr. 2013;40: 1741−1754.
94. Crawford DJ, Stuessy TF. Cryptic variation, molecular data, and the challenge of conserving plant diversity in oceanic archipelagos: the critical role of plant systematics. Korean J Pl Taxon. 2016;46: 129−148.