A comprehensive analysis integrating phenotypic assessment uncovering thornless cultivar lineages in Aralia elata.

Aralia elata is an Araliaceae woody plant species found in Northeastern Asia. To understand how genetic pools are distributed for A.elata clones, we were to analyze the population structure of A.elata cultivars and identify how these are correlated with thorn-related phenotype which determines the utility of A.elata. We found that the de novo assembled genome of 'Yeongchun' shared major genomic compartments with the public A.elata genome assembled from the wild-type from China. To identify the population structure of the 32 Korean and Japanese cultivars, we identified 44 SSR markers and revealed three main sub-clusters using ΔK analysis with one isolated cultivar. Machine-learning based clustering with thorn-related phenotype correlated moderately with population structure based on SSR analysis suggested multi-layered genetic regulation of thorn-related phenotypes. Thus, we revealed genetic lineage of A.elata and uncovered isolated cultivar which can provide new genetic material for further breeding.

Comparative chloroplast genome analyses of cultivated and wild Capsicum species shed light on evolution and phylogeny.

The chloroplast (cp.) genome, also known as plastome, plays crucial roles in plant survival, adaptation, and evolution. The stable genetic structure of cp. genomes provides an ideal system for investigating species evolution. We sequenced three complete cp. genome sequences of Capsicum species and analyzed them using sequences of various Capsicum species retrieved from the NCBI database. The cp. genome of Capsicum species maintains a well-preserved quadripartite structure consisting of two inverted repeats (IRs) flanked by a large single copy (LSC) region and a small single copy (SSC) region. The sizes of cp. genome sequences ranged from 156,583 bp (C. lycianthoides) to 157,390 bp (C.pubescens). A total of 127-132 unique genes, including 83-87 protein-coding, 36-37 tRNA, and eight rRNA genes, were predicted. Comparison of cp. genomes of 10 Capsicum species revealed high sequence similarity in genome-wide organization and gene arrangements. Fragments of trnT-UGU/trnL-UAA, ccsA, ndhD, rps12, and ycf1 were identified as variable regions, and nucleotide variability of LSC and SSC was higher than that of IR. Phylogenetic speciation analysis showed that the major domesticated C. annuum species were the most extensively divergent species and closely related to C. tovarii and C. frutescens. Analysis of divergent times suggested that a substantial range of speciation events started occurring ~ 25.79 million years ago (Mya). Overall, comparative analysis of cp. genomes of Capsicum species not only offers new insights into their genetic variation and phylogenetic relationships, but also lays a foundation for evolutionary history, genetic diversity, conservation, and biological breeding of Capsicum species.

Establishment of a UPLC-PDA/ESI-Q-TOF/MS-Based Approach for the Simultaneous Analysis of Multiple Phenolic Compounds in Amaranth (A. cruentus and A. tricolor).

We used ultraperformance liquid chromatography coupled with a photodiode-array detector and electrospray ionization quadrupole time-of-flight mass spectrometry (UPLC-PDA/ESI-Q-TOF/MS) to rapidly and accurately quantify 17 phenolic compounds. Then, we applied this method to the seed and leaf extracts of two Amaranthus species to identify and quantify phenolic compounds other than the 17 compounds mentioned above. Compounds were eluted within 30 min on a C18 column using a mobile phase (water and acetonitrile) containing 0.1% formic acid, and the specific wavelength and ion information of the compounds obtained by PDA and ESI-Q-TOF/MS were confirmed. The proposed method showed good linearity (r2 > 0.990). Limits of detection and quantification were less than 0.1 and 0.1 µg/mL, respectively. Intra- and interday precision were less than 2.4% and 1.8%, respectively. Analysis of amaranth seed and leaf extracts using the established method showed that the seeds contained high amounts of 2,4-dihydroxybenzoic acid and kaempferol, and leaves contained diverse phenolic compounds. In addition, six tentatively new phenolic compounds were identified. Moreover, seeds potentially contained 2,3-dihydroxybenzaldehyde, a beneficial bioactive compound. Thus, our method was an efficient approach for the qualitative and quantitative analysis of phenolic compounds, and could be used to investigate phenolic compounds in plants.

De novo transcriptome assembly and the identification of gene-associated single-nucleotide polymorphism markers in Asian and American ginseng roots.

