The objective of this study was to analyze the genetic relationships, using PCR-based ISSR markers, among 70 Indian flax (Linum usitatissimum L.) genotypes actively utilized in flax breeding programs. Twelve ISSR primers were used for the analysis yielding 136 loci, of which 87 were polymorphic. The average number of amplified loci and the average number of polymorphic loci per primer were 11.3 and 7.25, respectively, while the percent loci polymorphism ranged from 11.1 to 81.8 with an average of 63.9 across all the genotypes. The range of polymorphism information content scores was 0.03-0.49, with an average of 0.18. A dendrogram was generated based on the similarity matrix by the Unweighted Pair Group Method with Arithmetic Mean (UPGMA), wherein the flax genotypes were grouped in five clusters. The Jaccard's similarity coefficient among the genotypes ranged from 0.60 to 0.97. When the omega-3 alpha linolenic acid (ALA) contents of the individual genotypes were correlated with the clusters in the dendrogram, the high ALA containing genotypes were grouped in two clusters. This study identified SLS 50, Ayogi, and Sheetal to be the most diverse genotypes and suggested their use in breeding programs and for developing mapping populations.