INTRODUCTION
Job’s-tears (Coix lacryma-jobi L.) is an economically important annual crop mainly planted in Asian countries (Diao, 2017). It belongs to the Coix genus, Andropogoneae tribe within the family Poaceae, and is a tall grain-bearing tropical plant (Taylor, 1953; Zhou et al., 2010). Job’s-tears is also known as coix seed, tear grass, hato mugi, and adlay. The demand for Job’s-tears is rapidly increasing due to its medicinal functions as an ideal healthy food. It has been introduced and grown in almost all the tropical and subtropical areas in the world (Lim, 2012).
Inappropriate cultivation techniques coupled with the long-term use of landrace varieties have caused low seed quality and production. In addition, a few research studies about Job’s-tears have been conducted and led to a lack of new varieties with high yield and quality (Huang et al., 2009). Genetic resources are the most important part of a plant breeding program, and diversity analysis is a prerequisite for their more efficient management and utilization. Accurate identification of the genotype is very important during all the plant breeding steps, from initial parent selection to the final utilization of cultivars in production schemes (UPOV, 1991). Moreover, diversity analysis is an essential process to clearly identify the genetic relatedness of the available genetic resources. Modern objectives in plant breeding might be achieved by trait evaluation in genetic resources. Molecular marker methods for genotype description have proved useful, but these methods are expensive and need marker-linked trait information. Morphological characters must be recorded to select parents and their progenies, and they must always be used to describe and classify germplasm. Principle component
and cluster analysis are also useful tools to screen accessions (Karimi et al., 2009). To better conserve and utilize genetic resources, appropriate morphological variables should be carefully chosen (Giraldo et al., 2010). There are few studies that have reported the morphological analysis of Job’s-tears accessions. Ten major quantitative traits of 12 Job’s-tears germplasms from different regions were analyzed by path and principal component analysis (Li et al., 2010). The main components and cluster analysis of 25 Job’s-tears accessions were also conducted in 2013 (Wang et al., 2013). So far, studies on the genetic diversity of Job’s-tears have been conducted on accessions mainly originating from the Guangxi, Guizhou, and Yunnan Provinces in southern China (Li et al., 2001; Ma et al., 2010; Wang et al., 2015); however, no systematic study about the morphological diversity of Job’s-tears has been reported. Hence, the objectives were to analyze the genetic diversity of 94 Job’s-tears accessions from different regions in China and clarify the genetic relationships among the accessions to provide valuable information for the selection of excellent accessions as parents and utilize and conserve the Job’s-tears genetic resource.
MATERIALS AND METHODS
Plant material
The 94 Job’s-tears (Coix lacryma-jobi L.) accessions evaluated in this study are shown in Table 1. Accessions consisted of 40 cultivated and 54 wild accessions that were collected from different geographic regions in seven Chinese provinces (36 Guizhou, 5 Guangxi, 7 Sichaun, 7 Chongqing, 7 Hunan, 19 Yunnan, 4 Hubei, and 9 Qianxinan Institute of Agricultural Sciences accessions). All Job’s-tears accessions were collected and multiplied during 2013-2014.
Thirty seeds from each Job’s-tears accession were planted in seedling cups. Seedlings with 3-4 leaves were transplanted to the field at the Germplasm Resources Garden of Guizhou Academy of Agricultural Sciences (26°57’ N, 106°71’ E; 1074.3 m a.s.l.), Guizhou province, China, during April 2015. A completely randomized design (CRD) with three replicates was used. Each replicate consisted of 10 plants with 1 plant per hill. The row and plant spacing were 60 and 40 cm, respectively.
Morphological traits
Twelve different quantitative and qualitative traits were evaluated at harvest time. Five plants (normal growth, uniform performance, disease- and insect pest-free) from each accession were randomly selected for scoring. Seven quantitative characters, including stem node number (SNN), panicle branch number (PB), primer branch nodes (PBN), panicles per plant (PP), and grain number per plant (GNP) were evaluated. Plant height (PH) was measured from ground to tip of the main spike at physiological maturity. The 100-seed weight (SW) was measured with a scale sensitive to 0.1 g. Five qualitative characters, including total bract surface characteristics (TBSC), total bract texture (TBT), total bract shape (TBS), total bract color (TBC), and pericarp color (PC) were also recorded.
Data analysis
For the seven quantitative traits, descriptive statistics were computed for each accession with the Statistical Package for Social Science (SPSS 20.0; IBM, Armonk, New York, USA). The mean, maximum, minimum, standard deviation (SD), and coefficient of variation (CV) were calculated for each of the seven quantitative traits (Table 2). The five qualitative traits were scored according to Li et al. (2015): TBSC: smooth = 1, longitudinal convex stripes = 2; TBT: enamel = 1, crustaceous = 2; TBS: circular = 1, oval = 2, length circle = 3; TBC: white = 1, grayish-white = 2, blue-gray without dark stripes = 3, blue-gray with dark stripes = 4, yellow-white = 5, tawny without dark stripes = 6, tawny with dark stripes = 7, brown = 8, dark brown = 9; and PC: light yellow = 1, yellow = 2, brown = 3. Principal component analysis (PCA) was carried out using the original numerical data of quantitative characters and the assigned qualitative trait data. The Euclidean distance between Job’s-tears accessions, based on morphological characteristics, was calculated by the between-group linkage method in SPSS 20.0.
RESULTS
Quantitative trait analysis
Accessions Y35 (263.4 cm), Y41 (263.3 cm), Y63 (239.4 cm), and Y77 (239.8 cm) had the highest PH values, while Y16 (146.0 cm), Y30 (148.7 cm), and Y84 (59.6 cm) had the lowest PH (Table 2). The highest SNN values were observed in Y3 (13.0), Y23 (12.8), and Y41 (13.4), whereas the lowest SNN values were in Y16 (7.7) and Y39 (7.8). The highest PB was observed in Y25 (26.6), followed by Y36 (16.4) and Y74 (15.7), while the lowest PB was recorded in Y10 (5.0) and Y16 (5.0). The maximum PBN were found in Y35 (8.5), Y63 (8.5), and Y71 (8.7), whereas the minimum value was obtained in Y30 (2.7) and Y52 (1.7). The highest PP values were recorded in Y77 (120.4) and Y88 (117.0), while the lowest PP were in Y7 (19.9) and Y46 (19.9). The maximum GNP was observed in Y92 (316.7), followed by Y68 (302.8), Y3 (279.4), Y40 (273.5), and Y72 (266.9), while the lowest GNP values were in Y44 (51.4) and Y46 (45.0). The SW of all accessions ranged between 6.7 g and 25.3 g. The Y36 accession had the heaviest SW (25.3 g), followed by Y5, Y28, and Y60 with 19.7, 17.9, and 16.1 g, respectively. The lightest SW were observed in Y33 (6.8 g), Y38 (6.7 g), and Y49 (6.8 g).
Table 2 Mean of seven quantitative morphological traits of 94 Job’s-tears accessions.