We performed de novo transcriptome sequencing for Panax ginseng and Panax quinquefolius accessions using the 454 GS FLX Titanium System and discovered annotation-based genome-wide single-nucleotide polymorphism (SNPs) using next-generation ginseng transcriptome data without reference genome sequence. The comprehensive transcriptome characterization with the mature roots of four ginseng accessions generated 297,170 reads for 'Cheonryang' cultivar, 305,673 reads for 'Yunpoong' cultivar, 311,861 reads for the G03080 breeding line, and 308,313 reads for P. quinquefolius. In transcriptome assembly, the lengths of the sample read were 156.42 Mb for 'Cheonryang', 161.95 Mb for 'Yunpoong', 165.07 Mb for G03080 breeding line, and 166.48 Mb for P. quinquefolius. A total of 97 primer pairs were designed with the homozygous SNP presented in all four accessions. SNP genotyping using high-resolution melting (HRM) analysis was performed to validate the putative SNP markers of 97 primer pairs. Out of the 73 primer pairs, 73 primer pairs amplified the target sequence and 34 primer pairs showed polymorphic melting curves in samples from 11 P. ginseng cultivars and one P. quinquefolius accession. Among the 34 polymorphic HRM-SNP primers, four primers were useful to distinguish ginseng cultivars. In the present study, we demonstrated that de novo transcriptome assembly and mapping analyses are useful in providing four HRM-SNP primer pairs that reliably show a high degree of polymorphism among ginseng cultivars.

Novel microsatellite markers acquired from Rubus coreanus Miq. and cross-amplification in other Rubus species.

The Rubus genus consists of more than 600 species that are distributed globally. Only a few Rubus species, including raspberries and blueberries, have been domesticated. Genetic diversity within and between Rubus species is an important resource for breeding programs. We developed genomic microsatellite markers using an SSR-enriched R. coreanus library to study the diversity of the Rubus species. Microsatellite motifs were discovered in 546 of 646 unique clones, and a dinucleotide repeat was the most frequent (75.3%) type of repeat. From 97 microsatellite loci with reproducible amplicons, we acquired 29 polymorphic microsatellite markers in the Rubus coreanus collection. The transferability values ranged from 59.8% to 84% across six Rubus species, and Rubus parvifolius had the highest transferability value (84%). The average number of alleles and the polymorphism information content were 5.7 and 0.541, respectively, in the R. coreanus collection. The diversity index of R. coreanus was similar to the values reported for other Rubus species. A phylogenetic dendrogram based on SSR profiles revealed that seven Rubus species could be allocated to three groups, and that R. coreanus was genetically close to Rubus crataegifolius (mountain berry). These new microsatellite markers might prove useful in studies of the genetic diversity, population structure, and evolutionary relationships among Rubus species.

The complete chloroplast genome of Capsicum annuum var. glabriusculum using Illumina sequencing.

Chloroplast (cp) genome sequences provide a valuable source for DNA barcoding. Molecular phylogenetic studies have concentrated on DNA sequencing of conserved gene loci. However, this approach is time consuming and more difficult to implement when gene organization differs among species. Here we report the complete re-sequencing of the cp genome of Capsicum pepper (Capsicum annuum var. glabriusculum) using the Illumina platform. The total length of the cp genome is 156,817 bp with a 37.7% overall GC content. A pair of inverted repeats (IRs) of 50,284 bp were separated by a small single copy (SSC; 18,948 bp) and a large single copy (LSC; 87,446 bp). The number of cp genes in C. annuum var. glabriusculum is the same as that in other Capsicum species. Variations in the lengths of LSC; SSC and IR regions were the main contributors to the size variation in the cp genome of this species. A total of 125 simple sequence repeat (SSR) and 48 insertions or deletions variants were found by sequence alignment of Capsicum cp genome. These findings provide a foundation for further investigation of cp genome evolution in Capsicum and other higher plants.

Cross-amplification of Vicia sativa subsp. sativa microsatellites across 22 other Vicia species.

The temperate and herbaceous genus Vicia L. is a member of the legume tribe Fabeae of the subfamily Papilionoideae. The genus Vicia comprises 166 annual or perennial species distributed mainly in Europe, Asia, and North America, but also extending to the temperate regions of South America and tropical Africa. The use of simple sequence repeat (SSR) markers for Vicia species has not been investigated as extensively as for other crop species. In this study, we assessed the potential for cross-species amplification of cDNA microsatellite markers developed from common vetch (Vicia sativa subsp. sativa). For cross-species amplification of the SSRs, amplification was carried out with genomic DNA isolated from two to eight accessions of 22 different Vicia species. For individual species or subspecies, the transferability rates ranged from 33% for V. ervilia to 82% for V. sativa subsp. nigra with an average rate of 52.0%. Because the rate of successful SSR marker amplification generally correlates with genetic distance, these SSR markers are potentially useful for analyzing genetic relationships between or within Vicia species.

