Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Automated phenotyping for early vigour of field pea seedlings in controlled environment by colour imaging technology

  • Giao N. Nguyen ,

    Roles Conceptualization, Data curation, Formal analysis, Investigation, Methodology, Writing – original draft, Writing – review & editing

    giao.nguyen@ecodev.vic.gov.au

    Affiliation Australian Grains Genebank, Agriculture Victoria, Grains Innovation Park, Horsham, Victoria, Australia

  • Sally L. Norton,

    Roles Conceptualization, Data curation, Formal analysis, Funding acquisition, Methodology, Project administration, Supervision, Writing – original draft, Writing – review & editing

    Affiliation Australian Grains Genebank, Agriculture Victoria, Grains Innovation Park, Horsham, Victoria, Australia

  • Garry M. Rosewarne,

    Roles Conceptualization, Data curation, Funding acquisition, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliation Agriculture Victoria, Grains Innovation Park, Horsham, Victoria, Australia

  • Laura E. James,

    Roles Data curation, Investigation, Writing – review & editing

    Affiliation Agriculture Victoria, Grains Innovation Park, Horsham, Victoria, Australia

  • Anthony T. Slater

    Roles Conceptualization, Supervision, Writing – review & editing

    Affiliation Agriculture Victoria, AgriBio, Bundoora, Victoria, Australia

Correction

6 Dec 2018: The PLOS ONE Staff (2018) Correction: Automated phenotyping for early vigour of field pea seedlings in controlled environment by colour imaging technology. PLOS ONE 13(12): e0208972. https://doi.org/10.1371/journal.pone.0208972 View correction

Abstract

Early vigour of seedlings is a beneficial trait of field pea (Pisum sativum L.) that contributes to weed control, water use efficiency and is likely to contribute to yield under certain environments. Although breeding is considered the most effective approach to improve early vigour of field pea, the absence of a robust and high-throughput phenotyping tool to dissect this complex trait is currently a major obstacle of genetic improvement programs to address this issue. To develop this tool, separate trials on 44 genetically diverse field pea genotypes were conducted in the automated plant phenotyping platform of Plant Phenomics Victoria, Horsham and in the field, respectively. High correlation between estimated plant parameters derived from the automated phenotyping platform and important early vigour traits such as shoot biomass, leaf area and plant height indicated that the derived plant parameters can be used to predict vigour traits in field pea seedlings. Plant growth analysis demonstrated that the “broken-stick” model fitted well with the growth pattern of all field pea genotypes and can be used to determine the linear growth phase. Further analysis suggested that the estimated plant parameters collected at the linear growth phase can effectively differentiate early vigour across field pea genotypes. High correlation between normalised difference vegetation indices captured from the field trial and estimated shoot biomass and top-view area confirmed the consistent performance of early vigour field pea genotypes under controlled and field environments. Overall, our results demonstrated that this robust screening tool is highly applicable and will enable breeding programs to rapidly identify early vigour traits and utilise germplasm to contribute to the genetic improvement of field peas.

Introduction

Field pea (Pisum sativum L.) is a legume crop that is widely grown around the world with annual production of c. 11 million metric tonnes produced from 6.9 million hectares of cultivated land [1]. Australia is among the 10 largest field pea producing countries, where the crop accounts for 20% of pulse production in South Australia and Victoria, and is the second largest pulse crop grown in Western Australia and third in New South Wales [2]. The Australian field pea industry produces approximately 400,000 metric tonnes of grain annually. Of this, about 191,000 metric tonnes are exported with the market value of $A80 million [3]. Field pea production brings profit for growers as a cash crop and provides other benefits for the farming systems such as non-cereal crop rotation and biological nitrogen fixation. Field pea and other crop legumes annually contribute approximately 5–7 million metric tonnes of biologically fixed nitrogen to cultivated soil, saving farmers $US 8–12 billion on nitrogen fertilizer costs globally [4, 5].

Like many other agricultural crops, field pea production is critically affected by biotic and abiotic stresses such as weeds, drought and heat [6]. Competition from weeds is one of the major biotic constraints affecting field pea production, which can result in yield loss up to 25% [7]. Field pea is a very poor competitor against weeds compared to other crops due to its weak early vigour at the seedling stage [8, 9]. Globally, herbicides are widely used for weed control during pea cultivation, and although instantly effective, overuse of herbicides with similar active ingredients and modes of action is resulting in herbicide resistant weed biotypes [10] and increased production costs [11]. Herbicides can also potentially affect rhizobium and symbiotic nitrogen fixation with field pea, causing smaller positive impacts on subsequent crop rotations [12]. Likewise, drought and heat stresses cause particularly greater yield losses in field pea crops [13, 14]. These abiotic stresses can have critical effects if they occur during flowering and grain filling by affecting reproductive organs and pod setting, thus reducing seed number [15, 16].

Previous studies suggest that breeding for tolerant varieties is one of the most effective strategies to cope with biotic and abiotic stresses, and early vigour traits have been considered an important selection criterion by field pea breeders [6, 1719]. Although early vigour can be improved by using higher sowing rates and applying more nitrogen fertilizer, studies suggest that enhancing early vigour by genetic improvement is more effective and reliable [20]. Early vigour is the plant’s ability to establish quickly after sowing at the seedling stage and has been studied extensively in rice, wheat and other cereals [2123]. Genetic studies in wheat showed that seedlings of vigorous genotypes can produce biomass rapidly, tiller earlier, have more leaves and have greater water and nitrogen use efficiency [2426].

Although weed control can be managed by herbicides and other agronomical practices such as planting density, row spacing and orientation, use of vigorous genotypes with greater competitiveness is the most effective, non-chemical and environmentally friendly strategy [18, 27, 28]. Early vigour is also an important breeding trait for higher water use efficiency, especially in Mediterranean environments, since it minimizes soil water evaporation by boosting early vegetative ground cover [17, 29]. In water-limited environments, wheat genotypes with early vigour decreased water evaporation from the soil surface by reducing water loss by 90–110 mm and increased transpiration efficiency by 10% [30]. As a result, genotypes with early vigour have greater carbohydrate reserves before anthesis that can compensate for a photosynthesis reduction of up to 36% if drought occurs during the grain fill stage [31]. Early vigour is also an ideal trait of tropical crops grown in cold environments because it confers chilling tolerance [32]. Seedlings with high early vigour have a higher nitrogen uptake and photosynthetic nitrogen use efficiency [26, 33, 34]. Moreover, early vigour field pea varieties showed broader adaptation and yield maintenance under unfavorable growth conditions [3537]. Thus, there is a pressing need to develop early vigour field pea genotypes via breeding in response to biotic and abiotic stresses.

