Development of mapped simple sequence repeat markers from common bean (Phaseolus vulgaris L.) based on genome sequences of a Chinese landrace and diversity evaluation

Microsatellite or single sequence repeat (SSR) markers have been commonly used in genetic research in many crop species, including common bean (Phaseolus vulgaris L.). A limited number of existing SSR markers have been designed from high-throughput sequencing of the genome, warranting the exploitation of new SSR markers from genomic regions. In this paper, we sequenced total DNA from the genotype Hong Yundou with a 454-FLX pyrosequencer and found numerous SSR loci. Based on these, a large number of SSR markers were developed and 90 genomic-SSR markers with clear bands were tested for mapping and diversity detection. The new SSR markers proved to be highly polymorphic for molecular polymorphism, with an average polymorphism information content value of 0.44 in 131 Chinese genotypes and breeding lines, effective for distinguishing Andean and Mesoamerican genotypes. In addition, we integrated 85 primers of the 90 polymorphism markers into the bean map using an F2 segregating population derived from Hong Yundou crossed with Jingdou. The distribution of SSR markers among 11 chromosomes was not random and tended to cluster on the linkage map, with 14 new markers mapped on chromosome Pv01, whereas only four loci were located on chromosome Pv04. Overall, these new markers have potential for genetic mapping, genetic diversity studies and map-based cloning in common bean.     

基于SSR标记分析大豆疫霉根腐病抗源的遗传多样性

丰富的遗传多样性可为大豆育种提供宽阔的遗传基础,本研究基于35对SSR标记,对60份东北地区大豆疫霉根腐病抗性品种进行了遗传多样性分析,共检测到189个等位基因,平均每个位点等位变异数5.4个,多态性信息含量指数(PIC)为0.1550~0.8195,平均为0.6636;遗传相似系数的变异范围为0.31~0.74。利用5对高多态性SSR引物构建了60份抗性材料的指纹图谱,这5对SSR引物构建的指纹图谱可以将60份疫霉根腐病抗性材料逐一区分开。采用NTSYS2.10基于遗传距离的聚类分析,将60份抗性材料分为7个类群,其中78.33%的抗性品种(系)的遗传相似系数在0.45~0.74间,表明遗传差异相对较窄,品种间遗传多样性水平较低。聚类分析与群体遗传结构分析结果有部分重合,均反映出不同地区的抗性材料间存在一定的渗透和交流。     

Development of genomic simple sequence repeat markers in faba bean by next-generation sequencing

Faba bean (Vicia faba L.) is an important food legume crop with a huge genome. Development of genetic markers for faba bean is important to study diversity and for molecular breeding. In this study, we used Next Generation Sequencing (NGS) technology for the development of genomic simple sequence repeat (SSR) markers. A total of 14,027,500 sequence reads were obtained comprising 4,208 Mb. From these reads, 56,063 contigs were assembled (16,367 Mb) and 2138 SSRs were identified. Mono and dinucleotides were the most abundant, accounting for 57.5 % and 20.9 % of all SSR repeats, respectively. A total of 430 primer pairs were designed from contigs larger than 350 nucleotides and 50 primers pairs were tested for validation of SSR locus amplification. Nearly all (96 %) of the markers were found to produce clear amplicons and to be reproducible. Thirty-nine SSR markers were then applied to 46 faba bean accessions from worldwide origins, resulting in 161 alleles with 87.5 % polymorphism, and an average of 4.1 alleles per marker. Gene diversity (GD) of the markers ranged from 0 to 0.48 with an average of 0.27. Testing of the markers showed that they were useful in determining genetic relationships and population structure in faba bean accessions.     

Comparative evaluation of genetic diversity using RAPD,SSR and cytochrome P450 gene based markers with respect to calcium content in finger millet (<Emphasis Type="Italic">Eleusine coracana</Emphasis> L. Gaertn.)

