PhytoKeys 2I I: 33-44 (2022) A peer-rev iewed open-access journal 1] & doi: 10.3897/phytokeys.2 | 1.89452 46P h y toKe y https://phytokeys.pensoft.net Launched to accelerate biodiversity research Hydrangea marunoi (Hydrangeaceae), a new species from Osumi Peninsula, southern Japan Shuichiro Tagane', Shinji Fujii*, Shun K. Hirota’, Akiyo Naiki*, Tetsukazu Yahara® | The Kagoshima University Museum, Kagoshima University, 1-21-30 Korimoto, Kagoshima, 890-0065, Japan 2 Department of Field Ecology, University of Human Environments, Okazaki, Aichi, 444-3505, Japan 3 Field Science Center, Graduate School of Agricultural Science, Tohoku University, 232-3 Yomogida, Narukoonsen, Osaki, Miyagi, 989-6711, Japan 4 Iriomote Station, Tropical Biosphere Research Center, University of the Ryukyus, 870 Uehara, Taketomi-cho, Yaeyama-gun, Okinawa, 907-1541, Japan § Kyushu Open University, 744 Motooka, Fukuoka, 819-0395, Japan Corresponding author: Shuichiro Tagane (stagane29@gmail.com) Academic editor: Sandy Knapp | Received 24 June 2022 | Accepted 14 September 2022 | Published 11 October 2022 Citation: Tagane S, Fujii S, Hirota SK, Naiki A, Yahara T (2022) Hydrangea marunoi (Hydrangeaceae), a new species from Osumi Peninsula, southern Japan. PhytoKeys 211: 33-44. https://doi.org/10.3897/phytokeys.211.89452 Abstract Hydrangea marunoi Tagane & S. Fujii, from the Kimotsuki Mountains in the Ohsumi Peninsula, southern Japan, is described and illustrated. It is morphologically similar to H. alternifolia in having three-petaloid calyx lobes in marginal flowers, but is distinguished by the larger stamen number, and longer styles and seeds. Multiplex ISSR genotyping by sequencing (MIG-seq) demonstrated that the new species is mono- phyletic and closely related to H. amamiohsimensis and H. moellendorffii rather than H. alternifolia. Keywords Cardiandra, endemic species, flora, MIG-seq, plant taxonomy, Saxifragales Introduction Hydrangea L., including approximately 270 natural species (De Smet et al. 2015) and four widely cultivated species (Fulcher et al. 2016), is a well-known genus in Hydrangeaceae. Based on phylogenetic analysis, De Smet et al. (2015) proposed a broad circumscription of Hydrangea comprising approximately 200 species distributed across East and Southeast Asia and the Americas. Most Hydrangea species are shrubs or lianas. However, the species of Hydrangea L. sect. Cardiandra (Siebold & Zucc.) Y.De Smet Copyright Shuichiro Tagane et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 34 Shuichiro Tagane et al. / PhytoKeys 211: 33-44 (2022) & Samain are herbs that have been treated as members of the genus Cardiandra (Ohba 1985a, b, 2001; Wei and Bartholomew 2001; De Smet et al. 2015). In the current broad circumscription of Hydrangea, it is treated as a section of the genus Hydrangea, which additionally also includes eight groups previously known as genera i.e. Broussaisia Gaudich, Decumaria L., Deinanthe Maxim., Hydrangea s.str., and Pileostegia Hook. f. & Thomson, Platycrater Siebold & Zucc., and Schizophragma Siebold & Zucc. In the phylogeny by De Smet et al. (2015), Hydrangea sect. Cardiandra is monophyletic and a sister to the sect. Deinanthe, which comprises two known herbaceous species from China to Japan. In Hydrangea sect. Cardiandra, four species, H. alternifolia L., H. amamiohsimensis (Koidz.) Y. De Smet et Granados, H. moellendorffii Hance, and H. densifolia (C. E Wei) Y. De Smet & Granados are known from Japan, Taiwan and China (Ohba 1985a, b, 2001; Wei and Bartholomew 2001; Ohashi 2017), all but H. densifolia are known from Japan. During our floristic survey in Kagoshima Prefecture, southern Japan in 2021, we collected an unknown flowering species of Hydrangea sect Cardiandra. \t is similar to 7. alternifolia in appearance, but it differs from it in its habitat and some floral characters. To clarify the relationship between the unknown species and the other species of Hydrangea sect. Cardiandra in Japan, we examined the phylogenetic relationships of 52 samples of 1. alternifolia, H. amamiohsimensis, H. moellendorffii, and the unknown species using multiplex ISSR genotyping by sequencing (MIG- seq, Suyama and Matsuki 2015) and compared this with our observations of morphological characteristics. Multiplexed inter-simple sequence repeats (ISSR) genotyping by sequencing (MIG-seq) is a technique used to obtain many single nucleotide polymorphisms (SNPs) throughout a genome, which is valuable for determining molecular phylogenetic trees. It has been successfully applied to resolve the taxonomy of closely related taxa, including in Hydrangea (Hirota et al. 2022). Based on the phylogenetic hypotheses resulting from MIG-seq analysis and subsequent morphological observations, we describe Hydrangea marunoi, sp. nov. Materials and methods Morphological observation and assessment of conservation status To assess the novelty of the unknown species, we consulted the taxonomic literature (Ohba 1985a, b, 2001; Wei and Bartholomew 2001; Ohashi 2017) and herbarium specimens at FU, KAG, KAP, and TI, as well as the digitized specimen images of FKSE, TRPM, and those available at the Shimane Nature Museum of Mt. Sanbe avail- able on the web (Digital herbarium of Shimane University Faculty of Life and Envi- ronment Sciences http://tayousei.life.shimane-u.ac.jp/harbarium/). The conservation status was calculated following the IUCN Red List categories and criteria v3.1 (IUCN 2012) and IUCN guideline IUCN 2019). The Extent of Oc- currence (EOO) and Area of Occupancy (AOO) were calculated using the GeoCAT software (Bachman et al. 2011). Hydrangea marunoi (Hydrangeaceae), a new species from Japan 35 Taxon sampling for phylogenetic analysis To perform the phylogenetic analysis, 52 samples of the Hydrangea sect. Cardiandra were gathered from both our field surveys in Japan and herbarium specimens deposited at the Kagoshima University Museum (KAG): 10 1. marunoi samples, two H. amamiohsimensis samples, four H. moellendorffii samples, and 36 H. alternifolia samples (Suppl. material 1: Table S1; Fig. 1). Additionally, one sample of 1. bifida (Maxim) Y.De Smet & C.Granados of Hydrangea sect. Deinanthe was used as the outgroup. During the field survey, a small piece of leaf was cut, placed in a tea bag, and dried with silica gel in a zip-lock bag. — @® H. marunoi | ce’ m@ H. alternifolia (s.str.) eine Soo? ° OH. alternifolia (clade A) eh were @ H. amamiohsimensis A H. moellendorffii oe JPN10093 Os «> JPN10107 Shikoku Kyushu JPN1121 hs ) | (JPN10101, 10102) Mt. Takakuma (K1652-1654) a S £ iS © Q § Figure |. Collection localities of the four Hydrangea sect. Cardiandra species in this study. 36 Shuichiro Tagane et al. / PhytoKeys 211: 33-44 (2022) MIG-seq analysis Total DNA was extracted from dried leaves using the cetyl trimethylammonium bro- mide (CTAB) method (Doyle and Doyle 1990). For de novo SNP detection, MIG-seq (Suyama and Matsuki 2015) was performed according to the protocol described by Suyama et al. (2022). To prepare the MIG-seq library, a two-step PCR amplification process was performed: ISSR regions were amplified using the first PCR, and Illumina sequencing adaptors and indices were added to the first PCR products during the sec- ond PCR. Sequencing was performed on an Illumina MiSeq platform (Illumina, San Diego, CA, USA) using a MiSeq Reagent Kit v3 (150 cycles; Illumina). We skipped the sequencing of the first 17 bases of reads 1 and 2 (SSR primer regions and an- chors) using “DarkCycle.” Low-quality and extremely short reads containing adapter sequences were removed using Trimmomatic 0.39 (Bolger et al. 2014). The Stacks 2.60 pipeline (Rochette et al. 2019) was used for de novo SNP genotyping with the follow- ing parameters: the minimum depth of coverage required to create a stack (m) = 3, the maximum distance between stacks (M) = 2, and the maximum mismatches between loci when building the catalog (n) = 2. Three criteria were used for the SNP filtering. First, any SNP site where one of the two alleles had less than three counts was filtered out owing to the difficulty in distinguishing polymorphisms from sequencing errors when the minor allele count of SNPs is extremely low (Roesti et al. 2012). Second, SNPs with high heterozygosity (Ho 2 0.6) were removed because excess heterozygosity may have resulted from artifactual loci built from several paralogous genomic regions. Third, SNPs with a genotyping rate of < 30% were eliminated. Using the third crite- rion, the SNPs retained by at least 16 samples were included in the SNP dataset. Maximum likelihood phylogeny based on SNPs was inferred using the RAxML 8.2.10 software (Stamatakis 2014). We used a GIRCAT model with an ascertainment bias correction using the Lewis method and performed 1,000 replicates of parallelized tree search bootstrapping. Results Among the 17,753,114 raw reads (334,964 + 34,812 reads per sample) obtained, 13,254,044 reads (250,076 = 28,521 reads per sample) remained after quality control. After de novo SNP detection and filtering, 1875 loci and 4506 SNPs were identified. Hydrangea bifida (JPN4970) was removed from the SNP dataset because of its high proportion of missing data (0.982). The ten H. marunoi samples were monophyletic and formed two geographically defined groups; populations from Mt. Nokubi (K1658—1661) and that from the Oda River (JPN9950, JPN10103, K1633, K1637, K1638, KAG088891) (Fig. 2). Hydrangea marunoi was sister to a clade that included H. amamiohsimensis and H. moellendorffii. Hydrangea marunoi, H. amamiohsimensis, and H. moellendorffii were all supported as monophyletic by bootstrap values of 100%. Hydrangea alternifolia was sister to a clade that included these three species. Hydrangea marunoi (Hydrangeaceae), a new species from Japan 37 $s 22 Ah, g= | —~samples from Mt. Takakuma "% ae = s. & vy? Be Bi) & 5, = K “ea & te AF sf fe %, ro) &S es +