Early vigour is a polygenic trait that requires a large volume of high quality phenotypic data to dissect its genetic composition into smaller manageable and measurable components [38]. Conventionally, early vigour trait assessment involves manual methods such as visual scoring, and measuring leaf area, plant height, and shoot biomass [18, 39]. Although attainable, these methods are labor intensive, subjective and prone to human errors, and are not suitable for large scale trials. Therefore, robust and high-throughput phenotyping tools and platforms that can generate reliable and high quality phenotypic data for genomic selection have become the rate-limiting step in field pea’s genetic improvement [40].

Non-destructive phenotyping technology using sensors and cameras can offer high-throughput and reproducible screens of large scale trials as well as reliable, high quality data and dynamic growth analysis of crops [41]. This technology has also been recommended for studying early vigour for nitrogen use efficiency in agricultural crops [22]. The technology was built to detect and quantify the spectral reflectances resulting from the interaction between plant parts and electromagnetic radiation at different spectral regions such as visible (VIS, 400–700 nm), near infrared (NIR, 700–1000 nm) and short-wave infrared (SWIR, 1000–2500 nm) [41, 42]. Software and computer vision enable the analysis of these reflectances to derive digital plant objects that can be used as surrogates for plant architectural morphology, biomass, and grain yield [43]. Several automated plant phenotyping platforms that comprise of growth facilities, sensors and cameras are commercially available and have been successfully applied in crop research under controlled and field conditions [22, 4447]. For example, an automated high throughput phenotyping platform, PlantScreen (Photon Systems Instruments, Brno, Czech Republic) was used to assess the cold tolerance of field pea using digital colour imaging technology under controlled environments [48]. Similarly, Roth et al. [49] applied an aerial-based imaging phenotyping platform to estimate field pea biomass under field conditions. Vegetation indices such as normalized difference vegetation index (NDVI) derived from optical sensors has been used to analyse the growth of field pea and other crops under various field conditions [50, 51].

Here we report on the development of a high-throughput phenotyping method to evaluate early vigour of field pea in a controlled environment using an automated colour imaging technology and a comparative performance of the same genotypes under field conditions. The role of early seedling vigour in field pea breeding and production is also discussed.

Materials and methods

Plant material and experimental design

Forty-four genetically diverse field pea genotypes were used in these experiments to investigate early seedling vigour traits (S1 Table) [52]. Field pea seeds were carefully selected to ensure that seeds of the same genotype had similar size and shape to guarantee a similar level of germination. In the first experiment, field pea plants were grown in the greenhouse of Plant Phenomics Victoria, Horsham. Euro white pots (200 mm diameter x 190 mm deep, Garden City Plastics, Victoria, Australia) were filled by weight with 3.5 liters of potting mix consisting of 1,000 L legume mix (Biogro, SA), 1 kg Floranid 32, 1 kg Blue Macracote Coloniser Plus, 1 kg Nutricote N16, 1 kg Microplus trace element fertilizers, 225 g LibFer SP, 2 kg SaturAid, and 25 kg Lime. The pots were watered prior to sowing and placed on white saucers throughout the experiment to avoid water leaking on to the system. Three seeds were sown per pot and these were kept on rolling benches in the greenhouse of Plant Phenomics Victoria, Horsham. Each pot was thinned to one plant after seeds had germinated, approximately seven days after sowing (DAS), and blue wire cages were inserted into the pots to support plant growth. The colored cages facilitated differentiation of plant material from the support structure for imaging.

The first set of 352 plants (8 replicates per 44 genotypes) were loaded onto the fully automated plant phenotyping system of Plant Phenomics Victoria, Horsham, ten DAS and arranged in a randomised complete block design (RCBD). The automated plant phenotyping system is housed in a climate-controlled greenhouse and consists of conveyor belts, watering and weighing stations, and an imaging chamber with a Scannalyzer 3D imaging system (LemnaTec GmbH). A second set of 132 plants (44 genotypes, 3 replicates each) were grown on rolling benches for destructive harvest. The growth conditions in the greenhouse were controlled to 24°C during the day and 18°C during the night with a 12 h photoperiod. Enough water was applied automatically to maintain healthy plant growth during the experimental period and recorded into the system’s database (LemnaBase, LemnaTec GmbH).

In a second field experiment, all 44 pea genotypes were trialed in a RCBD design with three replicates during the 2016 winter–spring cropping season at the Plant Breeding Centre of Agriculture Victoria in Horsham, Victoria, Australia (36.74oS, 142.103oE; 133 m altitude). The experimental site has Vertosol heavy clay soil characteristics and a temperate climate with medium average annual rainfall of 450 mm [51]. Seeds were machine sown in plots (1 m width x 5 m length) at a density of 60 plants m-2. Fertilizer application and crop management for weed, pest and disease control were carried out in accordance with the standard practices in the area.

Image capture and processing

After loading onto the automated phenotyping system, plants were imaged daily by the Scannalyzer 3D plant-to-sensor imaging system which consists of two 28.8 megapixel red–green–blue (RGB) cameras (a side and a top camera), model Prosilica GT6600C (Allied Vision Technologies, Stadtroda, Germany). Side-view RGB images were acquired from three sides of the plant after consecutive rotations of 0, 120 and 240 degrees, and a top-view RGB image was taken from above the plant (Fig 1A and 1B). Captured images were automatically recorded in LemnaBase and analyzed by LemnaGrid software (LemnaTec GmbH).

thumbnail
Fig 1. Image acquisition and analysis by the Scannalyzer 3D, plant-to-sensor imaging system.