Genetic relationships among 52 Eleusine coracana (finger millet) genotypes collected from different districts of Uttarakhand were investigated by using randomly amplified polymorphic DNA (RAPD), simple sequence repeat (SSR) and cytochrome P450 gene based markers. A total of 18 RAPD primers, 10 SSR primers, and 10 pairs of cytochrome P450 gene based markers, respectively, revealed 49.4%, 50.2% and 58.7% polymorphism in 52 genotypes of E. coracana. Mean polymorphic information content (PIC) for each of these marker systems (0.351 for RAPD, 0.505 for SSR and 0.406 for cyt P450 gene based markers) suggested that all the marker systems were effective in determining polymorphisms. Pair-wise similarity index values ranged from 0.011 to 0.999 (RAPD), 0.010 to 0.999 (SSR) and 0.001 to 0.998 (cyt P450 gene based markers) and mean similarity index value of 0.505, 0.504 and 0.499, respectively. The dendrogram developed by RAPD, SSR and cytochrome P450 gene based primers analyses revealed that the genotypes are grouped in different clusters according to high calcium (300–450 mg/100 g), medium calcium (200–300 mg/100 g) and low calcium (100–200 mg/100 g). Mantel test employed for detection of goodness of fit established cophenetic correlation values above 0.95 for all the three marker systems. The dendrograms and principal coordinate analysis (PCA) plots derived from the binary data matrices of the three marker systems are highly concordant. High bootstrap values were obtained at major nodes of phenograms through WINBOOT software. Comparison of RAPD, SSR and cytochrome P450 gene based markers, in terms of the quality of data output, indicated that SSRs and cyt P450 gene based markers are particularly promising for the analysis of plant genome diversity. The genotypes of finger millet collected from different districts of Uttarakhand constitute a wide genetic base and clustered according to calcium contents. The identified genotypes could be used in breeding programmes and amajor input into conservation biology of cereal crops.     

用于绿豆种质资源遗传多样性分析的SSR及STS引物的筛选

目前能够用于绿豆(Vigna radiate)种质资源遗传多样性分析的PCR引物极其有限。通过12份农艺性状差异较大的绿豆种质对绿豆以及小豆(Vigna angularis)、豇豆(Vigna unguiculata)、菜豆(Phaseolus vulgaris)等近缘食用豆中的PCR引物进行筛选,结果表明41对绿豆SSR引物中能够有效扩增的有35对,6对有多态性;28对绿豆STS引物中有23对能够有效扩增,2对有多态性;8对小豆SSR引物能够有效扩增的有6对,但均无多态性;27对豇豆SSR引物能够有效扩增的有17对,1对有多态性;24对菜豆SSR引物能够有效扩增的有9对,1对有多态性。这些多态性引物的获得将有助于中国绿豆种质资源的遗传多样性分析。     

The use of microsatellites for detecting DNA polymorphism, genotype identification and genetic diversity in wheat

A set of 20 wheat microsatellite markers was used with 55 elite wheat genotypes to examine their utility (1) in detecting DNA polymorphism, (2)in the identifying genotypes and (3) in estimating genetic diversity among wheat genotypes. The 55 elite genotypes of wheat used in this study originated in 29 countries representing six continents. A total of 155 alleles were detected at 21 loci using the above microsatellite primer pairs (only 1 primer amplified 2 loci; all other primers amplified 1 locus each). Of the 20 primers amplifying 21 loci, 17 primers and their corresponding 18 loci were assigned to 13 different chromosomes (6 chromosomes of the A genome, 5 chromosomes of the B genome and 2 chromosomes of the D genome). The number of alleles per locus ranged from 1 to 13, with an average of 7.4 alleles per locus. The values of average polymorphic information content (PIC) and the marker index (MI) for these markers were estimated to be 0.71 and 0.70, respectively. The (GT)_n microsatellites were found to be the most polymorphic. The genetic similarity (GS) coefficient for all possible 1485 pairs of genotypes ranged from 0.05 to 0.88 with an average of 0.23. The dendrogram, prepared on the basis of similarity matrix using the UPGMA algorithm, delineated the above genotypes into two major clusters (I and II), each with two subclusters (Ia, Ib and IIa, IIb). One of these subclusters (Ib) consisted of a solitary genotype (E3111) from Portugal, so that it was unique and diverse with respect to all other genotypes belonging to cluster I and placed in subcluster Ia. Using a set of only 12 primer pairs, we were able to distinguish a maximum of 48 of the above 55 wheat genotypes. The results demonstrate the utility of microsatellite markers for detecting polymorphism leading to genotype identification and for estimating genetic diversity. Received: 15 May 1999 / Accepted: 27 July 1999     

Molecular characterization of <Emphasis Type="Italic">Vigna radiata</Emphasis> (L.) Wilczek genotypes based on nuclear ribosomal DNA and RAPD polymorphism

