Intraspecific Phylogenetic Relationships of Caryopteris incana in the Tsushima Islands, Japan, Using DNA Sequence Analysis

American Journal of Plant Sciences
Vol.06 No.14(2015), Article ID:59832,12 pages
10.4236/ajps.2015.614239

Intraspecific Phylogenetic Relationships of Caryopteris incana in the Tsushima Islands, Japan, Using DNA Sequence Analysis

Masaya Ando1, Hitoshi Watanabe2*, Kiyoshi Matsubara2, Akito Taniguchi3

1Graduate School of Horticulture, Chiba University, Matsudo, Japan

2Center for Environment, Health and Field Sciences, Chiba University, Kashiwa, Japan

3Faculty of Horticulture, Chiba University, Matsudo, Japan

Email: *hwatanabe@faculty.chiba-u.jp

Copyright © 2015 by authors and Scientific Research Publishing Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

Received 24 August 2015; accepted 20 September 2015; published 23 September 2015

ABSTRACT

Caryopteris incana is a perennial shrub distributed in the temperate zone of the East Asia. It is found in West Kyushu in Japan, where it is designated as an endangered species. Tsushima, Nagasaki, which experienced repeated connection and fragmentation between the Korean Peninsula and Japan, is an island on the route along which C. incana moved to Japan from continental Asia. We conducted field work and confirmed the genetic structure of populations using DNA sequence analysis to construct a detailed distribution map and clarify the intraspecific phylogenetic relationships of C. incana in Tsushima Island. We confirmed 72 populations in Tsushima. Using the leaves of individuals cultivated from seeds collected from each natural population, we analyzed the chloroplast and nuclear DNA sequence variations. Among the populations, sequence variations were confirmed in six regions of chloroplast DNA, and six haplotypes, including base substitutions, were distinguished. Two haplotypes were mainly divided at the border of the northern part of the southern island in Tsushima. One population in the northwestern part of the north island showed a haplotype derived from the southern part. This finding revealed that the distribution of C. incana had been artificially influenced. Several haplotypes were confirmed by sequence variations in the northern populations, but only one haplotype in the southern populations, suggesting that C. incana on the north island had separated early from the south island in Tsushima.

Keywords:

Caryopteris incana, Intraspecific Differentiation, Sequence Variation, Chloroplast DNA, ITS

1. Introduction

The distribution and variation of the existing organisms are the results of the accumulated effects of past population dynamics and evolutionary history. By obtaining effective characteristics via adaption to the local environment after expansion to a new location, the organisms will undergo speciation or intraspecific differentiation [1] . Biogeography is the study of the history of the organisms by considering the variations in their characteristics and environmental changes. The field focuses, in particular, on genetic variation in phylogeography, and has been used to examine various species worldwide in order to clarify the evolutionary history of organisms that have survived the glacial and interglacial periods over the past tens of thousands to millions of years [2] -[5] . The islands of Japan, located in East Asia, are a hotspot for immigration to, and emigration from, continental Asia, and, similar to the Mediterranean region, there has been repeated connection and fragmentation of lands over glacial and interglacial periods [6] -[8] .