Comparative Analysis of Polyphenolic Compounds in Different Amaranthus Species: Influence of Genotypes and Harvesting Year.

Amaranth is a nutritionally valuable crop, as it contains phenolic acids and flavonoids, yielding diverse plant secondary metabolites (PSMs) like phytosterol, tocopherols, and carotenoids. This study explored the variations in the contents of seventeen polyphenolic compounds within the leaves of one hundred twenty Amaranthus accessions representing nine Amaranthus species. The investigation entailed the analysis of phenolic content across nine Amaranthus species, specifically A. hypochondriacus, A. cruentus, A. caudatus, A. tricolor, A. dubius, A. blitum, A. crispus, A. hybridus, and A. viridis, utilizing ultra performance liquid chromatography with photodiode array detection (UPLC-PDA). The results revealed significant differences in polyphenolic compounds among accessions in which rutin content was predominant in all Amaranthus species in both 2018 and 2019. Among the nine Amaranthus species, the rutin content ranged from 95.72 ± 199.17 µg g-1 (A. dubius) to 1485.09 ± 679.51 µg g-1 (A. viridis) in 2018 and from 821.59 ± 709.95 µg g-1 (A. tricolor) to 3166.52 ± 1317.38 µg g-1 (A. hypochondriacus) in 2019. Correlation analysis revealed, significant positive correlations between rutin and kaempferol-3-O-ß-rutinoside (r = 0.93), benzoic acid and ferulic acid (r = 0.76), and benzoic acid and kaempferol-3-O-ß-rutinoside (r = 0.76), whereas gallic acid showed consistently negative correlations with each of the 16 phenolic compounds. Wide variations were identified among accessions and between plants grown in the two years. The nine species and one hundred twenty Amaranthus accessions were clustered into six groups based on their seventeen phenolic compounds in each year. These findings contribute to expanding our understanding of the phytochemical traits of accessions within nine Amaranthus species, which serve as valuable resources for Amaranthus component breeding and functional material development.

Development of 65 novel polymorphic cDNA-SSR markers in common vetch (Vicia sativa subsp. sativa) using next generation sequencing.

Vetch (Vicia sativa L.) is one of the most important annual forage legumes in the World due to its multiple uses (i.e., hay, grain, silage and green manure) and high nutritional value. However, detrimental cyanoalanine toxins in its plant parts including seeds and its vulnerability to hard winter conditions are currently reducing the agronomic values of vetch varieties. Moreover, the existence in the public domain of very few genomic resources, especially molecular markers, has further hampered breeding efforts. Polymorphic simple sequence repeat markers from transcript sequences (cDNA; simple sequence repeat [SSR]) were developed for Vicia sativa subsp. sativa. We found 3,811 SSR loci from 31,504 individual sequence reads, and 300 primer pairs were designed and synthesized. In total, 65 primer pairs were found to be consistently scorable when 32 accessions were tested. The numbers of alleles ranged from 2 to 19, frequency of major alleles per locus were 0.27-0.87, the genotype number was 2-19, the overall polymorphism information content (PIC) values were 0.20-0.86, and the observed and expected heterozygosity values were 0.00-0.41 and 0.264-0.852, respectively. These markers provide a useful tool for assessing genetic diversity, population structure, and positional cloning, facilitating vetch breeding programs.

Development and molecular characterization of 55 novel polymorphic cDNA-SSR markers in faba bean (Vicia faba L.) using 454 pyrosequencing.

Faba bean (Vicia faba L.) is a major food source and fodder legume, popularly known for its high content of seed-protein. Its role is critical in crop rotation, and for fixing nitrogen effectively. Polymorphic simple sequence repeat markers from transcript sequences (cDNA; simple sequence repeat [SSR]) were developed for faba bean (Vicia faba). We found that 1,729 SSR loci from 81,333 individual sequence reads and 240 primer pairs were designed and synthesized. In total, 55 primer pairs were found to be polymorphic and scorable consistently when screened in 32 accessions. The number of alleles ranged from 2 to 15, frequency of major alleles per locus varied from 0.17 to 0.91, the genotypes number ranged from 2 to 17, observed and expected heterozycosity values ranged from 0.00 to 0.44 and 0.17 to 0.89 and overall PIC values ranged from 0.16 to 0.88 respectively. These markers will be a useful tool for assessing the genetic diversity, understanding the population structure, and breeding patterns of faba bean.

