Introduction {#s1}
============

Genetic variability in natural plant populations holds the potential to deal with multiple biotic and abiotic stresses. The potential to select a superior line increases with genetic diversity, the discovery of which becomes an important tool in plant breeding. On depletion of genetic variability, plants are unable to cope with unfavourable environmental conditions or pathogens and pests. Diversity studies also facilitate the conservation and management aims of a particular plant species. For the effective use of genetic diversity in plant breeding, knowledge of its extent and distribution plays a crucial role. Considering its significance, a large number of studies have been performed to estimate genetic diversity employing various methodologies in multiple plant species. Assessment of genetic variability employing molecular markers has proved to be a keystone to understanding the genomic constitution, categorizing the genes responsible for important traits, the classification and conservation of genetic variation in plant germplasm and developing selective proliferation approaches for plant propagation.

Wheat production in developing countries moved from defective to surplus ([@PLV083C50]) during the Green Revolution ([@PLV083C26]). Being a good source of carbohydrate, protein, sugar, fat, fibre and minerals, it provides half of the energy requirements of the human population ([@PLV083C60]; [@PLV083C58]; [@PLV083C70]). A constantly rising population demands an increment in wheat production ([@PLV083C19]; [@PLV083C21]; [@PLV083C69]; [@PLV083C54]). As the crop already covers a wide agricultural area, there is a negligible possibility of area expansion ([@PLV083C77]; [@PLV083C10]; [@PLV083C12]). Wheat faces multiple demands including its growth under warmer conditions ([@PLV083C73]; [@PLV083C76]), fighting various diseases ([@PLV083C62]), reduced energy input for sustainable growth ([@PLV083C78]) and high nutritional quality ([@PLV083C58], [@PLV083C59]). Considering this situation, [@PLV083C36] has suggested a need for a 'second green revolution'. This second green revolution must place emphasis on the utilization of inherent resources and the thorough understanding of genetic diversity.

During the course of evolution, wheat gained sufficient genetic diversity along the road from einkorn to bread wheat. Today, however, its diversity is weakening due to repeated cultivation of landraces for specific characters, narrow adaptation, farmers\' varietal selection and the requirement of uniform varieties in industrial seed grain processing ([@PLV083C7]; [@PLV083C61]; [@PLV083C28]). Implementation of high-yielding commercial varieties played an important part in loss of genetic variation. This depletion has now encouraged the use of genetic resources in wheat breeding programmes.

Genetic diversity is crucial for adaptability and survival of wheat species against the threat of disease attack/onset ([@PLV083C25]). If all the individuals of a population are identical, they will behave similarly to a stress condition and potentially be equally unable to cope with the situation. Hence, it is beneficial to assess the genetic diversity at a particular level that may facilitate the efficient exploitation of the germplasm. Furthermore, in addition to the fact that genetic diversity plays a part in the development of high-yielding bread wheat varieties, issues like the spread of coeliac disease necessitate the development of new genetic variants of tetraploid wheat ([@PLV083C72]; [@PLV083C71]).

Polyploidy and genome evolution of wheat are also partially responsible for maintaining its genetic diversity. In a review, [@PLV083C75] shed light on several aspects of the genome duplication and divergence leading to the development of evolutionary genetic diversity. Polyploidy resulting from hybridization leads to gene duplication across the entire genome and thus underlies the emergence of genetic variation. The agriculturally important phenomenon of hybrid vigour in polyploids is a consequence of genetic variability. As wheat is a polyploid species, it is beneficial to include tetraploid and hexaploid varieties in genetic variability assessment programmes. Such assessment programmes are imperative for managing populations by identifying the breeding genotypes. For a long time, depiction of diversity was dependent on morphological characterization ([@PLV083C66]; [@PLV083C39]; [@PLV083C65]). But due to the influence of environmental conditions and changes during developmental stages, morphological traits are considered unreliable for diversity estimation, mainly for closely associated populations.

Momentous progress in molecular genetics benefitted our understanding of the wheat genome and provided approaches for breeding. With the expansion of novel technologies like molecular markers, researchers utilized a range of *Triticeae* species for genotypic identification ([@PLV083C34]). Molecular marker techniques vary from each other in data generation efficiency and the genome area covered in the study. Selection of the type of marker tool for a study depends on the target crop and the issue. For example, random amplified polymorphic DNA (RAPD) markers are known for their simplicity, cost efficiency, fast polymorphism assessment, no prior information of DNA sequences being required and extensive coverage of the intact genome being possible. However, due to low reproducibility of the RAPD system, expense of the amplified fragment length polymorphism (AFLP) marker system and requirement of prior information about DNA sequences in SSR analysis, another dominant marker system, inter-simple sequence repeats (ISSRs), was included in the study. Due to high annealing temperature and extended sequence in comparison to RAPD markers, ISSR primers can produce more reproducible and reliable band patterns. Inter-simple sequence repeat markers are employed for distinguishing DNA on the basis of single base variation or insertions and deletions, and are equivalent to the SSR system in reproducibility. These markers are widely implemented for DNA fingerprinting, identification of species association, genetic variability studies and for recognizing the geographic origin of different plant species along with their ploidy status ([@PLV083C74]; [@PLV083C32], [@PLV083C33]; [@PLV083C5]; [@PLV083C42]; [@PLV083C63]; [@PLV083C22]; [@PLV083C46]; [@PLV083C49]; [@PLV083C8]; [@PLV083C41]; [@PLV083C53]; [@PLV083C37]; [@PLV083C39]; [@PLV083C9]; [@PLV083C40]; [@PLV083C11]; [@PLV083C43]; [@PLV083C56]; [@PLV083C31]). Hence, in the present study, RAPD and ISSR were chosen among the various marker systems to yield the benefits of both the techniques, diminishing their drawbacks and increasing the credibility of our results.

India and Turkey play crucial roles in supporting food security through wheat production. India holds first and second place in wheat growing area and production, respectively. It has become a priority to replace the uniform high-yielding varieties spread during the Green Revolution with diverse high-quality varieties. Turkey is found to be the place of origin of both tetraploid and hexaploid wheat domestication ([@PLV083C29]; [@PLV083C44]; [@PLV083C17]) and India is known to be the centre of origin of some promising varieties. An assessment of genetic variability and association of tetraploid and hexaploid wheat varieties from the two developing countries would be of immense benefit to wheat improvement programmes.

Association and contrast among the wheat cultivars from different countries can provide a useful overview on the evolutionary record of the genotypes and, hence, can facilitate the reach of breeding improvement. Although a number of genetic similarity studies were conducted on diverse wheat germplasm using ISSR and RAPD marker systems, phylogenetic association of Indian and Turkish *Triticum* species has not been documented to date. The present study represents the first effort in this direction, its objectives being to gain a better understanding of the genetic association and population structure of Indian and Turkish wheat on the basis of both geographical origin and ploidy. The share of the genetic variations within and among populations was also revealed so that the information provided can be effectively used by scientists for the development of genetically diverse, promising and healthier wheat varieties.

Methods {#s2}
=======

Study materials {#s2a}
---------------

The object of the present diversity study was a collection of 95 Indian and Turkish wheat genotypes including tetraploid (*Triticum turgidum* ssp. *durum*) and hexaploid (*Triticum aestivum* L.) wheat cultivars (Table [1](#PLV083TB1){ref-type="table"}). Well-known varieties were chosen for the experiment to facilitate the use of results in future breeding programmes. Table 1.Name and ploidy of 95 Indian and Turkish wheat genotypes used in the study.Sl. no.Name of genotypeGenotype numberPloidyOrigin130_KR-8G16XIndia2AAI_2G26XIndia3AKAW_4006G36XIndia4AKDW_2997G44XIndia5DBW_14G56XIndia6DBW_39G66XIndia7DDK_1025G76XIndia8DT_132G84XIndia9GW_03-12G96XIndia10GW_03-2G106XIndia11GW_03-3G116XIndia12GW_03-4G126XIndia13GW_03-9G136XIndia14HD_2177G146XIndia15HD_2236G156XIndia16HD_2270G166XIndia17HD_2307G176XIndia18HD_2329G186XIndia19HD_2380G196XIndia20HD_2402G206XIndia21HD_2501G216XIndia22HD_2643G226XIndia23HD_2881G236XIndia24HUW_12G246XIndia25HUW_251G256XIndia26HUW_37G266XIndia27HUW_468G276XIndia28HUW_533G286XIndia29HUW_55G296XIndia30K_01006G306XIndia31K_0204G316XIndia32K_616G326XIndia33K_8020G336XIndia34K_86G346XIndia35K_88G356XIndia36K_911G366XIndia37KALYANSONAG376XIndia38KBD_65G384XIndia39KBD_821G394XIndia40KBD_921G404XIndia41KBD_922G414XIndia42KBD_925G424XIndia43KBD_9452G434XIndia44KBD_9915G444XIndia45KD_9851G454XIndia46KLP_306G466XIndia47KLP_307G476XIndia48KLPD_1106G484XIndia49NAW_1448G496XIndia50NIDW_295G504XIndia51NW_1076G516XIndia52NW_2036G526XIndia53PBW_550G536XIndia54RAJ_1482G546XIndia55RAJ_1555G554XIndia56RAJ_3072G566XIndia57RAJ_3077G576XIndia58RAJ_3777G586XIndia59RAJ_4027G596XIndia60RAJ_4037G606XIndia61RAJ_4120G616XIndia62RAJ_6560G624XIndia63RD_1008G634XIndia64RD_1063G644XIndia65RD_1093G654XIndia66RD_1097G664XIndia67SAW_327G676XIndia68SAW_337G686XIndia69SAW_94G696XIndia70SONALIKAG706XIndia71UP_2338G716XIndia72UP_2511G726XIndia73UP_2525G736XIndia74UP_2696G746XIndia75VEERIG756XIndia76VL_832G766XIndia77WR_1381G776XIndia78WR_1408G786XIndia79WR_1421G796XIndia80BAYRAKTAR 2000G806XTurkey81SEVALG816XTurkey82KENANBEYG826XTurkey83BEZOSTAJA 1G836XTurkey84GÜN_91G846XTurkey85KONYA_2002G856XTurkey86AKBUĞDAYG866XTurkey87GEREK_79G876XTurkey88KIRAÇ_66G886XTurkey89ESERG896XTurkey90SÖNMEZ 2001G906XTurkey91HARMANKAYA 99G916XTurkey92KINACI 97G926XTurkey93YÜREĞIR 89G936XTurkey94ALTAY 2000G946XTurkey95LÜTFIBEYG956XTurkey

Plant genomic DNA extraction {#s2b}
----------------------------

Two to three weeks grown seedlings were utilized for total wheat DNA extraction following the cetyltrimethylammonium bromide (CTAB) method ([@PLV083C16]) with some modifications. Initially, cells were disrupted and purified with 2 % CTAB buffer and 10 μL RNase A, respectively, followed by incubation at 65 °C. This was followed by protein extraction employing phenol : chloroform : isoamyl alcohol and finally the CTAB--DNA complex was precipitated with isopropanol. The DNA pellet was twice washed with 70 % ethanol, dried and ultimately, dissolved in 100 μL DNase--RNase-free water. Purified DNA quantity and quality were verified using spectrophotometry and 1 % agarose gel electrophoresis, respectively. The DNA samples were diluted to a concentration of 50 ng μL^−1^ as templates for polymerase chain reactions (PCRs).

Inter-simple sequence repeats analysis {#s2c}
--------------------------------------

Twenty-seven ISSR primers (Metabion) were examined for distinguishing the polymorphism patterns, and among those 10 primers showed positive outcomes (Table [2](#PLV083TB2){ref-type="table"}) against chosen wheat varieties. Every PCR mixture of 25 μL contained 2.5 μL of 10× *Taq* buffer containing ammonium sulfate (except ISSR F3 where KCl was used), 3 μL of 25 mM MgCl~2~, 0.4 μL of 25 mM dNTP, 0.5 μL of 10 μM primer, 1.5 units of *Taq* DNA Polymerase and 100 ng of template DNA. The two-step ISSR--PCR reactions were performed in a Eppendorf Master Cycler. The physical reaction conditions and the number of initial and final PCR cycles were optimized for each individual ISSR primer. Table 2.Characteristics and polymorphism revealed by ISSR primers for 95 wheat genotypes used in the study.ISSR primerSequenceMelting temperature (*T*~m~)Total number of bandsPolymorphic bandsPer cent polymorphism detectedISSR F35′-(AG)~8~ CG-3′56.08787.5ISSR F45′-(AG)~8~ TG-3′53.71212100ISSR F95′-(GAA)~5~-3′39.6131292.3ISSR M15′-(AGC)~6~ G-3′63.1121191.6ISSR M25′-(ACC)~6~ G-3′63.11414100ISSR M35′-(AGC)~6~ C-3′63.1171694.1ISSR M85′-(AC)~9~ G-3′56.7131292.3ISSR M95′-(AC)~8~ CG-3′56.01313100ISSR M125′-(GACAC)~4~-3′61.466100ISSR M175′-CAG (CA)~8-~3′56.78787.5Total11611094.8

Random amplified polymorphic DNA analysis {#s2d}
-----------------------------------------

For RAPD reactions, a total of 43 primers (MWG Biotech-AC) were screened for polymorphism using selected genotypes, and 10 primers were selected for the final reactions (Table [3](#PLV083TB3){ref-type="table"}). The PCR mixture contained 1.5 μL of 10× *Taq* buffer with ammonium sulfate, 2.5 μL of 25 mM MgCl~2~, 3 μL of 1 mM dNTP, 3 units *Taq* DNA polymerase (Thermoscientific), 1.5 μL of 5 μM RAPD primer and 50 ng of template DNA in a total volume of 15 μL. Polymerase chain reaction amplifications were carried out utilizing a Eppendorf Master Cycler with initial denaturation at 94 °C for 3 min, followed by repeated cycles of denaturation at 94 °C for 45 s, annealing as per the primer\'s melting temperature for 1 min and primer extension at 72 °C for 1 min. On completion of the repeated number of cycles, final extension was performed at 72 °C for 10 min. Table 3.Characteristics and polymorphism revealed by RAPD primers for 95 wheat genotypes used in the study.RAPD primerSequenceMelting temperature (*T*~m~)Total number of bandsPolymorphic bandsPer cent polymorphism detectedcRAPD15′-GAA ACG GGT G-3′326466.6cRAPD25′-GTG ACG TAG G-3′32121191.6RAPD B35′-GTG ACG TAG G-3′349777.7RAPD B45′-CTC ACC GTC C-3′346583.3RAPD B55′-GAC GGA TCA G-3′32111090.9RAPD B105′-CTA CTG CGC T-3′327685.7RAPD B135′-TTC AGG GTG G-3325360.0RAPD L25′-GTT TCG CTC C-3′328787.5RAPD L45′- AAG AGC CCG T-3′3211981.8RAPD L65′-CCC GTC AGC A-3′347571.4Total826781.7

Data analysis {#s2e}
-------------

All the ISSR- and RAPD-based PCRs were repeated three times for the identification of reproducible amplified bands. Amplified fragments were counted from smaller to larger size. A binary data matrix was obtained by scoring the gel as 1 and 0 to show the presence and absence of bands, respectively. Information capacity of the primers and polymorphism content of the genotypes were estimated by calculating the total number of bands and of polymorphic bands. The binary data matrix was used to obtain the similarity matrix depending on simple matching (SM) coefficient by Numerical Taxonomy and Multivariate Analysis System (NTSYS-PC) version 2.02e software ([@PLV083C55]). This similarity matrix was utilized in R software for constructing a combined dendrogram of RAPD and ISSR.

On the basis of SM coefficients, the similarity matrix was double centred using the DCENTER module of NTSYS-PC. Then eigen analyses were performed using the EIGEN module of NTSYS-PC to construct two-dimensional scatterplots by the R package. Scatterplots were drawn for the substantiation of the dendrograms and verification of genotypes clustering according to both ploidy and geographical origin.

To explain the population structure of Indian and Turkish wheat genotypes, analysis of molecular variance (AMOVA) was performed using GenAlEx 6.5 software ([@PLV083C47], [@PLV083C48]) with 1000 permutations. The programme was used for the determination of variance components and estimating the total variation within and among the populations.

Bayesian model-based clustering with assumed *K* populations was employed for genetically homogenous group estimation in Indian and Turkish wheat germplasm. A parameter of 50 000 burn-in period and 100 000 Markov Chain Monte Carlo replication, along with the admixture model and correlated allele frequencies, was used in STRUCTURE, version 2.3.4 ([@PLV083C51]; [@PLV083C23], [@PLV083C24]; [@PLV083C30]). A total of 10 independent runs were performed for each value of *K* (from 1 to 4 assumed) ([@PLV083C20]). For the determination of the best possible *K* value elucidating the genetically distinctive clusters in the data, the Structure Harvester v6.0 ([@PLV083C18]) programme was used implementing parameters described by [@PLV083C20].

Results {#s3}
=======

Genetic diversity {#s3a}
-----------------

Ninety-five Indian and Turkish wheat varieties were amplified using 43 and 27 ISSR and RAPD markers, respectively. The 10 most polymorphic ISSR and RAPD primers generated 116 and 82 genetic loci, respectively, with a total of 198 loci. Among ISSR primers, ISSR M3 generated the maximum number of polymorphic fragments (16) and cRAPD2 was the most prolific RAPD primer (11). In total, 94.8 and 81.7 % bands were found to be polymorphic among ISSR and RAPD markers. The average number of polymorphic bands per primer was 11.0 and 6.7 for ISSR and RAPD primers, respectively (Tables [2](#PLV083TB2){ref-type="table"} and [3](#PLV083TB3){ref-type="table"}). For both primer types, the main amplified region was in the range of 300--2000 bp (Figs [1](#PLV083F1){ref-type="fig"} and [2](#PLV083F2){ref-type="fig"}). Figure 1.The inter-simple sequence repeat M3 primer amplification profile of 95 Indian and Turkish wheat genotypes. Figure 2.Random amplified polymorphic DNA B5 primer amplification profile of 95 Indian and Turkish wheat genotypes.

Genetic relationships/association {#s3b}
---------------------------------

A Fan dendrogram of the combined RAPD and ISSR data showed clear groupings of genotypes on the basis of both ploidy and origin (Fig. [3](#PLV083F3){ref-type="fig"}). On combining both RAPD and ISSR data, individual errors of either marker system are reduced and combined the dendrogram provided a more robust overview of the relatedness of Indian and Turkish populations. Figure 3.Simple matching coefficient-based Fan dendrogram using NTSYS-PC and R software package of 95 Indian and Turkish wheat genotypes.

On the basis of ploidy, wheat varieties were divided into three clusters, containing 18 tetraploid and 77 hexaploid varieties. Among hexaploid varieties, two Indian genotypes, NW 2036 and RAJ 4027, were separated as outliers from the rest. However, all the hexaploid genotypes were basically divided into two groups, and six Turkish genotypes were clustered with the Indian Hexaploid group. Similarity coefficients among Indian hexaploid varieties ranged from 0.71 to 0.98 while among Turkish hexaploid varieties ranged from 0.42 to 0.95.

Furthermore, the molecular variance factor in both Indian and Turkish populations was compared as a further measure of genetic diversity. Results from AMOVA for geographical origin indicated 77 % genetic variation within populations, while the variation between the populations was 23 % (PhiPT = 0.232; *P* = 0.010). On the basis of ploidy, AMOVA detected higher genetic variation within tetraploid and hexaploid populations (92 %); however, the genetic variation between ploidies was only 8 % (PhiPT = 0.078; *P* = 0.010) (Table [4](#PLV083TB4){ref-type="table"}). Table 4.Analysis of molecular variance in Indian and Turkish wheat populations.Source of variationdfSquare sumVariance componentPercentageProbabilityGeographic origin Among populations1133.524.4623\<0.001 Within populations931370.77614.7477Ploidy Among populations154.211.328\<0.001 Within populations931450.0915.5992

Population structure {#s3c}
--------------------

Principal coordinate analysis (PCoA) serves as a platform to provide a spatial illustration of the comparative genetic distances between the individuals. It also assesses the robustness of the differentiation among the groups classified by the dendrogram ([@PLV083C35]). In our scatterplots, the first two principal components explained 17.6 and 10.7 % of the total variation, respectively. In accordance with the dendrogram, hexaploid individuals were clearly separated from tetraploid varieties by the first principal coordinate (Fig. [4](#PLV083F4){ref-type="fig"}A). Similarly, the second principal coordinate (10.71 % of total variation) divided the Turkish populations from the Indian ones (Fig. [4](#PLV083F4){ref-type="fig"}B). Figure 4.Principal coordinate analysis of 95 Indian and Turkish wheat genotypes based on (A) ploidy level of the genotypes and (B) geographic origin of the genotypes.

For population genetic structure analysis, Bayesian clustering modelling was executed in the STRUCTURE software using genotyping data generated by 177 RAPD and ISSR loci. As the clustering model presumes the underlying existence of *K* clusters, an Evano test was performed and yielded *K* = 3 as the highest log-likelihood. This means that 3 was the optimum number of subpopulations, indicating that the two major population groups actually represent three distinct clusters.

The analysis of structure according to the geographical origin was performed by setting the range of possible number of subpopulations (*K*) from 1 to 4. Indian and Turkish populations involved in this procedure showed separation from each other in accordance with clusters obtained in PCoA. At *K* = 3, wheat genotypes were divided into three clusters with two main populations, 1 and 2 (Fig. [5](#PLV083F5){ref-type="fig"}A) representing Indian and Turkish wheat gene pool, respectively. Red colour bars represent individuals belonging to the Indian wheat gene pool while those in green belong to the Turkish gene pool. The Indian wheat gene pool was again distributed into subclusters with blue representing the tetraploid wheat population (Fig. [5](#PLV083F5){ref-type="fig"}B). The Indian population consisted of 79 accessions, of which 72 % belonged to the first cluster, 4 % to the second cluster and 24 % to the third; whereas the Turkish group consisted of 16 accessions with 13, 86 and 1 % belonging to the first, second and third cluster, respectively (Table [5](#PLV083TB5){ref-type="table"}). Some of the Indian and Turkish hexaploid genotypes, including NW_2036, RAJ_4027, Bayraktar_2000, Seval, Gün_91, Konya_2002, showed admixture clustering (Fig. [5](#PLV083F5){ref-type="fig"}C). Within the first, second and third clusters, expected heterozygosity within individuals was found to be 0.18, 0.15 and 0.16, respectively. Table 5.Proportion of membership of each pre-defined population in each of the three clusters obtained from STRUCTURE analysis.Given populationInferred clustersNumber of individuals12310.7200.0390.2417920.1330.8570.00916 Figure 5.(A) Three clusters inferred from population STRUCTURE analysis; red zone consists of basically Indian varieties with blue zone representing Indian tetraploid subpopulation and green zone includes basically Turkish varieties. (B) For distinctive clusters, vertical coordinates denote membership coefficients and each vertical line along with the horizontal coordinate denotes individual genotypes. Numbers in brackets denote their main population group, India and Turkey. (C) Collection of genotypes on the basis of Q, which explains the proportion of every individual genome that belongs to two distinct clusters.

Discussion {#s4}
==========

The complex nature and huge size of the wheat genome pose serious challenges towards genetic means of increasing its production. Hence, furthering our understanding of the wheat genome utilizing a variety of analyses has assisted efforts towards the genetic improvement of modern cultivars. Our examination of the literature to date found no prior genotypic characterization of the Indian and Turkish wheat varieties, simultaneously using RAPD and ISSR markers ([@PLV083C34]). The present study constitutes the first attempt to better understand jointly the origin, evolution and molecular diversity of Indian and Turkish wheat varieties at different ploidy levels.

Evaluating genetic diversity in Indian and Turkish wheat {#s4a}
--------------------------------------------------------

Since their introduction, ISSR and RAPD markers have been broadly utilized for variability estimation of wheat genotypes. Several RAPD- and ISSR-based diversity studies including diploid, tetraploid and hexaploid wheat have been published ([@PLV083C13]; [@PLV083C42]; [@PLV083C52], [@PLV083C53]; [@PLV083C6]; [@PLV083C67]; [@PLV083C38]; [@PLV083C68]; [@PLV083C3]; [@PLV083C27]; [@PLV083C4]; [@PLV083C14]; [@PLV083C57]; [@PLV083C64]; [@PLV083C11]; [@PLV083C45]). Due to the possession of diverse (A, B, D) genomes of wheat, tetraploid and hexaploid varieties were involved in the study.

Although Indian and Turkish wheat germplasms were not simultaneously used earlier, the average RAPD- and ISSR primer-based polymorphism, 81.7 and 94.8 %, respectively, revealed in this study, were comparable with several prior diversity studies. The very first attempt made by other researchers among Indian tetraploid wheat varieties revealed high genetic variability in durum released cultivars (50.6 %) in comparison to landraces (44.8 %) ([@PLV083C52]). [@PLV083C67] found 79.6 % polymorphism among 27 tetra- and hexaploid Indian genotypes using RAPD markers. A detailed study on 96 commercial Indian wheat genotypes, including tetraploids, triticale and hexaploids, indicating 78.8 % polymorphism, revealed a narrow genetic base of tetraploid cultivars in comparison to hexaploids ([@PLV083C68]).

The similarity coefficient values range among Indian hexaploid varieties observed in our study (0.71--0.98) was found to be higher than that in a previous study by [@PLV083C27] (0.52--0.82) using RAPD markers. In the present work, the average count of polymorphic bands per primer was higher in the case of ISSR (11) compared with that in RAPD (6.7). These results were consistent with a previous study by [@PLV083C53] on Indian tetraploid wheat varieties. Although limited studies have been performed on diversity assessment of Turkish wheat, [@PLV083C1] assessed the genetic variability of 100 durum wheat varieties using RAPD markers and observed higher genetic diversity in landraces than in cultivars. [@PLV083C15] distinguished 16 Turkish bread wheat varieties using RAPD markers with product size in the range of 300--2800 bp, which was similar to our results.

Analysis of genetic relationships among wheat genotypes {#s4b}
-------------------------------------------------------

The dendrogram obtained in this study clearly clustered the genotypes according to their ploidy level, consistently with the evolution of wheat ([@PLV083C2]). Furthermore, Indian and Turkish varieties were grouped separately. The information revealed by the dendrogram highlighted the parentage association of the varieties. Varieties HD_2177 and HD_2329 grouped together in the dendrogram with 95 % similarity share three common parents, HD_1962, E_4870, K_65. HD_2402 also grouped with its parent variety HD_2236 and showed 96 % similarity. Varieties Raj_1482, Raj_3072 and Raj_3077 were clustered together. Within this cluster, Raj_1482 is the parent of Raj_3077, with 93 % similarity. HD_2307 and HD_2501, which were grouped separately from other 'HD' varieties, share as a common parent HD_2160 and consistently exhibited 93 % similarity. Not only hexaploids, but also some of the tetraploid varieties like AKDW_2997 were also allocated in the same subcluster with its parent Raj_1555 and showed 88 % similarity (Fig. [3](#PLV083F3){ref-type="fig"}).

Analysis of molecular variance results disclosed in the study were in agreement with the UPGMA clustering and supported a high level of diversity within-country samples. Although the variation between Indian and Turkish populations was lower in comparison to within-population variation, it was significant according to the partitioning value (*P* = 0.010) (Table [4](#PLV083TB4){ref-type="table"}). The results suggest that similarity association between the countries was affected by within-country inconsistencies of the varieties. This high variation within groups can be attributed to selective adaptation towards the growth conditions at the time of breeding.

Investigating the Indian and Turkish wheat population structure {#s4c}
---------------------------------------------------------------

Similar separation of Indian and Turkish wheat varieties was observed by PCoA on the basis of ploidy and geographical region. The outcomes of the two methods (cluster analysis and PCoA) were comparable. Both of them classified 95 wheat genotypes mainly into three clusters and offered similar alignment of the genotypes with a few negligible discrepancies. The groups attained were in agreement with the recognized geographical origin as well.

Population structure analyses indicated that wheat accessions can be efficiently categorized on the basis of both geographical origin and ploidy. Using the maximum membership probability in STRUCTURE, Indian and Turkish populations showed similar grouping to the UPGMA and PCoA clustering. The PCoA clustering divided Indian and Turkish populations basically into similar clusters to those produced by the Structure bar plot at *K* = 3. In PCoA, some of the Indian and Turkish varieties showed close association with each other, and similar varieties demonstrated admixture in Structure analysis confirming their relatedness within the diverse gene pool. In dendrogram also, these varieties showed distinct clustering with the main population groups. Closeness of some of the Turkish hexaploid genotypes with Indian hexaploid genotypes in PCoA was also supported by the population structure as well as the dendrogram. Similar and mutually supportive results from all the statistical analyses demonstrated the capability of RAPD and ISSR markers to distinguish Indian and Turkish wheat varieties efficiently.

Conclusions {#s5}
===========

Genetic diversity evaluation serves as a crucial platform in plant improvement. The present study provides a detailed understanding of the genetic association of Indian and Turkish hexaploid and tetraploid wheat. The Turkish hexaploid populations showed their closeness to Indian genotypes, confirming their alliance within the diverse gene pool. The present genetic diversity study of wheat material obtained from diverse regions will support breeders in expanding the genetic variation of breeding accessions and utilizing the studied wheat resources more effectively.

Sources of Funding {#s6}
==================

M.K.K. has been granted '2216 Research Fellowship for Foreign Citizens' by TÜBİTAK (The Scientific and Technological Research Council of Turkey) for performing the present research work. Also, a part of the research was funded by BAP (Turkish Scientific Research Project agency) under Project No. 14401106.

Contributions by the Authors {#s7}
============================

M.K.K. initiated and obtained the funding for the research work. M.K.K. and A.P. contributed in performing the research work and preparation of the manuscript. M.K.K., A.P., S.A.K. and Y.O. carried out all statistical analyses. All other authors have provided suggestions and guidance for the successful completion of the work. All authors read and approved the final manuscript.

Conflict of Interest Statement {#s8}
==============================

None declared.

We thank Enrique Lopez-Juez for language editing and Hall Cushman, Kermit Ritland and the anonymous reviewers for their valuable suggestions that have greatly improved the manuscript.

[^1]: **Associate Editor:** Kermit Ritland