Mungbean germplasm characterization, evaluation and improvement are fundamentally based on morpho-agronomic traits. The lack of break-through in mungbean production has been due to non-availability of genetic variability for high yield potential. Forty-four genotypes of mungbean [Vigna radiata (L.)Wilczek] were subjected to random amplified polymorphic DNA (RAPD) analysis to assess the genetic diversity and relationships among the genotypes. Multilocus genotyping by twelve RAPD primers generated 166 markers and detected an average of intraspecific variation amounting to 82% polymorphism in banding patterns. Dendrogram obtained from cluster analysis delineated all the 44 genotypes into six clusters. Higher values of Nei’s gene diversity (h) and Shannon information index (i) and genetic distance analysis validate existence of wide genetic diversity among mungbean genotypes tested. Besides internal transcribed spacer (ITS) length variations, single nucleotide polymorphisms (SNPs) and insertions/deletions (INDELS) were detected at number of sites in nuclear rDNA region and the sequences of representatives of each sub-cluster and all distinct genotypes have been submitted to NCBI database and assigned Gen accession numbers HQ 148136-148147. Multiple sequence alignment revealed further lineages of distinct genotypes with main RAPD clusters. The measures of relative genetic distances among the genotypes of mungbean did not completely correlate the geographical places of their development. The homogeneous phenotypic markers proved insufficient in exhibiting genetic divergence among mungbean genotypes studied. RMG-62, RMG-976, and NDM-56 have been identified as potential source of parents for crop improvement. RAPD primers, OPA-9 and OPA-2 as polymorphic genetic markers and number of pods/plant and number of seeds/plant as dependable phenotypic markers have been identified for improving yield potentials. This genetic diversity will be of significance in developing intraspecific crosses in mungbean crop improvement programme.     

Assessment of genetic diversity in Lepidium sativum (Chandrasur) a medicinal herb used in folklore remedies in India using RAPD

Lepidium sativum L. is a fast growing edible herb which belongs to family Brassicaceae. The seeds of L. sativum are aperient, diuretic, tonic, demulcent, carminative, galatogogue and emmenagogue. They have been used in the treatment of bacterial and fungal infections, as an aperient and also possess antibacterial and antifungal properties. The seeds of this plant possess rapid bone fracture healing ability. Despite of its diverse medicinal properties no molecular data for diversity analysis is available till date. During this study random amplified polymorphic DNA (RAPD) markers were used to detect genetic variations of L. sativum. Initially 50 decamer primers were screened, out of which only 32 primers showed reproducible fragments with easily recordable bands. A total of 414 reproducible and clear bands were distinguished across the selected primers and statistical analysis showed 361 polymorphic bands and 53 monomorphic bands. Cluster analysis of the genotypes based on UPGMA divided the 18 genotypes into two main clusters, with first cluster having only HCS-20 genotype of L. sativum and other having rest of all 17 genotypes. The dendrogram based on similarity matrix revealed 23–66% genetic relatedness among 18 genotypes. The results of the present study can be used for molecular breeding and improvement of L. sativum for various desired traits through hybridization in future.     

基于SRAP和TRAP标记的河北野生狗牙根种质遗传多样性分析

以SRAP和TRAP 2种标记技术对36份狗牙根材料的遗传多样性及亲缘关系进行了分析,其中包含34份河北省野生狗牙根种质资源。分别由238对SRAP和85对TRAP引物组合中筛选获得具有多态性的SRAP和TRAP引物组合各10对,PCR扩增总条带分别为186和161条,多态性条带156和132条,平均每对引物扩增出多态性条带各15.6和13.2条,多态性位点比率分别为83.4%和81.0%。2种标记合并进行聚类分析,所有供试的36份狗牙根材料遗传相似系数GS=0.519~0.983,平均为0.7。当GS=0.68时,可将36份供试材料分为4个类群。本研究结果表明河北野生狗牙根种质资源存在较丰富的遗传多样性,可为种质资源保护和选育优良狗牙根新品种提供科学依据。     

Genetic Diversity of Phaeoisariopsis griseola in Kenya as Revealed by AFLP and Group-specific Primers