A chromosome-level genome assembly of Korean mint (Agastache rugosa).

Agastache rugosa, also known as Korean mint, is a perennial plant from the Lamiaceae family that is traditionally used for various ailments and contains antioxidant and antibacterial phenolic compounds. Molecular breeding of A. rugosa can enhance secondary metabolite production and improve agricultural traits, but progress in this field has been delayed due to the lack of chromosome-scale genome information. Herein, we constructed a chromosome-level reference genome using Nanopore sequencing and Hi-C technology, resulting in a final genome assembly with a scaffold N50 of 52.15 Mbp and a total size of 410.67 Mbp. Nine pseudochromosomes accounted for 89.1% of the predicted genome. The BUSCO analysis indicated a high level of completeness in the assembly. Repeat annotation revealed 561,061 repeat elements, accounting for 61.65% of the genome, with Copia and Gypsy long terminal repeats being the most abundant. A total of 26,430 protein-coding genes were predicted, with an average length of 1,184 bp. The availability of this chromosome-scale genome will advance our understanding of A. rugosa's genetic makeup and its potential applications in various industries.

Volatile Compositions of Panax ginseng and Panax quinquifolium Grown for Different Cultivation Years.

The present study examined the volatile profiles of Panax ginseng (Korean ginseng) and Panax quinquefolium (American ginseng) grown for different cultivation years by using HS-SPME/GC-MS and determined the key discriminant volatile compounds by chemometric analysis including principal component analysis (PCA), hierarchical cluster analysis (HCA), and partial least squares-discrimination analysis (PLS-DA). Fifty-six compounds, including forty terpenes, eight alcohols, one alkane, one ketone, and one furan, were identified in the ginseng roots. The chemometric results identified two major clusters of American ginseng and Korean ginseng cultivars with distinct volatile compositions. The volatile compounds in fresh white ginseng roots were affected by the species, but the influence of different cultivation ages was ambiguous. The major volatile components of ginseng roots are terpenes, including monoterpenes and sesquiterpenes. In particular, panaginsene, ginsinsene, α-isocomene, and caryophyllene were predominant in Korean ginseng cultivars, whereas ß-farnesene levels were higher in American ginseng. The difference in volatile patterns between Panax ginseng and Panax quinquefolium could be attributed to the composition of sesquiterpenes such as ß-panaginsene, ginsinsene, caryophyllene, and ß-farnesene.

Pseudo-Chromosomal Genome Assembly in Combination with Comprehensive Transcriptome Analysis in Agaricus bisporus Strain KMCC00540 Reveals Mechanical Stimulus Responsive Genes Associated with Browning Effect.

Agaricus bisporus is one of the world's most popular edible mushrooms, including in South Korea. We performed de novo genome assembly with a South Korean white-colored cultivar of A. bisporus, KMCC00540. After generating a scaffold-level genomic sequence, we inferred chromosome-level assembly by genomic synteny analysis with the representative A. bisporus strains H97 and H39. The KMCC00540 genome had 13 pseudochromosomes comprising 33,030,236 bp mostly covering both strains. A comparative genomic analysis with cultivar H97 indicated that most genomic regions and annotated proteins were shared (over 90%), ensuring that our cultivar could be used as a representative genome. However, A. bisporus suffers from browning even from only a slight mechanical stimulus during transportation, which significantly lowers its commercial value. To identify which genes respond to a mechanical stimulus that induces browning, we performed a time-course transcriptome analysis based on the de novo assembled genome. Mechanical stimulus induces up-regulation in long fatty acid ligase activity-related genes, as well as melanin biosynthesis genes, especially at early time points. In summary, we assembled the chromosome-level genomic information on a Korean strain of A. bisporus and identified which genes respond to a mechanical stimulus, which provided key hints for improving the post-harvest biological control of A. bisporus.

Complete Chloroplast Genome of the Inverted Repeat-Lacking Species Vicia bungei and Development of Polymorphic Simple Sequence Repeat Markers.

