SRX23777384 SF5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602405 SF5 SF5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777385 SF6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602406 SF6 SF6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777386 SF7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602407 SF7 SF7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777387 SF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602408 SF8 SF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777388 ZF15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602409 ZF15 ZF15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777389 ZF2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602410 ZF2 ZF2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777390 ZF3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602411 ZF3 ZF3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777391 ZF4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602412 ZF4 ZF4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777392 CM21 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602413 CM21 CM21 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777393 ZF5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602414 ZF5 ZF5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777394 ZF6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602415 ZF6 ZF6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777395 ZF7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602416 ZF7 ZF7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777396 CM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602417 CM5 CM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777397 CM56 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602418 CM56 CM56 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777398 CM6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602419 CM6 CM6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777399 CM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602420 CM7 CM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777400 CF2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602421 CF2 CF2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777401 GF10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602422 GF10 GF10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777402 GF11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602423 GF11 GF11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777403 GF15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602424 GF15 GF15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777404 GM28 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602425 GM28 GM28 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777405 GM29 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602426 GM29 GM29 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777406 GM3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602427 GM3 GM3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777407 GM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602428 GM5 GM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777408 GM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602429 GM7 GM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777409 CF33 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602430 CF33 CF33 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777410 JF1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602431 JF1 JF1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777411 JF11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602432 JF11 JF11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777412 JM51 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602433 JM51 JM51 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777413 JM57 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602434 JM57 JM57 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777414 JM58 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602435 JM58 JM58 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777415 JM59 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602436 JM59 JM59 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777416 JM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602437 JM7 JM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777417 JM9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602438 JM9 JM9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777418 CF38 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602439 CF38 CF38 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777419 LF1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602440 LF1 LF1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777420 CF10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602441 CF10 CF10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777421 LM10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602442 LM10 LM10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777422 LM11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602443 LM11 LM11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777423 LM12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602444 LM12 LM12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777424 LM13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602445 LM13 LM13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777425 CF42 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602446 CF42 CF42 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777426 LM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602447 LM15 LM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777427 LM16 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602448 LM16 LM16 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777428 LM2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602449 LM2 LM2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777429 RM1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602450 RM1 RM1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777430 RM10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602451 RM10 RM10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777431 RM11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602452 RM11 RM11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777432 RM12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602453 RM12 RM12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777433 RM13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602454 RM13 RM13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777434 CF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602455 CF8 CF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777435 RM14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602456 RM14 RM14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777436 RM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602457 RM15 RM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777437 SF9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602458 SF9 SF9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777438 SM1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602459 SM1 SM1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777439 SM10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602460 SM10 SM10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777440 SM11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602461 SM11 SM11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777441 SM12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602462 SM12 SM12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777442 SM13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602463 SM13 SM13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777443 CM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602464 CM15 CM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777444 SM14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602465 SM14 SM14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777445 ZF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602466 ZF8 ZF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777446 ZF9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602467 ZF9 ZF9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777447 ZM1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602468 ZM1 ZM1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777448 ZM10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602469 ZM10 ZM10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777449 ZM11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602470 ZM11 ZM11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777450 ZM12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602471 ZM12 ZM12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777451 ZM13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602472 ZM13 ZM13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777452 CF17 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602473 CF17 CF17 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777453 GF16 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602474 GF16 GF16 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777454 GF18 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602475 GF18 GF18 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777455 GF19 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602476 GF19 GF19 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777456 GF2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602477 GF2 GF2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777457 GF20 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602478 GF20 GF20 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777458 GF21 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602479 GF21 GF21 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777459 GF22 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602480 GF22 GF22 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777460 CF23 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602481 CF23 CF23 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777461 JF16 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602482 JF16 JF16 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777462 JF17 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602483 JF17 JF17 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777463 JF19 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602484 JF19 JF19 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777464 JF22 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602485 JF22 JF22 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777465 JF31 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602486 JF31 JF31 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777466 JF47 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602487 JF47 JF47 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777467 JF52 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602488 JF52 JF52 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777468 JF60 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602489 JF60 JF60 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777469 LF10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602490 LF10 LF10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777470 LF11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602491 LF11 LF11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777471 LF12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602492 LF12 LF12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777472 LF13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602493 LF13 LF13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777473 LF14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602494 LF14 LF14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777474 LF15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602495 LF15 LF15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777475 LF2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602496 LF2 LF2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777476 LF3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602497 LF3 LF3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777477 CF14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602498 CF14 CF14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777478 CF4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602499 CF4 CF4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777479 LF5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602500 LF5 LF5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777480 LF6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602501 LF6 LF6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777481 LF7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602502 LF7 LF7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777482 LF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602503 LF8 LF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777483 LF9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602504 LF9 LF9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777484 LM1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602505 LM1 LM1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777485 RF1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602506 RF1 RF1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777486 RF10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602507 RF10 RF10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777487 RF11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602508 RF11 RF11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777488 RF12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602509 RF12 RF12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777489 RF13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602510 RF13 RF13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777490 RF14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602511 RF14 RF14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777491 RF15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602512 RF15 RF15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777492 RF2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602513 RF2 RF2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777493 CM1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602514 CM1 CM1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777494 SF1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602515 SF1 SF1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777495 SF10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602516 SF10 SF10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777496 SF11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602517 SF11 SF11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777497 SF12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602518 SF12 SF12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777498 SF13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602519 SF13 SF13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777499 SF14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602520 SF14 SF14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777500 SF15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602521 SF15 SF15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777501 SM9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602522 SM9 SM9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777502 CM18 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602523 CM18 CM18 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777503 ZF1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602524 ZF1 ZF1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777504 ZF10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602525 ZF10 ZF10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777505 ZF11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602526 ZF11 ZF11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777506 ZF12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602527 ZF12 ZF12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777507 ZF13 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602528 ZF13 ZF13 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777508 ZF14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602529 ZF14 ZF14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777509 ZM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602530 ZM7 ZM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777510 ZM8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602531 ZM8 ZM8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777511 ZM9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602532 ZM9 ZM9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777512 CM28 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602533 CM28 CM28 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777513 CM29 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602534 CM29 CM29 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777514 CM30 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602535 CM30 CM30 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777515 CM41 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602536 CM41 CM41 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777516 CM43 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602537 CM43 CM43 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777517 GM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602538 GM15 GM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777518 GM16 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604138 GM16 GM16 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777519 CF3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603191 CF3 CF3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777520 GM18 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604139 GM18 GM18 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777521 GM20 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602539 GM20 GM20 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777522 GM21 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603192 GM21 GM21 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777523 GM23 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602540 GM23 GM23 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777524 GM25 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602541 GM25 GM25 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777525 JM20 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604140 JM20 JM20 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777526 JM21 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602542 JM21 JM21 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777527 JM23 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602543 JM23 JM23 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777528 CF37 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602544 CF37 CF37 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777529 JM3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603193 JM3 JM3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777530 JM4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602545 JM4 JM4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777531 JM40 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602546 JM40 JM40 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777532 JM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602547 JM5 JM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777533 LM3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604633 LM3 LM3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777534 LM4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602550 LM4 LM4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777535 LM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602549 LM5 LM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777536 LM6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602548 LM6 LM6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777537 LM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603196 LM7 LM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777538 LM8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603195 LM8 LM8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777539 LM9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602551 LM9 LM9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777540 CF57 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602552 CF57 CF57 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777541 RM2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602553 RM2 RM2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777542 RM3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603194 RM3 RM3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777543 RM4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602554 RM4 RM4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777544 RM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602555 RM5 RM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777545 RM6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602556 RM6 RM6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777546 RM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604141 RM7 RM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777547 RM8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602557 RM8 RM8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777548 RM9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602558 RM9 RM9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777549 SM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602559 SM15 SM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777550 SM2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602560 SM2 SM2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777551 SM3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603198 SM3 SM3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777552 SM4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602561 SM4 SM4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777553 SM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602563 SM5 SM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777554 SM6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602562 SM6 SM6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777555 SM7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603197 SM7 SM7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777556 SM8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602564 SM8 SM8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777557 CM24 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604142 CM24 CM24 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777558 ZM14 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602565 ZM14 ZM14 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777559 ZM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602566 ZM15 ZM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777560 ZM2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602567 ZM2 ZM2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777561 ZM3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603199 ZM3 ZM3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777562 ZM4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602568 ZM4 ZM4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777563 ZM5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602569 ZM5 ZM5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777564 ZM6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602570 ZM6 ZM6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777565 GF24 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604453 GF24 GF24 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777566 GF4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602571 GF4 GF4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777567 GF5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602572 GF5 GF5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777568 GF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603201 GF8 GF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777569 GF9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602573 GF9 GF9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777570 GM1 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602574 GM1 GM1 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777571 GM10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602575 GM10 GM10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777572 GM11 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603202 GM11 GM11 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777573 CF36 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602577 CF36 CF36 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777574 JF61 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602576 JF61 JF61 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777575 JF62 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602578 JF62 JF62 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777576 JF63 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603200 JF63 JF63 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777577 JF65 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602579 JF65 JF65 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777578 JF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602580 JF8 JF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777579 JM10 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20604143 JM10 JM10 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777580 JM15 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602582 JM15 JM15 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777581 LF4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602581 LF4 LF4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777582 RF3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602583 RF3 RF3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777583 RF4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603204 RF4 RF4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777584 CF68 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602584 CF68 CF68 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777585 RF5 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602585 RF5 RF5 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777586 RF6 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602586 RF6 RF6 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777587 RF7 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20603203 RF7 RF7 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777588 RF8 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602587 RF8 RF8 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777589 RF9 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602588 RF9 RF9 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777590 SF2 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602589 SF2 SF2 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777591 SF3 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20617228 SF3 SF3 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777592 SF4 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602590 SF4 SF4 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000 SRX23777593 CM12 SRP492275 bp0 Genomic DNA of the samples were extracted from mature leaves using GeneMark Plant Genomic DNA Purification kit (GeneMark Technology Co., Ltd., Tainan, Taiwan). Quality and quantity of the extracted DNA were determined by Nanodrop UV spectrophotometer (Thermo) and stored in the 20C freezer.The adapters and primers were designed based on Shirasawa et al. (2016; 2020). Genomic DNA samples were digested with restriction enzyme PstI and MspI (FastDigest Restriction Enzymes, Thermo Scientific). AMPure XP beads (Beckman Coulter, Brea, CA, USA) were used for clean-up by removing the restriction enzyme. The adaptor ligation was conducted using the restriction enzyme digestion product, PstI adapter, MspI adapter, 10X T4 ligation buffer, and NxGen T4 DNA ligase (Lucigen Corporation) and incubated using a high precision temperature gradient reactor (TProfessional 96 Thermocycler, Biometra). Fragments size of 300800 bp of the ligated products were selected using AMPure XP beads (Beckman Coulter, Brea, CA, USA and followed by PCR amplification using primers with index and JMR hot start PCR mix kit (JMR Holdings, UK). The library quantities and qualities were evaluated using the Qubit dsDNA HS Assay Kit and qualitative analysis using the BiOptic Qsep400 and High Sensitivity (N1) Cartridge Kit. Libraries were sequenced using the Illumina Hiseq 4000 to perform 150 bp paired-end sequencing. The sequencing service was provided by Genomics BioSci & Tech Co Ltd. SRS20602591 CM12 CM12 RAD-Seq GENOMIC Restriction Digest Illumina HiSeq 4000