The Japanese Islands extend north and south; have wide climatic variation, from a southern subtropical zone to a north subarctic zone; and have maintained various endemic species in refugia where they avoided extinction during glacial periods. Further, there have been repeated invasions of organisms from the continent because of land connections with the continent during glacial periods [9] . After the formation of each land bridge, many areas were geographically isolated by the subsequent rising sea level; for instance, many endemic species inhabit the Yakushima Island, which is famous for the Yaku cedar. Among these regions, the islands of Tsushima, Nagasaki, have a mixed flora of species originating from the continent and from Japan as a result of the immigration and emigration of organisms with mainland Asia, because the Korean Peninsula and Tsushima were repeatedly connected and fragmentation by land [10] . Tsushima, Nagasaki, a continental island located northwest of Kyushu, extends 18 km from east to west and 82 km from north to south and has an area of approximately 700 km2 (Figure 1). It is approximately 50 km from the northwest of Tsushima to the Korean Peninsula, and approximately 132 km from the southeast of Tsushima to the Kyushu mainland. It has an oceanic climate with only small changes in temperature because of the Tsushima Current, but it is cool in the winter because of the northwest monsoon. Most of the islands consist of sedimentary rock rich in mud from the Tertiary, called the Taishu group. The islands are divided into the north and south, and a rias coast develops in the shoreline around the central region. The flora of Tsushima consists of the template plants from southern Japan, continental plant, and subtropical and tropical plants expanding their distribution to the north [9] [11] . Some continental plant species have expanded their distribution to Kyushu and West Japan through Tsushima, and the intraspecific and interspecific genetic differences between the Japanese Islands, the Korean Peninsula, and Tsushima, which is a halfway point between the two, are important phylogeographical indices when considered along with the geological history and paleontology of these areas [10] . Caryopteris incana (Thunb.) Miq. is one such species of continental plant.

Caryopteris is a genus of shrubs or subshrubs of the family Lamiaceae that is distributed in China, Mongolia, Tibet, Taiwan, Korea, and Japan. The terminal corymboid cymes have pedicellate ebracteate flowers, and

Figure 1. Map showing the location of the study area. (a) Geographical location of Tsushima, Japan. (b) Outline of Tsushima.

the corollas are bluish purple or pale green to yellowish. The leaves are strongly aromatic. The familial assignment and infrageneric classification of Caryopteris were determined based on their floral, fruit, and pollen morphology [12] and on phylogenetic analysis [13] . Caryopteris comprises seven species at present. Caryopteris incana is a perennial herb, shrub, or subshrub that is distributed in China, the Korean Peninsula, and Japan. Lavender-blue, cymose flowers in the axils of the opposite leaves appear naturally from late summer to fall and continue to flower for one to two months. C. incana is cultivated for use in flower arrangements and gardens. The hybrid cultivar C. × clandonensis (C. incana × C. mongolica Bunge) has a variety of superior properties. As a medicinal ingredient, all parts of C. incana include several kinds of phenylpropanoid glycosides, and it is used as a medicinal herb in China [14] -[18] . In addition, its essential oil extract has been reported to have an insecticidal effect [19] . Its distribution in the Japanese Islands is limited to West Kyushu, and it grows wild in the Kyushu mainland, the Goto Islands of Nagasaki, and the Koshiki Islands of Kagoshima, and, in particular, Tsushima is assumed to be a center of its distribution [20] [21] . C. incana grows wild mainly on sunny bare rock. This species may grow with Selaginella tamariscina, which likes similar environment. Populations of C. incana have decreased because of development and infrastructure maintenance around its natural habitat, and it has been designated as endangered species, listed as “vulnerable” in Japan [22] . However, detailed fieldwork related to C. incana in West Kyushu has not been conducted since 1988 [20] .

We performed this study along two purposes: to construct a detailed distribution map for Tsushima, which is a center of the distribution of C. incana, an endangered species; and to determine the genetic structure using the DNA sequence of populations in Tsushima and to compare this with the geographical structure. These results would offer useful information for determining the evolutionary history of C. incana in Tsushima and for devising a conservation plan.

2. Materials and Methods

2.1. Field Work

We investigated all of Tsushima based on the report by [20] to confirm six locations of exposed rock where C. incana would be expected to grow. In this article, when groups of individuals were separated by more than 2 km, we defined each as a different “population”. We recorded environmental data, such as the latitude/longitude, altitude, population area, and numbers of individual in each natural population. In addition, for the genetic investigation of each population, we collected seeds from mature individuals in all population. We performed the statistical analysis of environmental data by comparing with shore and inland populations using the Student’s t-test in SPSS.

2.2. Sampling and DNA Extraction