Background: Vicia bungei is an economically important forage crop in South Korea and China. Although detailed genetic and genomic data can improve population genetic studies, conservation efforts, and improved breeding of crops, few such data are available for Vicia species in general and none at all for V. bungei. Therefore, the main objectives of this study were to sequence, assemble, and annotate V. bungei chloroplast genome and to identify simple sequence repeats (SSRs) as polymorphic genetic markers. Results: The whole-genome sequence of V. bungei was generated using an Illumina MiSeq platform. De novo assembly of complete chloroplast genome sequences was performed for the low-coverage sequence using CLC Genome Assembler with a 200-600-bp overlap size. Vicia bungei chloroplast genome was 130,796-bp long. The genome lacked an inverted repeat unit and thus resembled those of species in the inverted repeat-lacking clade within Fabaceae. Genome annotation using Dual OrganellarGenoMe Annotator (DOGMA) identified 107 genes, comprising 75 protein-coding, 28 transfer RNA, and 4 ribosomal RNA genes. In total, 432 SSRs were detected in V. bungei chloroplast genome, including 64 mononucleotides, 14 dinucleotides, 5 trinucleotides, 4 tetranucleotides, 233 pentanucleotides, 90 hexanucleotides, and 14 complex repeated motifs. These were used to develop 232 novel chloroplast SSR markers, 39 of which were chosen at random to test amplification and genetic diversity in Vicia species (20 accessions from seven species). The unweighted pair group method with arithmetic mean cluster analysis identified seven clusters at the interspecies level and intraspecific differences within clusters. Conclusion: The complete chloroplast genome sequence of V. bungei was determined. This reference genome should facilitate chloroplast resequencing and future searches for additional genetic markers using population samples. The novel chloroplast genome resources and SSR markers will greatly contribute to the conservation of the genus Vicia and facilitate genetic and evolutionary studies of this genus and of other higher plants.

The complete mitochondrial genome of the poisonous mushroom Trichoderma cornu-damae (Hypocreaceae).

Trichoderma cornu-damae is a poisonous mushroom that contains trichothecene mycotoxins. The complete mitochondrial genome of this mushroom was determined using next-generation sequencing. This mitogenome is a circular molecule 94,608 bp in length with a GC content of 27.94% and contains 15 protein-coding genes, two rRNA genes (rnl and rns), and 25 tRNA genes. Phylogenetic analysis placed T. cornu-damae in the family Hypocreaceae group, which includes the genus Trichoderma. The mitogenome of T. cornu-damae will contribute to our understanding of the phylogeny, taxonomy, and population genetics of this mushroom.

Sequence information on simple sequence repeats and single nucleotide polymorphisms through transcriptome analysis of mungbean.

Mungbean (Vigna radiata (L.) Wilczek) is a unique species in its ability to fix atmospheric nitrogen, with early maturity, and relatively good drought resistance. We used 454 sequencing technology for transcriptome sequencing. A total of 150 159 and 142 993 reads produced 5 254 and 6 374 large contigs (≥ 500 bp) with an average length of 833 and 853 for Sunhwa and Jangan, respectively. Functional annotation to known sequences yielded 41.34% and 41.74% unigenes for Jangan and Sunhwa. A higher number of simple sequence repeat (SSR) motifs was identified in Jangan (1 630) compared with that of Sunhwa (1 334). A similar SSR distribution pattern was observed in both varieties. A total of 8 249 single nucleotide polymorphisms (SNPs) and indels with 2 098 high-confidence candidates were identified in the two mungbean varieties. The average distance between individual SNPs was approximately 860 bp. Our report demonstrates the utility of transcriptomic data for implementing a functional annotation and development of genetic markers. We also provide large resource sequence data for mungbean improvement programs.

Development of CAPS Markers for Evaluation of Genetic Diversity and Population Structure in the Germplasm of Button Mushroom (Agaricus bisporus).