(a) and (b) are raw images of a side-view and top-view (pea cv. Alma); (c) and (d) processed images showing the identification of the corresponding side-view and top-view objects (bright green plant). Estimated shoot biomass is the pixel sum of highlighted green objects in processed images (c, d). The height of the pink rectangle in (c) is the plant height. The perimeter enveloped by the pink line in (d) is the top-view convex hull.

https://doi.org/10.1371/journal.pone.0207788.g001

The region of interest consisting of the whole plant in raw images was separated from the background by LemnaGrid. In the subsequent steps, the image noise was removed from the region of interest and clear digital plant objects were determined (Fig 1C and 1D). The pixel sums of digital plant objects were generated by LemnaMiner software (LemnaTec GmbH) and subsequently used to estimate several morphological and physiological features of the plants (Table 1).

thumbnail
Table 1. List of field pea traits measured using digital colour imaging and conventional destructive methods.

https://doi.org/10.1371/journal.pone.0207788.t001

Manually destructive harvest

The second set of 132 field pea plants were destructively harvested at 25 DAS after being loaded onto the automated plant phenotyping platform and imaged the night before. Whole plants were weighed using a UniBloc electronic balance (Shimadzu, Melbourne, Australia) to determine fresh shoot biomass per pot (Table 1). The plant height of single plants was determined by measuring from the cut end from immediately above the soil to the tip of the main stalk (Table 1). All leaves from single plants were detached from stalks and leaf area was measured by a Portable Area Meter, model LI-3050A (LI-COR Inc., Lincoln, Nebraska, USA) (Table 1). The remaining 352 field pea plants were unloaded from the automated plant phenotyping platform and destructively harvested at 39 DAS. Fresh shoot biomass was determined as described above.

Normalized difference vegetation index (NDVI) measurements

Early vigour of pea genotypes grown in the field in the second experiment was assessed by a crop growth index NDVI derived from spectral reflectance measured by the Crop Circle sensing equipment (ACS-470; Holland Scientific Inc., Lincoln, NE, USA). NDVI was calculated using the formula from Rouse et al. [54]; (R760 –R670)/(R760+R670), where R670 and R760 are reflectance (R) at 670 nm (VIS region) and 760 nm (NIR region), respectively. Spectral reflectance signals were captured by scanning Crop Circle horizontally 0.75–0.90 m over the plant canopy at 52 DAS as described by Nguyen et al. [51].

Plant growth model and statistical analyses

Since biomass accumulation of cereal crops generally follows a nonlinear growth pattern [55], the “broken-stick” statistical model fitting two straight lines using regression split-line function of GENSTAT statistical software version 18.0 (VSN International Ltd, Hemel Hempstead, UK) was used to identify the linear growth phase of field pea plants as described by Kong et al. [56] and Kholová et al. [57].

Imaging-derived and manually measured data were checked for outliers by using boxplot function of GENSTAT statistical software and presented as means of eight replicates per genotype, with exception to the plants destructively harvested at 25 DAS as this data was a mean of three replicates. One-way analysis of variance (ANOVA) was performed to determine any varietal effects and linear regressions and Pearson’s correlation coefficients (r) were used to determine the relationship between estimated and measured plant traits by using R statistical software (version R-3.5.0) [58].

Results

Validation of nondestructive imaging phenotyping of growth indices

To validate the suitability of image analysis to predict the early vigour phenotype of field pea under controlled environments, we first analysed the estimated values captured through imaging against the measured values from destructive analysis of morphological and physiological parameters of 44 field pea genotypes (Table 1; Fig 1). The results showed that the estimated and measured traits are highly correlated for all 44 field pea genotypes (Fig 2). The most important estimated trait, estimated biomass (EB) is strongly correlated with measured traits such as measured biomass (MB) and leaf area (LA) with high Pearson’s correlation coefficients (r = 0.92 and 0.98, respectively; Fig 2). Similarly, two estimated traits top-view area (TVA) and top-view convex hull (TVCH) were also highly correlated with LA (r = 0.94 and 0.74; Fig 2). Other estimated traits such as estimated height (EH) and estimated water use efficiency (eWUE) also show high correlation with the corresponding measured traits (r = 0.95 and 0.92, respectively; Fig 2). Overall, these estimated and manually measured morphological and physiological parameters are highly correlated.

thumbnail
Fig 2. Relationship between estimated with measured morphological and physiological parameters of pea genotypes.

The cyan panels are the histograms of individual traits. Panels above and below the diagonal of each cyan panel are Pearson’s correlation coefficients (r) and bivariate scatter plots with trend lines, respectively. The asterisks indicate the statistically significant level (*** p<0.001). MB, measured shoot biomass; EB, estimated shoot biomass; LA, measured leaf area; TVA, top-view area; TVCH, top-view convex hull; MH, measured plant height; EH, estimated plant height; mWUE, measured water use efficiency; eWUE, estimated water use efficiency; sample number = 44.

https://doi.org/10.1371/journal.pone.0207788.g002

Dynamic growth analysis of field pea genotypes

Since early seedling vigour is strongly influenced by shoot biomass accumulation during the linear growth phase, we determined the earliest time point where estimated early vigour can be used to compare the performance of all field pea genotypes. Unlike conventionally destructive sampling methods, nondestructive digital imaging allows the calculation and observation of dynamic growth and shoot biomass accumulation of plants over time. Our data showed that the mean EB increased over the period from 11 to 39 DAS (Fig 3).

thumbnail
Fig 3. Dynamic growth of field pea genotypes.

Boxplots of the estimated shoot biomass accumulation over the growth period for 44 field pea genotypes. EB, estimated shoot biomass; DAS, days after sowing. In each box, the black line is the median; the blue squares are the mean; the light green dots are EB of individual field pea genotypes.

https://doi.org/10.1371/journal.pone.0207788.g003

These boxplots showed that the EB of 44 field pea genotypes could be separated into two distinct stages; the lag and the linear phases (Fig 3). Using the broken-stick statistical model, we identified the coordinates X which is the reference point of the days after sowing and Y, the estimated biomass at X of the breakpoints that separate the two growth phases for each genotype of field pea (Table 2). Data showed that the broken-stick model fitted well with the growth of all field pea varieties as indicated by the high adjusted coefficients of determination (R2 >0.99; Table 2). The slope of the regression after the breakpoint (slope 2) of all varieties exceeds that before the breakpoint (slope 1).

