Measurement of Arabidopsis thaliana Plant Traits Using the PHENOPSIS Phenotyping Platform

Materials and Reagents

1. Double-sided tape (Scotch 12 x 6.3 mm dispensed, permanent, clear)
   
2. Single-sided tape (Crystal clear tape Scotch 19 x 33 mm)
   
3. Template sheets for leaf scanning
   
4. Outer pots (APTE Society, http://www.apte.fr/)
   
5. Inner perforated pots (outer pots hand perforated with a Dremel 4000)
   
6. Pot labels (Point label, soft plastic 1.3 x 6 cm, BIER: http://www.fournitures-horticoles.com/magasin/catproduits.php?idgdf=9)
   
7. Wooden toothpicks
   
8. Microscope slide (Knittel Glass 76 x 26 mm, Starfrost)
   
9. Pencil
   
10. Fine permanent marker (Staedtler permanent, Lumacolor)
    
11. Cylindrical containers for soil drying (50 x 50 mm, Servilab, catalog number: 8668770 )
    
12. Paper bags for plant tissue drying (7 x 12 cm: http://www.beaumont-group.fr/produit/kraft-blanchi-frictionne-neutre-2/)
    
13. Scalpel (Swann-Morton, Carbon Steel Surgical blades)
    
14. Pincers (S Murray^TM Stainless Steel Watchmaker Forceps 11 cm, Fisher Scientific)
    
15. A. thaliana seeds
    
16. Agricultural soil (from Mauguio city, near Montpellier, France: N 43 37′ 01″, E 4 00′ 33″)
    
17. Compost (Neuhaus N2)
    
18. Clear nail polish (Maybelline New York, express Finish 40)
    
19. Ammonium dihydrogen phosphate (H₂PO₄NH₄) (VWR, catalog number: 21305.290 )
    
20. Potassium nitrate (KNO₃) (VWR, catalog number: 26869.291 )
    
21. Fe (E.D.D.H.A) (SEQUESTRENE 138 FE 100 SG)
    
22. Nitric acid (HNO₃) 52% (VWR, catalog number: 20412.362 )
    
23. Boric acid (H₃BO₃) (Merck, catalog number: 1.00165.0100 )
    
24. Manganese(II) sulfate monohydrate (MnSO₄·H₂O) (VWR, catalog number: 25303.233 )
    
25. Copper(II) sulfate pentahydrate (CuSO₄·5H₂O) (Sigma-Aldrich, catalog number: C2857-250G )
    Note: This product has been discontinued.
    
26. Zinc sulfate heptahydrate (ZnSO₄·7H₂O) (Avantor Performance Materials, J.T. Baker, catalog number: 4382-01 )
    
27. Ammonium molybdate tetrahydrate, (NH₄)₆Mo₇O₂₄∙4H₂O (Sigma-Aldrich, catalog number: A7302-100G )
    
28. Nutrient solution (see Recipes)
    
29. Microelement solution (see Recipes)
    

Equipment

1. PHENOPSIS phenotyping platform (Figure 1)
   1. Growth chamber
      
   2. Robotic arm (custom made by Optimalog: https://www.optimalog.com/phenopsis)
      
   3. Precision balance (Precisa, model: XB620C )
      
   4. Watering tubes
      
   5. RBG camera (Prosilica GC1600 (Allied Vision Technologies, model: PROSILICA GC 1600 ) with Computar Varifocal Megapixel M3Z1228C-MP (CBC, Computar, model: M3Z1228C-MP ), 12-36 mm, monture C zoom lens)
      
   6. Controlling computer (Dell)
      
2. Laptop computer for phenological stage notations (HP ProBook 650 G1 (HP Development, model: HP ProBook 650 G1 ), Intel Core i3-4000M Dual Core, 8 GB 1600 DDR3L 2DM)
   
3. Desktop scanner (Epson, model: Epson Perfection 4990 Photo )
   
4. Desktop computer for image analyses (HP Compaq Pro 6300 MT PC+ (HP Development, model: HP Compaq Pro 6300 Microtower ) Intel Core i5-3470 3.2 G 6 M HD 2500 CPU, GB DDR3-1600 DIMM (2 x 4 GB) RAM + 250 GB 7200 RPM 3.5 HDD, AMD Radeon HD 7450 DP)
   
5. Tablet monitor with a pen (Wacom, model: Cintiq 22HD )
   
6. Stereo microscope with camera attachment (Leitz DMRB, Manta G-201B camera)
   
7. Green LED lamp (Bulb LED GU10 Spot 1 W = 10 W green)
   
   
   Figure 1. PHENOPSIS phenotyping platform. A. Overview of the growth chamber; B-C. Robotic arm with RBG camera and watering tube visible; D. Precision balance; E. A screenshot from Optimalog software controlling the platform.
   

Software

1. Microsoft Excel
   
2. ImageJ v1.51 (https://imagej.nih.gov/ij/) (Schneider et al., 2012)
   
3. PHENOPSIS ImageJ macros (http://bioweb.supagro.inra.fr/phenopsis/MacroImageJ.php)
   
4. Application for automaton control:
   OPTIMA PLC software (https://optimalog.com/index.php?q=optimaplc_presentation)
   
5. Camera control: ProsilicaGigE (RGB Camera, Allied Vision); optionally: ThermaCAM Researcher (IR Camera, Flir), ImagingWin (Fluorescence camera, Walz)
   Note: All the software is adapted by Optimalog to be used with Phenopsis.
   
6. R v3.4.1 (https://www.r-project.org/)
   
7. Software for monitoring of environmental conditions:
   Campbell Logger Net (https://www.campbellsci.com/loggernet)
   

Procedure

1. Preparation of the experiment
   1. Design the experiment: experimental conditions (aerial conditions, photoperiod and water deficit treatments if considered), the duration of the experiment, tissue harvests, total number of plants in the platform and numbers of biological replicates. An example experiment could be: 21 A. thaliana genotypes grown in three watering regimes (well-watered, moderate and severe water deficit), 8 biological replicates per genotype per watering conditions. Using four experimental blocks, there will be 2 replicates per block (Figure 2). 
      
   2. Create a randomized experimental design, for example using a customized R script (Figure 2). Pot coordinate randomization is advisable to diminish the effect of the position in the growth chamber which might significantly affect the results, e.g., growth rate and flowering time. Block numbers can be assigned to different areas of the chamber to test for the pot position effect. Assign numbers to pots and enter the data to a spreadsheet.
      
      
      Figure 2. Schematic randomized experimental design in PHENOPSIS for a specific experiment with 21 A. thaliana genotypes, 3 watering regimes, and 8 biological replicates per genotype per watering conditions. Each point corresponds to one pot. Shapes correspond to watering regimes and colors to genotypes. Black lines divide the trays in PHENOPSIS growth chamber, and grey/white shading signify 4 blocks. Randomization reduces the effect of pot position in the growth chamber, and blocks can be used to investigate this effect.
      
      
   3. Prepare two series of pots (inner and outer) labeled with an adhesive tape. Inner pots are perforated with elongated holes, to enhance the soil drying process if necessary (Figure 3).
      
      
      Figure 3. Pots used in PHENOPSIS experiments. Inner pots (left) have extended holes for enhanced soil drying (visible as two lines on the surface of the pot). Holes are nearly identical in all inner pots, thus differences in soil drying dynamics are minimal. Outer pots (right) are not perforated. Source: http://bioweb.supagro.inra.fr/phenopsis/.
      
      
   4. Weigh the outer (P_out) and inner (P_in) pots. Enter the weights of the two empty pots into the spreadsheet.
      
   5. Prepare the substrate: Mix 50% agricultural soil with 50% compost.
      Note: This substrate is dedicated for water deficit experiments, as it provides optimal water retention for such purpose. Other substrates may be used for other types of experiments.
      
   6. Place the inner pot in the outer one and fill the inner pot with the mixed substrate, weigh the pots with the substrate (Total) and enter this weight into the spreadsheet. Collect and weigh one small soil sample (approx. 10 g) per every 10 pots (S^wet).
      
   7. Dry the soil samples at 80 °C for 3 days.
      
   8. Rinse the soil surface with 10 ml of nutrient solution (see Recipes) for better seed imbibition.
      
   9. Sow approximately 5 A. thaliana seeds onto the soil surface using a wet toothpick and put the pots in the growth chamber (20.5 °C and 70% relative air humidity) in the dark for two days.
      
   10. Transfer the pots into the PHENOPSIS growth-chamber and adjust the environmental conditions (vapor pressure deficit, photoperiod, temperature, PAR [max: 200 µmol]).
       
   11. Spray the soil surface in all pots with deionized water 3 times per day until germination.
       
   12. Measure the weights of soil samples (S^dry) and include them in the spreadsheet. Each soil sample will represent 10 pots. Calculate the initial SWC of the substrate as follows:
       
       
   13. Calculate the wet (PS^wet) and dry (PS^dry) weights of soil as follows:
       
       
   14. Calculate the target weights of pots in two combinations for desired SWCs–with single (Target_single) and double (Target_double) pot. SWC will be defined here as g H₂O g^-1 dry soil.
       
       If you plan to apply water deficit treatments, calculate target weights for both well-watered conditions and water deficit conditions.
       
       
2. Management of the experiment
   1. Declare the experiment in the database (PhenopsisDB, http://bioweb.supagro.inra.fr/phenopsis/). Include information about seed stocks, watering regimes and other details of the experiment.
      
   2. Note the germination date for each pot.
      
   3. After germination, adjust the watering conditions in the Optimalog software. This software controls functions of the robotic arm, including measurements of pot weight, watering and imaging, and cycles during the day. The watering can be programmed using a spreadsheet with calculated pot weights for desired conditions. Set 1-2 weighing/watering cycles per day. Set the target SWC to 0.35 g H₂O g^-1 dry soil.
      
      
      Figure 4. Example zenithal rosette image. A. Unprocessed image of A. thaliana rosette; B. Image processed using custom ImageJ macro that applies a color threshold to the picture. The resulting binary mask is created on the basis of green pixels. Source: http://bioweb.supagro.inra.fr/phenopsis/.
      
      
   4. Adjust the properties of zenithal imaging in prosilicaGigE software. Program the robot to take one picture per day after the light was turned on in the growth chamber (Figure 4).
      
   5. Perform phenological stage notations three times per week. Notations take ca. 2 h for 2 people and 504 pots, i.e., one in the chamber looking at the plant phenological stages and the other noting on the laptop. Critical developmental stages should be noted (Boyes et al., 2001), as well as leaf number for each plant. Notations are best stored in spreadsheets.
      
   6. At the 4-leaf stage determined individually for each pot, remove excessive plants, leaving one healthy plant per pot.
      
   7. If the experiment involves different SWC treatments, the target weight should be appropriately adjusted (Target_single) upon the appearance of first true leaves, and outer pots should be removed for the soil drying period. After reaching the desired SWC, reintroduce the outer pot and adjust the target weight to Target_double. This needs to be done individually for each pot.
      
   8. Harvest plant tissue at previously specified growth stages. If metabolomic and transcriptomic analyses are considered, it is advisable to harvest the tissue at the specific time of the day to avoid diurnal fluctuations. Night time harvests can be done under green light in order to avoid activation of photosynthetic processes and induction of light-dependent gene expression.
      
      
3. Biomass partitioning and leaf scans
   1. Remove the flowering stem and cut the rosette near the soil surface. Measure the rosette fresh weight (FW_R).
      
   2. Enclose the rosette between sheets of wet paper towel in the Petri dish, and place the Petri dish in the cold room for 24 h.
      
   3. Carefully remove water droplets from the rosette with a paper towel, and measure the water-saturated (turgid) weight of the rosette (TW_R).
      
   4. Excise all laminas from the rosette and stick them with double-sided tape to an A4 sheet of paper, according to their rank. The template for this sheet includes a calibration square for further analysis, and boxes for sample description (Figure 5). Measure the water-saturated weight of petioles after stripping the rosettes of laminas (TW_P), and stick the petioles to the same sheet of paper.
      
      
      Figure 5. Example leaf scan with indicated weights and final leaf numbers
      
      
   5. Count the final number of rosette (LN) and cauline leaves (L_CN), and write down these numbers.
      
   6. Scan the paper sheets with leaves and collect the laminas and petioles of each plant to separate bags.
      
   7. Place the bags at 80 °C for 3 days, and afterwards measure the dry weights of laminas (DW_L) and petioles (DW_P).
   
   
4. Stomatal density measurement
   1. Remove one leaf blade from the rosette (the largest leaf or a specific leaf rank). Cut it in half along the central nerve and stick both halves to a sheet of paper with double-sided tape, opposite sides up (Figure 6A).
      
   2. Cover the leaf blade with a thin layer of transparent nail polish. The nail polish needs to be not too viscous, and formation of air bubbles should be avoided.
      
   3. After 5 min of drying, stick one-sided tape to the nail polish layer, remove it carefully from the leaf blade, and stick it to the labeled microscope slide.
      
      
      Figure 6. Epidermal imprints. A. Leaf blades after peeling-off the nail polish; B. Example microphotograph of an epidermal imprint.
      

Data analysis

1. Leaf scans
   Analyze the leaf scans using a custom ImageJ macro (http://bioweb.supagro.inra.fr/phenopsis/MacroImageJ.php). The macro will prompt you to click on all the leaves and then calculate areas for each leaf. The sum of these areas for each plant is total leaf area (TLA).
   
   
2. Stomatal traits
   1. Place the slide under a stereo microscope. Take at least 3 good quality representative images from each epidermal imprint with 50-200 stomata visible at the area of the image (I_area) (e.g., Figure 6B).
      
   2. At each image count both the number of stomata (S_num) and the number of the pavement cells (Pav_num). Stomatal counting can be performed using a custom ImageJ macro available online (http://bioweb.supagro.inra.fr/phenopsis/MacroImageJ.php). The provided link also contains adequate documentation with clear step-by-step description of how to perform cell counting.
      
      
3. Rosette projections
   1. Collect zenithal images that are to be used in the analysis. For successful curve fitting, one image per 4 days is required.
      
   2. Extract projected rosette areas (PRA) using custom ImageJ macros (http://bioweb.supagro.inra.fr/phenopsis/MacroImageJ.php).
      
      
4. Calculated traits
   1. Calculate the trait values: rosette dry weight (DW_R), lamina/petiole ratio (DW_L/P), water content (WC), relative water content (RWC), leaf dry matter content (LDMC), and leaf mass per area (LMA):
      
      
   2. Calculate stomatal density (S_dens) and stomatal index (S_I):
      
      
   3. Estimate the dynamics of leaf formation (Vile et al., 2012) using LOESS or curve fitting (Figure 7A):
      
      where
      L_num is leaf number at the time of measurement,
      a is the number of days after the formation of first true leaves,
      a₀ is the time when LN/2 leaves developed,
      LN is the maximum leaf number,
      b is the maximum rate of leaf production.
      
   4. Estimate growth rate parameters (Bac-Molenaar et al., 2015) using curve fitting to projected rosette area data (Figure 7B). Depending on the complexity of the dataset, different functions can be used, such as simple exponential function:
      
      where
      A₀ is initial rosette size,
      r is the growth rate; or more complex sigmoidal function (Gompertz function):
      
      where
      A_max is the final rosette size, and b is the position along the time axis.
      
      
      Figure 7. Dynamic growth measurement of A. thaliana Col-0 accession under well-watered conditions (blue), and moderate (red) and severe (black) water deficit. A. Changes in rosette leaf number. Data are derived from phenological stage notations. Lines and shading correspond to LOESS and standard error respectively. B. Changes in projected rosette area. Lines correspond to Gompertz curves.
      
      
   All traits that can be measured using this protocol are listed in Table 1.
   
   Table 1. Traits measured in a basic PHENOPSIS experiment
   
   
   
   5. The effects of the genotype (G), environment (E) and G x E interaction can be calculated using two-way analysis of variance (ANOVA). Multivariate effects may be analyzed by multivariate ANOVA (MANOVA). For these statistical methods, additional data analysis software is needed.
      In A. thaliana accessions, the genotype effect tends to be strong, what makes direct observation of plant responses difficult. It is then advisable to, in addition to the analysis of absolute trait values, calculate and analyze the response ratios to different conditions as it was done by Rymaszewski et al. (2017). A response ratio is a logarithmized ratio of mean trait value under treatment conditions and mean trait value under control conditions.
      
   6. Working on multiple genotypes and multiple traits, it is often beneficial to reduce the dimensionality using such methods as principal component analysis (PCA). These methods can also be utilized to visualize and better understand the data structure prior to the in-depth analysis.
      

Recipes

Note: Add deionized water to the required volume.

1. Nutrient solution (150 L)
   

   H2PO4NH4	60 ml of 28.76 g/L stock
   KNO3	120 ml of 88.45 g/L stock
   Microelement solution	60 ml
   Fe (E.D.D.H.A)	7.2 g

   Adjust pH to 5.2-5.8 with HNO₃
   
2. Microelement solution (5 L)
   

   H3BO3	1.85 g
   MnSO4∙H2O	1.90 g
   CuSO4∙5H2O	0.03 g
   ZnSO4∙7H2O	1.44 g
   (NH4)6Mo7O24∙4H2O	0.06 g

Acknowledgments

Access to the PHENOPSIS platform was permitted thanks to funds from the European Plant Phenotyping Network (EPPN, grant agreement No. 284443) funded by the FP7 Research Infrastructures Programme of the European Union. This research was also supported by a PRELUDIUM (2012/05/N/NZ9/01396) grant from the National Science Centre, Poland, awarded to W.R. PHENOPSIS developments are supported by FEDER-FSE-IEJ 2014-2020 Languedoc-Roussillon APSEVIR and PHENOPSIS 2.0 projects. Authors declare no conflict of interest.

References

1. Bac-Molenaar, J. A., Granier, C., Keurentjes, J. J. and Vreugdenhil, D. (2016). Genome-wide association mapping of time-dependent growth responses to moderate drought stress in Arabidopsis. Plant Cell Environ 39(1): 88-102.
   
2. Bac-Molenaar, J. A., Vreugdenhil, D., Granier, C. and Keurentjes, J. J. (2015). Genome-wide association mapping of growth dynamics detects time-specific and general quantitative trait loci. J Exp Bot 66(18): 5567-5580.
   
3. Baerenfaller, K., Massonnet, C., Walsh, S., Baginsky, S., Buhlmann, P., Hennig, L., Hirsch-Hoffmann, M., Howell, K. A., Kahlau, S., Radziejwoski, A., Russenberger, D., Rutishauser, D., Small, I., Stekhoven, D., Sulpice, R., Svozil, J., Wuyts, N., Stitt, M., Hilson, P., Granier, C. and Gruissem, W. (2012). Systems-based analysis of Arabidopsis leaf growth reveals adaptation to water deficit. Mol Syst Biol 8: 606.
   
4. Boyes, D. C., Zayed, A. M., Ascenzi, R., McCaskill, A. J., Hoffman, N. E., Davis, K. R. and Gorlach, J. (2001). Growth stage-based phenotypic analysis of Arabidopsis: a model for high throughput functional genomics in plants. Plant Cell 13(7): 1499-1510.
   
5. Fiorani, F. and Schurr, U. (2013). Future scenarios for plant phenotyping. Annu Rev Plant Biol 64: 267-291.
   
6. Granier, C., Aguirrezabal, L., Chenu, K., Cookson, S. J., Dauzat, M., Hamard, P., Thioux, J. J., Rolland, G., Bouchier-Combaud, S., Lebaudy, A., Muller, B., Simonneau, T. and Tardieu, F. (2006). PHENOPSIS, an automated platform for reproducible phenotyping of plant responses to soil water deficit in Arabidopsis thaliana permitted the identification of an accession with low sensitivity to soil water deficit. New Phytol 169(3): 623-635.
   
7. Granier, C. and Vile, D. (2014). Phenotyping and beyond: modelling the relationships between traits. Curr Opin Plant Biol 18: 96-102.
   
8. Rymaszewski, W., Vile, D., Bediee, A., Dauzat, M., Luchaire, N., Kamrowska, D., Granier, C. and Hennig, J. (2017). Stress-related gene expression reflects morphophysiological responses to water deficit. Plant Physiol 174(3): 1913-1930.
   
9. Schneider, C. A., Rasband, W. S. and Eliceiri, K. W. (2012). NIH Image to ImageJ: 25 years of image analysis. Nat Methods 9(7): 671-675.
   
10. Tisné, S., Schmalenbach, I., Reymond, M., Dauzat, M., Pervent, M., Vile, D. and Granier, C. (2010). Keep on growing under drought: genetic and developmental bases of the response of rosette area using a recombinant inbred line population. Plant Cell Environ 33(11): 1875-1887.
    
11. Vile, D., Pervent, M., Belluau, M., Vasseur, F., Bresson, J., Muller, B., Granier, C. and Simonneau, T. (2012). Arabidopsis growth under prolonged high temperature and water deficit: independent or interactive effects? Plant Cell Environ 35(4): 702-718.