Non-destructive phenotyping for early seedling vigor in direct-seeded rice
Plant Methods volume 16, Article number: 127 (2020)
Early seedling vigor is an essential trait of direct-seeded rice. It helps the seedlings to compete with weeds for water and nutrient availability, and contributes to better seedling establishment during the initial phase of crop growth. Seedling vigor is a complex trait, and phenotyping by a destructive method limits the improvement of this trait through traditional breeding. Hence, a non-invasive, rapid, and precise image-based phenotyping technique is developed to increase the possibility to improve early seedling vigor through breeding in rice and other field crops.
To establish and assess the methodology using free-source software, early seedling vigor was estimated from images captured with a digital SLR camera in a non-destructive way. Here, the legitimacy and strength of the method have been proved through screening seven diverse rice cultivars varying for early seedling vigor. In the regression analysis, whole-plant area (WPA) estimated by destructive-flatbed scanner (WPAs) and non-destructive imaging (WPAi) approaches was strongly related (R2 > 83%) and suggested that WPAi can be adapted in place of destructive methods to estimate seedling vigor. In addition, this study has identified a set of new geometric traits (convex hull and top view area) for screening breeding lines for early seedling vigor in rice, which decreased the time by 80% and halved the cost of labor in data observation.
The method demonstrated here is affordable and easy to establish as a phenotypic platform. It is suitable for most glasshouses/net houses for characterizing genotypes to understand the plasticity of shoots under a given environment at the seedling stage. The methodology explained in this experiment has been proven to be practical and suggested as a technique for researchers involved in direct-seeded rice. Consequently, it will help in the simultaneous screening of genotypes in large numbers, the identification of donors, and in gaining information on the genetic basis of the trait to design a breeding program for direct-seeded rice.
The benefits of decreasing the water footprint along with less labor use and an increase in the cost–benefit ratio have led rice farmers to shift their puddled-transplanted rice to direct-seeded rice (DSR). For successful crop establishment under DSR, rapid uniform emergence and accumulation of biomass in the early phase of crop growth are the key factors . Thus, understanding the spatio-temporal changes in shoot biomass in the early phase of the crop by imaging would help to differentiate lines for vigor and provide insight into the physiology of rice seedlings under direct-seeded conditions [2, 3]. Therefore, developing an automated non-destructive screening method for an essential agronomic trait would enhance the productivity of rice under direct-seeded conditions.
Non-destructive phenotyping techniques are the key factors for screening and developing suitable rice genotypes for the target environment in a brief period [4, 5]. Several phenotyping methods have been optimized for screening biotic and abiotic stresses . The absence of a suitable non-destructive-based high-throughput phenotyping system has restricted the exploitation of agronomically important traits in rice. Accumulation of biomass in the early phase of crop growth is necessary under DSR and is considered an important parameter for seedling vigor [1, 7]. Therefore, a genotype with early seedling vigor has significance in smothering the effect of weed competitiveness and water use efficiency to maintain the sustainability of rice production in rainfed and DSR conditions. The key limiting factor in developing rice varieties for direct-seeded conditions with early seedling vigor is the non-availability of a suitable non-destructive phenotyping screening technique to select genotypes against weed competitiveness with a consistent result. As biomass/shoot weight was observed to be closely related to early vigor measured at 14 and 28 days after sowing (DAS) , screening genotypes by measuring biomass with a non-destructive method has added advantages over other methods.
The early phase of a crop is more fragile and dynamic in response to the environment, and the complex nature of the trait poses a problem in phenotyping for early seedling vigor . Further, phenotyping by destructive sampling and collecting of seeds from the same individuals, in the case of segregating generations, would be an additional impossible target. Recent advances in genomic technologies have changed the way of breeding programs by generating more genotypic data. Nevertheless, the same breeding programs failed to achieve the objective when the translation of such data failed to identify genotypes with the desirable trait . Therefore, to overcome this bottleneck and to use those genotypic data in an efficient way, non-destructive phenotyping with precision is highly valued. In recent years, several reports have presented the advantages of non-destructive phenotyping by imaging techniques using near-infrared reflectance and spectral imaging using fluorescence and thermal wavelengths [10,11,12]. The application of image-based phenotyping is picking up in the area of field crops to understand the complex traits that are highly influenced by the environment. Several screening protocols and pipelines for data analysis were developed for some of the intricate stress factors such as salinity, nitrogen, and water deficiency, and nodal root angle in barley, rice, and sorghum, etc. [3, 9, 13,14,15,16,17,18,19,20]. Conversely, imaging techniques for early seedling vigor have not been standardized in rice. Most of the protocols developed by imaging are automated and require high-end facilities. Irrespective of the countries possessing those high-end automated facilities, they are not affordable to all researchers.
The existing field-based screening methodology to estimate seedling vigor is based on harvesting samples over relative time , which is labor-intensive. As the early phase of seedlings is dynamic in nature and in collecting data to estimate growth analysis, biomass or leaf area index from a subset of the population would not provide reasonable information. However, whole-plant area (WPA) is associated with seedling vigor, but the destructive method cannot capture the actual area. Therefore, a phenotyping screening protocol needs to be developed in a cost-effective way that is easy to handle, less labor-intensive, suitable for screening year-round, and amenable to integrating those phenotypic data with genotypic data generated from forward genetic tools such as genome-wide association mapping, linkage mapping, and gene sequencing. In addition, a non-destructive image-based phenotyping protocol should be flexible and experimentally verified by comparing it with existing field-based techniques. On the other hand, high-throughput phenotyping integrated with imaging techniques would be more flexible to capture the dynamic changes taking place in plants over a time interval. This would decrease genotype x environment interaction and several parameters such as compactness, leaf rolling, and drying related to abiotic stress and leaf damage due to pests and diseases would be measured seamlessly .
In the present study, we focused on establishing a non-destructive phenotyping protocol to estimate early seedling vigor in rice using images. Seven rice (Oryza sativa L.) genotypes of improved and traditional lines were grown in a pot under normal conditions without stagnation of water. Growth rates and related agro-morphological traits of those genotypes were measured by proposed non-destructive image-based and conventional destructive harvests to test the protocol efficiency, reproducibility, and ability to differentiate vigorous genotypes.
Plant growth and partitioning of biomass among genotypes at 14 and 28 DAS
Significant differences were observed among the seven genotypes for traits studied at 14 and 28 DAS. On the 14th day after sowing, 16 traits exhibited significant differences among the 28 traits studied, while 19 traits showed significant differences between genotypes at 28 DAS (Table 1). However, traits observed after manual sampling such as shoot length and leaf number per plant were found to be significant across the two dates of observation (Tables 1, 2). Third-leaf width, eccentricity, convex hull, caliper length, whole-plant area by destructive-flatbed scanner (WPAs), whole-plant area by non-destructive imaging (WPAi), top view area, and compactness exhibited strong significant differences across the two dates of observation (Tables 3, 4) and these differences were captured well by the images from all seven genotypes grown under net house conditions. This suggests that observing growth parameters by imaging could capture subtle differences across genotypes, which is not possible in the traditional way of screening genotypes by destructive sampling. Among the seven genotypes, GM-217 and Vandana registered higher values for growth parameters such as shoot length, leaf number, WPAi, convex hull, and compactness at 14 DAS. On the other hand, LB-46 and Varshadhan were observed to have maximum leaf number, WPAi, and convex hull at 28 DAS. Traits such as root length, shoot dry weight, root dry weight, tiller number, and stem weight displayed a significant difference at 28 DAS, whereas these traits displayed non-significant growth at 14 DAS.
Comparison between destructive sampling and automated image phenotyping (WPAi) for seedling vigor
Seedling vigor was generally predicted conventionally through growth analysis. Growth analysis was calculated as the absolute growth rate (AGR), crop growth rate (CGR), and relative growth rate (RGR), which are destructive methods (explained in detail under the Methods section (Method-2 & -3). All three of these growth parameters at 14 and 28 DAS were plotted against WPAi (non-destructive method) using regression curves to find out the association and contribution in variation for WPAi. In the present experiment, WPAi is considered a principal parameter from image-based phenotyping, since WPA is the target trait measured through automated image analysis. AGR was calculated from shoot and root length, while CGR and RGR were calculated based on the dry weight of shoots with tillers and roots. The AGR for shoot length and root length plotted against the AGR of WPAi exhibited no relationship with the AGR of WPAi with nearly zero regression (≤ 0.04) (Fig. 1a, b). On the other hand, CGR of shoot dry weight with tillers and root dry weight exhibited strong and positive correlation with CGR of WPAi (Fig. 2a, b). The CGR of shoot dry weight with tillers explained 74.26% of the variation and the CGR of root dry weight explained 45.20% of the variation of CGR-WPAi. Similarly, RGR of shoot dry weight with tillers and root dry weight showed a positive relationship with the RGR of WPAi (Fig. 3a, b). However, only the RGR of shoot dry weight had a positive correlation and explained 13.80% of the variation of the RGR of WPAi. Thus, the CGR of shoot and root dry weight and RGR of shoot dry weight had a positive relationship with CGR of WPAi and RGR of WPAi, respectively. Another method of destructive sampling was WPAs, which was a more precise method than the phenotypic measurement for WPAi estimation. Both WPAi and WPAs were plotted on the graph, where WPAs was plotted on the X-axis against WPAi on the Y-axis as a dependent variable. The correlation of WPAs at both 14 and 28 DAS was strong and positive with WPAi. The regression showed that WPAs explained 83.11% of the variation at 14 DAS and 87.33% of the variation at 28 DAS of WPAi (Fig. 4). This was the strongest relationship among all the traits observed.
Relationship between WPAi and other seedling traits
The relationship of seedling phenotypic traits, individual leaf traits, stem area, and WPAs with WPAi was calculated to understand the correlation and percentage of variation contributed by each trait (morphological traits and geometric traits) toward WPAi. As WPAi was used as a trait in the identification of early seedling vigor in plants, it was considered as a primary parameter in the image-based phenotyping method.
Relationship of seedling phenotypic traits and geometric traits with WPAi
Both positive and negative correlation was observed between seedling traits and WPAi. At 28 DAS, WPAi had exhibited a strong positive relationship with morphological traits such as the number of leaves (r = 0.95, p < 0.01), shoot and root dry weight, and tiller number (Table 5). In general, no strong negative association was observed other than with some specific leaves at 28 DAS (presented in the following section). Similarly, geometric traits that related to the size of the plant such as caliper length (r = 0.88, p < 0.05), convex hull, and top view area had a strong and positive correlation with WPAi at 28 DAS. Parallel to the data observed at 28 DAS, traits observed at 14 DAS also exhibited a similar trend of relationship (Table 5). The contribution of variation by shoot length explained 90.4% of the WPAi variation at 14 DAS, while it had a negative contribution of 3.6% at 28 DAS. At 28 DAS, 91.20% of the variation in WPAi was explained by the number of leaves (Table 5).
Geometric traits such as caliper length, eccentricity, convex hull, and top view area explained 48.3%, 27.3%, 25.2%, and 1.9% of the variation, respectively, at 14 DAS in a negative direction. Conversely, at 28 DAS, they explained the variation positively (Table 5). Overall, seedling phenotypic traits and geometric traits were highly correlated with WPAi and thus can be predictable using linear regression.
Relationship of the different leaves with WPAi
The relationship between WPAi and individual leaf length, width, and area at 14 and 28 DAS was estimated using a linear regression curve (Table 5). There were three leaves per plant at 14 DAS and seven leaves at 28 DAS in all the genotypes. Some genotypes produced a fourth leaf at 14 DAS and some produced an eighth leaf at 28 DAS. Therefore, the fourth leaf (at 14 DAS) and eighth leaf (at 28 DAS) were excluded from the analysis. Both positive and negative regression were observed between WPAi and individual leaf length. Among all the leaves, third leaf length had a positive (0.89, p < 0.01) correlation with WPAi, which explained 77.64% of the WPAi variation at 14 DAS. Similarly, length of the fifth-seventh leaf at 28 DAS had a positive association and explained > 40% of the variation (Table 5).
With reference to leaf width, no strong relationship between WPAi and leaf width was observed at both 14 and 28 DAS, although the highest variation was explained by seventh leaf width (35.33%) at 28 DAS (Table 5). In the case of leaf area, the second-third leaves had a positive relationship with WPAi at 14 DAS and explained 51.29% of the WPAi variation. For 28 DAS, the fifth-seventh leaves had a positive relationship with WPAi, which was similar to the results obtained for leaf length. However, the leaf area of the fifth and sixth leaves had a strong correlation, which explained the extent of variation (77.79% and 53.25%, respectively) of WPAi. For leaf weight, all leaves had a positive correlation (Table 5). Among them, the second and third leaves of 14-day-old seedlings explained variation of more than 40%. On the other hand, the leaf weight of 28-day-old seedlings at different levels had a weak correlation with WPAi. Of these, the sixth and seventh leaves explained variation of > 15%.
Relationship of stem area with WPAi and stage-specific traits with WPAi
The greenness in the stem also helps in photosynthesis, which contributes to the overall growth of seedlings. Hence, stem area was also measured to find out the relationship with WPAi. The relationship between stem area and WPAi was positive (0.84, p < 0.01) at 14 DAS and negative (0.12ns) at 28 DAS. At 14 DAS, a strong relationship was observed, which explained 61.33% of the variation, while at 28 DAS the correlation was negative (1.59%) and very weak. Some of the unique traits in correlation to seedlings were also measured, which were seedling age-specific. Traits such as seed weight with mesocotyl, seed weight, and mesocotyl length were measured at 14 DAS and bulk eighth leaf (terminal) weight, stem weight, and tiller number were measured at 28 DAS. The traits that were measured at 14 DAS were not available to measure at 28 DAS and vice versa. The traits seed weight with mesocotyl, seed weight, and mesocotyl length were negative and had a weak correlation with WPAi at 14 DAS, whereas bulk eighth leaf weight, stem weight, and tiller number showed a positive relationship with WPAi at 28 DAS. Single stem weight and tiller number had a strong relationship with WPAi and explained 82.30% and 79.80% of the variation of WPAi, respectively.
Relationship between morphological traits and geometric traits
Understanding the trait association between morphological and geometric traits observed by image analysis would be helpful in the identification of surrogate traits in the absence of an automated non-destructive imaging system. At 28 DAS, number of leaves per plant expressed a strong positive association with top view area (0.89, p < 0.01), WPAi (0.95), caliper length (0.91, p < 0.01), convex hull (0.83), root dry weight (0.81, p < 0.05), shoot dry weight (0.93, p < 0.01), tiller number (0.97, p < 0.01), and stem weight (0.88, p < 0.01). The RGR-related trait compactness derived from the differences in top view area and convex hull exhibited a strong positive association with leaf width irrespective of all levels and a strong negative association with shoot length, eccentricity, convex hull, and caliper length on both dates of observation. Further, it has a weak negative association with leaf number (− 0.20 at 14 DAS and − 0.28 at 28 DAS). The parameter convex hull displays the degree of leaves spreading that helps to cover the ground. The number of leaves showed a strong positive association with convex hull on both observation dates (0.55 at 14 DAS and 0.83 at 28 DAS) and a negative association with leaf width at all levels.
Identification of genotypes with high seedling vigor
Generally, to identify genotypes with high seedling vigor, AGR, CGR, and RGR were used. Among these, CGR was commonly used to find vigorous genotypes. The CGR, RGR, and AGR of genotypes were compared to the WPAi of the respective growth rates. The CGR of shoot weight and root weight, AGR of shoot length and root length, and RGR of root dry weight and shoot dry weight were plotted against the CGR of WPAi, AGR of WPAi, and RGR of WPAi, respectively. The highest growth rate was observed in LB-46 (CGR: 0.02 (shoot dry weight) and 0.01 (root dry weight); RGR: 0.15 (shoot dry weight) and 0.19 (root dry weight)), considered as a genotype with high seedling vigor, followed by GM-217 and Varshadhan, based on the destructive method (CGR and RGR) for shoot and root dry weight. Meanwhile, AGR showed that ARC10656 and AC38399 were superior to other genotypes.
Image analysis by the non-destructive way has identified LB-46 (CGR of WPAi: 212.36; RGR of WPAi: 0.09) as a genotype with high seedling vigor, followed by Varshadhan. In both methods, LB-46 is common. Thus, WPAi is certainly an alternative to the destructive method (Fig. 4). On the basis of imaging and scanning methods (WPAi and WPAs) at 28 DAS, the highest WPAi was gained in order as LB-46 (4068 mm2), Varshadhan, and GM-217, while LB-46, GM-217, AC38399, and Varshadhan were judged as top genotypes by the WPAs method (Table 4) (Fig. 4). Overall at 14 DAS, both WPAi and WPAs identified GM-217 (1498 mm2 (WPAi)/1264 mm2 (WPAs)) and Vandana as top contenders, followed by LB-46 in WPAi and AC38399 in WPAs as the next best genotypes (Table 3).
Grouping pattern of genotypes and association between variables
Principal component analysis (PCA) was employed for 29 traits observed at 14 DAS, for which it has explained 98.99% of the variation by PC1 and 0.96% by PC2 (Fig. 5a). On the basis of magnitudes of loadings/eigenvalues, nine highly variable traits (first and third leaf area, caliper length, convex hull, eccentricity, stem area, top view area, WPAi, WPAs) were identified. Similarly, PCA was performed for 44 traits at 28 DAS, which governs 99.49% of the variation on the PC1 axis and 0.47% on the PC2 axis (Fig. 5b). Out of 44 traits on the basis of magnitude of PCA, fifth, sixth, and seventh leaf area; sixth and seventh leaf length; convex hull; stem area; top view area; WPAi; and WPAs were identified as highly variable traits. Among them, convex hull, stem area, top view area, WPAi, and WPAs were highly variable common traits between the two dates of observation.
The cultivar-by-trait biplots (Fig. 5a, b) were analyzed for both 14 and 28 DAS. They showed a strong relationship between WPAi and WPAs at both 14 and 28 DAS. On both dates, biplot classifies the traits into two groups, based on their association with whole-plant area (WPAi and WPAs). At 14 DAS, traits such as width, area, and dry weight of first, second, and third leaves; leaf length of first and second leaves; and root length were clustered together. Conversely, traits that related to the image that contributes directly to WPA such as shoot length, caliper length, eccentricity, convex hull, third leaf length, stem thickness, shoot dry weight, stem area, top view area, and leaf number were grouped together. The vector of compactness stayed away from both groups. A similar trend was also observed at 28 DAS. Traits such as leaf length, width, weight, and area of 1–4; leaf dry weight of 5–8; stem area; root length; and compactness were grouped together. On the other hand, geometric and morphological traits that contributed to WPA such as shoot length; seventh leaf length, width, and weight; shoot length; shoot dry weight; stem weight; caliper length; convex hull; eccentricity; leaf length, width, and area of fifth and sixth leaf; top view area; and leaf number were clustered together. As the seedling grows, the association of stem area and compactness was clustered with traits that had a minimum role in estimating WPAi and WPAs. In agreement with the preceding section, genotypes LB-46, GM-217, and Varshadhan exhibited the highest magnitude toward geometric traits on both dates of observation.
Rationale of early seedling vigor trait in rice
As DSR saves water and labor [1, 22], it is becoming popular across all fronts of rice ecosystems by covering 29 million ha of the rice area in Asia (21% of the total rice area) . Good crop establishment is a major challenge in the DSR system. Germinating seeds and/or seedlings under DSR are exposed directly to an array of stresses [1, 7]. Thus, strong and early seedling vigor is an imperative breeding trait for DSR varieties. Rice cultivars with high early seedling vigor decrease crop yield loss due to weeds (16% yield loss), water , and nutrient by exceeding the growth of weeds . In addition, they help to achieve rapid and uniform emergence from the field. Genetic improvement for early vigor has been practiced in rice for better crop establishment. Caton et al.  reported that early vigor was a highly repeatable trait among rice cultivars. Further, 87% grain yield variation in rice was accounted for by vegetative vigor (2 weeks after sowing) in a comparison between weed and weed-free field conditions . Therefore, the identification and use of suitable donors with relevant traits associated with early seedling vigor and variability available in genotypes are essential. This helps in selecting traits and designing plants for early seedling vigor. However, the introgression of seedling vigor into any recipient parent or selection of lines with seedling vigor in segregating generations could not be achieved with destructive sampling. This limits breeding for seedling vigor traits. The absence of a non-destructive method to estimate seedling vigor hampers the whole experiment and creates hurdles in exploiting early seedling vigor in rice. Thus, a robust automated non-destructive image phenotyping technique will help to overcome these limitations in the area of early seedling vigor. Seedling vigor has several component traits: biomass accumulation, canopy coverage, plant height, etc. These components are traditionally recorded visually and in a destructive way [28, 29]. To establish a relation of those traits with seedling vigor and to judge the robustness of WPAi, an automated non-destructive phenotyping technique in rice was developed and the method established with genetically diverse rice genotypes. Phenotyping using RGB imaging has been designed for field crops, often for abiotic stresses (drought, salinity, and cold) [15, 30, 31], but, surprisingly, rare attempts were made to screen traits such as early seedling vigor.
The necessity of automated image-based phenotyping for seedling vigor
Early seedling vigor is a polygenic trait, and it requires measurements of phenotypic data of component traits for genetic dissection into smaller manageable and measurable components . Conventionally, early seedling vigor assessment involves manual visual scoring, leaf area measurement, shoot biomass measurement, etc. [28, 33]. Manual methods are labor-intensive, in particular, the measurements are prone to human error, manual data management, and data keeping, and may not be suitable for handling a large number of samples. In rice, 2-week-old seedlings are small and delicate, and often lead to error. Therefore, robust and automated phenotyping platforms that can capture high-quality and reliable phenotypic data would be error-free and straightforward to handle. Image-based phenotyping offers several advantages over destructive methods, in which digital color images are used to quantify phenotype [13, 14, 34]. In our experiment, we established a phenotypic platform with available resources. It saves nearly 80% of the time (660 s were required per sample of five plants with four persons for observing data by destructive sampling method vis-à-vis two persons with 152 s of proposed imaging protocol) and ~ 50% on the cost of labor. A population developed for early seedling vigor would segregate for the early seedling vigor component traits and differ at the genic level. Therefore, destructive sampling for early seedling vigor populations might lead to a loss in variation and deviation from a normal probability distribution and Hardy–Weinberg law with biased results. Therefore, image-based phenotyping would overcome those constraints to achieve precise phenotyping with better reproducibility.
Destructive versus non-destructive method of phenotyping
Data observed by RGB imaging have identified subtle differences between genotypes at both dates (14 and 28 DAS) of observation, while the same could not be achieved with the traditional way of measurement observed with respect to 14 DAS (Tables 1, 2, 3, 4). This suggests that phenotyping by imaging would be a better technique to find differences in the early stage of growth, which is dynamic and delicate for manual handling. Using destructive sampling, based on a higher growth rate, genotypes with seedling vigor were identified and were compared and selected through WPAi. This comparison was made to understand the potential and accuracy of the measurements obtained from RGB imaging. Among the destructive methods, CGR and CGR-WPAi for shoots (R2 > 74%) and roots (R2 > 45%) were strongly related and RGR of shoot and root dry weight with RGR-WPAi had weak similarity in the selection of genotypes, whereas AGR and AGR-WPAi had no similarity. Thus, the high magnitude of association for seedling vigor estimated from the biomass-based destructive method by CGR can be replaced with the CGR-WPAi technique, which was our primary experimental objective. To find out the fidelity of the WPAi method, the traditional way of sampling was adopted to estimate whole-plant area (WPAs) by scanning individual leaf blades and stems. In the regression analysis, WPAs and WPAi were strongly related (R2 > 83%) on both (14 and 28 DAS) observation dates. The relationship between them was very high and WPAi can also be used in place of WPAs (destructive) to estimate seedling vigor. Similarly, Hairmansis et al.  and Campbell et al.  have identified tolerant rice genotypes under salinity by capturing morphological and physiological responses by processing RGB images in a non-destructive way. They have proved the successful introgression of imaging techniques in high-throughput phenotyping. Further, they have explained the reliability of imaging techniques by the linear relationship between fresh weight and projected image area. However, the dry weight of samples would always be preferred over the fresh weight of samples to avoid variability in moisture content among the samples and genotypes. Therefore, to comprehend the fidelity of imaging techniques, the relationship between WPAi and shoot dry biomass was studied. A strong correlation existed between WPAi and shoot dry weight at 28 DAS (r = 92, p < 0.01; R2 = 86%) and medium association at 14 DAS (r = 0.31; R2 = 9.6%). It is optimistic that accuracy might increase with the age of the seedlings with more biomass. Therefore, shoot area by images would be a good surrogate for estimating seedling vigor using shoot biomass up to 4 weeks of age (Fig. 6a, b) and, further, to determine the relationship between WPAi and root dry weight. WPAi of shoot and root dry weight was correlated and we could find a strong association between them at 28 DAS (r = 86, p < 0.01; R2 = 74%) and 14 DAS (r = 0.58; R2 = 34%). This suggests that WPA obtained by RGB imaging would also be useful for understanding the growth rate of below-ground parts of rice seedlings in relative time .
Grouping of traits across DAS to determine genotypes with higher variability
PCA was carried out for all seven genotypes to identify trends among the genotypes and the traits responsible for the source of variation for seedling vigor. Nine highly variable traits on day 14 and 10 traits on day 28 govern the highest variation among the seven genotypes. Across the two observation dates, the traits WPAi, WPAs, maximum leaf growth at 28 DAS (6th and 7th leaf area), convex hull, and top view area of both observation dates delivered the highest variation among 73 seedling traits (29 traits at 14 DAS and 44 traits at 28 DAS) (Fig. 7a, b). Thus, for the image-based measured parameters of whole-plant area, convex hull captured the maximum variability (R2 = 0.25 (14 DAS), R2 = 0.79 (28 DAS)), which has maximum variation toward WPAi and is considered as an important trait for the selection of genotypes for seedling vigor. Thus, these traits were highly variable and contributed to the selection of genotypes for early seedling vigor and are thought to be useful in improving seedling vigor in rice through automated image phenotyping. Among the traits studied at two different times, the traits measured at 28 DAS contributed much to differentiate genotypes. Therefore, for the study involved in the identification of vigorous genotypes at the seedling stage, the parameters observed at 28 DAS by image-based phenotyping would be adequate. This would help in circumventing destructive sampling, thus saving resources, time, and labor involved in observing data at two different times. The present methodology in combination with the automated handling system would efficiently screen the dynamic responses of breeding lines in limited time. Further, this will help us to understand the mechanisms involved in enhancing the growth rate and genes to design a breeding program.
Magnitude of component traits in determining vigorous genotypes
The digital color image data were used for estimating WPAi and had a relationship with the component traits of seedling vigor, such as third leaf length (77%), third leaf area (51.29%), and stem area (61.33%) at 14 DAS, which explained the highest variation for WPAi. Whereas, at 28 DAS, fifth and sixth leaf area (77.79% and 53.25%, respectively), stem weight (82.30%), and tiller number (79.80%) had a higher magnitude of regression for WPAi. Leaf area at both 14 and 28 DAS seems to be a major determinant to predict genotype performance; in particular, the recent fully expanded leaves were found to have a higher contribution. This is in line with Hairmansis et al.  and Nguyen et al. , who reported that leaf traits would be the best predictor in identifying the performance of field pea and rice, respectively. This confirms that the physiologically active and fully expanded leaves could be used as representations to predict early vigor in rice. Shoot length at 14 and 28 DAS had a different level of contribution toward WPAi. Variation in shoot length at 14 DAS (0.95, p < 0.01) had a greater contribution toward WPAi than shoot length at 28 DAS (− 0.19 ns) and vice versa in the case of geometric traits observed by image analysis (caliper length, eccentricity, convex hull, and top view area). This might be due to the variation for shoot length expressed by genotypes at 14 DAS that was significantly different at 1% (Tables 1, 2), while at 28 DAS it had significance at 5%.
Understanding the association between traits would help to identify target genotypes with the surrogate traits. Length, width, and area of the leaf at each level made a different contribution toward WPAi in determining vigorous genotypes. Third leaf length at 14 DAS had a positive relationship with WPAi and explained 77.64% of WPAi variation in comparison with early formed leaves (first and second leaf). The weak association of first and second leaf length with WPAi might be due to their tiny nature and their decreased visibility in the image. On the other hand, the increase in leaf area of early formed leaves (first and second) had a positive correlation with root length. In the present experiment, leaf number had a strong positive association with biomass. The rise in the number of leaves in the early stage of seedling growth would be due to an increase in tiller number; ultimately, that would increase biomass. Further, an increase in leaf number would have a cascade positive effect on top view area, caliper length, and convex hull. The enhanced convex hull area increases ground cover by the canopy, which will have ample access to water and fertilizer. The existence of temporal and genotypic differences in canopy cover is considered an important trait for distinguishing genotypes with early vigor . Good ground cover increases the weed smothering effect, provides ample access to fertilizer and water , and decreases loss of water by covering the ground. However, the increase in leaf number had a negative association with compactness and leaf width. This is in line with Richards : the high rate of canopy cover was associated with an increased rate of tillers rather than leaf expansion. According to the ideal plant type concept, the plant should be compact in nature with few productive tillers and broad leaves. The ideal plant type concept would be more suitable for an ecosystem in which there is no weed competitiveness. The same concept would not be appropriate for an ecosystem in which weed competitiveness is a regular phenomenon, as in DSR. Therefore, a genotype with early seedling vigor accumulating high biomass and having enhanced convex hull and evenly spaced leaves with minimum compactness would be more suitable. A compact plant type would always have overlapped leaves with a 45-degree angle to avoid shade effects for more light interference. Genotypes with a better convex hull and evenly spaced narrow leaves with decreased leaf width and decreased shade effects for the contemporary leaves for proper light interference to improve photosynthesis would be preferred. Therefore, these geometric traits might play a big role in studying the architecture of the plant.
At 14 DAS, GM-217 and Vandana attained maximum growth with more leaves and increased shoot length and biomass. The high biomass and leaf number had increased the WPAi, top view area, eccentricity, and convex hull. This had allowed the genotypes to secure the top position while their growth rate from 14 to 28 DAS slowed. Genotypes LB-46 and Varshadhan had gained a pronounced growth rate by increasing tiller number. The increase in tiller number eventually increases leaf number, biomass, WPAi, top view area, convex hull, caliper length, and leaf area of terminal leaves. This helps the genotypes to be more vigorous at 28 DAS. The increase in top view area and convex hull enhances ground cover. Thus, the enhanced soil cover improves the weed smothering effect. Further, the increase in the size of the canopy (caliper length) and leaf area augments ground cover. These types of geometric trait data are difficult to generate through manual systems and are time-consuming .
A comparative study was done between the destructive and non-destructive methods to identify genotypes having high seedling vigor. Overall, LB-46 was found to be a common genotype across the methods used to judge the genotypes, with the highest seedling vigor. Varshadhan and AC38399 were identified as the next best genotypes by WPAi, while the traditional method identified GM-217 and Varshadhan as the next best performers. The variability in the position of genotypes between WPAi and the traditional method might be due to the handling of different plant samples of the same variety at a relative time. In addition, their growth rate and manual handling of 2-week-old small seedlings would cause some errors. Therefore, the image analysis technique was found to be a very effective determinant of seedling vigor without human interference. These variations were easily captured through WPAi and chances of human error could be decreased in such cases. Further, the traditional way of estimating vigor in the case of often and highly cross-pollinated crops would not be more reliable. Therefore, WPAi has quite a few advantages and can be used at any stage of seedlings across different crops. On the other hand, we found some minor differences in the ranking of genotypes at 14 and 28 DAS between WPAi and WPAs. These differences in the ranking of genotypes between imaging and scanning might be due to the overlapping of leaves during imaging and some unexposed area of droopy leaves.
Seven rice (Oryza sativa L.) genotypes of improved and traditional lines, LB-46, GM-217, AC38399, ARC10656, Vandana, Sabita (NC492), and Varshadhan, were used in this study. Varshadhan and Sabita were developed at the International Rice Research Institute (IRRI), Philippines, and in Chinsurah, West Bengal, India, respectively, for the semi-deep ecosystem. Vandana was developed for upland conditions by ICAR-National Rice Research Institute (NRRI). ARC10656 and GM-217 belong to traditional rice collections of Assam and Tamil Nadu, respectively. LB-46 was the progeny of Oryza sativa x O. nivara, developed by NRRI, Cuttack.
Seed selection and sterilization
Seeds of all the genotypes were grown and harvested in the wet season of 2017 and packed separately for drying. All the genotypes were sorted by uniform seed size and underwent heat treatment to break seed dormancy. The seeds were kept in a hot-air oven at 50 °C for 45 h. Later, the seeds were surface-sterilized with 75% ethanol for 1 min. These seeds were further sterilized with 2.5% sodium hypochlorite for 20 min and washed five times with sterile distilled water to remove any traces of sterilizing agent.
Growing conditions and experimental design
Five seeds of each genotype were sown into a pot (white color, 20 cm height × 15 cm diameter) containing 2.5 kg of clayey loam without any external fertilizer. Each genotype was raised in five pots with five biological replicates. The spacing between plants was maintained at 20 × 15 cm to simulate the recommended spacing for direct-seeded conditions. Therefore, pots were arranged in such a way to have 20 cm between rows and 15 cm within rows. After one week, the seedlings were thinned to a single seedling per pot by maintaining sufficient moisture with a proper drainage hole at the base. The experiment was conducted in the net house at NRRI (20°27′09″ N, 85°55′57″ E, 26 masl), Cuttack, during March 2018. The plants received 13 h/11 h of day/night cycle. The average temperature in the net house was 33 °C in the day and 23 °C at night, with an average light intensity of ~ 1200 µmol m−2 s−1 during the observation period.
Image capture and analysis
At 14 and 28 DAS, images of five biological replicates per genotype were captured, and the same set of plants was used to measure the morphometric data to estimate growth parameters on the same day. Three different techniques were followed to measure the morphometric data of the whole plant: (1) non-destructive imaging and (2) two versions of destructive growth analysis.
Method-1: Non-destructive imaging
In the first method, data were observed non-destructively on a whole-plant basis using a 12-megapixel Nikon camera (RGB images) at a distance of 1.5 m. To have a uniform background, the potted plant was placed over a raised platform having a dark background behind and over the platform. High-intensity artificial light was used to obtain a uniform background effect and strong wind or airflow was checked to avoid movement of leaves during imaging. A known scale indicator/ruler was placed above and adjacent to the pot to calculate the leaf/whole-plant area with proper labeling. Three colored images per plant were captured from the top of the plant and two from either side of the plant at 90° (Fig. 8). A uniform distance between plant and camera, camera setting, and background light were maintained throughout the imaging process.
The recorded images were processed using Image J software. In brief, the respective image was cropped to remove any excess area on all four sides by retaining the potted plants. To assess the green portion of the plant, the cropped image was used to separate the plants from the imaging background by selecting the color threshold. To identify the plant as a whole and to remove any further noise, all portions of the plant were highlighted using hue, saturation, and brightness (Fig. 9). Using the known indicator, the whole-plant area/desired portion of the plant was converted from pixels to square millimeters. The summed area of all three images (top and two side views) was used to estimate the whole-plant area (WPAi) and expressed in square millimeters. Additionally, geometric traits such as convex hull, compactness, caliper length, and eccentricity were determined from the top view image (captured from the top of the plant) as described by Neilson et al. . The observations collected from imaging were used to calculate relative (RGRi), absolute (AGRi), and crop (CGRi) growth rate of the plants. The growth rate was calculated using the average of final and initial WPAi, simulating the way the growth rate was calculated for the destructive method.
where A1 and A2 are the whole-plant area at times t1 and t2, respectively, i = image-based, loge = natural logarithm, and P = spacing (m2).
Method-2 and -3: Destructive growth analysis
The plants used for imaging were uprooted and the roots were washed to record morphometric traits. This was followed by the same plants being used for the destructive method. The leaves were meticulously cut from the stem and placed between the paper pages of a clean notebook to maintain their shape. Later, the leaves and stem were scanned using a flatbed scanner with a ruler to calculate the leaf/whole-plant area using Image J software. The scanned images were used to estimate the area of each leaf and stem following the steps adapted to convert pixels to square millimeters as described earlier in the section on image capture and analysis. Traits such as tiller number, number of leaves, stem thickness (mm), shoot length (mm), root length (mm), shoot dry weight (g), root dry weight (g), specific leaf weight (g), and stem weight (g) were manually observed, while specific leaf length (mm), width (mm), and area (mm2) and stem area (mm2) were measured from scanned images of the flatbed scanner by destructive sampling. The summed area of all leaves and stems observed from the scanned images was used to estimate the whole-plant area (WPAs). In addition, the area of a single leaf blade calculated from the scanned image was compared with the biomass of the respective leaf and WPAi to understand the relationship and percentage of variation contributed by them toward WPAi.
All three methods were compared to assess the fidelity of the data obtained through the imaging process. To assess the growth rate of seedlings, absolute growth rate, relative growth rate, and crop growth rate were calculated accordingly:
where h1 and h2 are plant height at times t1 and t2, respectively, and m = manual method.
RGR was determined by using the dry weight of periodical observations and represented as mg g−1 day−1.
where W1 and W2 are plant dry weights at times t1 and t2, respectively.
CGR was calculated by measuring plant dry weight at a regular interval of time divided by land area and represented as g m−2 day−1:
where W1 and W2 are plant dry weights at times t1 and t2, respectively, and P = spacing (m2).
Linear regression was estimated between WPAi and seedling traits using MS Office Excel 2016. Principal component analysis (PCA) was performed with 29 traits at 14 DAS and 44 traits at 28 DAS to estimate the variability among genotypes and traits. Biplot figures explain the variances of the variables and correlation between the variables through vectors and similarity between genotypes in the multivariate space based on the nature of growth rate [41, 42]. These analyses were performed using Windostat 7.5 software.
In the present experiment, the non-destructive-based imaging technique captured the dynamic responses of plants in the early stage and revealed significant differences across genotypes. Consistency in the ranking of genotypes across different methods and a strong correlation between morphological and image-associated traits confirm the reliability and reproducibility of the proposed method. The proposed imaging technique and the identified geometric traits convex hull and top view area were found to be significant in determining promising genotypes for early seedling vigor during the initial phase of plant establishment. Further, the method saves resources, time, and labor by validating that the parameters observed at 28 DAS are adequate in identifying vigorous genotypes. This has raised confidence that imaging techniques have the potential to identify and differentiate small differences that are considered as phenotypically identical and difficult to distinguish by unidentifiable traits through conventional measurements. To capitalize on the present non-destructive imaging technique as a high-throughput to handle more samples in a given time, the method needs to be fully automated. This would enable integrating the platform as a tool with the forward genetics approach in the identification of QTLs/genes for the traits identified under this system for possible future improvements in the study area of direct-seeded rice.
Availability of data and materials
The datasets during and/or analyzed during the current study are available from the corresponding author upon request.
absolute growth rate
crop growth rate
days after sowing
principal component analysis
relative growth rate
- WPAi :
Mahender A, Anandan A, Pradhan SKK. Early seedling vigour, an imperative trait for direct-seeded rice: an overview on physio-morphological parameters and molecular markers. Planta. 2015;241:1027–50.
Humplík JF, Lazár D, Husičková A, Spíchal L. Automated phenotyping of plant shoots using imaging methods for analysis of plant stress responses: a review. Plant Methods. 2015;11:29.
Crowell S, Falcão AX, Shah A, Wilson Z, Greenberg AJ, McCouch SR. High-resolution inflorescence phenotyping using a novel image-analysis pipeline, PANorama. Plant Physiol. 2014;165:479–95.
Yang W, Feng H, Zhang X, Zhang J, Doonan JH, Batchelor WD, et al. Crop phenomics and high-throughput phenotyping: past decades, current challenges and future perspectives. Mol Plant. 2020. https://doi.org/10.1016/j.molp.2020.01.008.
Sun D, Cen H, Weng H, Wan L, Abdalla A, El-Manawy AI, et al. Using hyperspectral analysis as a potential high throughput phenotyping tool in GWAS for protein content of rice quality. Plant Methods. 2019;15:54.
Rahaman M, Chen D, Gillani Z, Klukas C, Chen M. Advanced phenotyping and phenotype data analysis for the study of plant growth and development. Front Plant Sci. 2015;6:619.
Anandan A, Anumalla M, Pradhan SK, Ali J. Population structure, diversity and trait association analysis in rice (Oryza sativa L) germplasm for early seedling vigor (ESV) using trait linked SSR markers. PLoS ONE. 2016;11:e0152406.
Anandan A, Meher J, Sah RP, Samantaray S, Parameswaran C, Panneerselvam P, et al. Enhancing Input Use Efficiency in Direct-seeded Rice with Classical and Molecular Breeding. In: Rice reseach for enhancing productivity, profitability and climate resilience, 2018, p. 73–89.
Neilson EH, Edwards AM, Blomstedt CK, Berger B, Møller BL, Gleadow RM. Utilization of a high-throughput shoot imaging system to examine the dynamic phenotypic responses of a C4 cereal crop plant to nitrogen and water deficiency over time. J Exp Bot. 2015;66:1817–32.
Winterhalter L, Mistele B, Jampatong S, Schmidhalter U. High throughput phenotyping of canopy water mass and canopy temperature in well-watered and drought stressed tropical maize hybrids in the vegetative stage. Eur J Agron. 2011;35:22–32.
Romano G, Zia S, Spreer W, Sanchez C, Cairns J, Araus JL, et al. Use of thermography for high throughput phenotyping of tropical maize adaptation in water stress. Comput Electron Agric. 2011;79:67–74.
Harbinson J, Prinzenberg AE, Kruijer W, Aarts MGM. High throughput screening with chlorophyll fluorescence imaging and its use in crop improvement. Curr Opin Biotechnol. 2012;23:221–6.
Rajendran K, Tester M, Roy SJ. Quantifying the three main components of salinity tolerance in cereals. Plant Cell Environ. 2009;32:237–49.
Golzarian MR, Frick RA, Rajendran K, Berger B, Roy S, Tester M, et al. Accurate inference of shoot biomass from high-throughput images of cereal plants. Plant Methods. 2011;7:2.
Hairmansis A, Berger B, Tester M, Roy SJ. Image-based phenotyping for non-destructive screening of different salinity tolerance traits in rice. Rice. 2014;7:16.
Joshi DC, Singh V, Hunt C, Mace E, van Oosterom E, Sulman R, et al. Development of a phenotyping platform for high throughput screening of nodal root angle in sorghum. Plant Methods. 2017;13:56.
Wu D, Guo Z, Ye J, Feng H, Liu J, Chen G, et al. Combining high-throughput micro-CT-RGB phenotyping and genome-wide association study to dissect the genetic architecture of tiller growth in rice. J Exp Bot. 2019;70:545–61.
Atkinson JA, Wingen LU, Griffiths M, Pound MP, Gaju O, Foulkes MJ, et al. Phenotyping pipeline reveals major seedling root growth QTL in hexaploid wheat. J Exp Bot. 2015;66:2283–92.
Narisetti N, Henke M, Seiler C, Shi R, Junker A, Altmann T, et al. Semi-automated root image analysis (saRIA). Sci Rep. 2019;9:1–10.
Turner SD, Ellison SL, Senalik DA, Simon PW, Spalding EP, Miller ND. An automated image analysis pipeline enables genetic studies of shoot and root morphology in carrot (Daucus carota L.). Front Plant Sci. 2018;9:1703.
Richards RA, Rebetzke GJ, Watt M, Condon AGT, Spielmeyer W, Dolferus R. Breeding for improved water productivity in temperate cereals: phenotyping, quantitative trait loci, markers and the selection environment. Funct Plant Biol. 2010;37:85–97.
Richards M, Sander BO. Alternate wetting and drying in irrigated rice. Climate-Smart Agriculture Practice Brief. Copenhagen, Denmark: CGIAR Research Program on Climate Change, Agriculture and Food Security (CCAFS), Copenhagen, Denmark. http://www.ccafs.cgiar.org. 2014.
Pandey S, Velasco L. Economics of direct seeding in Asia: patterns of adoption and research priorities. In: Direct seeding: Research strategies, opportunities. International Rice Research Institute, Los Baños, Philippines; 2002. p. 3–14.
Rao AN, Johnson DE, Sivaprasad B, Ladha JK, Mortimer AM. Weed management in direct-seeded rice. Adv Agron. 2007;93:153–255.
Ali J, Jewel ZA, Mahender A, Anandan A, Hernandez J, Li Z. Molecular genetics and breeding for nutrient use efficiency in rice. Int J Mol Sci. 2018;1:27.
Caton BP, Cope AE, Mortimer M. Growth traits of diverse rice cultivars under severe competition: implications for screening for competitiveness. Field Crops Res. 2003;83:157–72.
Zhao DL, Atlin GN, Bastiaans L, Spiertz JHJ. Developing selection protocols for weed competitiveness in aerobic rice. Field Crops Res. 2006;97:272–85.
McDonald GK. Competitiveness against grass weeds in field pea genotypes. Weed Res. 2003;43:48–58.
Armstrong EL, Pate JS. The field pea crop in SW Australia I Patterns of growth, biomass production and photosynthetic performance in genotypes of contrasting morphology. Aust J Agric Res. 1994;45:1347–62.
Honsdorf N, March TJ, Berger B, Tester M, Pillen K. High-throughput phenotyping to detect drought tolerance QTL in wild barley introgression lines. PLoS ONE. 2014;9:e97047.
Humplík JF, Lazár D, Fürst T, Husičková A, Hýbl M, Spíchal L. Automated integrative high-throughput phenotyping of plant shoots: a case study of the cold-tolerance of pea (Pisum sativum L). Plant Methods. 2015;11:20.
Rebolledo MC, Dingkuhn M, Courtois B, Gibon Y, Clément-Vidal A, Cruz DF, et al. Phenotypic and genetic dissection of component traits for early vigour in rice using plant growth modelling, sugar content analyses and association mapping. J Exp Bot. 2015;66:5555–66.
Moot DJ. Harvest index variability within and between field pea (Pisum sativum L.) crops. Lincoln University; PhD Thesis. Lincoln University, Canterbury, New Zealand, 1993.
Harris BN, Sadras VO, Tester M. A water-centred framework to assess the effects of salinity on the growth and yield of wheat and barley. Plant Soil. 2010;336:377–89.
Campbell MT, Knecht AC, Berger B, Brien CJ, Wang D, Walia H. Integrating image-based phenomics and association analysis to dissect the genetic architecture of temporal salinity responses in rice. Plant Physiol. 2015;168:1476–89.
Kumar R, Kumawat N, Kumar S, Kumar R, Kumar M, Sah RP, et al. Direct seeded rice: research strategies and opportunities for water and weed management. ORYZA-An Int J Rice. 2016;53:354–65.
Nguyen GN, Norton SL, Rosewarne GM, James LE, Slater AT. Automated phenotyping for early vigour of field pea seedlings in controlled environment by colour imaging technology. PLoS ONE. 2018;13:e0207788.
Walter A, Liebisch F, Hund A. Plant phenotyping: from bean weighing to image analysis. Plant Methods. 2015;11:14.
Poorter H, Bühler J, van Dusschoten D, Climent J, Postma JA. Pot size matters: a meta-analysis of the effects of rooting volume on plant growth. Funct Plant Biol. 2012;39:839–50.
Richards RA. Crop improvement for temperate Australia: future opportunities. Field Crops Res. 1991;26:141–69.
Anandan A, Rajiv G, Ramarao A, Prakash M. Internode elongation pattern and differential response of rice genotypes to varying levels of flood water. Funct Plant Biol. 2012;39:137–45.
Anandan A, Pradhan SK, Das SK, Behera L, Sangeetha G. Differential responses of rice genotypes and physiological mechanism under prolonged deepwater flooding. Field Crops Res. 2015;172:153–63.
This work was supported by the ICAR-National Rice Research Institute (NRRI). We gratefully acknowledge the support extended by the director, NRRI, to conduct this experiment and special thanks go to the administrative staff of our institute who supported the conduct of this experiment in official terms.
The authors would like to thank and acknowledge the Bill & Melinda Gates Foundation (BMGF) for providing a research sub-grant to Jauhar Ali for the Green Super Rice Project under ID OPP1130530. We would also like to thank the ICAR-National Rice Research Institute, Cuttack, Odisha, for providing funds to Annamalai Anandan for the institutional project. Funders did not have any role in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Anandan, A., Mahender, A., Sah, R.P. et al. Non-destructive phenotyping for early seedling vigor in direct-seeded rice. Plant Methods 16, 127 (2020). https://doi.org/10.1186/s13007-020-00666-6