thumbnail
Table 2. Regression parameters as determined by the split-line linear regression model of 44 field pea genotypesa.

https://doi.org/10.1371/journal.pone.0207788.t002

Pearson’s correlation analysis between MB, and parameters of the broken-stick model showed that MB was highly negatively correlated with X coordinate, while it was highly positively correlated with Y coordinate, slope 1 and slope 2 (Fig 4).

thumbnail
Fig 4. Relationship between early vigour of field pea genotypes and parameters of the “broken-stick” model.

The cyan panels are the histograms of measured biomass and parameters of the broken-stick model. Panels above and below the diagonal of each cyan panel are Pearson’s correlation coefficients and bivariate scatter plots with trend lines, respectively. The asterisks indicate the statistically significant level (*** p<0.001). MB, measured shoot biomass; X and Y coordinates, the reference points of the days after sowing and the estimated biomass at X; Slope 1, the coefficient before the break point; Slope 2, the coefficient after the break point; sample number = 44.

https://doi.org/10.1371/journal.pone.0207788.g004

Data also showed that X coordinates of several varieties were between 21–23 DAS; such as Alma, Dunn, and Whero, whereas many other varieties had their X coordinates at later dates over 26 DAS; Bluey, King, Maki, Mukta, PBA Pearl, PBA Twilight, Sturt, and Yarrum (Table 2). The latest X coordinates of several field pea genotypes was approximately 26.3 DAS, suggesting that any growth points after this date fell into linear growth phase and only EB values from this point forward should be used for the comparison of early vigour across 44 field pea varieties (Table 2). For consistency, we used the estimated morphological and physiological values collected at 27 DAS hereafter to compare the performance of field pea genotypes in the following sections.

Assessment of early vigour traits of field pea genotypes

To determine how well the estimated traits correlate with early seedling vigour of field pea, we compared the MB harvested at 39 DAS, a time point lying in the linear growth phase, against estimated morphological and physiological values of 44 pea genotypes at 27 DAS (Table 3). Overall, the performance of all varieties estimated morphological and physiological values per genotype were relatively consistent with MB (Table 3, Fig 2). However, there was significant variation among estimated traits, with the most consistent traits relative to MB being EB, TVA, eWUE and to a lesser extent for TVCH, EH, and RGR, while TVCOM was the least consistent trait (Table 3, Fig 2). For example, varieties Alma, Laura, Cressy Blue, and Cooke are the most vigorous genotypes, whereas, Santi, PBA Oura, Mukta, and Yarrum are the least vigorous genotypes in respect of MB (Table 3). However, Laura has a significantly smaller EB, TVA, EH, and eWUE compared to Cressy Blue.

thumbnail
Table 3. Comparative performance of pea genotypes for early vigoura.

https://doi.org/10.1371/journal.pone.0207788.t003

As expected, conventionally leafed field pea genotypes had higher biomass accumulation and TVCOM than semi-leafless lines, indicating their stronger early vigour. Nevertheless, it’s noteworthy to mention that some semi-leafless varieties such as Glenroy and Maitland also had strong early vigour as indicated by high biomass accumulation (Table 3).

Correlation between normalized difference vegetation index captured in the field with estimated shoot biomass and top-view area in controlled environment

To compare the performance of pea genotypes in the field and controlled environment, we identified the association between NDVI values collected from the field trial and EB and TVA derived from imaging-based phenotyping in the Plant Phenomics Victoria, Horsham (Fig 5). Our data showed that NDVI values captured by Crop Circle are highly correlated with EB and TVA for all pea genotypes with correlation coefficients (r) of 0.7 and 0.75, respectively (Fig 5). This suggests that field pea genotypes selected for early vigour using this imaging method from the greenhouse are likely to show their early vigour under field conditions.

thumbnail
Fig 5. Comparative performance of 44 field pea genotypes under field and controlled environment.

The cyan panels are the histograms of estimated biomass (EB) and top-view area (TVA) collected at 27 DAS in Plant Phenomics Victoria, Horsham and the vegetation indices (NDVI) measured by Crop Circle at 52 DAS in the field. Panels above and below the diagonal of each cyan panel are Pearson’s correlation coefficients and bivariate scatter plots with trend lines, respectively. The asterisks indicate the statistically significant level (*** p<0.001). Sample number = 44.

https://doi.org/10.1371/journal.pone.0207788.g005

Discussion

The overall aim of this research was to develop a high-throughput and reliable phenotyping method to rapidly assess early seedling vigour of field pea using digital colour imaging technology in a controlled environment. This is a crucial trait that contributes to biological weed control [27], drought tolerance and water use efficiency [17, 29, 30], and improved carbohydrate reserves [31] and nutrient uptake [26, 34, 59] in crops. Genetic improvement for early vigour has been proposed for many agricultural crops such as rice [38, 60], wheat [24, 61] and maize [62]. Therefore, developing new genotypes with improved early seedling vigour is a priority for field pea breeding programs [6, 17]. However, early seedling vigour is a complex trait [38] and the lack of a robust, high-throughput, and reliable phenotyping method that is powerful enough to dissect the component traits for genetic improvement is a bottleneck [63, 64]. The development and validation of high throughput phenotyping techniques, such as this one, will help to remove this limitation.

Seedling vigour can be broken down into several components including biomass accumulation, canopy coverage, and plant height. These components are generally evaluated visually and destructively [18, 35, 39]. The recent emergence of digital colour imaging techniques for plant phenotyping has provided opportunities to non-destructively evaluate morphological and physiological features of agricultural crops such as field pea, wheat, maize, barley and rice under various growing conditions [43, 48, 6567]. To identify a robust and high-throughput plant phenotyping method applicable to screen field pea for early vigour traits, an experiment on genetically diverse field pea genotypes was conducted on the automated plant phenotyping platform of Plant Phenomics Victoria, Horsham over the winter–spring cropping season of 2016. The digital colour imaging data acquired by the plant phenotyping platform showed that estimated traits such as EB, TVA, EH, are highly correlated with shoot biomass accumulation, leaf area coverage and plant height. This confirms that they can be used as representations of these traits to predict the early vigour of field pea. These high correlations between digitally estimated and conventionally measured morphological and physiological traits were also observed by other studies using a similar automated imaging phenotyping system [55, 68]. The data also demonstrated that the early vigour across field pea genotypes can be quantitatively assessed as early as 27 DAS without destructive analysis of the plants, that can contribute to a cost-effective screen of field pea plants. Economical rapid assessment of less than 30 days will also allow for the rapid assessment of larger populations of accessions from genebanks for gene discovery, or even breeding populations for the development of varieties with a combination of desirable traits including early vigour.