Agaricus bisporus is a globally cultivated mushroom with high economic value. Despite its widespread cultivation, commercial button mushroom strains have little genetic diversity and discrimination of strains for identification and breeding purposes is challenging. Molecular markers suitable for diversity analyses of germplasms with similar genotypes and discrimination between accessions are needed to support the development of new varieties. To develop cleaved amplified polymorphic sequences (CAPs) markers, single nucleotide polymorphism (SNP) mining was performed based on the A. bisporus genome and resequencing data. A total of 70 sets of CAPs markers were developed and applied to 41 A. bisporus accessions for diversity, multivariate, and population structure analyses. Of the 70 SNPs, 62.85% (44/70) were transitions (G/A or C/T) and 37.15% (26/70) were transversions (A/C, A/T, C/G, or G/T). The number of alleles per locus was 1 or 2 (average = 1.9), and expected heterozygosity and gene diversity were 0.0-0.499 (mean = 0.265) and 0.0-0.9367 (mean = 0.3599), respectively. Multivariate and cluster analyses of accessions produced similar groups, with F-statistic values of 0.134 and 0.153 for distance-based and model-based groups, respectively. A minimum set of 10 markers optimized for accession identification were selected based on high index of genetic diversity (GD, range 0.299-0.499) and major allele frequency (MAF, range 0.524-0.817). The CAPS markers can be used to evaluate genetic diversity and population structure and will facilitate the management of emerging genetic resources.

Identification of Vicia Species Native to South Korea Using Molecular and Morphological Characteristics.

Recently, within the Fabaceae family, the Vicia genus has been recognized for its vital role in sustainable agriculture. Vicia species are economically important grain and forage crops. However, the presence of complex morphological characteristics makes identification and recognition of native species difficult. In this study, the possibility of using DNA barcoding regions (ITS2, matK, and rbcL) to distinguish among 19 Vicia taxa (59 accessions) found in South Korea was evaluated. The sequence alignment analysis revealed considerable nucleotide diversity (π) between the loci, in which ITS2 showed the highest mean interspecific distance, whereas there was no intraspecific variability among the barcode regions in 12 of the 19 taxa. Phylogenetic analysis of combined barcoding regions revealed well-resolved phylogeny with the highest species level discrimination. Combinations of barcode loci were also used in classification at the subgenera and section levels. The results revealed that the combined barcoding regions can be used effectively to differentiate the following species: Vicia angustifolia var. segetilis, Vicia bungei, Vicia villosa, Vicia cracca, Vicia dasycarpa, Vicia hirsuta, Vicia tetrasperma, Vicia amurensis, Vicia hirticalycina, and Vicia chosenensis. However, it is difficult to differentiate the species of Vicia unijuga, Vicia unijuga var. kaussanensis, Vicia linearifolia, Vicia unijuga f. angustifolia, Vicia nipponica, Vicia amoena, Vicia venosa var. cuspidata, Vicia pseudo-orobus, and Vicia japonica with the tested barcode regions. These species come under sect. Vicilla and are found to be closely related or species that have recently undergone speciation; thus, it has limitation to distinguish with recommended barcodes. Hence, to differentiate the unclassified species, 39 morphological characteristics were investigated, in which 16 useful characteristics were selected for efficient classification. Finally, the 16 selected morphological useful traits efficiently differentiated all the Vicia species. In conclusion, a combination of barcoding loci together with morphological characteristics of this study efficiently discriminated all the Korean Vicia species.

Complete mitochondrial genome sequence of Gymnopilus junonius.

Gymnopilus junonius, a well-known poisonous mushroom, is distributed worldwide. It contains a hallucinogenic alkaloid psilocybin and several other bioactive compounds. The mitochondrial genome, a circular DNA molecule of 161,145 bp, comprises 15 protein-coding genes, 24 transfer RNA genes, and 2 ribosomal RNA genes. The guanine-cytosine content was 31.56%. Based on the mitochondrial genome sequence, a phylogenetic tree was constructed to demonstrate the phylogenetic relationship. In this study, the phylogenetic positions of G. junonius and its related genera were determined.

Interspecific variation of seed morphological and micro-morphological traits in the genus Vicia (Fabaceae).

Seed macro- and micro-morphology were analyzed to evaluate their capacity to discriminate species in the genus Vicia (Fabaceae). To assess the interspecific variation of the taxa in the genus Vicia, 41 accessions were obtained from the USDA-ARS germplasm collection in the USA and 19 accessions were collected from Korea. Seed morphological characteristics such as shape, color, mottling, finish, length, width, diameter, hilum shape, hilum color, hilum length, and lens distance from the hilum were examined under a stereomicroscope. Testa texture characteristics such as testa pattern, papillae type, density, height, ribbing, surface deposits, and peaks topped with wax were examined under scanning electron microscopy. Various gross-morphological traits of seeds of Vicia species have been analyzed and compared. The present study revealed significant variation in testa traits. Testa were papillose and papillose with mounds, the latter being observed only in Vicia lathyroides. The present study revealed 20 key traits that could be used to diagnose Vicia species and classify them.