Phytotaxa 373 (1): 053–070 http://www.mapress.com/j/pt/ Copyright © 2018 Magnolia Press ISSN 1179-3155 (print edition) Article PHYTOTAXA ISSN 1179-3163 (online edition) https://doi.org/10.11646/phytotaxa.373.1.2 Molecular and cytological evidences denied the immediate-hybrid hypothesis for Saxifraga yuparensis (sect. Bronchiales, Saxifragaceae) endemic to Mt. Yubari in Hokkaido, northern Japan SAYA TAMURA1,2, TOMOKO FUKUDA3, ELENA A. PIMENOVA4, EKATERINA A. PETRUNENKO4, PAVEL V. KRESTOV4, SVETLANA N. BONDARCHUK5, OLGA A. CHERNYAGINA6, YOSHIHISA SUYAMA7, YOSHIHIRO TSUNAMOTO7, AYUMU MATSUO7, HAYATO TSUBOI8, HIDEKI TAKAHASHI9, KEN SATO9, YOKO NISHIKAWA10, TAKASHI SHIMAMURA10, HIROKO FUJITA11 & KOH NAKAMURA11 1 Division of Environmental Resources, Graduate School of Agriculture, Hokkaido University, Sapporo 060-8589, Japan Japan Wildlife Research Center, Tokyo 130-8606, Japan 3 College of Liberal Arts and Sciences, Mie University, Tsu 514-8507, Japan 4 Botanical Garden-Institute, Far Eastern Branch of the Russian Academy of Sciences, Makovskii Str. 142, Vladivostok, 690024, Russia 5 Sikhote-Alin State Nature Biosphere Reserve, Ministry of Nature Resources and Environment, Terney 692150, Russia 6 Kamchatka Branch of Pacific Institute of Geography, Far Eastern Branch of the Russian Academy of Sciences, PetropavlovskKamchatsky 683000, Russia 7 Field Science Center, Graduate School of Agricultural Science, Tohoku University, Naruko-onsen, Osaki 989-6711, Miyagi, Japan 8 Hakuba-Goryu Alpine Botanical Garden, Hakuba 399-9211, Japan 9 The Hokkaido University Museum, Hokkaido University, Sapporo 060-0810, Japan 10 Institute of Environmental Sciences, Hokkaido Research Organization, Sapporo 060-0819, Japan 11 Botanic Garden, Field Science Center for Northern Biosphere, Hokkaido University, Sapporo 060-0003, Japan; kohnakamur@gmail.com 2 Abstract An alpine plant Saxifraga yuparensis is endemic to a scree consisting of greenschist of Mt. Yubari in Hokkaido, Japan and it has been proposed as an immediate hybrid derived from two species of the same section Bronchiales based on morphological intermediacy: namely S. nishidae, a diploid species endemic to a nearby cliff composed of greenschist and tetraploid S. rebunshirensis comparatively broadly distributed in Japan and Russian Far East. Saxifraga yuparensis is red-listed and it is crucial for conservation planning to clarify whether this is an immediate hybrid and lacks a unique gene pool. The immediate-hybrid hypothesis was tested by molecular and cytological data. In nuclear ribosomal and chloroplast DNA trees based on maximum parsimony and Bayesian criteria, S. yuparensis and S. rebunshirensis formed a clade with several other congeners while S. nishidae formed another distinct clade. Genome-wide SNP data clearly separated these three species in principal coordinate space, placing S. yuparensis not in-between of S. rebunshirensis and S. nishidae. Chromosome observation indicated that S. yuparensis is tetraploid, not triploid directly derived from diploid-tetraploid crossing. Additionally, observation of herbarium specimens revealed that leaf apex shape of S. yuparensis fell within the variation of S. rebunshirensis. These results indicate that S. yuparensis is not an immediate hybrid of S. rebunshirensis and S. nishidae but a distinct lineage and an extremely narrow endemic species, that deserves for intensive conservation. Keywords: Alpine plant, Bronchiales, diabasic plant, genome-wide SNP, hybrid, Japan, MIG-seq, Russia, Saxifraga, tetraploid Introduction The genus Saxifraga L. (1753: 398) (Saxifragaceae, e.g., Soltis et al. 2001) encompasses more than 440 species and is composed of at least 13 sections and nine subsections (Tkach et al. 2015). One of the Saxifraga sections, sect. Bronchiales DeChaine E.G. (2014: 27) was subdivided from sect. Trachyphyllum Koch W.D.J. (1836: 270) based on molecular phylogenetic analyses (DeChaine et al. 2013, DeChaine 2014). This section has approximately 20 species (DeChaine et al. 2013, DeChaine 2014) mostly in northern high latitudes, and among them, three species are distributed in Hokkaido, northern Japan: S. yuparensis Nosaka (1974: 149), S. rebunshirensis Sipliv. (1971: 155), and S. nishidae Accepted by Nicodemo Passalacqua: 7 Sept. 2018; published: 23 Oct. 2018 Licensed under a Creative Commons Attribution License http://creativecommons.org/licenses/by/3.0 53 FIGURE 1. Saxifraga yuparensis (A, D & J 26 July 2016, G 19 July 2017, Gama-iwa at Mt. Yubari), S. rebunshirensis (B, E, H & K, 10 June 2016, Rebun Island), and S. nishidae (C, F, I & L, 26 July 2016, Gama-iwa at Mt. Yubari). A, B & C Habitat. D, E & F Habit. G, H & I Spots on petal. J, K & L Leaf. Scale bars are 1 cm for D–F and 5 mm for G–L. The figure is available in color online. Miyabe & Kudô (1917: 170) (Fig. 1 & 2). Saxifraga rebunshirensis is broadly distributed in Japan (Hokkaido and high mountains in Honshu) and Russian Far East (Sakhalin and the Kuril Islands), while S. yuparensis is endemic to a scree consisting of greenschist at Mt. Yubari in Hokkaido and S. nishidae is endemic to a cliff composed of greenschist at Mt. Yubari and a cliff of schalstein at Mt. Ashibetsu, Hokkaido (Shimizu 1983, Charkevicz 1989, Sato 2007, Nakagawa & Sato 2015). Mountains composed of diabasic rocks (including greenschist, shalestein, etc.) are known to house plants 54 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. adapted to infertile soils and fragile geological features, as is the case in several endemics in Mt. Yubari (Watanabe 1971, Nosaka 1974). Diabasic plants are typically narrow endemics and grow on unstable soils (Watanabe 1971), and thereby include many endangered species. The number of individuals and patches of S. yuparensis are decreasing due to collapse of rocks and collection for horticultural purposes, and it is designated as an endangered species: CR (critically endangered) at the national level (Ministry of the environment of Japan 2017) and En (Endangered) at the provincial level (Hokkaido 2001). The diagnostic characters of these three species are the color of petal spots, presence/absence of a conspicuous claw of petal base, and the shape of leaf apex. Saxifraga yuparensis has yellow petal spots, petal claws, and slightly tricuspidate leaves; S. rebunshirensis has yellow and red (irregularly only yellow) petal spots, claw-less petals, and cuspidate leaves; S. nishidae has yellow and red petal spots, claw-less petals, and tricuspidate leaves (Miyabe & Kudô 1917, Hara 1952, Murata 1961, Nosaka 1974, Nosaka 1980, Shimizu 1983, Ohba 1999, Nakagawa & Sato 2015). However, several authors considered that S. yuparensis is an immediate hybrid derived from crossing between S. rebunshirensis and S. nishidae based on morphological intermediacy in the shape of leaf apex (Toyokuni 1988, Nishikawa et al.1992, Iwatsuki & Kato 1994, Umezawa 2004, Shimizu et al. 2014, Takahashi 2015a). Additionally, seeds of S. yuparensis rarely germinate (Shimizu 1983) and this suggests that the species is an infertile triploid hybrid, because S. nishidae is diploid (2n = 26, Funamoto & Nakamura 1993) while S. rebunshirensis is tetraploid (2n = 48 using samples with no locality information, Sakai, 1935; 2n = 50 with samples from Nagano of central Japan, Funamoto & Nakamura 1990), although the chromosome number of S. yuparensis has not been reported. The hybridorigin hypothesis, however, has not been tested by molecular analyses. On the other hand, S. yuparensis is sometimes treated as a synonym or variety of S. rebunshirensis (Shimizu 1983, DeChaine 2014, Takahashi 2015a, Okuyama 2016; note the last one is the latest Japanese flora). For planning the conservation of S. yuparensis, it is crucial to clarify whether this is an immediate hybrid and lacks a unique gene pool. In a preceding molecular phylogenetic study of Saxifraga, the species relationships were partially revealed, but the relationship among the three species remained unclear (DeChaine et al. 2013). The aim of this study is to elucidate the species relationships and test the immediate-hybridity of S. yuparensis using molecular analyses and cytological observation. In molecular analyses, DNA sequencing and genome-wide single nucleotide polymorphism (SNP) detection were employed. Collaterally, the diagnostic character of leaf apex shape was investigated using specimens of the three species to reevaluate the morphological intermediacy of S. yuparensis. FIGURE 2. Maps of East Asia (A) showing the species range of Saxifraga rebunshirensis and of Hokkaido (B) showing the range of S. yuparensis (red; square), S. rebunshirensis (blue; circle), and S. nishidae (green; triangle). (the map image from Google Earth). The figure is available in color online. MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 55 Material & Methods Molecular analyses Taxon sampling Molecular phylogenetic analyses were conducted using samples of our collection and DNA sequence data from GenBank (Table 1). Saxifraga yuparensis grows on a rocky scree of Mt. Yubari and has only two patches (0.3 x 0.2 m and 0.5 x 1.5 m), being located only ca. 5 m away from each other. One sample from each of the two patches (two samples in total) was collected to avoid collecting the same clones because this species can propagate vegetatively (authors’ observation). DeChaine et al. (2013) sequenced a sample of “S. nishidae” (RBGE-E00295524) but we identified this specimen as S. yuparensis based on leaf apex morphology and found that the specimen was confused and mislabeled with genuine S. nishidae (RBGE-E00295525) collected on the same day on Mt. Yubari. The data of “S. nishidae” (RBGE-E00295524) was used for S. yuparensis. Saxifraga nishidae has more than 100 patches on a nearby rocky cliff of greenschist of Mt. Yubari. Six samples were collected from different patches with at least 5 m interval. In Mt. Yubari, S. rebunshirensis (Ken Sato 85.0207 in SAPS, Appendix) had been collected more than 30 years ago from Rosoku-iwa, that is several kilometers away from the rocky place where S. yuparensis and S. nishidae grow; but this time we could not collect the species there. We did not conduct DNA extraction/PCR with the herbarium specimen because the specimen was old and we were afraid to damage the valuable specimen in vain. Instead, two samples each from two localities in Hokkaido, Rebun Island and Mt. Nishi-Kumaneshiri, were used. Also, two plants from Mt. Hakuba in central Japan were hired. GenBank data for S. rebunshirensis from Moneron Island to the south of Sakhalin (UBC-V164570) were also used. To elucidate the phylogenetic positions of these three species, we incorporated allied 16 species of sect. Bronchiales (DeChaine et al. 2013) in the analyses. In Russian Far East, we collected S. ascoldica Sipliv. (1971: 156) in Primorsky Krai and S. cherlerioides D.Don (1822: 382) and S. funstonii (Small) Fedde (1905: 613) in Kamchatka. For the other species, DNA sequence data were obtained from GenBank. Distribution of the 19 species covered the whole distribution of sect. Bronchiales. For outgroups, two species of sect. Gymnopera D.Don (1822: 343) and one species of sect. Cymbalaria Griseb. (1843: 336) were selected (Table 1), following the result of Tkach et al. (2015). DNA extraction and sequencing Total genomic DNA was extracted from fresh leaves using the cetyltrimethyl ammonium bromide (CTAB) method (Doyle & Doyle, 1987) with some modifications. Preceding phylogenetic studies on Saxifraga employed internal transcribed spacer (ITS) of nuclear ribosomal DNA (nrDNA) and trnL-F intergenic spacer of chloroplast DNA (cpDNA) (DeChaine et al. 2013, Tkach et al. 2015). We used primers following these preceding studies: trnL(UAA) and trnF(GAA) for trnL-F region (Taberlet et al. 1991), and 17SE and 26SE for the entire ITS1, 5.8S and ITS2 region (Sun et al. 1994). PCR was performed in 25 μl total volume with the following reagents: about 10 ng of genomic DNA, 1 unit of Taq DNA polymerase master mix (Ampliqon, Rødovre, Denmark), 0.4 μM of each primer, and 4 % DMSO. After an initial heating step at 94°C for 3 min, samples were incubated for 25 cycles of 94°C for 1 min, 60°C (ITS) or 52°C (trnL–F) for 1 min, 72°C for 2 min, with final extension at 72°C for 5 min. The cycle sequencing reaction was carried out with a Big Dye Terminator Cycle Sequencing Kit ver. 3.1 (Applied Biosystems, Foster City, CA, USA) with the same primers used in the PCR. Direct sequencing was performed on an ABI Prism 3130 DNA analyzer (Applied Biosystems). The sequence data were deposited in DDBJ (DNA Data Bank of Japan) database (2017, Table 1). Phylogenetic analyses DNA sequences were aligned using ClustalX ver. 1.8 (Thompson et al. 1997) implemented in BioEdit ver. 7.2.5 (Hall 1999) and then manually adjusted. Phylogenetic analyses were based on a maximum parsimony (MP) criterion using PAUP* ver. 4.0b10 (Swofford 2002) and a Bayesian approach using MrBayes ver. 3.2 (Ronquist et al. 2012). In the Bayesian phylogenetic analysis, MrModeltest 2.3 (Nylander 2008) was used to estimate the appropriate evolutionary model of nucleotide substitution based on the Akaike Information Criterion (AIC). Based on the model selected, two separate runs of Metropolis coupled Markov chain Monte Carlo (MCMCMC) analyses were performed, each with a random starting tree and four chains (one cold and three heated). The MCMCMC length was one million generations, and the chain was sampled every one hundredth generation from the cold chain. The mixing and convergence of the MCMC chains of the two runs was assessed by inspection of the trace plots of parameters using Tracer ver. 1.6 (Rambaut et al. 2013). The first 250 sample trees (2.5% of the total 10,000 sample trees) were discarded as burn-in. After the burn-in, the effective sample sizes (ESS) of all parameters were more than 500, indicating that 56 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. MOLECULAR AND CYTOLOGICAL EVIDENCES TABLE 1. Samples used for ITS and cpDNA sequencing analyses. N is the number of sequenced individuals. Phytotaxa 373 (1) © 2018 Magnolia Press • 57 Taxon Sect. Bronchiales S. anadyrensis Losinsk. S. ascoldica Sipliv. S. austromontana Wiegand Distribution* N Voucher ITS accession no. trnL-F accession no. RU: Magadan to mouth of Lena River RU: Ussuri NA: Cascade Mountains, Rocky Mountains 2 2 11 KF196358 LC314071, LC314072 KF196332–KF196347 KF196374 LC314090, LC314091 KF196398–KF196409 S. bronchialis L. RU: northern RU to Pacific Coast 2 154288, 78524 (ALA) 042206, 042207(SAPT) 619568, 750642, 750643, 750647, 758235, 779011, 780148 (RMH), 361429, 369601 (WTU), 22800, 22801 (WWB) 283068, 282875 (WTU) KF196360, KF196361 S. caulescens Sipliv. S. cherlerioides D.Don RU: Siberia; Mongolia RU: Sakhalin, Kuril Islands, Okhotsk, NA: Alaska, Yukon 1 8 00258082 (RBGE) KF196325 358172, 385184, 385185 (WTU), 042200– KF196363, KF196364, 042205 (SAPT) KF196366, LC314087 S. codyana Zhmylev S. derbekii Sipliv. S. funstonii (Small) Fedde NA: Alaska, Yukon RU: Far East near Sea of Okhotsk RU: Taimyr to Chukotka, Kamchatka, NA: Alaska, Yukon, British Columbia 1 1 10 212762 (UBC) 164312 (UBC) 121070, 156966, 156968 (ALA), 356796 (WTU), 22807–22809 (WWB), 042197042199(SAPT) KF19 6396, KF196397 KF196375 KF196376–KF196378, LC314096, LC314098, LC314100, LC314105– 314107 KF196379 KF196380 KF196393–KF196395, KF196410– KF196412, KF196417, LC314099, LC314101, LC314999, LC316189 S. kruhsiana Fisch. ex Ser. S. nishidae Miyabe & Kudô RU: Kamchatka, Okhotsk Japan: Hokkaido; Yubari Mountains, Ashibetsu Mountains RU: Chukotka to Okhotsk Japan: Hokkaido, high mountains in Honshu, RU: Sakhalin, Kuril Islands RU: Urals to Chukotka, Kamchatka RU: Angara to Okhotsk, Zea, Lk. Baikal NA: British Columbia Haida Gwaii NA: Alaska, Canada, Greenland NA: Olympic Mountains, Cascade Mountains Japan: Hokkaido; Yubari Mountains 2 6 107779 (ALA) 042191–042196(SAPT) 1 7 129124 (ALA) 164570 (UBC), 042185–042190(SAPT) 2 1 3 4 3 107682, 168369 (ALA) 154286 (ALA) 214354, 214965 (UBC), 22805 (WWB) 363418 (RMH), 22802–22804 (WWB) 284404 (WTU), 159757 (UBC), 22806 (WWB) 00295524 (RBGE), 042183–042184 (SAPT) S. omolojensis A.P.Khokhr. S. rebunshirensis (Engl. & Irmsch.) Sipliv. S. spinulosa Adams S. stelleriana Merk. ex Ser. S. taylorii Calder & Savile S. tricuspidata Rottb. S. vespertina (Small) Fedde S. yuparensis Nosaka Sect. Cymbalaria – S. hederacea L. – S. cymbalaria L. Sect. Gymnopera – S. spathularis Brot. * NA, North America; RU, Russia. 3 KF196327 KF196367 KF196326, KF196328– KF196331, KF196362, KF196365, LC314070, LC314088, LC314089 KF196356, KF196357 KF196381, KF196382 LC314081–LC314086 LC314108–LC314113 KF196368 KF196323, LC314073– 314078 KF196354, KF196355 KF196359 KF196351–KF196353 KF196339–KF196342 KF196348–KF196350 KF196384 KF196385, LC314092–314095, LC314102, LC314103 KF196386, KF196387 KF196388 KF196389, KF196390 KF196413–KF196416 KF196391, KF196392 KF196324, LC314079, KF196383, LC314097, LC314104 LC314080 1 1 – – AF261182 – – AF374777 1 – AJ233858 AF374785 the analyses sampled the posterior distributions of each parameter satisfactorily, and the values of Average Standard Deviation of Split Frequency (ASDSF) were below 0.008. The 50% majority rule consensus tree and Bayesian posterior probabilities (PP) of all the post-burn-in trees was generated using FigTree ver. 1.4.2 (Rambaut 2014). In the MP phylogenetic analysis, indels were treated as missing data. The characters were treated as unordered, and the character transformations were weighted equally. The branch collapse option was set to collapse at a minimum length of zero. A heuristic parsimony search was performed with 1000 replicates of random additions of sequences with ACCTRAN character optimization, tree bisection-reconnection (TBR) branch swapping, and MULTREES and STEEPEST DESCENT options on. Statistical support for each clade was assessed by bootstrap analysis (Felsenstein 1985). One thousand replicates of heuristic searches, with the TBR branch swapping switched on and MULTREES options off, were performed to calculate bootstrap percentages (BP). Genome-wide SNP analysis To test the immediate hybridity and species relationships based on multiple independent nuclear loci, genome-wide SNP detection was conducted with the focal three species. The samples of our own collection used in the sequencing analysis were employed. Preceding studies (Rutledge et al. 2015, Karimi et al. 2016) based on genome-wide SNP data (observed and simulated data) successfully indicated that hybrids were plotted at the intermediate between the parental species in a principal coordinate analysis (PCA) (e.g., F1 and F2 individuals occurred at intermediate between the two parental species, and backcross individuals occurred between the F1 and the one parental species; Rutledge et al. 2015). For SNP detection, we employed a method MIG-seq [i.e., multiplexed inter simple sequence repeats (ISSRs) Genotyping by sequencing; Suyama & Matsuki 2015]. MIG-seq is a type of reduced representation sequencing, where ISSRs are amplified by multiplex PCR. Thereby the method is more applicable to low-quality genomic DNA than RADseq. Most of SNPs obtained by MIG-seq are putatively selectively neutral because these are immediately adjacent to two repeat regions (Takahashi et al. 2016). Highly reduced representation libraries were constructed following Suyama & Matsuki (2015) and sequenced on an Illumina MiSeq sequencer (Illumina, San Diego, CA, USA). Both ends of the fragments were read by paired-end sequencing. Sequences of reads were quality-filtered using FASTX-Toolkit (http://hannonlab.cshl.edu/fastx_toolkit/), with the settings of q = 30 and p = 40. The quality-filtered reads were then bundled together to form putative loci (stacks) with the software Stacks v1.15 (Catchen et al. 2013) using the ‘ustacks’ option, with the settings of minimum depth of coverage required (m) = 20, maximum distance between stacks (M) = 2, and the deleveraging (d) and removal (r) algorithms enabled. Using the ‘cstacks’ option, a catalog was created for all possible loci and alleles with the parameter ‘number of allowed mismatches between samples (n)’ = 4. The stacks created by ‘ustacks’ were then searched against the catalogue using the ‘sstacks’ option. The data set of all the samples was considered as a single population and SNPs were retrieved using the ‘populations’ option, with the ‘write_single_snp’ option to select only the first SNP per locus to avoid the linkage between SNPs. Only loci that were shared among at least 75 % of the samples were considered. If a locus was not recovered in a sample, the genotype of the sample was treated as missing data. PCA of SNP genotypes was conducted using GENODIVE (Meirmans & Van Tienderen 2004) with the default parameter settings. Chromosome observation For chromosome observation, we used four individual of Saxifraga yuparensis from Mt. Yubari (Saya Tamura 992, 993, 994, 995, SAPT), four individuals of S. rebunshirensis from Mt. Nishi-kumaneshiri (Saya Tamura 434, 436, SAPT) and Rebun Island (Saya Tamura 328, 329, SAPT) of Hokkaido, and four individuals of S. nishidae from Mt. Yubari (Saya Tamura 996, 997, 998, 999, SAPT). Plant materials were collected in the field and cultivated in Botanic Garden, Hokkaido University. The methods for the observation of somatic chromosome numbers followed Fukuda et al. (2007, 2016) with some modifications. Root tips were pretreated in 0.002 M 8-hydroxyquinoline solution for 10 hours at 4°C, and fixed in Farmer’s solution (glacial acetic acid: 99% ethanol = 1 : 3) for more than 16 hours at 4°C. Then the root tips were macerated in 1N HCl for 10 minutes at 60°C, stained with 1 % aceto-orcein for 10 minutes, and squashed. Chromosomes were observed under a binocular Zeiss Axio Imager A1 microscope (Carl Zeiss SpA, Italy), and photomicrographs were taken using a Anyty 3R-DKMC01 camera (3R solution corporation, Japan). Morphological observation A diagnostic character of leaf apex shape was observed in the three species. Typically, S. rebunshirensis has cuspidate leaves while S. nishidae has tricuspidate leaves; and S. yuparensis has slightly tricuspidate leaves (Shimizu 1983, Ohba 58 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. 1999, Wakabayashi 2001, Shimizu 2014, Okuyama 2016). Herbarium specimens of the three species (Appendix) were examined in the herbaria of Hokkaido University Museum (SAPS), Botanic Garden, Hokkaido University (SAPT) and the National Museum of Nature and Science (TNS) to elucidate the variability of the diagnostic character especially within S. rebunshirensis, that has large species range. Results Phylogenetic relationships based on ITS The aligned length of the ITS data was 822 bp. In S. rebunshirensis, five sequence types were recognized among the seven samples, and S. nishidae had two sequence types in the six samples. The two samples of S. yuparensis of our collection (SAPT-042183, 042184) showed one type of sequence with a double-peak signal (adenin and guanin) at one site (position 423 in the registered sequences) in the sequence chromatogram, while all the samples of S. rebunshirensis and S. nishidae had only adenin but not guanin at the same site. This site was coded as R in the two S. yuparensis samples. In the Bayesian analysis, the AIC selected the SYM+I+G model. The 50% majority rule consensus tree of all the post-burn-in trees is depicted (Fig 3). In the MP analysis, 201 nucleotide substitutions were found in 57 variable sites and 114 sites were parsimony informative among them. 15,104 equally parsimonious trees of 314 steps were obtained with a consistency index (CI) = 0.78, a retention index (RI) = 0.91, and a resealed consistency index (RC) = 0.71. The topology of the strict consensus tree was the same as that of the Bayesian tree, and therefore BPs of the MP analysis are plotted on the Bayesian tree (Fig 3). In the following, only clades with PP ≥ 0.95 and/or BP ≥ 70% were considered adequately supported. The phylogenetic tree had moderate resolution and supported monophyly for several species, including S. nishidae (PP = 1.00 / BP = 100). On the other hand, monophyly was not supported for each S. yuparensis and S. rebunshirensis; the two species were included in a clade (1.00 / 55.4) with a portion of samples (three samples) of S. funstonii and S. caulescens Sipliv. (1971: 151). This clade was nested within a clade with the other five samples of S. funstonii and one S. codyana (1.00 / 66.3). Thereby, S. yuparensis plus S. rebunshirensis were phylogenetically distinct from S. nishidae. Phylogenetic relationships based on trnL-F The aligned length of the trnL-F data was 446 bp. All the samples of S. rebunshirensis and S. yuparensis had the same sequence, that was shared with other 11 species in the same clade (below). In the six samples of S. nishidae, only one type of sequence was recovered. In the Bayesian analysis, the AIC selected the GTR model. The 50% majority rule consensus tree of all the post-burn-in trees is depicted (Fig 4). In the MP analysis, 28 nucleotide substitutions were found in 12 variable sites and 13 sites were parsimony informative among them. 574,318 equally parsimonious trees of 30 steps were obtained with CI = 1. The topology of the strict consensus tree was the same as that of the Bayesian tree, and therefore BPs of the MP analysis are plotted on the Bayesian tree (Fig 4). While the phylogenetic tree largely had low resolution, but nevertheless S. nishidae formed a well-supported clade (1.00 / 97.2), and this clade was sister to a portion of samples (three samples) of S. austromontana (0.97 / 63.0). In contrast, each S. yuparensis and S. rebunshirensis was not monophyletic; these two species were recovered in the same polyphyletic clade (0.98 / 61.1) with other 11 species. Species relationships based on genome-wide SNPs The resultant data set comprised 699 SNP loci. The result of the PCA based on the SNP genotyping data is shown (Fig. 5). The first and second principal components explained 51.6 % and 11.5 % of the variances. The SNP genotyping data clearly separated the three species in the principal coordinate space; the plots of S. nishidae were separated from S. yuparensis and S. rebunshirensis along the first principal coordinate, and those of S. yuparensis and S. rebunshirensis were distinguished along the second principal coordinate. Note that S. yuparensis was not placed in-between of S. rebunshirensis and S. nishidae. MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 59 FIGURE 3. Bayesian majority-rule consensus tree of Saxifraga sect. Bronchiales based on ITS. The numerals on branches are Bayesian posterior probabilities (PP: upper) and bootstrap percentages (BP: lower) in the MP analysis. The figure is available in color online. 60 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. FIGURE 4. Bayesian majority-rule consensus tree of Saxifraga sect. Bronchiales based on trnL-F. The numerals on branches are Bayesian posterior probabilities (PP: upper) and bootstrap percentages (BP: lower) in the MP analysis. The figure is available in color online. Chromosome number The chromosome number of S. yuparensis was first determined as 2n=48 (Fig 6–A & D). In S. rebunshirensis the chromosome number was 2n=48 in all the three individuals examined (Fig 6–B & E) and this is the same as reported by Sakai (1935) but different from the number reported by Funamoto & Nakamura (1990) (2n = 50). The chromosome number of S. nishidae was 2n=26 in all the four individuals (Fig 6–C & F), and that was the same as reported by Funamoto & Nakamura (1993). MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 61 FIGURE 5. Plots of principal coordinate analysis (PCA) based on genome-wide SNPs of Saxifraga yuparensis (red; square), S. rebunshirensis (blue; circle), and S. nishidae (green; triangle). The figure is available in color online. FIGURE 6. Photomicrographs of somatic metaphase chromosomes of Saxifraga yuparensis (A), S. rebunshirensis (B) and S. nishidae (C). D, E and F are explanatory drawings of A, B and C, respectively. Scale bar represents 5 μm. Variability of leaf apex morphology All the 14 specimens examined for S. yuparensis had slightly tricuspidate leaf apex except for leaves on flowering stems (Fig 7–A & D, Appendix). In S. rebunshirensis, 199 out of 213 specimens examined had cuspidate leaves (Fig 7–B & E, Appendix) but 14 specimens collected from Hokkaido, Honshu, and Sakhalin had slightly tricuspidate leaves on vegetative stems (Fig 7–F, Appendix). On the other hand, all the 26 specimens of S. nishidae had conspicuously tricuspidate leaf apex except for ones on flowering stems (Fig 7–C & G, Appendix), and the specimens of S. nishidae were clearly distinguished from the other two species. 62 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. FIGURE 7. Specimens of Saxifraga yuparensis (A & D), S. rebunshirensis (B, E & F), and S. nishidae (C & G). Scale bars are 5 mm for D–G. The figure is available in color online. MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 63 Discussion In the molecular phylogenetic analyses based on ITS and cpDNA, S. yuparensis was not monophyletic and recovered in the same clade with S. rebunshirensis plus several other species (S. funstonii and S. caulescens in ITS; 11 species in trnL-F). On the other hand, S. nishidae formed a well-supported monophyletic clade in both ITS and cpDNA, and this clade was distinct from the clades including S. yuparensis and S. rebunshirensis. If S. yuparensis were a hybrid derivative of the two species, two distinct sequences from each of the S. yuparensis samples would be recovered in both the clades of S. rebunshirensis and S. nishidae in the ITS phylogeny, although it cannot be excluded that the ITS allele of one parental species (here S. rebunshirensis) became dominant by concerted evolution via unequal crossingover and gene conversion through time (given that S. yuparensis is not an immediate but old hybrid). However, the PCA based on the SNP genotyping data clearly separated the three species and S. yuparensis samples were not placed in-between of S. rebunshirensis and S. nishidae. This result, being based on multiple independent nuclear loci, also negate the hybrid hypothesis of S. yuparensis. The result also likely negates the scenario that ITS allele of S. rebunshirensis became dominant by concerted evolution in S. yuparensis. Chromosome observation revealed that S. yuparensis is tetraploid. If it were an immediate hybrid derived from diploid-tetraploid crossing, it would be triploid. Thus, the chromosome data negate, at least, the immediate hybridity. Considered in a comprehensive way, all the results based on ITS and cpDNA sequences, genome-wide SNPs, and chromosome count negate the hypothesis that S. yuparensis is derived from hybridization between S. rebunshirensis and S. nishidae. Concerning the infertility of S. yuparensis, Shimizu (1983) reported that the seeds of S. yuparensis did not germinate, and we have not observed mature fruits in the wild population. The infertility of S. yuparensis is probably caused by self-incompatibility, not by triploidy. Self-incompatibility in Saxifraga was reported for species of the same section Bronchiales (Goertzen 1996). Currently there are two patches of S. yuparensis at Mt. Yubari but these patches were originally one patch, that were subsequently fragmented by rock debris (authors’ observation), and the two patches likely contain very closely related plants such as those of the same sibling that are incompatible. Concerning the taxonomic treatment of S. yuparensis, the SNP analysis suggests that S. yuparensis is a lineage distinct from S. rebunshirensis. Note, however, that S. rebunshirensis itself has taxonomic problems and is sometimes treated as a subspecies or variety of S. cherlerioides, S. bronchialis L. (1753: 400), or S. funstonii (Engler et al. 1919, Hara 1937, Hara 1952, Toyokuni 1975, Takahashi 2015b). The present phylogenetic analyses are not conclusive about the species relationships because monophyly was not supported for each of these species and because the trees lacked necessary resolution. The SNP analysis did not include the other species treated in the phylogenetic analyses, specifically S. funstonii and S. caulescens clustered with S. yuparensis and S. rebunshirensis in the ITS tree (and in the cpDNA tree with the other nine species, although it showed lower resolution). Full revision of their taxonomy needs more intensive analysis of the allied species using genome-wide SNPs, but it is beyond the scope of the present study (i.e., test on the hybrid hypothesis). When we follow the idea to treat S. rebunshirensis as an independent species (Siplivinsky 1971; DeChaine 2014), what is the diagnostic character to separate it from S. yuparensis? Slightly tricuspidate leaf is a character based on which S. yuparensis has been distinguished from S. rebunshirensis. In the morphological observations of specimens, it was revealed that S. rebunshirensis from Hokkaido, Honshu, and Sakhalin sometimes had slightly tricuspidate leaves. The leaf apex shape of S. yuparensis fell within the variation of S. rebunshirensis and thereby is not a diagnostic character. In addition, Barkalov (2009) conducted morphological comparison between S. yuparensis from Mt. Yubari and S. rebunshirensis samples from the southern Kuril Islands, where S. yuparensis is not reported and which should be treated as S. rebunshirensis (Charkevics 1989; Takahashi 2015b); he recognized no morphological difference. According to literature, another diagnostic character of S. yuparensis and S. rebunshirensis is a conspicuous claw of petal base in the former species, that is lacked in the latter species (Nosaka 1980; Shimizu et al. 2014), while Barkalov (2009) did not mention the character. Although it was difficult to examine the character with dried herbarium specimens, by observing living plants in wild and botanic gardens, we also recognized the difference in the two species (authors’ observations, Appendix). In conclusion, S. yuparensis is not an immediate hybrid of S. rebunshirensis and S. nishidae but a distinct lineage, as was revealed by the DNA sequence, genome-wide SNP, and cytological evidences. This is an extremely narrow endemic species, that is morphologically distinguished from the congeners in Hokkaido by having a claw of petal base, and deserves for intensive conservation. 64 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. Acknowledgement We thank Seiichiro Miyamato, Minori Somada, Kenji Murakami, and Shogo Kikuta for assistance in field sampling and Hiroyuki Sato for help in herbarium work. This study was supported in part by Grants-in-Aid for Scientific Research, KAKENHI to K.N. (16K18596) and a research grant from the Mitsui & Co. Environment Fund to K.N. (R15–0067). Reference Barkalov, V.Y. (2009) Flora of the Kuril Islands. Dalnauka, Vladivostok, 468 pp. Catchen, J., Hohenlohe, P., Bassham, S., Amores, A. & Cresko, W. (2013) Stacks: an analysis tool set for population genomics. Molecular Ecology 22 (11): 3124–3140. https://doi.org/10.1111/mec.12354 Charkevicz, S.S. (1989) Saxifragaceae. In: Charkevicz, S.S. (Ed.) Plantae vasculares orientis extremi sovietici 4. NAUK, Moscow, pp. 122–190. DDBJ (2017) DNA Data Bank of Japan, National Institute of Genetics, Japan, Mishima. Available from: http://www.ddbj.nig.ac.jp/ (accessed 19 August 2017) DeChaine, E.G., Anderson, S.A., McNew, J.M. & Wendling, B.M. (2013) On the evolutionary and biogeographic history of Saxifraga sect. Trachyphyllum (Gaud.) Koch (Saxifragaceae Juss.). PloS one 8: e69814. https://doi.org/10.1371/journal.pone.0069814 DeChaine, E.G. (2014) Introducing the Spotted Saxifrages: Saxifraga Sect. Bronchiales, Sect. Nov. (Saxifragaceae). Rhodora 116: 25– 40. https://doi.org/10.3119/13-09 Don, D. (1822) Monograph of the genus Saxifraga. Transactions of the Linnean Society of London 13: 341–452. https://doi.org/10.1111/j.1095-8339.1821.tb00066.x Doyle, J.J. & Doyle, J.L. (1987) A rapid DNA isolation procedure for small quantities of fresh leaf tissue. Phytochemistry Bulltein 19: 11–15. Engler, H.G.A. & Irmscher, E. (1919) Saxifragaceae-Saxifraga. In: Engler, H.G.A. (Ed.) Das Pflanzenreich 69 (IV. 117. II) Wilhelm Engelmann, Leipzig, pp. 449–709. Fedde, F. (1905) Novorum generum, specierum, varietatum, formarumque Siphonogamarum index. Just’s Botanischer Jahresbericht 33: 613. Felsenstein, J. (1985) Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783–791. https://doi.org/10.1111/j.1558-5646.1985.tb00420.x Fukuda, T., Naiki, A. & Nagamasu, H. (2007) Karyotypic analysis of Skimmia japonica (Rutaceae) and related species. Journal of Plant Research 120 (1): 113–121. https://doi.org/10.1007/s10265-006-0033-4 Fukuda, T., Andreeva, E., Taran, A., Takahashi, H. & Ikeda, H. (2016) Chromosome number of Micranthes nelsoniana (D. Don) Small (Saxifragaceae) in northeast Asia. Caryologia 69 (4): 325–329. Funamoto, T. & Nakamura, T. (1990) Notes on somatic chromosome numbers in Japanese Saxifraga (1). Chromosome Information Service 49: 4–6. Funamoto, T. & Nakamura, T. (1993) Notes on somatic chromosome numbers in Japanese Saxifraga (2), Karyotypes of four alpine Saxifraga species. La Kromosomo 70: 2412–2416. Goertzen, L.R. (1996) Genetic diversity and origin of two queen charlotte islands plants: Senecio newcombei and Saxifraga taylori. Master Thesis, University of British Columbia, Vancouver. 62 pp. Grisebach, A.H.R. (1843) Spicilegium florae rumelicae et bithynicae exhibens synopsin plantarum quas in aest., volume 1. Vieweg & Son, Brunsvigae, 407 pp. 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. Hara, H. (1937) Observations ad plantas asia orientalis XIII. Journal of Japanese Botany 13: 174. Hara, H. (1952) Contributions to the study of variations in the Japanese plants closely related to those of Europe or North America. Part 1. Journal of the Faculty of Science University of Tokyo, Section III (Botany) 6: 29–96. Hokkaido (2001) Hokkaido red list (plants). Hokkaido Government, Sapporo, 10 pp. Available from: http://www.pref.hokkaido.lg.jp/ks/ skn/grp/03/redlist1.pdf (accessed 20 August 2017) MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 65 Iwatsuki, K. & Kato, T. (1994) Red data plants—a pictorial of Japanese flora facing extinction. Takarajimasha, Tokyo, 208 pp. Karimi, K., Strucken, E., Moghaddar, N., Ferdosi, M., Esmailizadeh, A. & Gondro, C. (2016) Local and global patterns of admixture and population structure in Iranian native cattle. BMC Geneticsgenetics 17 (1): 108. https://doi.org/10.1186/s12863-016-0416-z Koch, W.D.J. (1836) Synopsis florae Germanicae et Helveticae. Wilmans, Frankfurt, 844 pp. Linnaeus, C. (1753) Species Plantarum 1. Laurentii Salvii, Holmiae, 560 pp. Meirmans, P.G., Van Tienderen, P.H. (2004) GENOTYPE and GENODIVE: two programs for the analysis of genetic diversity of asexual organisms. Molecular Ecology Resources 4 (4): 792–794. https://doi.org/10.1111/j.1471-8286.2004.00770.x Ministry of the Environment, Government of Japan (2017) Ministry of the Environment red list 2017. Tokyo, Ministry of the Environment, Government of Japan. Available from: http://www.env.go.jp/press/files/jp/105449.pdf (accessed 20 August 2017) Miyabe, K. & Kudô, Y. (1917) Materials for a flora of Hokkaido VII. Sapporo Natural History Society 6: 170–171. Murata, G. (1961) Saxifragaceae. In: Kitamura, S. & Murata, G. (Eds.) Colored illustrations of herbaceous plants of Japan (Choripetalae), Hoikusha, Osaka, pp. 136–157. Nakagawa, H. & Sato, K. (2015) A new Locality of Saxifraga nishidae (Saxifragaceae) from Hokkaido, Japan. Journal of Japanese Botany 90: 356–358. Nishikawa, Y., Miyaki, M. & Hori, S. (1992) Community formations of Saxifraga yuparensis and S. nishidae in the rocky ridge region of Mt. Yubari. Report of Hokkaido Institute of Environmental Sciences 19: 61–66. Nosaka, S. (1974) The Phanerogam Flora of Mt. Yûpari, Prov. Ishikari, Hokkaido, Japan. Journal of the Faculty of Science, Hokkaido University Series 5 Botany 9: 55–300. Nosaka, S. (1980) The plants of genus Saxifraga in Hokkaido. Hoppousansou 1: 2–10. Nylander, J.A.A. (2008) MrModeltest 2.3. Department of Systematic Zoology, Uppsala University. Uppsala. Ohba, H. (1999) Saxifragaceae. In: Satake, Y., Kitamura, S., Tominari, T., Ohwi, J. & Watari, S. (Eds.) Wild flowers of Japan II. Heibonsha, Tokyo, pp. 153–172. Okuyama, Y. (2016) Saxifragaceae. In: Ohashi, H., Kadota, Y., Kihara, H., Murata, J. & Yonekura, K. (Eds.) Wild flowers of Japan II revised new edition. Heibonsha, Tokyo, pp. 210–212. Rambaut, A., Drummond, A.J. & Suchard, M. (2013) Tracer v1.6: MCMC Trace Analysis Package. Institute of Evolutionary Biology, University of Edinburgh, UK. Available from: http://tree.bio.ed.ac.uk/software/tracer/ (accessed 20 August 2017) Rambaut, A. (2014) FigTree 1.4.2 software, Institute of Evolutionary Biology, University of Edinburgh. Available from: http://tree.bio. ed.ac.uk/software/figtree/ (accessed 20 August 2017) Ronquist, F., Teslenko, M., van der Mark, P., Ayres, D.L., Darling, A., Höhna, S., Larget, B., Liu, L., Suchard, M.A. & Huelsenbeck, J.P. (2012) MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Systematic Biology 61: 539–542. https://doi.org/10.1093/sysbio/sys029 Rutledge, L.Y., Devillard, S., Boone, J.Q., Hohenlohe, P.A. & White, B.N. (2015) RAD sequencing and genomic simulations resolve hybrid origins within North American Canis. Biology Letters 11 (7): 20150303. https://doi.org/10.1098/rsbl.2015.0303 Sakai, K. (1935) Studies on the chromosome number in alpine plants II. Genetics 9: 226–230. Sato, K. (2007) Geobotanical study on the alpine vegetation of Hokkaido, Japan. Hokkaido University Press, Sapporo, 688 pp. Shimizu, T. (1983) The new alpine flora of Japan in color II. Hoikusha, Tokyo, 331 pp. Shimizu, T., Kadota, Y. & Kihara, H. (2014) Alpine flowers of Japan. Yama-kei Publishers, Tokyo, 512 pp. Siplivinsky, V. (1971) Generis Saxifraga L. species asiaticae e sectione Trachyphyllum Gaud. Novitates Systematicae Plantarum Vascularium 8: 147–158. Small, J.K. (1905) Saxifragaceae. In: Britton, N.L. (Ed.) North American flora. New York Botanical Garden, New York, pp. 81–158. Soltis, D.E., Kuzoff, R.K., Mort, M.E., Zanis, M., Fishbein, M., Hufford, L., Koontz, J. & Arroyo, M.K. (2001) Elucidating deep-level phylogenetic relationships in Saxifragaceae using sequences for six chloroplastic and nuclear DNA regions. Annals of the Missouri Botanical Garden 88: 669–693. https://doi.org/10.2307/3298639 Sun, Y., Skinner, D.Z., Liang, G.H. & Hulbert, S.H. (1994) Phylogenetic analysis of Sorghum and related taxa using internal transcribed spacers of nuclear ribosomal DNA. Theoretical and Applied Genetics 89: 26–32. https://doi.org/10.1007/BF00226978 Suyama, Y., Matsuki, Y. (2015) MIG-seq: an effective PCR-based method for genome-wide single-nucleotide polymorphism genotyping using the next-generation sequencing platform. Scientific Reports 5: 16963. https://doi.org/10.1038/srep16963 66 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. Swofford, D.L. (2002) PAUP*: phylogenetic analysis using parsimony, version 4.0b10. Sinauer Associates, Sunderland. Available from: http://paup.sc.fsu.edu/ (accessed 20 August 2017) Taberlet, P., Gielly, L., Pautou, G. & Bouvet, J. (1991) Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Molecular Biology 17: 1105–1109. https://doi.org/10.1007/BF00037152 Takahashi, H. (2015a) Saxifraga yuparensis, Saxifraga nishidae. In: Yawara, T., Fujii, S., Ito, M. & Ebihara, A. (Eds.) Red data plants. Yama-kei Publishers, Tokyo, pp. 314. Takahashi, H. (2015b) Plants of the Kuril Islands. Hokkaido University Press, Sapporo, pp. 509. Takahashi, Y., Suyama, Y., Matsuki, Y., Funayama, R., Nakayama, K., & Kawata, M. (2016) Lack of genetic variation prevents adaptation at the geographic range margin in a damselfly. Molecular ecology 25: 4450–4460. https://doi.org/10.1111/mec.13782 Tkach, N., Röser, M., Miehe, G., Muellner-Riehl, A.N., Ebersbach, J., Favre, A. & Hoffmann, M.H. (2015) Molecular phylogenetics, morphology and a revised classification of the complex genus Saxifraga (Saxifragaceae). Taxon 64: 1159–1187. https://doi.org/10.12705/646.4 Thompson, J.D., Gibson, T.J., Plewniak, F., Jeanmougin, F. & Higgins, D.G. (1997) The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research 25: 4876–4882. https://doi.org/10.1093/nar/25.24.4876 Toyokuni, H. (1975) Index plantarum in regionibus alpinis Hokkaidoensibus sponte crescentium 2. Journal of the Asahikawa University 3: 158. Toyokuni, H. (1988) Alpine flowers of Japan. Yama-kei Publishers, Tokyo, 719 pp. Umezawa, S. (2004) Flora of the Yubari & Hidaka Range, Hokkaido. Hokkaido Shimbun Press, Sapporo, 239 pp. Wakabayashi, M. (2001) Saxifraga L. In: Iwatsuki, K. Boufford, David, E. & Ohba, H. (Eds.) Flora of Japan. Angiospermae, Dicotyledoneae, Archichlamydeae (b) vol IIb. Koudansha, Tokyo, pp. 54–55. Watanabe, T. (1971) Phytogeographical studies of the alpine plants (vascular plants) on the Hidaka-Yubari ranges, Hokkaido. Memories of the National Science Museum 4: 95–126. MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 67 Appendix Herbarium specimens used for morphological observations. Specimens that were observed for petal shape are marked with daggers (†), and ones of Saxifraga rebunshirensis that had slightly tricuspidate leaves are marked with asterisks (*). Superscript letters (A-F) indicate scientific names (synonyms) or Japanese names on the labels, as shown in the footnote. Saxifraga yuparensis: JAPAN. Hokkaido: Mt. Yubari Gamaiwa, 22-26 July 1928, Koidzumi Hideo 16910 (TNS). Mt. Yubari, July 1934, Ohwi Jizaburo 224482 (TNS). Mt. Yûpari Gama-iwa, 1 August 1967, Mikio Tohyama s.n. (SAPS). Mt. Yupari Gama-iwa,17 August 1973, Shiro Nosaka 031584, 031585 (SAPS). Hokkaido: Mt. Yubari Gamaiwa, 6 August 1991, Ken Sato 91.1708, 91.1328, 91.1329 (SAPS). Mt. Yubari Gama-iwa, 26 July 2016, Saya Tamura 582†, 586† (SAPT). Mt. Yubari Gama-iwa 19 July 2017, Saya Tamura 992†, 993†, 994†, 995† (SAPT). S. rebunshirensis: RUSSIA. Sakhalin: S. Saghallen Mt. Tosso, 23 July 1927, Hiratsuka N., Nagai M. & Iwadare S. s.n.A * (SAPS). Mt. Tosso, July 1927, Sugawara s.n.A (TNS). Kusu river, 18 July 1928, Okamae Tsutomu 252162F (TNS). Shirutoru,13 July 1929, Haga T. 13939A (SAPS). Mt. Tosso, 13 July 1929, Haga T. 13936A (SAPS). Shirutoru, 14 July 1929, Haga T. 13937A (SAPS). Mt. Tosso, July 1929, 20855F (TNS). Kasipo, 1 August 1929, 21329F*(TNS). Sakhalin, July 1930, Shigezou Sugawara 12105C* (SAPT). Shirutori Mt. Ishiyama, 20 July 1935, Shigezou Sugawara 12106C (SAPT). Hoe, 23 August 1935, Nakaji Masayoshi s.n.A (TNS). South of the Peninsula of Schmidt, 16 August 2001, Fukuda T. 2211E* (SAPS). Peninsula of Schmidt, 14 August 2001, Barkarov V. Y. 10812D* (SAPS). Eastern coast of the Peninsula of Schmidt, 14 August 2001, Fukuda T. 2117E (SAPS). South of the Peninsula of Schmidt, 16 August 2001, Fukuda T. 2198E, 2211E (SAPS). Tikhaya village, 28 August 2001, Fukuda T. 2483E (SAPS). Mt. Zhdanko, 2 August 2002, Hideki Takahashi 30172C (SAPS). Mt. Ostrovskaya, 3 August 2002, Noriyuki Fujii 01133C (SAPS). Mt. Changa, 7 August 2002, Hideki Takahashi 30264C, 30309C (SAPS). Central Sakhalin Mt. Chamga, 17 August 2006, Ken Sato, 06.0201 (SAPS). Central Sakhalin Mt. Vaida, 20 August 2006, Ken Sato 06.0316 (SAPS). Mt. Calcareous, 11 July 2007, Takahashi H., Kawahara T., Kitamura K. & Taran A., 33.763E (SAPS). RUSSIA. Kuril Islands: Shikotan, 3 July 1930, Tôno S. s.n.C (SAPT). JAPAN. Hokkaido: Mt. Rishiri, 6 July 1914, s.n.A (TNS). Rebun Island, 10 July 1914, Tamaki Seiichi 26106F (TNS). Mt. Ashibetsu Ishikari, 30 July 1915, Koidzumi Hideo 67488A (TNS). Kamifurano, 8 September 1915, Koidzumi Hideo 74161F*, 74162F (TNS). Mt. Kamifurano, 13 July 1917, Koidzumi Hideo 76906F (TNS). Mt. Rishiri, 5 August 1922, Hiratsuka Naohide s.n.A (TNS). Kamikawamura, 24-27 August 1926, Koidzumi Hideo 13157F (TNS). Mt. Chubetsuhondake, 19 August 1927, Koidzumi Hideo 15046F (TNS). Mt. Nipesotsu, 25-26 August 1927, Koidzumi Hideo 15339F (TNS). Mt. Yubari, 25 July 1928, Koidzumi Hideo 63705F (TNS). Mt. Rishiri, 23 July 1929, Shigeo Akiyama s.n.A (SAPS). Mt. Rishiri, July 1929, Akiyama S. s.n.A (SAPS). Shakotan, 3 July 1930, Tôno S. 010397C (SAPT). Mt. Rishiri, 13 August 1952, Furuse Miyoshi 130160F (TNS). Mt. Jozankei-Tengu, 3 July 1955, Toyokuni H. s.n.C (SAPS). Mt. Rishiri, 29 July 1956, Haginiwa Joju JH046164C (TNS). Mt. Taisetsu, 11 August 1956, Toyokuni H. 027749C (SAPS). Mt. Rishiri, 3 August 1957, Satake Yoshisuke & Ito Eiko 266166C, 273271F (TNS). Mt. Shari, 31 July 1959, Kawashiro Z. s.n.B (SAPS). Mt. Rausu Shiretoko, 11 July 1965, Haga S.& Saitou R. 308044B (TNS). Mt. Rishiri, 13 August 1965, 498 (TNS). Mt. Rhishiri, 13 August 1965, Andoh Yukio 347F (TNS). Mt. Rhishiri, 13 August 1965, Harajima Hisao 889F (TNS). Mt. Rishiri, 13 August 1965, Masuda Shinya 487F, 498F, 711B* (TNS). Mt. Zyozankeitengu,15 June 1969, Ken Sato 69.0101B (SAPS). Mt. Zyozankei-Tengu, 18 June 1970, Ken Sato 70.2802F, 70.2801F (SAPS). Mt.Zyozankei-Tengu, 24 June 1970, Watanabe S. W153F, 154F, 155F, 156F (SAPS). Mt. Rausu, 6 July 1970, Togashi M. 273107F (TNS). Mt. Rausu, 9 september 1971, Sato K. & Nishikawa T. 71.5062F (SAPS). Mt. Yûpari at Rôsoku-iwa, 6 August 1973, Shiro Nosaka 031581C*, 031582C* (SAPS). Mt. Rausu, 14 August 1979, Sato K. & Samejima J. 79.5258F (SAPS). Mt. Kumano, 31 August 1979, Ken Sato 79.3115F, 79.3093F (SAPS). Mt. Kumano, 17 July 1981, Ken Sato 81.0135F, 81.0136F (SAPS). Pankechiebu, 24 June 1982, Ken Sato 82.0517F, 82.0516F (SAPS). Mt. Ashibetsu, 10 July 1982, Ken Sato 82.0841C (SAPS). Mt. Kumano, 18 August 1982, Ken Sato 82.0337F (SAPS). Rebun Island, 13 July 1983, Ken Sato 83.0202F (SAPS). Mt. Kumano, 2 August 1983, Tsuchida T. 010396F (SAPT). Mt. Yubari Rosoku-iwa, 22 August 1985, Ken Sato 85.0207F (SAPS). Hoheikyo, 11 July 1986, Matsuji Hara 6470C (SAPS). Mt. Nishi-kumaneshiri-dake, 12 August 1987, Takahashi H., Kawabata K., Kushibiki E., Kikuzawa Y.& Kudo T. 7872 C, 7885C (SAPS). Mt. Tengu-dake, 24 July 1988, Takahashi H., Hara M., Sato K. & Umezawa S. 8548C (SAPS). Mt. Rausu, 7 July 1991, Takita K. 4766B (SAPT). Mt. Nipesotsu, 4 August 1991, Takita K. 5076B* (SAPT). Kamishihoro town Mt. Nishi-kumaneshiri, 27 July 1994, Ken Sato 94.0786, 94.0834 (SAPS). Mt. Rishiri Kutsugata-peak, 11 August 1994, Ken Sato 94.0990 (SAPS). Rebun Island Momo-iwa, 10 June 2016, Saya Tamura 328†, 329†(SAPT). Mt. NishiKumaneshiri cultivated in the Botanic Garden, Hokkaido University, 1 July 2016, Saya Tamura 434†, 436†(SAPT). JAPAN.Honshu: Mt. Shirouma, 26 August 1902, Yabe Yoshisada 24445F (TNS). Mt. Yatsugadake Koshu, 30 August 68 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL. 1902, Hayata B. 18763A (SAPS). Mt. Yatsugadake Shinshu, August 1902, Hayata B. 18761, 18762, 18764F (SAPS). Mt. Amida, 25 July 1903, Takeda H. s.n.A (SAPS). Mt. Togakushiura, 19 July 1904, Takeda H. s.n.A (SAPS). Mt. Yatsugadake Shinshu, 21 July 1904, Takeda H. s.n.A (SAPS). Mt. Shirouma Shinshu, 16 August 1905, Takeda F. s.n. A (TNS). Mt. Shirouma Shinshu, 11 August 1906, Miura Michiya s.n.A (TNS). Mt. Yatsugatake, August 1906, Sakurai Hanzaburo 6635F (TNS). Mt. Yatsugatake Shinano, 26 July 1908, Iishiba Eikichi s.n.C (TNS). Mt. Ohrenge Echigo, 10 August 1908, Nakamura Masao 94366F* (TNS). Mt. Shirouma Shinano, 15-16 August 1908, Shimura H. 138822F (TNS). Mt. Yatsugatake Shinshu, 13 August 1910, Koidzumi Hideo 67489A (TNS). Mt. Hariki Shinano, 25 July 1913, Yamakawa S. 264A (TNS). Mt. Shirouma, 4 August 1913, Yazawa Yonesaburo s.n.A* (TNS). Mt. Senjo Nagano, 9 August 1913, Takeda H. 3663F (TNS). Mt. Senjogatake Shinshu, 10 August 1916, Koidzumi Hideo 67492A (TNS). Mt. Yatsugatake Shinshu, August 1916, Koidzumi Hideo 67491F (TNS). Mt. Hotaka Hida, 1 September 1916, Koidzumi Hideo 26206F (TNS). Mt. Shirouma Shinano, 27 July 1917, Matsushima M. 37204C (TNS). Mt. Shirouma, 28 July 1920, Yazawa Yonesaburo s.n.A (TNS). Mt. Shirouma Shinano, 26-30 August 1920, Koidzumi Hideo 785F (TNS). Mt. Shirouma Shinano, 21 July 1921, Koidzumi Hideo 45051F (TNS). Mt. Kaisuruga Shinano, 27 July 1921, Koidzumi Hideo 2035F (TNS). Mt. Tsuba-Kuro Shinano, 22 July 1921, Koidzumi Hideo 1909F, 67367A (TNS). Mt. Senjo Shinano, 25 July 1922, Koidzumi Hideo 3516F (TNS). Mt. Senjo Shinano, 26 July 1922, Koidzumi Hideo 73067F (TNS). Mt. Shiomi Shinano, 28 July 1922, Koidzumi Hideo 3652F (TNS). Mt. Kogoochi Shinano, 29 July 1922, Koidzumi Hideo 3717F (TNS). Shimoina Shinano, 30 July 1922, Koidzumi Hideo 3725F (TNS). Mt. Arakawadake Shinshu, 31 July 1922, Koidzumi Hideo 67490F (TNS). Mt. Kitadake Koshu, 7 August 1922, Koidzumi Hideo 3882F (TNS). Mt. Yatsugadake Shinano, 28 August 1922, Koidzumi Hideo 4486F (TNS). Mt. Nokogiri Shinano, 29 July 1923, Koidzumi Hideo 5309F (TNS). Mt. Yatsugatake Nagano, 29 July 1923, Koidzumi Hideo 91244F (TNS). Mt. Arakawa Shinshu, 31 July 1923, Koidzumi Hideo 67494A (TNS). Kitaadzumigun Shinano, 6 August 1923, Koidzumi Hideo 48821F (TNS). Mt. Daitenjo Shinano, 6 August 1923, Koidzumi Hideo 64254F (TNS). Mt. Yari Shinano, 7 August 1923, Koidzumi Hideo 5087F, 64257F (TNS). Mt. Kitahotaka Shinano, 7 August 1923, Koidzumi Hideo 48387F (TNS). Mt. Ohtenjo Shinano, 8 August 1923, Koidzumi Hideo 48274F (TNS). Mt. Harinoki Shinano, 18 August 1923, Koidzumi Hideo 5734F (TNS). Mt. Yari Hida, 24 August 1923, Koidzumi Hideo 5954F (TNS). Mt. Yatsugatake Gongendake Shinano, 30 August 1923, Koidzumi Hideo 6005F (TNS). Mt. Ushikubi Shinano, August 1923, Koidzumi Hideo 5043F (TNS). Mt. Kurodake Kamishinkawagun, 1923, Koidzumi Hideo 5859F (TNS). Mt. Maekogouchidake Shinano, 1923, Koidzumi Hideo 44023F (TNS). Mt. Tsugamura Suruga, 22 July 1924, Koidzumi Hideo 496871F (TNS). Mt. Shiomi, July 1924, Yazawa Yonesaburo s.n. A* (TNS). Mt. Tsurugi Yetstu, 8 August 1924, Koidzumi Hideo 7343F (TNS). Mt. Tateyama Yetstu, 9 August 1924, Koidzumi Hideo 7890F (TNS). Mt. Akaishi range Kai-Suruga, 2-9 July 1925, Koidzumi Hideo 9562F, 61484F, 51496F (TNS). Mt. Shirouma Shinano, July 1925, Koidzumi Hideo 10929F (TNS). Mt. Akaishi range, 8 August 1925, Koidzumi Hideo 9563F (TNS). Mt. Senjo Kai, 29 August 1925, Ohwi Jizaburo 224189 (TNS). Mt. Tateyama Yetsiu, July 1926, Koidzumi Hideo, 12136F (TNS). Mt. Yatsugatake, 18 July 1927, s.n.A (TNS). Mt. Yatsugatake Shinano, 7 July 1928, Koidzumi Hideo 88652F (TNS). Mt. Yatsugadake Shinshu, 28 July 1928, Inagaki K. s.n.A (SAPS). Mt. Kisokoma Shinano, 4 August 1928, Koidzumi Hideo 1044F (TNS). Mt. Hotaka Shinano, 2 September 1929, Koidzumi Hideo 19029F, 19030F, 19032F*, 55500F (TNS). Mt. Norikura Shinano, 27 July 1930, Koidzumi Hideo & Yokouchi I. 87496F (TNS). Mt. Ioh Nagano, July 1930, Maekawa Fumio 85567B (TNS). Mt. Shiomi Shinano, 3 August 1930, Muramatsu Kouki 3453 (TNS). Mt. Kanenashi Shinano, 23 July 1931, Koidzumi Hideo 45346F (TNS). Mt. Arakawa Suruga, 7 August 1931, Hashimoto 44088F (TNS). Mt. Akaishi Suruga, 8 August 1931, Suzuki S. s.n.A (TNS). Sannomine Echizen, 22 July 1933, Tashiro Yoshitaro 42366F (TNS). Mt. Kisokoma, 22 August 1934, Koidzumi Hideo 105682F (TNS). Mt. Kisomae, 13 August 1935, Koidzumi Hideo 105787F (TNS). Mt. Shishidake Ecchu Tateyama, 28 July 1937, s.n.C (TNS). Mt. Senjo Kai, 1 August 1938, Okuyama Shunki 4192C (TNS). Mt. Shimizu Etchu, 2 August 1938, Hayashi, 69396F (TNS). Mt. Higashi Suruga, 1 August 1948, Hirano Hideo. 79776F (TNS). Mt. Yatsugatake Shinano-Kai, 27 July 1949, Haginiwa Joju JH009022C, JH009023C, JH046160A (TNS). Mt. Higashi Suruga, 2 August 1949, Hirano Hideo. 32F (TNS). Mt. Kamikouchi Shinano, 26 July 1951, Muramatsu Kouki 3052F (TNS). Mt. Higashi Shinano, 29 July 1951, Muramatsu Kouki 3324D (TNS). Mt. Senjo, August 1952, Iwatsuki Z. s.n.F (TNS). Mt. Yatsugatake Shinano, 12 July 1953, Haginiwa Joju JH046163C (TNS). Mt. Yatsugatake Shinano, 28 July 1953, Ueda Minoru 104289F (TNS). Mt. Akusawa Suruga, 23 August 1953, Haginiwa Joju JH009021A (TNS). Kisawa village, 3 August 1956, Muramatsu Kouki 2603D (TNS). Yamadera-oku Yamagata, 10 June 1957, Kato Nobuo 133441C (TNS). Mt. Akaishi Shinano, 4 August 1958, Ibaragi Yoichi 139447C (TNS). Mt. Harinoki Shinano, 30 July 1959, Okuyama Shunki & Utsumi H. 15618C (TNS). Mt. Harinoki Shinano, 30 July 1959, Okuyama Shunki 15675F (TNS). Sanpusi-touge Nagano, 4 August 1962, Kogure Kazuo 152954B (TNS). Mt. Kitadake Yamanashi, 21 August 1962, Ishizawa 8240F (TNS). Mt. Kitadake Kai, 6 August 1963, Nakamura Takehisa 594C (TNS). Mt. Kitadake Kai, 7 August 1963, Ohtake Iwao 387, 388B (TNS). Mt. Kitadake Kai, 7 August 1963, Masuda Shinya 2395, 2424, 2455, 2530B MOLECULAR AND CYTOLOGICAL EVIDENCES Phytotaxa 373 (1) © 2018 Magnolia Press • 69 (TNS). Mt. Higashi-Arakawa Shizuoka, 1 September 1963, Iwatsuki K.& Koyama K. 153E (TNS). Mt. Kogouchi Shinano, 16 August 1964, Masuda Shinya 3435B (TNS). Otomiyama-touge Echigo, 3 August 1965, Okuyama Shunki 24735C (TNS). Mt. Senjo Shinano, 20 August 1966, Masuda Shinya 5289B (TNS). Mt. Shiomi to Mt. Kita-arakawa in Minami Alps National Park, 1 August 1967, Konta F. & Takahashi H. 511B (TNS). Mt. Yatsu Shinano, 20 August 1968, Okuyama Shunki 12444C (TNS). Asahigoya-Mt. Yukikura Toyama, 29 August 1969, Satomi N.& Mochizuki R. 256064C (TNS). Sanpukutouge-Eboshi Nagano, 7 July 1971, Kojima Yu 340B (TNS). Mt. Amakazari, 29 July 1975, Nomura Kouhei 39383B (TNS). Mt. Amakazari Niigata, 14 August 1978, Sakuma Yoshihisa s.n.F (TNS). Mt. Amakazari Niigata, 12 August 1980, Sakuma Yoshihisa s.n.F (TNS). Mt. Tsurugi Toyama, 1 September 1980, Takai Yumi 949B (TNS). Mt. Senmai South alps, 31 August 1984, Konta Fumihiro. 15688B (TNS). Mt. Hakuba cultivated at Hakuba Goryu Alpine Botanical Garden, 16 June 2016, Saya Tamura 334†, 335†(SAPT). Mt. Kamikouchi Suruga, Koidzumi Hideo 6882F (TNS). Mt. Kitahotaka Shinshu, Koidzumi Hideo 67470F (TNS). S. nishidae: JAPAN. Hokkaido: Holotype (designated by Miyabe & Kudo 1917)—Mt. Yubari at Gamaiwa, 7 August 1916, Shozo Nishida 037669 (SAPS!). Mt. Yubari, 3 August 1921, Takeda & Tatewaki 016363, 016364, 016365 (SAPS). Mt. Yubari Gamaiwa, 22-26 July 1928, Koidzumi Hideo 16912, 16913 (TNS). Mt. Yubari, 24 July 1928, Koidzumi Hideo 16911 (TNS). Mt. Yubari, 1 August 1933, Sugimoto Junichi 45613 (TNS). Mt. Yupari Gama-iwa, 18 August 1957, Shiro Nosaka 031624 (SAPS). Mt. Yubari, 1 August 1987, Hideki Takahashi 016362 (SAPS). Mt. Yubari, 6 August 1991, Ken Sato 91.1327, 91.1707 (SAPS). Mt. Yubari, 22 July 1995, Takita K. 6312 (SAPT). Mt. Yubari, 3 August 2014, Ken Sato 14.200 (SAPS). Mt. Ashibetsu meotoiwa, 22 July 2014, Nakagawa 14019, 14020 (SAPS). Mt. Yubari Gama-iwa, 26 July 2016, Saya Tamura 589, 594, 610, 615, 620, 627 (SAPT). Mt. Yubari Gamaiwa 19 July 2017, Saya Tamura 996†, 997†, 998†, 999† (SAPT). A, Saxifraga bronchialis L.; B, Saxifraga bronchialis L. subsp. funstonii (Small) Hulten var. rebunshirensis Hara; C, Saxifraga cherlerioides D.Don var. rebunshirensis Hara; D, Saxifraga rebunshirensis (Engl. Et Irmsh.) Sipl.; E, Saxifraga sp.; F, Shikotanso (Japanese name) 70 • Phytotaxa 373 (1) © 2018 Magnolia Press TAMURA ET AL.