We planted the seeds from the natural populations and grew them under the same conditions in a greenhouse. We transplanted seeds into seven bowls and managed the plants in a non-temperate/climate controlled environment from June to July after planting in April 2013. The growth medium consisted of red soil: peat moss: perlite = 7:2:1 without basal fertilizer. Fresh young leaves were gathered from each individual and were kept at −80˚C until DNA extraction. Genomic DNA was extracted by the modified cetyltrimethylammonium bromide (CTAB) method [23] [24] . We adjusted extraction DNA to a concentration of 100 ng/µl by a spectrum altimeter and used polymerase chain reaction (PCR). We used a sample of one or two individuals for the chloroplast DNA sequence analysis, and a bulk sample, which mixed DNA of the same concentration that we extracted and adjusted from 10 individuals for the nuclear DNA sequence analysis.

2.3. DNA Sequence Analyses

To investigate the chloroplast DNA sequence, We used three regions that were registered in Genbank in C. incana that were suitable for the phylogenetic analysis and were more subtle than variation within the genus; matK [25] [26] , trnL-trnF [27] [28] , and rpl32-trnL [29] [30] . The Genbank accession numbers of matK, trnL- trnF, and rpl32-trnL are AF315295, JF301359, and JQ669280 respectively. Therefore, we chose eight intergenic spacer sequences that were determined to be suitable for the classification of closely related species by [29] : trnQ-rps16, atpI-atpH, ndhF-rpl32, petL-psbE, psbD-trnT, psbJ-petA, rps16-trnK, and trnV-ndhC. On the other hand, we used the ITS [31] -[33] of the Genbank registration sequence for the nuclear DNA sequence. The Genbank accession number of the ITS is EF508064. Each primer was prepared from the original paper or sequence information newly in Genbank registration regions, and universal primers were used in the remaining regions (Table 1). The PCR solution was adjusted according to the instruction of the KAPA extra Taq kit (NIPPON Genetics). All reactions were performed with the following program using a Veriti Thermal Cycler (Life technologies): initial denaturation for 2 min at 95˚C; 35 cycles of 20 s at 95˚C, 15 s of annealing at 50˚C for all the chloroplast regions and 56˚C for the ITS region, and 1 or 2 min at 68˚C; a final extension for 2 min at 68 ˚C; and kept at 4˚C until further processing. The amplification was confirmed by agarose gel electrophoresis. PCR products were purified using ExoSAP-IT (GE), which reacted for 30 min at 37˚C, 15 min at 80˚C. Sequencing reactions were performed using a BigdyeTM Terminator v3.1 Cycle Sequencing Kit (Life technologies), and gel filtration was performed for each sample. Sequences were analyzed in a 3500 Genetic Analyzer (Life technologies).

2.4. Phylogenetic Analyses

The base sequences of each population were aligned using the BioEdit software (version 7.2.5) and compared with the registration sequences in the Genbank registration region [34] . Phylogenetic trees were constructed by the neighbor-joining method using the Kimura 2-parameter model in MEGA 6 [35] . The reliability of the topology was assessed with the bootstrap analysis by 1000 replications. Tripora divaricate (Maxim.) P.D. Cantino, which is a related genus of Caryopteris, was included as an outgroup in the phylogenetic analysis. The network figure was drawn using the SplitsTree 4.0 software package [36] .

Table 1. Primer names and sequences for the amplification and cycle sequencing of chloroplast and nuclear DNA.

aK = G or T, R = A or G, Y = C or T.

3. Results and Discussion

3.1. Local Environments

From our field work, we confirmed 72 natural populations in 111 locations throughout Tsushima that were suitable for the growth of C. incana (Figure 2). In agreement with the report of [20] , it seems to be widely distributed in each place over the islands. We confirmed that C. incana grew locally on open rocky sites, such as roadsides, mountains, and shorelines. In such places, it grew wild while avoiding competition with plants at some environments, such as gaps in a natural forest, artificial rocky places by infrastructure maintenance, gaps in concrete surfaces, and bare rock places facing the sea. There were few individuals of C. incana on rocky places covered by creeping vines. Natural populations were not confirmed on rock surfaces along the shore exposed to strong winds, and there were many individuals in inlets. Thus, natural populations of C. incana tended to be locally discontinuous. In addition, populations were confirmed near the slope faces for greening along roadside, and it is possible that their presence at such locations was affected by human activities, such as planting.