Water use efficiency is one beneficial trait in field pea breeding for water-limited environments [19]. It can be defined as a unit of crop production gain per unit of water supply, and increased WUE is generally associated with higher grain yield [11]. The plant phenotyping platform used here can capture digital imaging data to assess early vigour traits and precisely record the amount of water supplied to each pot, which has also been reported by Ge et al. [69]. This is one advanced feature of an automated plant phenotyping method in comparison to conventional methods, where the former allows precise monitoring and documenting of the supplied water amount. Therefore, eWUE can be directly determined by the estimated biomass and known water usage for each field pea genotype. Armstrong et al [37] used a conventional phenotyping method to study WUE in field peas and discovered that conventionally leafed and vigorous varieties, such as Dundale and Wirrega, had higher WUEs than the most vigorous semi-leafless variety Dinkum, which is in line with findings from the current research. This again confirms the feasibility of the automated imaging phenotyping system for WUE studies. However, it’s important to note that estimated traits do not fully account for the variation observed in measured traits when compared to other imaging methods due to overlapping leaves that cannot be differentiated in the captured images (Fig 1). This limitation of the phenotyping method by imagery has been discussed extensively elsewhere [70]. Some estimated traits are more highly correlated with measured traits than others. These variations suggest that multiple estimated traits should be taken into consideration when comparing early vigour of field pea genotypes to ensure the selection of the best performers.

Most of the semi-leafless field pea genotypes have reduced internode length and this genetic combination is favoured by field pea breeding programs due to improved lodging tolerance and greater air flow through the canopy that results in reduced disease pressure [71]. However, the trait combination of semi-leafless and short internode length has a direct impact on early seedling vigour, leaf area index, and ground cover [35]. Interestingly, some semi-leafless varieties such as Glenroy and Maitland had very strong early seedling vigour. Thus, such early vigour semi-leafless genotypes identified in the present study might possess multiple desirable, inheritable attributes, and should be recommended for selection as parents in field pea breeding programs.

Although greenhouse screening may be a cost-effective method to select the best performers for breeding programs, the screening outcomes need to be further validated under field conditions, as crop performance is highly influenced by environmental conditions such as soil, water availability, and temperature [17, 72, 73]. In the present study, NDVI values as a proxy of early vigour collected from the field trial using Crop Circle were well associated with EB and TVA, confirming that selecting field pea genotypes using the automated imaging phenotyping method were reasonably reliable under field conditions. The stronger correlation between NDVI and TVA was probably due to both metrics being derived from 2D observations downward from sensors to the top of the plants and canopies, respectively; whereas, EB was calculated from 3D observations of plants. Digital cameras have been successfully applied to study early vigour traits and ground cover percentage of field pea and other crops [35, 67, 74]. Therefore, it is possible that similar sensory or imagery systems mounted on ground- or aerial-vehicles could be applied to render the high-throughput phenotyping of pea under field conditions [22, 44, 46, 49, 7577] and this should be included in future investigations.

Generally, early seedling vigour directly contributes to crop establishment, biotic and abiotic stress tolerance and finally, grain yield at harvest [19]. However, no direct association between early seedling vigour with grain yield of crops grown under favourable conditions has been reported [24, 78]. In contrast, there is mounting evidence of a relationship between early seedling vigour and grain yield of field pea under adverse growing conditions such as limited rainfall environments [37]. Early seedling vigour of some crops including field pea is not a trait suitable for all environmental conditions, but it could be a valuable phenotypic trait for targeted environments [50]. In sub-tropical environments where crops rely on water availability stored in the soil profile, early seedling vigour might even be a harmful trait, since plants exhaust reserved water quickly resulting in terminal drought at the reproductive phase. In temperate and Mediterranean type climates such as southern Australia, cropping systems mainly depend on seasonal rainfall where early vigour enhances ground cover, minimizing water losses through run-off and evaporation, and facilitating crop growth [11, 14].

Armstrong and Pate [36] reported on the trialing of six contrasting morphological field pea genotypes at three different rainfall locations in Western Australia. Their results showed that Wirrega, a conventionally leafed and early vigorous variety, outperformed other varieties in the dry regions, which was attributed to better ground cover, long main stem, and larger green area index; whereas, its performance was much poorer than the others under cool and high rainfall environments [36]. In the present study, no grain yield data of the pea genotypes was collected, though it might be noteworthy to investigate the contribution of early vigor to grain yield in future studies, particularly under adverse growing conditions. Apparently, weak vigour genotypes are undesirable for any growth conditions and too vigorous genotypes are also not ideal for some environments. Therefore, genotypes with certain degrees of early vigour should be selected so that they can perform best in targeted environments.

Conclusions

We have developed a high-throughput digital image phenotyping method to assess early seedling vigour of genetically diverse field pea genotypes using an automated plant phenotyping platform. Our results have demonstrated that the imaging method is fully capable of detecting variations of early vigour traits of field pea varieties under controlled environments, and this has been further validated for their comparative performance under field conditions. Therefore, this robust screening method is highly applicable and will enable breeding programs to rapidly identify early vigour traits and utilise germplasm to contribute to the genetic improvement of field peas. To our best knowledge, this is a very first method detailing the application of high-throughput automated imaging phenotyping technology to assess early vigour of field pea under a controlled environment.

Supporting information

S1 Table. Origin of 44 genetically diverse field pea varieties used in this study.

https://doi.org/10.1371/journal.pone.0207788.s001

(DOCX)

Acknowledgments