PH: Plant height; SNN: stem node number; PB: panicle branch number; PBN: primer branch nodes; PP: panicles per plant; GNP: grain number per plant; SW: 100-seed weight.
There were significant genetic variations in all the quantitative traits of cultivated and wild accessions (Table 3). For cultivated species, the PH mean value was 198. 6 cm and ranged between 59.6 and 263.4 cm, SD was 32.4 cm, and CV was 16.3%. For the wild species, PH mean value was 200.9 cm and ranged between 148.7 and 263.3 cm, SD was 25.4 cm, and CV was 12.6%. Thus, the variation of cultivated species PH was greater than for wild species. Although the maximum value of cultivated species was similar to that of wild plants, Y35 (263.4 cm) and Y41 (263.3 cm), the minimum was quite different at 56.87 and 150 cm, respectively. Wild species mean variations were higher than for cultivated species in PB, PP, and SW, while SNN and PBN were not significantly different. For GNP, the mean of wild species was higher than for cultivated species at 165.7 and 135.2, respectively.
Qualitative trait analysis
Qualitative traits of Job’s-tears are shown in Tables 4 and 5. The bract surface characteristic in 80 accessions was smooth (85.1% of the total), and the remaining 14 accessions had longitudinal convex stripes (14.9%). There were 40 accessions with crustaceous bract texture (42.6% of the total), whereas 51 accessions had an enamel texture (57.4%). As for bract shape, five accessions had a circular shape (5.3%), 54 accessions were oval-shaped (57.4%), and 35 accessions length circle shape (37.3%). Regarding TBC, there were 27 white accessions, 6 grayish-white, 3 blue-gray without dark stripes, 3 blue-gray with dark stripes, 1 yellow-white, 4 tawny without dark stripes, 6 tawny with dark stripes, 19 brown , and 25 dark brown that accounted for 28.6%, 6.4%, 3.2%, 3.2%, 1.1%, 4.3%, 6.4%, 20.2%, and 26.6%, respectively. With respect to PC, there were 48 light yellow accessions (51.1%), 12 yellow (12.8%), and 34 brown (36.1%).
Table 4 Qualitative morphological traits of 94 Job’s-tears accessions.