There were 39 populations near shores at around 20 m from the sea, and the average altitude in these locations was 5.6 ± 1.4 m above the high tide line (Table 2). On the other hand, there were 33 populations on mountains or small inlands, and the average altitude of these locations was 47.7 ± 6.3 m, covering a wider range than shore populations. There were 22 populations on the south island, at an average altitude of 38.2 ± 8.8 m, and 50 populations on the north island, at an average altitude of 19.0 ± 3.6 m. It was thought that these values were due to higher altitudes and a higher ratio of inland populations on the south island than on the north island, but there were few populations in the center of either island. In addition, the number of individuals in populations along

Figure 2. Geographic distribution of C. incana in tsushima. Black circles indicate each population. Numbers label different populations. Boxed labels indicate populations reported by [20] .

Table 2. The latitude, longitude, altitude, area, and number of individuals of study populations in C. incana.

aover 3000 m2. bover 500 individuals.

the shore tended to be lower than in inland populations (120.5 ± 39.3 < 172.9 ± 36.5, F (2.591) = 0.112, P = 0.268; no significant difference). Because the places where it could grow were more limited along the shore than inland and were surrounded by trees and the sea, it seemed that the shore populations were less likely to expand their distribution in the future. Particularly, isolated populations with few individuals were presumed to have lower fitness from inbreeding depression [37] [38] or Allee effect [39] .

3.2. Sequence Variations

Sequence variations among populations in six of the 11 regions of amplified chloroplast DNA were confirmed using primers prepared from the Genbank registration sequence or universal primers. The intergenic spacer sequence of trnL-trnF region distinguished H2 from other haplotypes by sequence variation in the number of repetitions (Table 3). Similarly, the sequence of the rpl32-trnL region distinguished H1 from other haplotypes by such variation. The sequence of the psbD-trnT region distinguished H1, H2, and H3 from other haplotypes by sequence variations in the number of repetitions and an insertion-deletion (indel) of 9 bases. The sequence of the trnQ-rps16 region distinguished H2 from other haplotypes by sequence variations in the number of repetitions and an indel of 13 bases, and distinguished H1, H2, and H4 from other haplotypes by three substitution sites. The sequence of the rps16-trnK region distinguished H1 and H2 from other haplotypes by sequence variations

Table 3.Polymorphic sites and chloroplast DNA haplotypes based on sequences of seven non-coding regions in C. incana.

in the number of repetitions and an indel of 13 bases, and distinguished haplotypes besides H3 and H4 by four substitution sites. The sequence of H1 was different from other haplotypes in multiple sites. H4, H5 and H6 differed from H3 by sequence variations in a substitution site and the number of repetitions.

The numbers of polymorphic sites in each sequence region indicated a trend that was similar to the results of the numbers of polymorphic sites of each region according to [29] . Several substitution sites between Genbank registration sequences and those of study populations were confirmed in the matK, trnL-trnF and rpl32-trnL regions. However, there were no indications of substitution among the study populations in these sites. Because these regions were used for wider phylogenetic analyses, such as interspecific comparisons, than the regions with substitution among the study populations, the distribution of C. incana in Tsushima may have been relatively recent. On the other hand, the sequence of the ITS region in the nuclear DNA did not have clearly polymorphic sites among study populations.

Two haplotypes were shown mainly in the regions that confirmed sequence variation between populations (Figure 3). The border between H1 and H3 was located in the northern part of the south island (Figure 3), and sequence variations between the haplotypes were present at several sites (Table 3). The shortest straight line distance between populations of the difference haplotypes was approximately 3 km across the northwestern border and approximately 13 km across the northeastern border on the south island. Interestingly, H1, which was located on the south island, was found in a population in the northwest of the north island. This was a large population on the slope along a roadside near public facilities. Because of this local environment, we speculated that