Genetic diversity of 50 Phaeoisariopsis griseola isolates collected from different agroecological zones in Kenya was studied using group‐specific primers and amplified fragment length polymorphism (AFLP) markers. Group‐specific primers differentiated the isolates into Andean and Mesoamerican groups, corresponding to the two common‐bean gene pools. Significant polymorphisms were observed with all the AFLP primer combinations used, reflecting a wide genetic diversity in the P. griseola population. A total of 207 fingerprints was generated, of which 178 were polymorphic. Cluster analysis of the polymorphic bands also separated the isolates into the two groups defined by group‐specific primers. All the isolates examined were grouped into three virulence populations; Andean, Afro‐Andean and Mesoamerican, and their genetic diversity measured. On average, greater diversity (91%) was detected within populations than between populations (9%). The genetic distance between Andean and Mesoamerican populations was higher (D = 0.0269) than between Andean and Afro‐Andean (D = 0.0095). The wide genetic diversity reported here has significant implications in breeding for resistance to angular leaf spot and should be taken into consideration when screening and deploying resistant bean genotypes.     

Qualitative and quantitative characterization of RAPD variation among snap bean (Phaseolus vulgaris) genotypes

Ten snap bean (Phaseolus vulgaris) genotypes were screened for polymorphism with 400 RAPD (random amplified polymorphic DNA) primers. Polymorphic RAPDs were scored and classified into three categories based on ethidium bromide staining intensity. An average of 5.19 RAPD bands were scored per primer for the 364 primers that gave scorable amplification products. An average of 2.15 polymorphic RAPDs were detected per primer. The results show that primer screening may reduce the number of RAPD reactions required for the analysis of genetic relationships among snap-bean genotypes by over 60%. Based on the analysis of the distribution of RAPD amplification, the same number of polymorphic RAPDs were amplified from different genotypes for all RAPD band intensity levels. A comparison of RAPD band amplification frequency among genotypes for the three categories of bands classified by amplification strength revealed a measurable difference in the frequencies of RAPDs classified as faint (weakly amplifying) compared to RAPD bands classified as bold (strongly amplifying) indicating a possible scoring error due to the underscoring of faint bands. Correlation analysis showed that RAPD bands amplified by the same primer are not more closely correlated then RAPD bands amplified by different primers but are more highly correlated then expected by chance. Pairwise comparisons of RAPD bands indicate that the distribution of RAPD amplification among genotypes will be a useful criterion for establishing RAPD band identity. For the average pairwise comparison of genotypes, 50% of primers tested and 15.8% of all scored RAPDs detected polymorphism. Based on RAPD data Nei's average gene diversity at a locus was 0.158 based on all scorable RAPD bands and 0.388 if only polymorphic RAPD loci were considered. RAPD-derived 1 relationships among genotypes are reported for the ten genotypes included in this study. The data presented here demonstrate that many informative, polymorphic RAPDs can be found among snap bean cultivars. These RAPDs may be useful for the unique identification of bean varieties, the organization of bean germplasm, and applications of molecular markers to bean breeding.     

CONSTRUCTION OF A GENETIC LINKAGE MAP FOR PORPHYRA HAITANENSIS (BANGIALES,RHODOPHYTA) BASED ON SEQUENCE‐RELATED AMPLIFIED POLYMORPHISM AND SIMPLE SEQUENCE REPEAT MARKERS1

Molecular markers and molecular genetic maps are prerequisites for molecular breeding in any plant species. A comprehensive genetic linkage map for cultivated Porphyra haitanensis T. J. Chang et B. F. Zheng has not yet been developed. In this study, 157 double haploid (DH) lines [derived from a YSIII (wildtype) × RTPM (red‐type artificial pigmentation mutant) cross] were used as a mapping population in P. haitanensis. A total of 60 pairs of sequence‐related amplified polymorphism (SRAP) primers and 39 pairs of simple sequence repeat (SSR) primers were used to detect polymorphisms between the two parents. Fifteen SRAP and 16 SSR polymorphic primer pairs were selected to analyze the DH population. A linkage genetic map comprising 67 SRAP markers and 20 SSR markers in five linkage groups, with a total length of 830.6 cM and an average of 10.13 cM between markers, was constructed. The markers were distributed evenly in all linkage groups without clustering. The linkage groups comprised 12–23 markers ranging in length from 134.2 to 197.3 cM. The estimated genome length of P. haitanensis was 942.4 cM, with 88.1% coverage. This is the first report of a comprehensive genetic map in P. haitanensis. The map presented here will provide a basis for the development of high‐density genetic linkage maps and lay the foundation for molecular breeding work in P. haitanensis.     

Polymorphism analysis in advanced mutant population of oat (<Emphasis Type="Italic">Avena sativa</Emphasis> L.) using ISSR markers