TBSC: Total bract surface characteristics; TBT: total bract texture; TBS; total bract shape: TBC: total bract color; FT: fruit type; PC: pericarp color.
Correlation coefficient analysis
The Pearson correlation coefficient revealed significant correlations among some variables measured in Job’s-tears accessions (Table 6). The PC, TBC, TBS, TBT, and TBSC were positively correlated as well as PH, PBN, and SNN. The PH, PP, and GNP were also positively correlated one with the other. The TBSC and TBT were correlated with PB, PP and GNP. The PB was also positively correlated with TBC and TBS, whereas GNP was correlated with SW.
Principal component analysis
To fully reflect the various factors that played a principal role in the comprehensive indicators, PCA was carried out on seven quantitative traits and five qualitative traits. The Eigenvalues, contribution rate, and accumulative contribution rate were also gathered (Table 7). According to the 85% accumulative contribution rate standard, most agronomic trait information can be summarized in the first seven PCs of the experiment. The accumulative contribution rate accounted for 87.31% of the total variation.
Table 6 Simple correlation matrix for 12 morphological traits of 94 Job’s-tears accessions.

*, **Significant at p = 0.05 and p = 0.01, respectively.
PH: Plant height; SNN: stem node number; PB: panicle branch number; PBN: primer branch nodes; PP: panicles per plant; GNP: grain number per plant; SW: 100-seed weight; TBSC: total bract surface characteristics; TBT: total bract texture; TBS: total bract shape; TBC: total bract color; PC: pericarp color.
Table 7 Eigenvalues, proportion of variance, and morphological traits that contributed to first seven principal components (PCs).

The TBSC, TBT, and TBS were highly loaded in PC1. This indicated that PC1 reflected the main factor of involucre plant parts with an eigenvalue of 3.60, and the contribution rate was 30.03% of the total variation of the studied samples. In PC2, the eigenvalue was 2.15 and the contribution rate was 17.92% of the total morphological variation; among the accessions, these values were mainly explained by PH, SNN, and PBN. This indicated that PC2 was reflected as the main factor of plant type. In PC3, the eigenvalue was 1.44 and the contribution rate was 12.03% of the total variation; GNP and PP was more loaded in PC3. This indicated that PC3 reflected the main factor of yield. In PC4, the eigenvalue was 1.05 and the contribution rate contributed 8.74% of the total morphological variation in these accessions in only SW. The PC5 was 7.69% of the total variation with TBS highly loaded. The PC6 and PC7 accounted for 6.23% and 4.68% and mainly loaded PC and TBC, respectively. In general, for the 12 morphological traits studied, PC1 and PC2 constituted 47.5% of the total morphological variation with most seed-related traits and vegetative traits.
Cluster analysis
According to the cluster analysis results (Figure 1) and the mean value of each group (Table 8), tested Job’s-tears accessions were classified in different groups with significant differences in the morphological characteristics. Based on the genetic distance of 10.5, 94 Job’s-tears accessions were grouped into seven major clusters.
Table 8 Arithmetic mean for traits of each group.