Figure 3. Geographic distribution of chloroplast DNA haplotypes detected in C. incana.

this population was formed from soil used for construction that mixed in seeds of C. incana from some southern population, or planting after slope face spray constructions along a roadside. This suggested that there had been human influence on the distribution of C. incana in Tsushima. In addition, we confirmed the variation in the eastern side of the north island and parts of the central region.

Unlike chloroplast DNA, haplotype in the ITS regions of nuclear DNA did not present genetic a clear genetic structure in the sequence of substitution parts indicating the hetero model. If the sequences in other nuclear DNA regions beside the ITS region were similar among populations with different chloroplast DNA haplotypes, it could be that there is cross-pollination among all populations in Tsushima. In this case, it would be suggested that the present distribution of haplotypes was the result of pollination among populations from throughout the islands after the distribution and fragmentation of populations. Studies using other nuclear DNA regions or microsatellites would be necessary to determine whether this hypothesis is correct. If the genetic structure of the nuclear DNA in Tsushima was confirmed in regions other than the ITS region in accordance with the speculated cross-pollination among populations, the one northern population with a southern chloroplast DNA haplotype might cause introgression to occur by crossing with neighboring populations in the future.

3.3. Phylogenetic Analyses

Because of the limited variation in the intergenic regions where sequence variation can easily occur in chloroplast DNA, we cannot construct reliable family trees among haplotype of each population (Figure 4). Thus, the distribution was more likely to have developed relatively recently, so comparison to sequence variations in other natural habitat areas is necessary. The haplotype network in Figure 5 showed that the main northern haplotype (H3) was genetically different from the southern haplotype (H1). H4, H5, and H6 were speculated to have derived from H3 because they only vary by one base substitution and one site of repetition from H3. H4, which is distributed in the northern part of the south island, was suggested to have expanded after mutation from H3. Because H5 and H6 were confirmed in small populations that seemed to be unlikely to have seed dispersal to or from the shore, they were speculated to be derived by genetic drift from H3 after the present distribution was formed. In contrast to these haplotypes, the haplotype that expanded along the southern part of the south island did not have any derived haplotypes, so it was suggested that the distribution on the north island occurred earlier than on the south island. H2 had a very different sequence from H1 or H3. These populations may have represented

Figure 4. Neighbor-joining tree of chloroplast DNA haplotypes based on the sequences of seven non-coding regions in C. incana. Numbers below the branches indicate the bootstrap values.

Figure 5. Haplotype network among all Tsushima haplotypes in C. incana. Haplotype numbers correspond to Table 3. Small circles in black indicate a substitution.

early original variation because they did not seem to have been affected by human influence from their neighboring environments.

The different distribution of haplotypes in the northern and the southern islands in Tsushima suggested that the distribution might have expanded to north from glacial refugia on the south island, or that C. incana colonized the north and south islands in Tsushima at different times during the repeated northward and southward movement across the land bridge through glacial and interglacial cycles. In the former situation, there would be more differentiated haplotypes on the south island; however, the southern populations had only one haplotype, unlike the northern population. Furthermore, it seemed that the latter situation was likely because the southern haplotype was genetically different from northern haplotypes, but these data could not identify when the northern and southern populations had colonized Tsushima. To estimate these times, we would have to compare these populations with those on the Goto Islands and Kyushu mainland. Then, we might be able to estimate when C. incana expanded its distribution to Japan by studying wild individuals from the Korean Peninsula.

4. Conclusion

In conclusion, we confirmed the present natural environment of Caryopteris incana, which was an endangered species in Japan, and constructed its distribution map in Tsushima. In addition, we investigated the genetic structure among the populations in Tsushima based on chloroplast DNA sequences. Then, we confirmed that C. incana colonized the north and south islands at different times. However, we did not find a clear genetic structure from the nuclear DNA sequence. A population on the north island with the southern haplotype suggested that there might have been some human influence on the distribution of C. incana in Tsushima, and that this population might be able to cause introgression by crossing with neighborhood populations in the future. Additional comparisons of the genetic structure among natural populations outside Tsushima would offer useful information on the evolutionary history of C. incana.

