search for




 

Assessment of Genetic Relationship among Rhynchostylis Species based on Inter-Simple Sequence Repeat (ISSR) Markers
Plant Breed. Biotech. 2024;12:69-81
Published online July 17, 2024
© 2024 Korean Society of Breeding Science.

Juthaporn Saengprajak1*, Jirapa Phetsom1, Aphidech Sangdee1, Porntip Atichart1, Sununta Chuncher1, Piyada Theerakulpisut2, Arnusorn Saengprajak3, and Sudarat Thanonkaew4

1Department of Biology, Faculty of Science, Mahasarakham University, Kantarawichai District, Maha Sarakham 44150, Thailand
2Department of Biology, Faculty of Science, Khon Kaen University, Khon Kaen 40002, Thailand
3Department of Physics, Faculty of Science, Mahasarakham University, Kantarawichai District, Maha Sarakham 44150, Thailand
4Walairukkhawet Botanical Research Institute, Mahasarakham University, Kantarawichai District, Maha Sarakham 44150, Thailand
Corresponding author: Juthaporn Seangprajak
TEL. +66-8-6862-6808
E-mail. juthaporn.s@msu.ac.th
Received April 15, 2024; Revised June 7, 2024; Accepted June 25, 2024.
This is an Open-Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (http://creativecommons.org/licenses/by-nc/4.0) which permits unrestricted non-commercial use, distribution, and reproduction in any medium, provided the original work is properly cited.
Abstract
Abstract The genus Rhynchostylis contains important commercial orchids in Thailand with high diversity, but limited information is available regarding the genetic diversity of these orchids. Rhynchostylis has a short blooming period, which makes more difficult to distinguish between variations based only on their appearances. This study evaluated the genetic variation among six accessions of Rhynchostylis, along with one Vanda and one Aerides accession collected from different locations in the northeast of Thailand using 16 polymorphic inter-simple sequence repeat (ISSR) markers. The ISSR markers comprised 83 alleles with sizes ranging from 208 to 1,223 bp. The average number of alleles per locus was 5.19, with a standard deviation of 3.49. The average percentage of polymorphic bands was 95.21%. Polymorphism information content (PIC) values ranged between 0.05 and 0.44, with a mean average of 0.21. The calculated genetic similarity coefficients ranged from 0.383 to 0.914, suggesting a high level of genetic diversity among all the samples. UPGMA dendrograms were created using genetic similarity coefficients and divided into three main clusters. Cluster I contained four closely related R. gigantea orchids. Cluster II comprised two accessions, R. gigantea var. vivaphandhul and V. coelestis, while Cluster III contained two accessions of R. retusa and A. houlettiana. Based on ISSR data, the genetic similarities among the 8 orchid accessions do not correlate with flower color phenotypes and sampling locations, except for Cluster I. The results suggest that ISSR markers can effectively assess the genetic information of both wild and cultivated orchid resources. This study provides useful information for further development of novel markers specific to orchid varieties and for assisting the success of orchid breeding programs through the selection of parent plants.
Keywords : DNA fingerprinting, Endemic orchid, Epiphytic orchid, Genetic variation, ISSR, Rhynchostylis
Introduction

Rhynchostylis orchids, belonging to the Orchidaceae family, are tropical plants that produce stunning waxy flowers arranged in long and elegant inflorescences. These plants are native to Southeast Asia, Myanmar, and Thailand and grow on trees as epiphytes. Rhynchostylis orchids are characterized by their long, dark green leaves and a distinctive scent that is common across the genus, regardless of the blooming season of district species. There are about 1,300 species and 180 to 190 different genera of orchids in Thailand including major tropical orchids such as Vanda, Rhynchostylis, Ascocentrum, and Aerides. These orchids make significant contributions to the Thai orchid industry, particularly as cut flowers and potted plants (Thammasiri 2016). Currently, 54% of the orchids produced in Thailand are exported, while the remaining 46% are consumed in the domestic market (Thammasiri 2014).

Rhynchostylis, commonly referred to as the foxtail orchid (like Aerides), derives its generic name from the Greek rhyncho, beak and stylis, column describing the beaked column exhibited in this genus. The World Checklist of Vascular Plants (WCSP) currently recognizes four species as Rhynchostylis gigantea (Lindl.) Ridl., Rhynchostylis retusa (L.) Blume, Rhynchostylis rieferi (Higgins) and Rhynchostylis cymifera Yohannan, J. Mathew & Szlach. Two subspecies of Rhynchostylis gigantea; Rhynchostylis gigantea subsp. gigantea and Rhynchostylis gigantea subsp. violacea (Lindl.) Christenson are also recognized (Teoh 2022). Two species of this genus grow in Thailand, namely R. gigantea and R. retusa. The flower colors of R. gigantea range from white to white with red-purple spots, red-purple, and orange. This genus is popularly grown as a potted plant and used as a parent for producing interspecific hybrids (Thammasiri 2016).

In Thailand, Rhynchostylis orchids show large diversity. Diverse genetic characteristics are observed among closely related species, particularly in the distinctive floral forms of Rhynchostylis retusa. However, surveys on Rhynchostylis and closely related species encounter difficulty in distinguishing between specific species based on only the characteristics of their appearance (Sumaythachotphong et al. 2020).

Phenotyping, biochemistry, and DNA diversity are all utilized in assessing the diversity of genetic resources but methods for evaluating diversity using phenotype and biochemical characteristics are not always reliable due to environmental influences, labor intensity, and associated costs(McNally et al. 2009). Assessing genetic diversity through DNA analysis is the approach commonly utilized due to its repeatability, stability, and reliability. Several DNA-based molecular marker techniques have been developed to analyze genetic diversity including the analysis of restriction fragment length polymorphism (RFLP) using restriction enzymes, random amplified polymorphic DNA (RAPD), amplified fragment length polymorphism (AFLP), simple sequence repeats (SSR), inter-simple sequence repeats (ISSR), and single nucleotide polymorphism (SNP) detection via polymerase chain reaction (PCR) (Jeung et al. 2005; Navinder Saini et al. 2004). The reproducibility of RAPD and AFLP techniques presents limitations, and the markers they generate are dominant, leading to constraints in assessing relationships among genetic resources (Amiteye 2021). A single 16 to 25 bp primer is used in ISSR, targeting identical regions among microsatellites. These primers include a nucleotide anchoring motif that targets microsatellite regions comprising eight repeating dinucleotide units, six repeating trinucleotide units, or multiple repeating tetra or penta-nucleotide units (Jing-Yuan et al. 2018). Recently, ISSR markers have emerged as an alternative system with the reliability and advantages of microsatellites (SSRs) due to their high polymorphism, rapidity, simplicity, reproducibility, and cost-effectiveness (Pradeep Reddy et al. 2002). This technique has been widely used in studies on orchids for genetic fingerprinting (Pongsrila et al. 2014), species or variety identification (Sharma et al. 2013), genetic mapping (Gholami et al. 2021), genetic diversity (Lal et al. 2023), evolution and molecular ecology (Suetsugu et al. 2022), and hybridization analysis (Fajardo et al. 2014).

Thai orchids are good examples, showing how an ornamental crop can develop with different varieties in its natural habitat to become a major crop. The development of the Thai orchid industry has taken a long time, but it is now a high-income business and an important part of the national economy (Sumaythachotphong et al. 2020). This study focused on R. gigantea, which has grown wild for over a century in Wat Pa Mancha Khiri, Mancha Khiri District, Khon Kaen Province, Thailand. Over 4,000 Chang Kra wild orchids bloom in full splendor on the branches of about 280 tree species every year during January and February (Thongtham 2012). Rhynchostylis has a short blooming a period which makes more difficult to distinguish variations among varieties based only on their appearance. Being considered endangered inhabitants, Rhynchostylis species such as R. gigantea and R. retusa are reported as endemic in Southeast Asian countries, including Thailand, Laos PDR, and Myanmar. Among these species, R. gigantea is the most well-known, with the highest demand compared to other native orchid species (Anuttato et al. 2017). As such, these orchids have been the subject of conventional breeding projects aiming to establish the genetics of their economically desirable traits (Sumaythachotphong et al. 2020).

Eight orchid plants were selected for study: Rhynchostylis gigantea (Lindl.) Ridl., R. gigantea var. harrisonianum Holtt., R. gigantea (white with red-blotches flower), R. gigantea var. rubrum Sagarik, R. gigantea var. vivaphandhul, R. retusa (L.) Blume., Vanda coelestis (Rchb.f.) Motes and Aerides houlettiana Rchb. f. These species were chosen because they are endemic to Thailand, easily cultivated, widely available, and popular(Julsrigival et al. 2013). Notably, Rhynchostylis subjected to continuously decreasing genetic variability and thus listed in the Appendix-II of CITES (Jariyajirawattana et al. 2020).

In this study, we have focused on investigating the genetic relationship within the species R. gigantea, including wild varieties and commercially favored hybrid cultivars, and comparing them with closely related species V. coelestis and A. houlettiana, utilizing ISSR molecular markers. The results can serve as a basis for implementing strategies in the future aimed at conserving, sustainably utilizing, and hybridizing these orchids.

Materials and Methods

Plant materials

This study was carried out in the Laboratory of Molecular Genetics, Department of Biology, Faculty of Science, Mahasarakham University, Thailand. Six Rhynchostylis accessions, one Vanda genotype and one Aerides accession were surveyed and collected from various locations in Northeast Thailand, including both wild and cultivated areas (Table 1). The species R. gigantea was studied at Wat Pa Mancha Khiri, Mancha Khiri District, Khon Kaen Province, Thailand in its natural habitat. The orchids grow naturally on tamarind trees and identified plants were collected. The collection locations of R. gigantea and the other seven orchid plants, along with GPS coordinates and location altitudes presented in Fig. 1. Eight orchid samples were collected and grown in the greenhouse of the Department of Biology, Faculty of Science, Mahasarakham University. The visual appearance of the specimens as well as their corresponding names in the collection, general geographical distribution, and specific origin are presented in Table 1, Figs. 1 and 2.

Table 1 . Locations of Rhynchostylis sp., V. coelestis and A. houlettiana accessions collected from Northeast Thailand.

S. No.LocationGPS
Coordinates
AltitudeSpecies collected
1Wat Pa Mancha Khiri, Mancha Khiri District, Khon Kaen Province16°07'07.2"N 102°32'00.8"E163.16 mR. gigantea (Lindl.) Ridl.

2Agricultural Technology Park, Faculty of Agriculture, Khon Kaen University, Khon Kaen Province16°27'52.9"N 102°48'47.0"E176.11 mR. gigantea var. harrisonianum Holtt.

3Agricultural Technology Park, Faculty of Agriculture, Khon Kaen University, Khon Kaen Province16°27'52.9"N 102°48'47.0"E176.11 mR. gigantea (white witd red-blotches flower)

4Mancha Khiri District, Khon Kaen Province16°07'28.6"N 102°32'20.9"E163.16 mR. gigantea var. rubrum Sagarik

5Keung Subdistrict, Muang District, Maha Sarakham Province16°12'28.3"N 103°15'55.2"E143.98 mR. gigantea var. vivaphandhul

6Agricultural Technology Park, Faculty of Agriculture, Khon Kaen University, Khon Kaen Province16°27'52.9"N 102°48'47.0"E176.11 mR. retusa (L.) Blume

7Kut bak Subdistrict, Kut bak District, Sakon Nakhon Province17°04'18.5"N 103°49'04.9"E204.92 mV. coelestis (Rchb.f.) Motes

8Suranaree Subdistrict, Muang District, Nakhon Ratchasima Province14°52'13"N 102°00'22"E227.03 mAerides houlettiana Rchb. f.

*Species identification was based on appearance (Sumaythachotphong et al. 2020).


Figure 1. Sample collection sites. Here, 1: Wat Pa Mancha Khiri, Mancha Khiri, Khon Kaen Province; 2, 3, 6: Khon Kaen University, Khon Kaen Province; 4: Mancha Khiri, Khon Kaen Province; 5: Muang, Mahasarakham Province; 7: Kut bak, Sakon Nakhon Province; 8: Muang, Nakhon Ratchasima Province.
Figure 2. Flowers of R. gigantea (A), R. gigantea var. harrisonianum (B), R. gigantea (white with red-blotches flower) (C), R. gigantea var. rubrum (D), R. gigantea var. vivaphandhul (E), R. retusa (F), V. coelestis (G), and A. houlettiana (H).

DNA extraction

A small piece (2 cm×2 cm) of fresh leaf from each sample was used for DNA extraction and ground into a fine powder using liquid nitrogen. Genomic DNA was extracted from the homogenized mixture using a commercially available kit (PureDireX Genomic DNA Isolation Kit (Plants), Thermo Fisher Scientific Co., Ltd.), following the manufacturer's instructions meticulously.

Quantification of the extracted genomic DNA was performed by running the dissolved DNA in 1.0% agarose gel in 1X TBE (Tris-Borate-EDTA) buffer alongside uncut λ DNA of known concentration. To prepare for ISSR analysis, the genomic DNA was diluted to a concentration of 50 ng⋅μL-1 using nuclease-free H2O (Invitrogen™, Thermo Fisher Scientific Co., Ltd.), maintaining optimal conditions for the experiment.

ISSR-PCR fingerprinting

Sixteen ISSR primers (University of British Columbia (UBC), Biotechnology Laboratory, Vancouver, BC, Canada) were used for PCR amplification to study the polymorphism and banding patterns (Table 2). Polymerase chain reaction (PCR) was conducted in thermal cyclers (Biometra TAdvanced, Bio-Active Co., Ltd.). The ISSR-PCR reaction was prepared as a 20 µl reaction mixture containing 1 µl template DNA (50 ng⋅μL-1), 16.55 µl nuclease-free H2O, 0.8 µl MgCl2 (50 mM), 1X PCR buffer, 0.4 µl dNTPs (10 mM), 0.05 µl ISSR primer (100 pmol⋅μL-1) and 0.2 µl Taq DNA Polymerase (5 U⋅μL-1). The amplification program was 4 minutes preheating and initial denaturation at 94°C, then 40 cycles of 1 minutes at 94°C denaturation step, 1 minutes primer annealing at 51°C/52°C/53°C/54°C/55°C depending on the primer used, and 2 minutes extension at 72°C, with a final extension of 10 minutes at 72°C. The PCR products were kept at 4°C until gel electrophoresis. The amplified DNA fragments were then fractionated by electrophoresis on 1.0% (w/v) agarose gel with 1X TBE buffer. Electrophoresis was performed under a constant voltage of 100 V for 60-90 minutes. The size of each DNA band produced was estimated using a 100 bp DNA ladder RTU (Bio-Helix Co, Ltd). The gel electrophoresis results were visualized using Visafe Green Gel Stain (Vivantis®, Malaysia), and photographed under an ultraviolet transilluminator. Sequences and optimized annealing temperatures of the 16 ISSR primers used in this study are presented in Table 2.

Table 2 . List of ISSR (Inter-simple sequence repeat) primers used to assess genetic diversity in this study.

No.PrimerSequence (5’ - 3’)Annealingtemperature (°C)References
1PCP-1GAC GAC GAC GAC GAC55Pant et al. (2024)
2PCP-2 AGG AGG AGG AGG AGG AGG55Pant et al. (2024)
3PCP-3 GTG CGT GCG TGC GTG C55Pant et al. (2024)
4UBC 810 GAG AGA GAG AGA GAG AT53Tran et al. (2022).
5UBC 811GAG AGA GAG AGA GAG AC53Pathak et al. (2022)
6UBC 813 CTC TCT CTC TCT CTC TT53Fajardo et al. (2014)
7UBC 815CTC TCT CTC TCT CTC TG 53Itsuji et al. (2015)
8UBC 817CAC ACA CAC ACA CAC AA 52Itsuji et al. (2015)
9UBC 819 GTG TGT GTG TGT GTG TA 52Itsuji et al. (2015)
10UBC 820 GTG TGT GTG TGT GTG TC54Tikendra et al. (2021)
11UBC 827 ACA CAC ACA CAC ACA CG54Tikendra et al. (2021)
12UBC 834 AGA GAG AGA GAG AGA GT 52Shukla et al. (2017)
13UBC 835 AGA GAG AGA GAG AGA GYC52Pant et al. (2024)
Pant et al. (2024)
14UBC 848 CAC ACA CAC ACA CAC ARG52Pant et al. (2024)
15UBC 888 BDB CAC ACA CAC ACA CA51Pant et al. (2024)
16UBC 890CTC TCT CTC TCT CTV HV51Pant et al. (2024)

UBC = University of British Columbia primer

B, D, H, R, V, Y = IUB code



Data scoring and analysis

The DNA profiling data generated by 16 ISSR primers for the eight investigated genotypes, including six Rhynchostylis accessions, one Vanda genotype, and one Aerides accession, were used for statistical analysis. Distinctive, independent, and reproducible amplified bands/alleles were scored as present (1) or absent (0) across all orchid genotypes for each primer. The data were entered into a binary matrix and subsequently analyzed using the computer package NTSYS-pc (ver. 2.1) (Rohlf 1992).

The polymorphism information content (PIC) value was calculated for each ISSR locus using the formula developed by Roldan-Ruiz et al. (2000) PIC = 2fi(1- fi), where fi is the frequency of marker bands which were present and 1-fi is the frequency of marker bands which were absent. The dendrogram was created by the UPGMA (Unweighted Pairgroup Method with Arithmetical Averages) method in the SAHN (Sequential, Hierarchical, Agglomerative, and Nested Clustering) module of NTSYS-pc 2.1 software.

Results

Banding patterns and ISSR analysis

All plant samples were collected from various locations in the northeast of Thailand (Table 1). In this study, all extracted genomic DNA showed sufficiently good quality as PCR templates with ISSR primers. Sixteen ISSR primers examined displayed the appropriate banding pattern with DNA samples from six Rhynchostylis specimens, one Vanda accession, and one Aerides specimen, and were subsequently used for further analysis (Table 2). The sequences of ISSR primers were longer and bound to the conserved regions between SSR regions, meaning that this method had superior characteristics such as simplicity and high reproducibility. The ISSR method is more effective for distinguishing genotypes with high genetic similarity due to the high mutation rate commonly seen in ISSR loci (Verma et al. 2017).

The PCR amplification with 16 ISSR primers yielded 83 bands ranging from 208 to 1,223 bp across six Rhynchostylis accessions, one Vanda accession, and one Aerides specimen. Representative images of band profiles generated with two ISSR primer pairs in all orchid accessions are shown in Figs. 3A and 3B and Figs. 4A and 4B. Among the 16 ISSR primers, band numbers varied from 1 to 13, averaging 5.19 bands per primer, with a standard deviation of 3.49. Primer PCP-3 displayed the maximum 13 bands, while primer UBC 848 displayed the minimum 1 band, and 80 DNA fragments (96.39%) were polymorphic. The average percentage of polymorphic bands was recorded as 95.21%. The PIC values of these 16 primers varied from 0.05 (UBC 848) to 0.44 (PCP-3) (Table 3), with a mean average of 0.21.

Table 3 . Band variation and PIC values of the sixteen ISSR markers used in this study.

S.NoPrimerTotal number of bandsNumber of polymorphic bandsNumber of monomorphic bandsPIC valuePolymorphism percentageProduct size (bp)
1PCP-1111100.40100404-1155
2PCP-2 7700.29100540-1213
3PCP-3 131300.44100279-927
4UBC 810 4400.18100550-1150
5UBC 8115500.22100275-879
6UBC 813 4400.18100540-1180
7UBC 8155500.22100220-1223
8UBC 8173210.1466.67721-960
9UBC 819 2200.09100796-960
10UBC 820 2200.09100560-776
11UBC 827 5500.22100618-1091
12UBC 834 2200.09100355-678
13UBC 835 6600.25100425-1060
14UBC 848 1100.05100556-953
15UBC 888 10910.3790245-1126
16UBC 8903210.1466.67208-822

Total83803
Mean5.19±3.495.000.190.2195.21

Figure 3. DNA amplification results with ISSR primers UBC835 (A) and PCP-3 (B). Lane M represents 100 bp DNA ladder; Lane 1, R. gigantea; 2, R. gigantea var. harrisonianum; 3, R. gigantea (white with red-blotches flower); 4, R. gigantea var. rubrum; 5, R. gigantea var. vivaphandhul; 6, V. coelestis; 7, R. retusa and 8, A. houlettiana.
Figure 4. Band profile generated by ISSR primers UBC815 (A) and PCP-1 (B). Lane M represents 100 bp DNA ladder; Lane 1, R. gigantea; 2, R. gigantea var. harrisonianum; 3, R. gigantea (white with red-blotches flower); 4, R. gigantea var. rubrum; 5, R. gigantea var. vivaphandhul; 6, V. coelestis; 7, R. retusa and 8, A. houlettiana.

Genetic diversity

The results showed that the 16 ISSR primers generated 83 bands ranging from 208 to 1,223 bp. In the DNA fragments, the average polymorphic DNA bands were up to 95.21%. This was important information for the molecular study and showed that the primers provided more detailed information about the genetic makeup of the individuals being analyzed. The PIC value is an important parameter showing the effectiveness of specific primers in genomic characterization. In this study, the PIC values of the 16 ISSR primers varied from 0.05 (UBC 848) to 0.44 (PCP-3) (Table 3), with a mean average of 0.21. Therefore, all these primers were deemed suitable for genetic diversity study according to the classification of Botstein et al. (1980).

Genetic similarity and clustering analysis

A UPGMA dendrogram derived from the 16 ISSR data was generated to explore the genetic relationships among the eight orchid accessions. The calculated genetic similarity coefficients showed a large variation between samples from 0.383 to 0.914 (full data is presented in Table 4). The largest genetic similarity coefficient was 0.914 between accessions R. gigantea var. harrisonianum (Khon Kaen University, Khon Kaen Province) and R. gigantea var. rubrum (Mancha Khiri, Khon Kaen Province), while the lowest was 0.383 between accessions R. gigantea var. harrisonianum (Khon Kaen University, Khon Kaen Province) and R. retusa (Khon Kaen University, Khon Kaen Province). These results suggested that a high level of genetic diversity was exhibited among all the samples.

Table 4 . Genetic identity matrix among the six Rhynchostylis specimens, one Vanda accession, and one Aerides specimen based on the analysis of sixteen ISSR markers.

S1S2S3S4S5S6S7S8
S11.000
S20.5681.000
S30.7900.6791.000
S40.6300.9140.7161.000
S50.5800.5430.6170.5801.000
S60.5680.4810.5560.5190.7651.000
S70.5680.3830.5310.4200.5680.6051.000
S80.5560.3950.5430.4320.5800.5190.6671.000

Note: S1, R. gigantea; S2, R. gigantea var. harrisonianum; S3, R. gigantea (white with red-blotches flower); S4, R. gigantea var. rubrum; S5, R. gigantea var. vivaphandhul; S6, V. coelestis; S7, R. retusa and S8, A. houlettiana, respectively.



The result from the ISSR phylogenetic analysis was correlated with the morphological characteristics. The accessions were divided into three main clusters (I, II, and III) (Fig. 5). Cluster I comprised of four samples from Khon Kaen Province; R. gigantea, R. gigantea (white with red-blotches flower), R. gigantea var. harrisonianum, and R. gigantea var. rubrum. Cluster II included two accessions of Rhynchostylis and Vanda from Mahasarakham and Sakon Nakhon Provinces as R. gigantea var. vivaphandhul and V. coelestis, respectively while Cluster III comprised two samples including R. retusa and A. houlettiana from Khon Kaen and Nakhon Ratchasima Provinces, respectively. Accessions from the same location or nearby were grouped within the same cluster.

Figure 5. Unweighted pair-group method with arithmetic average clustering for six Rhynchostylis accessions, one Vanda and one Aerides accession based on sixteen ISSR markers.
Discussion

Sixteen ISSR markers were used to assess the genetic diversity of six Rhynchostylis accessions, one Vanda genotype and one Aerides specimen. Each marker distinguished only a few accessions but using a multiple marker set simultaneously allowed for the assessment of genetic diversity across various genetic resources. The ISSR primer sequences were longer and bound to conserved regions between the SSR regions, demonstrate superior characteristics in simplicity and high reproducibility. ISSR is more effective for distinguishing genotypes with high genetic similarity due to the high mutation rate commonly seen in ISSR loci (Verma et al. 2017).

In this study, ISSR primers generated 83 bands, with average polymorphic DNA bands up to 95.21%. Results showed that the primers provided more detailed information about the genetic makeup of the individuals being analyzed. The PIC value is an important parameter showing the effectiveness of specific primers in genomic characterization. In this study, the PIC values ranged from 0.05 to 0.44, with a mean average of 0.21. With the highest PIC value of 0.44, the primer PCP-3 was determined to be appropriate for genetic diversity research. According to the classification of Botstein et al. (1980), PIC values ≥ 0.5 are deemed to provide very high information, 0.5> PIC ≥ 0.25 indicate moderate information, and PIC <0.25 suggest little information. Therefore, the moderate PIC values observed in this study may be attributed to the requirement for additional marker types and/or a greater number of markers. Indeed, the PIC value of primer UBC 848 was considered to be extremely low, suggesting the need to increase the number of primers used or choose other DNA markers such as RAPD, SSR, or SNP to gain more genetic information.

ISSR data were used to calculate genetic similarity coefficients, with results indicating significant diversity among the samples, ranging from 0.383 to 0.914. The highest values were observed among the accessions of R. gigantea var. harrisonianum (pure-white flower) and R. gigantea var. rubrum (solid-red flower). The high genetic similarity between the distinct color variants of both R. gigantea suggests a close genetic relationship, likely due to recent divergence or hybridization (Bock et al. 2023). In contrast, the lowest value between R. gigantea var. harrisonianum (pure-white flower) and R. retusa (white with pink spots) indicates their significant genetic divergence. The results showed that genetic variations within R. gigantea do not correlate with the phenotypic flower color germplasm present.

UPGMA dendrograms of eight related orchids were created using genetic similarity coefficients and NTSYS-pc 2.1 software. The clustering pattern observed in the dendrogram suggests a close genetic relationship among the orchid specimens, with Cluster I, comprising only four accessions of R. gigantea, indicating a high genetic similarity among specimens, particularly those from Khon Kaen Province. The results support previous findings regarding the high genetic relationship within R. gigantea populations, likely due to their close geographic proximity and the potential for gene flow among individuals within the same region (Hamlin et al. 2020). The presence of both R. gigantea var. vivaphandhul and V. coelestis in Cluster II suggests a moderate level of genetic similarity (0.481) between these two accessions, despite them from different location. Similarly, Cluster III, comprised of R. retusa and A. houlettiana specimens from different provinces, suggests genetic divergence between these species, likely influenced by geographic isolation and genetic diversity (Andriamihaja et al. 2021).

An increasing number of studies have used DNA markers to investigate the genetic relationships between orchid species. Parab et al. (2008) investigated the level of genetic variation among populations of A. maculosum, an epiphytic orchid from Goa, India, using RAPD and ISSR markers. A report on Dendrobium orchids concluded that RAPD markers outperformed SSR markers in phylogenetic analysis, as they more accurately correlated with the descriptive morphological characteristics (Basavaraj et al. 2020). Tran et al. (2022) utilized a combination of 10 RAPD and 10 ISSR markers to evaluate the genetic diversity of 20 jewel orchids that did not correspond to the taxa previously studied. However, the RAPD marker has a significant drawback due to its lack of reproducibility and susceptibility to experimental conditions (Konzen et al. 2017).

ISSR markers have been utilized for assessing genetic diversity and other purposes. Pathak et al. (2022) studied the regeneration competence of an endangered orchid, Vanda cristata Wall. ex Lindl. using leaf explants with five ISSR primers. Results showed a high degree of genetic stability to obtain true-to-type plantlets. The protocol developed in this study may be useful for conserving V. cristata. Lal et al. (2023) analyzed the genetic variability among natural populations of two endemic orchids, A. multiflora and R. retusa in Uttarakhand, India using ISSR markers.

Sixteen ISSR markers were selected for this study based on the primers utilized in earlier studies. The PCR results showed high polymorphism and excellent repeatability of the allele pattern. The study of genetic relationships across eight orchid specimens categorized them based on their sampling locations and flower color phenotypes. Our study provided useful information of DNA profile to assess the genetic relationships among six accessions of Rhynchostylis, along with one Vanda and one Aerides accession. The DNA profiles were utilized to assess the genetic information of orchid germplasm and to preserve wild and cultivated resources. Based on this study, further analysis can lead to the development of novel markers specific to orchid varieties. The results can also be utilized to select parent plants and contribute to the success of orchid breeding programs.

Acknowledgments

This study was financially supported by Mahasarakham University, Thailand. The authors thank the Department of Biology, Faculty of Science, Mahasarakham University for the facilities provided.

References
  1. Amiteye S. 2021. Basic concepts and methodologies of DNA marker systems in plant molecular breeding. Heliyon 7(10): 1-20.
    Pubmed KoreaMed CrossRef
  2. Andriamihaja CF, Ramarosandratana AV, Grisoni M, Jeannoda VH, Besse P. 2021. Drivers of population divergence and species differentiation in a recent group of indigenous orchids (Vanilla spp.) in Madagascar. Ecol Evol 11(6): 2681-2700.
    Pubmed KoreaMed CrossRef
  3. Anuttato S, Boonruangrod R, Kongsamai B, Chanprame S. 2017. Morphological characterization of wild Rhynchostylis gigantea in Thailand. J ISSAAS 23(2): 20-32.
  4. Basavaraj B, Nagesha N, Jadeyegowda M. 2020. Molecular characterization of Dendrobium orchid species from western Ghat region of Karnataka using RAPD and SSR Markers. Intl J Curr Microbiol Appl Sci 9(1): 2157-2169.
    CrossRef
  5. Bock DG, Cai Z, Elphinstone C, Gonzalez-Segovia E, Hirabayashi K, Huang K, et al. 2023. Genomics of plant speciation. Plant Commun 4(5): 1-22.
    Pubmed KoreaMed CrossRef
  6. Botstein D, White RL, Skolnick M, Davis RW. 1980. Construction of a genetic linkage map in man using restriction fragment length polymorphisms. Am J Hum Genet 32(3): 314-331.
  7. Fajardo CG, de Almeida Vieira F, Molina WF. 2014. Interspecific genetic analysis of orchids in Brazil using molecular markers. Plant Syst Evol 300: 1825-1832.
    CrossRef
  8. Gholami S, Vafaee Y, Nazari F, Ghorbani A. 2021. Molecular characterization of endangered Iranian terrestrial orchids using ISSR markers and association with floral and tuber-related phenotypic traits. Physiol Mol Biol Plants 27: 53-68.
    Pubmed KoreaMed CrossRef
  9. Hamlin JAP, Hibbins MS, Moyle LC. 2020. Assessing biological factors affecting postspeciation introgression. Evol Lett 4(2): 137-154.
    Pubmed KoreaMed CrossRef
  10. Itsuji M, Takahashi C, Tatarenko IV, Kondo K. 2015. ISSR and ITS comparisons among Pogonia japonica, P. minor and P.'Miyamatokisou'nomen nudum. Chromosom Bot 10(3): 109-118.
    CrossRef
  11. Jariyajirawattana D, Bunnag S. 2020. The effects of vitrification method on the seed survival and genetic stability of Rhynchostylis gigantea (Rindl.) Lidl. Asia Pac J Sci Technol 25(03): 1-9.
  12. Jeung J, Hwang H, Moon H, Jena K. 2005. Fingerprinting temperate japonica and tropical indica rice genotypes by comparative analysis of DNA markers. Euphytica 146: 239-251.
    CrossRef
  13. Jing-Yuan X, Yan Z, Ze Y, Gang W, Guo-Yong X, Min-Jian Q. 2018. Molecular diversity analysis of Tetradium ruticarpum (WuZhuYu) in China based on inter-primer binding site (iPBS) markers and inter-simple sequence repeat (ISSR) markers. Chin J Nat Med 16(1): 1-9.
    Pubmed CrossRef
  14. Julsrigival J, Songsak T, Kirdmanee C, Chansakaow S. 2013. Determination of volatile constituents of Thai fragrant orchids by gas chromatography-mass spectrometry with solid-phase microextraction. J Nat Sci 12(1): 43-57.
    CrossRef
  15. Konzen ER, Peron R, Ito MA, Brondani GE, TSA S. 2017. Molecular identification of bamboo genera and species based on RAPD-RFLP markers.
    CrossRef
  16. Lal A, Pant M, Pant G, Palni L, Kumar A, Kumar G. 2023. ISSR marker assisted genetic diversity assessment in natural populations of two endemic orchids Aerides multiflora and Rhynchostylis retusa from Uttarakhand, India. S Afr J Bot 157: 151-158.
    CrossRef
  17. McNally KL, Childs KL, Bohnert R, Davidson RM, Zhao K, Ulat VJ, et al. 2009. Genomewide SNP variation reveals relationships among landraces and modern varieties of rice. Proc Natl Acad Sci U S A 106(30): 12273-12278.
    Pubmed KoreaMed CrossRef
  18. Navinder Saini NS, Neelu Jain NJ, Sunita Jain SJ, Jain R. 2004. Assessment of genetic diversity within and among Basmati and non-Basmati rice varieties using AFLP, ISSR and SSR markers. Euphytica 140: 133-146.
    CrossRef
  19. Pant M, Pant K, Chauhan A, Negi AS, Nainwal P, Srivastava A, et al. 2024. Molecular marker-assisted genetic diversity analysis in soybean cultivars from Himalayan region of Uttarakhand, India. Ecol Genet Genom 31: 100237.
    CrossRef
  20. Parab G, Krishnan S. 2008. Assessment of genetic variation among populations of Rhynchostylis retusa, an epiphytic orchid from Goa, India using ISSR and RAPD markers. Indian J Biotechnol 7: 133-139.
  21. Pathak P, Kumari A, Thakur B. 2022. Regeneration competence of an endangered orchid, Vanda cristata Wall. ex Lindl. using leaf explants: A study in vitro. S Afr J Bot 151: 1018-1024.
    CrossRef
  22. Pongsrila P, Chomnawang P, Srisamoot N, Wiriyaampaiwong P. 2014. DNA finger printing analysis of fourteen species of orchids using ISSR. KKU Res J 19: 210-216.
  23. Pradeep Reddy M, Sarla N, Siddiq E. 2002. Inter simple sequence repeat (ISSR) polymorphism and its application in plant breeding. Euphytica 128: 9-17.
    CrossRef
  24. Rohlf FJ. 1992. NTSYS-pc: numerical taxonomy and multivariate analysis system: Applied Biostatistics.
  25. Roldan-Ruiz I, Dendauw J, Van Bockstaele E, Depicker A, De Loose M. 2000. AFLP markers reveal high polymorphic rates in ryegrasses (Lolium spp.). Mol Breed 6: 125-134.
    CrossRef
  26. Sharma SK, Kumaria S, Tandon P, Rao SR. 2013. Assessment of genetic variation and identification of species-specific ISSR markers in five species of Cymbidium (Orchidaceae). J Plant Biochem Biotechnol 22: 250-255.
    CrossRef
  27. Shukla SP, Sharma A. 2017. In vitro seed germination, proliferation, and ISSR marker-based clonal fidelity analysis of Shorea tumbuggaia Roxb.: An endangered and high trade medicinal tree of Eastern Ghats. In Vitro Cell Dev Biol Plant 53: 200-208.
    CrossRef
  28. Suetsugu K, Hirota SK, Nakato N, Suyama Y, Serizawa S. 2022. Morphological, ecological, and molecular phylogenetic approaches reveal species boundaries and evolutionary history of Goodyeracrassifolia (Orchidaceae, Orchidoideae) and its closely related taxa. PhytoKeys 212: 111-134.
    Pubmed KoreaMed CrossRef
  29. Sumaythachotphong N, Bundithya W, Potapohn N. 2020. Flower color inheritance of Rhynchostylis gigantea (Lindl.) Ridl. Thai J Agric Sci 53(4): 178-191.
  30. Teoh ES. 2022. Orchid Species from Himalaya and Southeast Asia Vol. 3 (R-Z). Springer.
    CrossRef
  31. Thammasiri K. 2014. Presented at II International Orchid Symposium.
    CrossRef
  32. Thammasiri K. 2016. Thai orchid genetic resources and their improvement. Horticulturae 2(3): 9.
    CrossRef
  33. Thongtham N. 2012. In "Bangkok Post". Tailand.
  34. Tikendra L, Potshangbam AM, Dey A, Devi TR, Sahoo MR, Nongdam P. 2021. RAPD, ISSR, and SCoT markers based genetic stability assessment of micropropagated Dendrobium fimbriatum Lindl. var. oculatum Hk. f.-an important endangered orchid. Physiol Mol Biol Plants 27: 341-357.
    Pubmed KoreaMed CrossRef
  35. Tran TKP, Pham MH, Trinh TH, Widiarsih S. 2022. Investigation of the genetic diversity of jewel orchid in Vietnam using RAPD and ISSR markers. Biodiversitas 23(9): 4816-4825.
    CrossRef
  36. Verma KS, ul Haq S, Kachhwaha S, Kothari S. 2017. RAPD and ISSR marker assessment of genetic diversity in Citrullus colocynthis (L.) Schrad: a unique source of germplasm highly adapted to drought and high-temperature stress. 3 Biotech 7: 1-24.
    Pubmed KoreaMed CrossRef


March 2024, 12
Full Text(PDF) Free

Cited By Articles
  • CrossRef (0)

Social Network Service
Services
  • Science Central