PH: Plant height; SNN: stem node number; PB: panicle branch number; PBN: primer branch nodes; PP: panicles per plant; GNP: grain number per plant; SW: 100-seed weight.
Cluster I contained 21 accessions from Guizhou (7), Yunnan (6), Sichuan (3), Qianxinan Institute (3), Hebei (1), and Guangxi (1). The Y22, Y85, Y83, Y76, Y88, Y23, Y25, Y17, Y18, Y63, Y74, Y35, Y41, Y51, Y64, Y90, Y58, Y67, Y87, Y73, and Y69 accessions were classified in this cluster. The main features were PH between 188.8 and 263.4 cm, SNN between 9 and 14.4, PB between 5.5 and 25.5, PBN between 3.8 and 8.5, PP between 42.2 and 117.0, GNP between 136.0 and 227.8, and SW between 7.4 and 12.6 g. This group had relatively high PH, SNN, PB, and PP values (Table 8).
Cluster Ⅱ contained only Y77. It was a cultivated variety collected in Maguan in Yunan province, which had the highest PH (239.8 cm), SNN (11.2), and PP (120.4) and higher SW (13.4 g).
Cluster Ⅲ consisted of 4 accessions from Guizhou (1), Yunnan (1), Sichuan (3), Hebei (1), and Chongqing (1). The Y11, Y52, Y86, and Y17 accessions were classified in this cluster. The main features were PH between 155.4 and 195.0 cm, SNN between 8.3 and 12.0, PB between 6.4 and 7.7, PBN between 1.7 and 4.4, PP between 19.9 and 39.5, GNP between 237.0 and 266.0, and SW between 8.1 and 15.2 g. This group had relatively low PH, PB, PBN, and PP, but it had relatively high GNP and SW.
Cluster Ⅳ had 10 accessions from Guizhou (2), Yunnan (3), Sichuan (2), Qianxinan Institute (2), and Chongqing (1). This group had the highest PB (11.5) and GNP (270.5), whereas SW (8.8) was the lowest.
Cluster Ⅴ only contained Y10 and Y18 collected in two different adjacent counties (Changshun and Puding) in Guizhou province. This group had the highest PBN (6.5). Although this group had lower GNP, with only 78.9, it had the heaviest SW (15.7 g).
Cluster Ⅵ had 55 accessions, and this group was collected in Guizhou (24), Yunnan (7), Sichuan (2), Qianxinan Institute (4), Hebei (2), Guangxi (4), Hunan (7), and Chongqing (5). All the characteristics of this group were at a relatively stable intermediate level.
Cluster Ⅶ contained only Y84, which was a cultivated variety collected in Quejing in Yunnan province; it had the shortest PH (59.6 cm), the lowest SNN (7.4), PBN (1.8), PP (25.6), and GNP (55.5). In addition, PB (7.7) and SW (9.1) were also relatively lower.
DISCUSSION
Germplasm resource collections are an important step to breed improved crops (Nelson et al., 2011; Andini et al., 2013). Excavating the excellent Job’s-tears germplasm was the most important task to improve the species. The results of this study showed that both cultivated and wild accessions have large variations in PB, PBN, PP, and GNP. These accessions have abundant diversity and a great range of optional resources for breeding. The cultivated variety Y36 had 25.3 g SW more than previous findings of 12.5-21.5 g by Wang et al. (2013); it could be used as excellent gene resources to improve varieties with a good combination of agronomic materials such as higher PB and medium PH. The Y92, Y68, Y3, Y40, and Y72 varieties could be used as materials for a high yield variety with more GNP. The cultivated variety Y30 and wild variety Y84, with lower PH, could be used as excellent gene resources for new dwarf varieties. The Y77 accession had higher PH, SNN, PB, GNP, SW, and the highest PP. It was an excellent resource with better comprehensive traits.
Grain yield was closely associated with GNP at harvest (Wang et al., 2013). The analysis of simple correlations among traits revealed that accessions with high PH, SNN, PBN, and PP had high GNP, and positive correlation among accessions suggested that these traits could be used as selection criteria to evaluate accessions. The PP was also positively and significantly related to GNP, indicating that increasing PP increased GNP. The PP is a good trait for selecting excellent accessions. There were no quantitative traits significantly correlated to SW except GNP, but it was a negative correlation. Generally, when GNP was higher, SW was lighter; therefore, we cannot judge if an accession is good or not only using SW. The qualitative traits had a positive or negative correlation with one another. Some quantitative traits were significantly correlated one another; for example, PB was significantly correlated with TBSC, TBT, TBS, TBC, the PP, and GNP. Positive correlations were recorded for TBSC and TBT. However, the reasons for such a correlation between quantitative traits and qualitative traits are still unclear and need to be studied further.
It is more difficult to analyze the multi-index problems with a number of indicators related one to another. The PCA can simplify multi-index analysis by converting original and more related indices into a new index. Previous studies show that PCA was an effective method for comprehensive crop evaluation (Wang et al., 2013). In the present study, 12 traits of the 94 Coix accessions were reduced to seven main components using the PCA method and with the accumulative contribution rate up to 87.31%. The PC1, PC2, and PC3 reflected the effective main factor of involucre, plant type, and yield, respectively. The PC4, PC5, PC6, and PC7 reflected the main factor of SW, TBS, PC, and TBC, respectively. Therefore, 12 trait indices were reduced to seven comprehensive indicators used to represent the original variables, simplify the data, and reveal the relationship between the variables. They can also provide a favorable scientific basis for parent selection in breeding programs (Li et al., 2015). It was also found that PB is one of the most important yield indicators, but it was summarized in PC1 instead of PC2. Therefore, it was necessary to combine the original data and use dialectical analysis methods to remove the apparent phenomenon when we evaluated the quality of germplasm resources by PCA. In addition, breeders must take advantage of the main factors and expand the different traits of breeding materials to accelerate breeding programs for new varieties.
Genetic diversity analysis of germplasms using morphological traits is an initial step for crop improvement (Julia et al., 2016; Peratoner et al., 2016; Loumerem and Alercia, 2016). The variations in morphological traits can be used to classify materials in different groups. The shape, color, and texture of seeds are important to classify Coix species (Schaaffhausen, 1952; Rao and Nirmala, 2010). The present study combined quantitative and quantitative traits. The 94 accessions were grouped into seven clusters. Cluster I contained 21 accessions from Guizhou (7), Yunnan (6), Sichuan (3), Qianxinan Institute (3), Hebei (1), and Guangxi (1). The accessions in this group had different plant characteristics that are good for breeding material. Cluster II contained only Y77, which had the highest PH SNN, PP, and SW: this accession can be used as a high-yielding variety and parental material to improve a superior dwarf variety. Cluster Ⅲ contained 4 accessions, which can be used as excellent resources to promote dwarf, anti-lodging, and high yield varieties. Cluster Ⅳ had the highest PB and PP and can be used as excellent parental material to improve crop yield. Cluster Ⅴ contained only Y10 and Y18 collected from two different adjacent locations, Changshun and Puding in Guizhou province; this indicated that accessions can cluster with geographic locations. Cluster Ⅵ was the largest group consisting of 55 accessions. Cluster Ⅶ contained only Y84 that had the shortest PH (59.6 cm). It can be used as excellent parent material to cultivate new dwarf and anti-lodging varieties.
Clustering results showed that the accessions collected from different areas could be grouped, such as Y11, Y52, Y86, and Y7 collected from Guizhou, Yunnan, Sichuan, Hebei, and Chongqing, respectively. Some accessions collected from the same area might not be grouped, such as Y35 and Y36 in Cluster Ⅰ and cluster Ⅵ, respectively, which were collected from Ceheng in Guizhou province. This indicated that genetic differences were not influenced by geographic differences, and accessions collected from the same area can be very different as to genetic variation (Li et al., 2010). Genetic differences were based on variety type (wild or cultivated) and qualitative traits were inconsistent with the results of clustering by geographic distribution (Xi et al., 2016). This is because agronomic traits are easily affected by environmental conditions and cultural practices (Bruschi et al., 2003; Wang et al., 2013).
CONCLUSIONS
Morphological diversity analysis in the present study showed high variations among the materials according to ANOVA and simple correlations and multivariate analysis. Significant and positive correlations were found between grain number per plant and among other yield-related attributes. The reasons for such a correlation between quantitative and qualitative traits are still unclear and need further study. However, the information about the correlation of characteristics can be used for indirect selection and as breeding criteria. The results also indicated that these accessions have the potential to improve yield in the Job’s-tears breeding program. Principal component analysis showed that the observed variations were mainly caused by traits such as total bract surface characteristics, total bract texture, total bract shape, plant height, stem node number, and primer branch nodes, and this indicated that their contribution is important to discriminate accessions. Cluster analysis grouped the 94 Job’s-tears accessions into seven clusters. This indicated high diversity for most of the traits, and demonstrated that genetic differences cannot only be based on geographic differences, variety type (wild or cultivated), and qualitative traits. Since morphological traits are easily affected by environmental conditions, the genetic relationship cannot only be reflected by the similarity between morphological characteristics.