Present investigation was carried out to evaluate genetic diversity among 38 M6 population of oat cv. JO-1. To validate the observed morpho-physiological variations, these lines were analyzed with 21 ISSR primers. A total of 132 loci were amplified by these 21 ISSR markers and 116 loci were found to be polymorphic (87.87 %). The genetic similarity coefficient values among 39 oat genotypes based on ISSR analysis ranged from 0.305 to 0.957. The cluster analysis divided the oat genotypes into two groups. Mutants JMO 81 and JMO 82 were found to be most divergent, hence can be used as parents in breeding program for the development of superior cultivars.     

Genomic diversity among sorghum genotypes with resistance to sorghum shoot fly, Atherigona soccata

Host plant resistance is one of the important components for management of sorghum shoot fly, Atherigona soccata. The levels of resistance in cultivated germplasm are low to moderate, and therefore, it is important to identify sorghum genotypes with diverse mechanisms of resistance based on physico-chemical and or molecular markers. We assessed the genetic diversity of 15 sorghum genotypes with different levels of resistance/susceptibility to shoot fly, A. soccata using 93 sorghum simple sequence repeat (SSR) primer pairs and simultaneously characterized for 15 morpho-biochemical traits associated with shoot fly resistance. Of these 93 SSR primer pairs, amplification products from 79, thought to correspond to single-copy loci distributed across all ten sorghum chromosome pairs, showed good polymorphism across the 15 sorghum genotypes. The polymorphic information content (PIC) values of these 79 SSR markers ranged from 0.06 to 0.86. The Principal Coordinate Analyses (PCoA) and cluster analyses based on dissimilarity matrices derived from SSR based allelic variation (Neighbor-Joining distance) and variation in 15 morpho-biochemical traits (based on Gower??s distance), revealed grouping of most susceptible genotypes in single cluster. The improved breeding lines grouped with resistant or susceptible genotypes, based on shared pedigree. Based on these results, three resistant accessions viz., IS 1054, IS 1057 and IS 4664 were found diverse to IS 18551, which is widely used as shoot fly resistance donor. These diverse sources, after further characterization for resistance mechanisms, can be used in breeding programs for improving shoot fly resistance.     

Genetic diversity and population structure of <Emphasis Type="Italic">Salvia miltiorrhiza</Emphasis> Bge in China revealed by ISSR and SRAP

Salvia miltiorrhiza Bge is a traditional Chinese medicinal herb used as an important drug to cure cardiovascular diseases. In this work, inter simple sequence repeats (ISSR) and sequence related amplified polymorphism (SRAP) markers, were applied to assess the level and pattern of genetic diversity in five important cultivated populations of S. miltiorrhiza. Among these populations, 120 bands were amplified by 5 ISSR primers, of which all were polymorphic, and 110 polymorphic bands (90.16%) were observed in 122 bands amplified by 6 SRAP primers. A high levels of genetic diversity at the species level was detected with Hs = 0.1951, 0.1927 respectively. Analysis of molecular variance revealed that a greater proportion of total genetic variation existed within populations (86.64 and 84.83% respectively) rather than among populations (13.36 and 15.17% respectively). Cluster analysis divided the five populations into two groups. The genetic relationships among populations have low correlation with their geographical distribution (Mantel test; r = 0.4870 and 0.5740 respectively). The study indicated that both ISSR and SRAP markers were effective and reliable for assessing the degree of genetic variation of S. miltiorrhiza. Our results suggested that random collecting, preserving and planting seeds without deliberate selection might be an efficient way to conserve genetic resources of medicinal plants. Their effective use was also discussed on the further breeding.     

Developing and characterising Ricinus communis SSR markers by data mining of whole-genome sequences

Ricinus communis is a versatile industrial oil crop that is cultivated worldwide. Genetic improvement and marker-assisted breeding of castor bean have been slowed owing to the lack of abundant and efficient molecular markers. As co-dominant markers, simple sequence repeats (SSRs) are useful for genetic evaluation and molecular breeding. The recently released whole-genome sequence of castor bean provides useful genomic resources for developing markers on a genome-wide scale. In the present study, the distribution and frequency of microsatellites in the castor bean genome were characterised and numerous SSR markers were developed using genomic data mining. In total, 18,647 SSR loci at a density of one SSR per 18.89 Kb in the castor bean genome sequence (representing approximately 352.27 Mb) were identified. Dinucleotide repeats were the most frequently observed microsatellites, although the AAT repeat motif was also prevalent. Using six cultivars as screening samples, 670 polymorphic SSR markers from 1,435 primer pairs (46.7 %) were developed. Trinucleotide motif loci contained a higher proportion of polymorphisms (48.5 %) than dinucleotide motif loci (39.2 %). The polymorphism level in the SSR loci was positively correlated with the increasing number of repeat units in the microsatellites. The phylogenetic relationship among 32 varieties was evaluated using the developed SSR markers. Cultivars developed at the same institute clustered together, suggesting that these cultivars have a narrow genetic background. The large number of SSR markers developed in this study will be useful for genetic mapping and for breeding improved castor-oil plants. These markers will also facilitate genetic and genomic studies of Euphorbiaceae.     