Cite this paper

MasayaAndo,HitoshiWatanabe,KiyoshiMatsubara,AkitoTaniguchi, (2015) Intraspecific Phylogenetic Relationships of Caryopteris incana in the Tsushima Islands, Japan, Using DNA Sequence Analysis. American Journal of Plant Sciences,06,2361-2373. doi: 10.4236/ajps.2015.614239

References

  1. 1. Charlesworth, D. and Charlesworth, B. (1987) Inbreeding Depression and Its Evolutionary Consequences. Annual Review of Ecology and Systematics, 18, 237-268.
    http://dx.doi.org/10.1146/annurev.es.18.110187.001321

  2. 2. Keller, L.F. and Waller, D.M. (2002) Inbreeding Effects in Wild Populations. Trends in Ecology and Evolution, 17, 230-241.
    http://dx.doi.org/10.1016/S0169-5347(02)02489-8

  3. 3. Stephens, P.A., Sutherland, W.J. and Freckleton, R.P. (1999) What Is the Allee Effect? Oikos, 87, 185-190.
    http://dx.doi.org/10.2307/3547011

  4. 4. Huson, D.H. (1998) Splitstree: Analyzing and Visualizing Evolutionary Data. Bioinformatics, 14, 68-73.
    http://dx.doi.org/10.1093/bioinformatics/14.1.68

  5. 5. Tamura, K., Stecher, G., Peterson, D., Filipski, A. and Kumar, S. (2013) MEGA6: Molecular Evolutionary Genetics Analysis Version 6.0. Molecular Biology and Evolution, 30, 2725-2729.
    http://dx.doi.org/10.1093/molbev/mst197

  6. 6. Hall, T.A. (1999) BioEdit: A User-Friendly Biological Sequence Alignment Editor and Analysis Program for Windows 95/98/NT. Nucleic Acids Symposium Series, 41, 95-98.

  7. 7. Huang, M., Crawford, D.J., Freudenstein, J.V. and Cantino, P.D. (2008) Systematics of Trichostema (Lamiaceae): Evidence from ITS, ndhF, and Morphology. Systematic Botany, 33, 437-446.
    http://dx.doi.org/10.1600/036364408784571554

  8. 8. Steane, D.A., Scotland, R.W., Mabberley, D.J. and Olmstead, R.G. (1999) Molecular Systematics of Clerodendrum (Lamiaceae): ITS Sequences and Total Evidence. American Journal of Botany, 86, 98-107.
    http://dx.doi.org/10.2307/2656958

  9. 9. White, T.J., Bruns, T., Lee, S. and Taylor, J. (1990) Amplification and Direct Sequencing of Fungal Ribosomal RNA Genes for Phylogenetics. In: Innis, M.A., Gelfand, D.H., Shinsky, J.J. and White, T.J., Eds., PCR Protocols: A Guide to Methods and Applications, Academic Press, Inc., New York, 315-322.
    http://dx.doi.org/10.1016/b978-0-12-372180-8.50042-1

  10. 10. Drew, B.T. and Sytsma, K.J. (2012) Phylogenetics, Biogeography, and Staminal Evolution in the Tribe Mentheae (Lamiaceae). American Journal of Botany, 99, 933-953.
    http://dx.doi.org/10.3732/ajb.1100549

  11. 11. Shaw, J., Lickey, E.B., Schilling, E.E. and Small, R.L. (2007) Comparison of Whole Chloroplast Genome Sequences to Choose Noncoding Regions for Phylogenetic Studies in Angiosperms: The Tortoise and the Hare III. American Journal of Botany, 94, 275-288.
    http://dx.doi.org/10.3732/ajb.94.3.275

  12. 12. Drew, B.T. and Sytsma, K.J. (2011) Testing the Monophyly and Placement of Lepechinia in the Tribe Mentheae (Lamiaceae). Systematic Botany, 34, 1038-1049.
    http://dx.doi.org/10.1600/036364411X605047

  13. 13. Taberlet, P., Gielly, L., Pautou, G. and Bouvet, J. (1991) Universal Primers for Amplification of Three Non-Coding Regions of Chloroplast DNA. Plant Molecular Biology, 17, 1105-1109.
    http://dx.doi.org/10.1007/BF00037152

  14. 14. Shi, S., Du Y., Boufford, D.E., Gong, X., Huang, Y., He, H. and Zhong, Y. (2003) Phylogenetic Position of Schnabelia, a Genus Endemic to China: Evidence from Sequences of cpDNA matK Gene and nrDNA ITS Regions. Chinese Sci- ence Bulletin, 48, 1576-1580.
    http://dx.doi.org/10.1360/03wc0081

  15. 15. Hilu, K.W. and Liang, H.P. (1997) The matK Gene: Sequence Variation and Application in Plant Systematics. American Journal of Botany, 84, 830-839.
    http://dx.doi.org/10.2307/2445819

  16. 16. Lassner, M.W., Peterson, P. and Yoder, J.I. (1989) Simultaneous Amplification of Multiple DNA Fragments by Polymerase Chain Reaction in the Analysis of Transgenic Plants and Their Progeny. Plant Molecular Biology Reporter, 7, 116-128.
    http://dx.doi.org/10.1007/BF02669627

  17. 17. Doyle, J.J. and Doyle, J.L. (1987) A Rapid DNA Isolation Procedure for Small Quantities of Fresh Leaf Tissue. Phytochemical Bulletin, 19, 11-15.

  18. 18. Environment Agency of Japan (2000) Threatened Wildlife of Japan, Red Data Book 2nd Edition. Japan Wildlife Research Center, 8, 521.

  19. 19. Itow, S., Kim, C.S., Kim, M.H., Kawasato, H. and Oh, J.G. (1993) Biogeography and Ecology of Caryopteris incana Maxim. (Verbenaceae) in Regions of Japan-Korea Strait. Bulletin of the Faculty of Liberal Arts, Nagasaki University, Natural Science, 34, 45-50.

  20. 20. Itow, S. and Kawasato, H. (1988) The Distribution and Ecology of Caryopteris incana Maxim. (Verbenaceae) in Western Kyushu, Japan. Hikobia, 10, 135-143.

  21. 21. Chu, S.S., Liu, Q.Z., Zhou, L., Du, S.S. and Liu, Z.L. (2011) Chemical Composition and Toxic Activity of Essential Oil of Caryopteris incana against Sitophilus zeamais. African Journal of Biotechnology, 10, 8476-8480.

  22. 22. Zhao, D.P., Matsunami, K. and Otsuka, H. (2009) Iridoid Glucoside, (3R)-oct-1-en-3-ol Glycosides, and Phenylethanoid from the Aerial Parts of Caryopteris incana. Journal of Natural Medicines, 63, 241-247.
    http://dx.doi.org/10.1007/s11418-009-0317-9

  23. 23. Li, J. and Wang, Y. (2004) Synthesis of Trisaccharide of Incanoside from Caryopteris incana. Synthetic Communications, 34, 515-522.
    http://dx.doi.org/10.1081/SCC-120027292

  24. 24. Gao, J., Igalashi, K. and Nukina, M. (2000) Three New Phenylethanoid Glycosides from Caryopteris incana and Their Antioxidative Activity. Chemical and Pharmaceutical Bulletin, 48, 1075-1078.
    http://dx.doi.org/10.1248/cpb.48.1075

  25. 25. Gao, J., Igalashi, K. and Nukina, M. (1999) Radical Scavenging Activity of Phenylpropanoid Glycosides in Caryopteris incana. Bioscience, Biotechnology, and Biochemistry, 63, 983-988.
    http://dx.doi.org/10.1271/bbb.63.983

  26. 26. Gao, J. and Han, G. (1997) Cytotoxic Abietane Diterpenoids from Caryopteris incana. Phytochemistry, 44, 759-561.
    http://dx.doi.org/10.1016/S0031-9422(96)00609-7

  27. 27. Cantino, P.D., Wagstaff, S.J. and Olmstead, R.G. (1999) Caryopteris (Lamiaceae) and the Conflict between Phylogenetic and Pragmatic Considerations in Botanical Nomenclature. Systematic Botany, 23, 369-386.
    http://dx.doi.org/10.2307/2419511

  28. 28. Abu-Asab, M.S., Cantino, P.D., Nowicke, J.W. and Sang, T. (1993) Systematic Implications of Pollen Morphology in Caryopteris (Labiatae). Systematic Botany, 18, 502-515.
    http://dx.doi.org/10.2307/2419422

  29. 29. Nakanishi, H. (1996) Plant Species with Northbound Distribution in Western-Kyushu, Japan: Definition, Composition and Origin. Acta phytotaxonomica et Geobotanica, 47, 113-124.

  30. 30. Itow, S. (1997) A Review of Phyto- and Vegetation Geography in the Japan-Korea Strait Region. Bulletin of the Faculty of Liberal Arts, Nagasaki University, Natural Science, 38, 25-51.

  31. 31. Nakanishi, H. (2010) Distribution and Ecology of Islet Biased Plants in northern Kyushu, Japan. Vegetation Science, 27, 1-9.

  32. 32. Aoki, K., Suzuki, T., Hsu, T.W. and Murakami, N. (2004) Phylogeography of the Component Species of Broad-Leaved Evergreen Forests in Japan, Based on Chloroplast DNA Variation. Journal of Plant Research, 117, 77-94.
    http://dx.doi.org/10.1007/s10265-003-0132-4

  33. 33. Bai, W.N., Liao, W.J. and Zhang, D.Y. (2010) Nuclear and Chloroplast DNA Phylogeography Reveal Two Refuge Areas with Asymmetrical Gene Flow in a Temperate Walnut Tree from East Asia. New Phytologist, 188, 892-901.
    http://dx.doi.org/10.1111/j.1469-8137.2010.03407.x

  34. 34. Qiu, Y.X., Fu, C.X. and Comes, H.P. (2011) Plant Molecular Phylogeography in China and Adjacent Regions: Tracing the Genetic Imprints of Quaternary Climate and Environmental Change in the World’s Most Diverse Temperate Flora. Molecular Phylogenetics and Evolution, 59, 225-244.
    http://dx.doi.org/10.1016/j.ympev.2011.01.012

  35. 35. Comes, H.P. and Kadereit, J.W. (1998) The Effect of Quaternary Climatic Changes on Plant Distribution and Evolution. Trends in Plant Science, 3, 1360-1385.
    http://dx.doi.org/10.1016/S1360-1385(98)01327-2

  36. 36. Medail, F. and Diadema, K. (2009) Glacial Refugia Influence Plant Diversity Patterns in the Mediterranean Basin. Journal of Biogeography, 36, 1333-1345.
    http://dx.doi.org/10.1111/j.1365-2699.2008.02051.x

  37. 37. Hewitt, G.M. (1996) Some Genetic Consequences of Ice Ages, and Their Role in Divergence and Speciation. Biological Journal of the Linnean Society, 58, 247-276.
    http://dx.doi.org/10.1111/j.1095-8312.1996.tb01434.x

  38. 38. Avise, J.C., Arnold J., Ball, R.M., Bermingham, E., Lamb, T., Neigel, J.E., Reeb, C.A. and Saunders, N.C. (1987) Intraspecific Phylogeography: The Mitochondrial DNA BRIDGe between Population Genetics and Systematics. Annual Reviews of Ecology and Systematics, 18, 489-522.
    http://dx.doi.org/10.1146/annurev.es.18.110187.002421

  39. 39. Darwin, C. (1859) On the Origins of Species by Means of Natural Selection. John Murray, London.

NOTES

*Corresponding author.