Special thanks are given to our late colleague, Kevin Murray, for his significant contribution in conducting the greenhouse experiments. Authors wish to thank George Alcock and Nicole Sawyer for assistance in conducting greenhouse experiments; Debra Partington and Matthew Hayes for critical suggestions and comments on crop growth analysis.

References

  1. 1. FAOSTAT. 2014 [30 August 2018]. Available from: http://faostat.fao.org
  2. 2. ABARES. Australian crop report. February 2018, Canberra: Australian Bureau of Agricultural and Resource Economics and Sciences, 2018.
  3. 3. Pulse Australia. Field Pea 2018. Available from: http://www.pulseaus.com.au/growing-pulses/bmp/field-pea. [verified 30 August 2018].
  4. 4. Foyer CH, Lam H-M, Nguyen HT, Siddique KHM, Varshney RK, Colmer TD, et al. Neglecting legumes has compromised human health and sustainable food production. Nat Plants. 2016;2:16112. pmid:28221372
  5. 5. Herridge DF, Peoples MB, Boddey RM. Global inputs of biological nitrogen fixation in agricultural systems. Plant Soil. 2008;311(1):1–18.
  6. 6. Siddique KHM, Erskine W, Hobson K, Knights EJ, Leonforte A, Khan TN, et al. Cool-season grain legume improvement in Australia—use of genetic resources. Crop Pasture Sci. 2013;64(4):347–60. https://doi.org/10.1071/CP13071.
  7. 7. Pandey A, Prasad K, Singh P, Singh R. Comparative yield loss assessment and crop-weed association in major winter crops of Mid Hills of NW Himalayas. Indian J Weed Sci. 1998;30(1 & 2):54–7.
  8. 8. Spies JM, Warkentin TD, Shirtliffe SJ. Variation in field pea (Pisum sativum) cultivars for basal branching and weed competition. Weed Sci. 2011;59(2):218–23.
  9. 9. Lemerle D, Verbeek B, Coombes N. Losses in grain yield of winter crops from Lolium rigidum competition depend on crop species, cultivar and season. Weed Res. 1995;35(6):503–9.
  10. 10. Bajwa AA, Akhter MJ, Iqbal N, Peerzada AM, Hanif Z, Manalil S, et al. Biology and management of Avena fatua and Avena ludoviciana: two noxious weed species of agro-ecosystems. Environ Sci Pollut Res. 2017;24(24):19465–79. pmid:28766148
  11. 11. Siddique KHM, Johansen C, Turner NC, Jeuffroy M-H, Hashem A, Sakar D, et al. Innovations in agronomy for food legumes. A review. Agron Sustain Dev. 2012;32(1):45–64.
  12. 12. Drew EA, Gupta VVSR, Roget DK. Herbicide use, productivity, and nitrogen fixation in field pea (Pisum sativum). Aust J Agric Res. 2007;58(12):1204–14. https://doi.org/10.1071/AR06394.
  13. 13. Jha UC, Bohra A, Parida SK, Jha R. Integrated “omics” approaches to sustain global productivity of major grain legumes under heat stress. Plant Breeding. 2017;136(4):437–59.
  14. 14. Blessing CH, Mariette A, Kaloki P, Bramley H. Profligate and conservative: water use strategies in grain legumes. J Exp Bot. 2018;69(3):349–69. pmid:29370385
  15. 15. Jeuffroy M, Lecoeur J, Roche R. The seed number. In: Munier-Jolain N, Biarnès V, Chaillet I, Lecoeur J, Jeuffroy M, editors. Physiology of the pea crop. Enfield, NH: CRC Press; 2010. p. 104–31.
  16. 16. Guilioni L, Wéry J, Lecoeur J. High temperature and water deficit may reduce seed number in field pea purely by decreasing plant growth rate. Funct Plant Biol. 2003;30(11):1151–64. https://doi.org/10.1071/FP03105.
  17. 17. Sadras VO, Lake L, Chenu K, McMurray LS, Leonforte A. Water and thermal regimes for field pea in Australia and their implications for breeding. Crop Pasture Sci. 2012;63(1):33–44. https://doi.org/10.1071/CP11321.
  18. 18. McDonald GK. Competitiveness against grass weeds in field pea genotypes. Weed Res. 2003;43(1):48–58.
  19. 19. Siddique KHM, Regan KL, Tennant D, Thomson BD. Water use and water use efficiency of cool season grain legumes in low rainfall Mediterranean-type environments. Eur J Agron. 2001;15(4):267–80. https://doi.org/10.1016/S1161-0301(01)00106-X.
  20. 20. Sloane DHG, Gill GS, McDonald GK. The impact of agronomic manipulation of early vigour in wheat on growth and yield in South Australia. Aust J Agric Res. 2004;55(6):645–54. https://doi.org/10.1071/AR03170.
  21. 21. Mahender A, Anandan A, Pradhan SK. Early seedling vigour, an imperative trait for direct-seeded rice: an overview on physio-morphological parameters and molecular markers. Planta. 2015;241(5):1027–50. pmid:25805338
  22. 22. Nguyen GN, Kant S. Improving nitrogen use efficiency in plants: effective phenotyping in conjunction with agronomic and genetic approaches. Funct Plant Biol. 2018;45(6):606–19. https://doi.org/10.1071/FP17266.
  23. 23. Presterl T, Ouzunova M, Schmidt W, Möller EM, Röber FK, Knaak C, et al. Quantitative trait loci for early plant vigour of maize grown in chilly environments. Theor Appl Genet. 2007;114(6):1059–70. pmid:17340099
  24. 24. Rebetzke GJ, Richards RA. Genetic improvement of early vigour in wheat. Aust J Agric Res. 1999;50(3):291–302. https://doi.org/10.1071/A98125.
  25. 25. Richards RA, Lukacs Z. Seedling vigour in wheat—sources of variation for genetic and agronomic improvement. Aust J Agric Res. 2002;53(1):41–50. https://doi.org/10.1071/AR00147.
  26. 26. Liao M, Fillery IRP, Palta JA. Early vigorous growth is a major factor influencing nitrogen uptake in wheat. Funct Plant Biol. 2004;31(2):121–9. http://dx.doi.org/10.1071/FP03060.
  27. 27. Bajwa AA, Walsh M, Chauhan BS. Weed management using crop competition in Australia. Crop Prot. 2017;95:8–13. http://dx.doi.org/10.1016/j.cropro.2016.08.021.
  28. 28. Sardana V, Mahajan G, Jabran K, Chauhan BS. Role of competition in managing weeds: An introduction to the special issue. Crop Prot. 2017;95:1–7. http://dx.doi.org/10.1016/j.cropro.2016.09.011.
  29. 29. Condon AG, Richards RA, Rebetzke GJ, Farquhar GD. Breeding for high water-use efficiency. J Exp Bot. 2004;55(407):2447–60. pmid:15475373
  30. 30. Richards RA. Current and emerging environmental challenges in Australian agriculture: the role of plant breeding. Aust J Agric Res. 2002;53(8):881–92. https://doi.org/10.1071/AR02052.
  31. 31. Palta JA, Kobata T, Turner NC, Fillery IR. Remobilization of carbon and nitrogen in wheat as influenced by postanthesis water deficits. Crop Sci. 1994;34(1):118–24.
  32. 32. Farooq M, Aziz T, Wahid A, Lee D-J, Siddique KHM. Chilling tolerance in maize: agronomic and physiological approaches. Crop Pasture Sci. 2009;60(6):501–16. https://doi.org/10.1071/CP08427.
  33. 33. Pang J, Palta JA, Rebetzke GJ, Milroy SP. Wheat genotypes with high early vigour accumulate more nitrogen and have higher photosynthetic nitrogen use efficiency during early growth. Funct Plant Biol. 2014;41(2):215–22. http://dx.doi.org/10.1071/FP13143.
  34. 34. Liao M, Palta JA, Fillery IRP. Root characteristics of vigorous wheat improve early nitrogen uptake. Aust J Agric Res. 2006;57(10):1097–107. http://dx.doi.org/10.1071/AR05439.
  35. 35. Armstrong E, Pate J. The field pea crop in S.W. Australia. I. Patterns of growth, biomass production and photosynthetic performance in genotypes of contrasting morphology. Aust J Agric Res. 1994;45(7):1347–62. https://doi.org/10.1071/AR9941347.
  36. 36. Armstrong E, Pate J. The field pea crop in S.W. Australia. II. Effects of contrasting morphology and environment on reproductive performance. Aust J Agric Res. 1994;45(7):1363–78. https://doi.org/10.1071/AR9941363.
  37. 37. Armstrong E, Pate J, Tennant D. The field pea crop in South Western Australia—Patterns of water use and root growth in genotypes of contrasting morphology and growth habit. Funct Plant Biol. 1994;21(4):517–32. https://doi.org/10.1071/PP9940517.
  38. 38. 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(18):5555–66. pmid:26022255
  39. 39. Moot DJ. Harvest index variability within and between field pea (Pisum sativum L.) crops. Canterbury, New Zealand: Lincoln University; 1993.
  40. 40. Araus JL, Cairns JE. Field high-throughput phenotyping: the new crop breeding frontier. Trends Plant Sci. 2013;19(0):52–61. http://dx.doi.org/10.1016/j.tplants.2013.09.008.
  41. 41. Fahlgren N, Gehan MA, Baxter I. Lights, camera, action: high-throughput plant phenotyping is ready for a close-up. Curr Opin Plant Biol. 2015;24:93–9. http://dx.doi.org/10.1016/j.pbi.2015.02.006. pmid:25733069
  42. 42. Mulla DJ. Twenty five years of remote sensing in precision agriculture: Key advances and remaining knowledge gaps. Biosyst Eng. 2013;114(4):358–71. https://doi.org/10.1016/j.biosystemseng.2012.08.009.
  43. 43. Chen D, Neumann K, Friedel S, Kilian B, Chen M, Altmann T, et al. Dissecting the phenotypic components of crop plant growth and drought responses based on high-throughput image analysis. Plant Cell. 2014;26(12):4636–55. PMC4311194. pmid:25501589
  44. 44. Virlet N, Sabermanesh K, Sadeghi-Tehran P, Hawkesford MJ. Field Scanalyzer: An automated robotic field phenotyping platform for detailed crop monitoring. Funct Plant Biol. 2017;44(1):143–53. https://doi.org/10.1071/FP16163.
  45. 45. Deery D, Jimenez-Berni J, Jones H, Sirault X, Furbank R. Proximal remote sensing buggies and potential applications for field-based phenotyping. Agronomy. 2014;4(3):349–79.
  46. 46. Shi Y, Thomasson JA, Murray SC, Pugh NA, Rooney WL, Shafian S, et al. Unmanned aerial vehicles for high-throughput phenotyping and agronomic research. PLoS ONE. 2016;11(7):e0159781. pmid:27472222
  47. 47. 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(1):29. pmid:25904970
  48. 48. 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. pmid:25798184.
  49. 49. Roth L, Streit B. Predicting cover crop biomass by lightweight UAS-based RGB and NIR photography: an applied photogrammetric approach. Precision Agriculture. 2018;19(1):93–114.
  50. 50. Sadras VO, Lake L, Leonforte A, McMurray LS, Paull JG. Screening field pea for adaptation to water and heat stress: Associations between yield, crop growth rate and seed abortion. Field Crops Res. 2013;150:63–73. http://doi.org/10.1016/j.fcr.2013.05.023.
  51. 51. Nguyen GN, Panozzo J, Spangenberg G, Kant S. Phenotyping approaches to evaluate nitrogen-use efficiency related traits of diverse wheat varieties under field conditions. Crop Pasture Sci. 2016;67(11):1139–48. https://doi.org/10.1071/CP16173.
  52. 52. Armstrong E, McCaffery D, Carpenter D, Regan K, Harries M, Pritchard I, editors. Field Pea Focus 2003; Griffith, NSW, Australia, 1–3 October 2003 and Northam, WA, 17–18 September 2003: Pulse Australia, Edgecliffe, NSW, Australia).
  53. 53. Hoffmann WA, Poorter H. Avoiding bias in calculations of relative growth rate. Ann Bot. 2002;90(1):37–42. pmid:12125771
  54. 54. Rouse J, Haas J, Schell J, Deering D. Monitoring vegetation systems in the Great Plains with ERTS. In: Freden S, Mercanti E, Becker M, editors. NASA SP-351, Third Earth Resources Technology Satellite-1 Symposium Vol I: Technical presentations. NASA: Washington, DC1974. p. 309–17.
  55. 55. 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(7):1817–32. pmid:25697789
  56. 56. Kong H, Palta JA, Siddique KHM, Stefanova K, Xiong YC, Turner NC. Photosynthesis is reduced, and seeds fail to set and fill at similar soil water contents in grass pea (Lathyrus sativus L.) subjected to terminal drought. J Agron Crop Sci. 2015;201(4):241–52.
  57. 57. Kholová J, Zindy P, Malayee S, Baddam R, Murugesan T, Kaliamoorthy S, et al. Component traits of plant water use are modulated by vapour pressure deficit in pearl millet (Pennisetum glaucum (L.) R.Br.). Funct Plant Biol. 2016;43(5):423–37. https://doi.org/10.1071/FP15115.
  58. 58. R Core Team. R: A language and environment for statistical computing. https://www.r-project.org/. 2017.
  59. 59. Ryan PR, Liao M, Delhaize E, Rebetzke GJ, Weligama C, Spielmeyer W, et al. Early vigour improves phosphate uptake in wheat. J Exp Bot. 2015;66(22):7089–100. pmid:26320241
  60. 60. Rebolledo MC, Dingkuhn M, Péré P, McNally KL, Luquet D. Developmental dynamics and early growth vigour in rice. I. Relationship between development rate (1/Phyllochron) and growth. J Agron Crop Sci. 2012;198(5):374–84.
  61. 61. Rebetzke GJ, Richards RA, Holland JB. Population extremes for assessing trait value and correlated response of genetically complex traits. Field Crops Res. 2017;201:122–32. http://dx.doi.org/10.1016/j.fcr.2016.10.019.
  62. 62. Jha UC, Bohra A, Jha R. Breeding approaches and genomics technologies to increase crop yield under low-temperature stress. Plant Cell Rep. 2017;36(1):1–35. pmid:27878342
  63. 63. Furbank RT, Tester M. Phenomics–technologies to relieve the phenotyping bottleneck. Trends Plant Sci. 2011;16(12):635–44. https://doi.org/10.1016/j.tplants.2011.09.005. pmid:22074787
  64. 64. Rebolledo MC, Luquet D, Courtois B, Henry A, Soulié J-C, Rouan L, et al. Can early vigour occur in combination with drought tolerance and efficient water use in rice genotypes? Funct Plant Biol. 2013;40(6):582–94. https://doi.org/10.1071/FP12312.
  65. 65. Klukas C, Chen D, Pape J-M. Integrated analysis platform: an open-source information system for high-throughput plant phenotyping. Plant Physiol. 2014;165(2):506–18. pmid:24760818
  66. 66. Campbell MT, Du Q, Liu K, Brien CJ, Berger B, Zhang C, et al. A comprehensive image-based phenomic analysis reveals the complex genetic architecture of shoot growth dynamics in rice (Oryza sativa). Plant Genome. 2017;10(2):1–14. pmid:28724075
  67. 67. Armoniené R, Odilbekov F, Vivekanand V, Chawade A. Affordable imaging lab for noninvasive analysis of biomass and early vigour in cereal crops. Biomed Res Int. 2018;2018:9. https://doi.org/10.1155/2018/5713158.
  68. 68. 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(5):e97047. pmid:24823485
  69. 69. Ge Y, Bai G, Stoerger V, Schnable JC. Temporal dynamics of maize plant growth, water use, and leaf water content using automated high throughput RGB and hyperspectral imaging. Comput Electron Agric. 2016;127:625–32. https://doi.org/10.1016/j.compag.2016.07.028.
  70. 70. Li L, Zhang Q, Huang D. A review of imaging techniques for plant phenotyping. Sensors. 2014;14(11):20078. pmid:25347588
  71. 71. Annicchiarico P, Filippi L. A field pea ideotype for organic systems of northern Italy. J Crop Improv. 2007;20(1–2):193–203.
  72. 72. McMurray LS, Davidson JA, Lines MD, Leonforte A, Salam MU. Combining management and breeding advances to improve field pea (Pisum sativum L.) grain yields under changing climatic conditions in south-eastern Australia. Euphytica. 2011;180(1):69–88.
  73. 73. Nguyen GN, Joshi S, Kant S. Water availability and nitrogen use in plants: effects, interaction and underlying molecular mechanisms. In: Hossain MA, Kamiya T, Burritt D, Tran L-SP, Fujiwara T, editors. Plant macro-nutrient use efficiency: molecular and genomic perspectives in crop plants. 1st ed. San Diego, CA, USA: Academic Press; 2017. p. 233–43.
  74. 74. Duan T, Chapman SC, Holland E, Rebetzke GJ, Guo Y, Zheng B. Dynamic quantification of canopy structure to characterize early plant vigour in wheat genotypes. J Exp Bot. 2016;67(15):4523–34. pmid:27312669
  75. 75. Slater AT, Cogan NO, Rodoni BC, Daetwyler HD, Hayes BJ, Caruana B, et al. Breeding differently—the digital revolution: high-throughput phenotyping and genotyping. Potato Res. 2018:1–16. https://doi.org/10.1007/s11540-018-9388-x.
  76. 76. Araus JL, Kefauver SC. Breeding to adapt agriculture to climate change: affordable phenotyping solutions. Curr Opin Plant Biol. 2018. https://doi.org/10.1016/j.pbi.2018.05.003.
  77. 77. Kefauver SC, Vicente R, Vergara-Díaz O, Fernandez-Gallego JA, Kerfal S, Lopez A, et al. Comparative UAV and field phenotyping to assess yield and nitrogen use efficiency in hybrid and conventional barley. Front Plant Sci. 2017;8(1733). pmid:29067032
  78. 78. TeKrony DM, Egli DB. Relationship of seed vigor to crop yield: A review. Crop Sci. 1991;31(3):816–22.