Development and application of 15 novel polymorphic microsatellite markers for sect. Paeonia (Paeonia L.)

Herbaceous peony, which belongs to sect. Paeonia, Paeonia L., Paeoniaceae, is not only a famous traditional flower in China; it is also highly regarded as an ornamental in Europe and the USA. On account of the abundance in germplasm resources and wide distribution ranges, herbaceous peony can be divided into three cultivar groups, the Chinese Peony Cultivar Group, the European Peony Cultivar Group and the Hybrid Peony Cultivar Group. However, most studies on genetic relationships in Paeonia are limited to the first group. This study used 15 polymorphic microsatellite markers that were developed by magnetic bead enrichment to explore the genetic diversity of 89 genotypes collected from Beijing, the USA and Canada and introduced to China. From 145 allelic loci that were detected from 15 primers, 140 (96.6%) were polymorphic. The number of allelic loci ranged from 4 to 16 with an average of 9.3, and the polymorphic index content ranged from 0.362 to 0.825 with an average of 0.678. Both cluster analysis and principal component analysis based on the 15 polymorphic primer pairs could distinguish the species from cultivars although the cluster analysis was better able to reflect the genetic relationships between all individuals. Some new insights about Paeonia L. are suggested. Furthermore, the 15 SSR markers have great significance as they will allow for furthering the studies on genetic diversity, genetic relationships and even the molecular marker-assisted breeding of Paeonia L.     

New gene-derived simple sequence repeat markers for common bean (Phaseolus vulgaris L.)

Common bean is an important and diverse crop legume with several wild relatives that are all part of the Phaseoleae tribe of tropical crop legumes. Sequence databases have been a good source of sequences to mine for simple sequence repeats (SSRs). The objective of this research was to evaluate 14 sequence collections from common bean for SSRs and to evaluate the diversity of the polymorphic microsatellites derived from these collections. SSRs were found in 10 of the GenBank sequence collections with an average of 11.3% of sequences containing microsatellite motifs. The most common motifs were based on tri- and dinucleotides. In a marker development programme, primers were designed for 125 microsatellites which were tested on a panel of 18 common bean genotypes. The markers were named as part of the bean microsatellite-database (BMd) series, and the average polymorphism information content was 0.404 for polymorphic markers and predicted well the genepool structure of common beans and the status of the wild and cultivated accessions that were included in the study. Therefore, the BMd series of microsatellites is useful for multiple studies of genetic relatedness and as anchor markers in future mapping of wide crosses in the species.     

Association analysis,genetic diversity and structure analysis of tobacco based on AFLP markers

Knowledge in the area of genetic diversity could aid in providing useful information in the selection of material for breeding such as hybridization programs and quantitative trait loci mapping. To this end, 50 Nicotiana tabacum genotypes were genotyped with 21 primer combination of amplified fragment length polymorphism (AFLP). A total of 480 unambiguous DNA fragments and 373 polymorphic bands were produced with an average of 17.76 per primer combination. Also, the results revealed high polymorphic rate varing from 52.63 to 92.59 %, demonstrating that AFLP technique utilized in this research can be a powerful and valuable tool in the breeding program of N. tabacum. Cluster analysis based on complete linkage method using Jaccard’s genetic distance, grouped the 50 tobacco genotypes into eight clusters including three relatively big clusters, one cluster including Golden gift, Burly 7022 and Burly Kreuzung, one cluster consisting of two individuals (Pereg234, R9) and three single-member clusters (Pennbel69, Coker176 and Budisher Burley E), Recent genotypes showed high genetic distance from other genotypes belonging to cluster I and II. Association analysis between seven important traits and AFLP markers were performed using four statistical models. The results revealed the model containing both the factors, population structure (Q) and general similarity in genetic background arising from shared kinship (K), reduces false positive associations between markers and phenotype. According to the results nine markers were determined that could be considered to be the most interesting candidates for further studies.