Phénotypage basé sur l’image à haut débit pour déterminer les réponses morphologiques et physiologiques à des stress uniques et combinés chez la pomme de terre

Les effets potentiels du changement climatique, y compris l’augmentation de l’intensité et de la fréquence des vagues de chaleur, des inondations et des sécheresses, ont des impacts négatifs sur les cultures^{en croissance 1}. Il est important de comprendre l’influence du changement climatique sur la variabilité des cultures et les fluctuations qui en résultent dans la production agricole annuelle². Avec l’augmentation de la population et de la demande alimentaire, le maintien du rendement des plantes cultivées est un défi, il est donc urgent de trouver des cultures résistantes au climat pour la sélection ^3,4. La pomme de terre (Solanum tuberosum L.) est l’une des cultures vivrières essentielles qui contribue à la sécurité alimentaire mondiale en raison de sa valeur nutritionnelle élevée et de l’efficacité accrue de l’utilisation de l’eau. Cependant, la réduction de la croissance et du rendement dans des conditions défavorables est un problème majeur, en particulier chez les variétés sensibles ^5,6. De nombreuses études ont souligné l’importance d’étudier d’autres approches pour maintenir la productivité des cultures de pommes de terre, y compris les pratiques agricoles, la recherche de génotypes tolérants et la compréhension de l’impact du stress sur le développement et^{le rendement7,8,9}, qui est également très demandé par les producteurs (ou agriculteurs) européens de pommes de ^terre10.

Les plateformes de phénotypage automatisées, y compris le phénotypage basé sur l’image, permettent des analyses quantitatives de la structure et de la fonction des plantes qui sont essentielles pour sélectionner les caractères d’intérêt^{pertinents 11,12}. Le phénotypage à haut débit est une technique avancée et non invasive permettant de déterminer divers traits morphologiques et physiologiques d’intérêt de manière reproductible et rapide ¹³. Bien que le phénotype reflète des différences génotypiques en relation avec les effets environnementaux, la comparaison des plantes dans des conditions contrôlées avec un stress permet de relier les informations de phénotypage étendues à une condition spécifique (stress)¹⁴. Le phénotypage basé sur l’image est essentiel pour décrire la variabilité phénotypique, et il est également capable de cribler un ensemble de traits tout au long du développement des plantes, quelle que soit la taille de la population¹⁵. Par exemple, la mesure des traits morphologiques, y compris la forme, la taille et l’indice de couleur des feuilles à l’aide de capteurs d’imagerie rouge-vert-bleu (RVB), est utilisée pour déterminer la croissance et le développement des plantes. De plus, les mesures des traits physiologiques, y compris les performances photosynthétiques, la température du couvert forestier et la réflectance des feuilles, sont quantifiées à l’aide de plusieurs types de capteurs, tels que la fluorescence de la chlorophylle, l’infrarouge thermique (IR) et l’imagerie hyperspectrale¹⁶. Des études récentes dans des environnements contrôlés ont montré le potentiel de l’utilisation du phénotypage basé sur l’image pour évaluer différents mécanismes et réponses physiologiques des plantes soumises à des stress abiotiques tels que la chaleur chez la pomme de terre¹⁷, la sécheresse dans l’orge¹⁸, le riz¹⁹ et les traitements combinés de sécheresse et de chaleur chez le blé²⁰. Même si l’étude des réponses des plantes à de multiples interactions de stress est complexe, les résultats révèlent de nouvelles perspectives dans la compréhension des mécanismes des plantes pour faire face aux changements rapides des conditions climatiques²¹.

Les réponses physiologiques et morphologiques des plantes sont directement influencées par les conditions de stress abiotique (température élevée, déficit hydrique et inondations), ce qui entraîne une réduction du rendement²². Même si les pommes de terre ont une efficacité d’utilisation de l’eau élevée par rapport à d’autres cultures, le déficit hydrique affecte négativement la quantité et la qualité du rendement en raison de l’architecture racinaire peu profonde⁵. En fonction de l’intensité et de la durée du niveau de sécheresse, l’indice de surface foliaire est réduit, et le retard de croissance de la canopée avec inhibition de la formation de nouvelles feuilles est prononcé au cours des stades ultérieurs du stress, conduisant à une réduction du taux de photosynthèse²³. Le niveau seuil d’eau est critique en cas d’excès d’eau ou de périodes de sécheresse prolongées, ce qui a un effet négatif sur la croissance des plantes et le développement des tubercules en raison de la limitation de l’oxygène, de la diminution de la conductivité hydraulique des racines et de la restriction des échanges gazeux^24,25. De plus, les pommes de terre sont sensibles aux températures élevées où les températures supérieures aux niveaux optimaux entraînent un retard dans l’initiation, la croissance et les taux d’assimilation des tubercules²⁶. Lorsque les stress apparaissent en combinaison, les régulations biochimiques et les réponses physiologiques diffèrent des réponses individuelles au stress, soulignant la nécessité d’étudier les réponses des plantes aux combinaisons de stress²⁷. Les stress combinés peuvent entraîner (encore plus) de réductions sévères de la croissance des plantes et des effets déterminants sur les caractères liés à la reproduction²⁸. L’impact de la combinaison de stress dépend de la dominance de chaque stress sur les autres, ce qui entraîne une réponse accrue ou supprimée de la plante (p. ex., la sécheresse entraîne généralement la fermeture des stomates alors que les stomates sont ouverts pour permettre le refroidissement de la surface des feuilles en cas de stress thermique). Cependant, la recherche combinée sur le stress est encore émergente et d’autres investigations sont nécessaires pour mieux comprendre la régulation complexe qui médie les réponses des plantes dans ces conditions²⁹. Ainsi, cette étude vise à mettre en évidence et à recommander un protocole de phénotypage utilisant plusieurs capteurs d’imagerie qui peut être adapté pour évaluer les réponses morpho-physiologiques et comprendre les mécanismes sous-jacents de la performance globale de la pomme de terre dans le cadre de traitements de stress simples et combinés. Comme on l’avait supposé, la combinaison de plusieurs capteurs d’imagerie s’est avérée être un outil précieux pour caractériser les stratégies précoces et ultérieures lors de la réponse au stress des plantes. L’optimisation du protocole de phénotypage basé sur l’image sera un outil interactif permettant aux chercheurs et aux sélectionneurs de plantes de trouver des caractères d’intérêt pour la tolérance au stress abiotique.

1. Préparation du matériel végétal et conditions de croissance

1. Transpiquer in vitro des boutures de pommes de terre (Solanum tuberosum L., cv. Lady Rosetta) issues d’une culture tissulaire dans des pots de 250 ml.
2. Remplissez les pots avec du substrat Klasmman 2 entièrement saturé et conservez-les dans la chambre de culture dans des conditions de faible luminosité pendant 1 semaine.
3. Ajustez les conditions d’éclairage au niveau de la canopée à 160 μmol·^m-2·^s-1 avec une combinaison de 25 % de lumière blanche et de 35 % d’infrarouge à l’aide d’un posemètre.
4. Repiquez les plantes après 10 jours de culture des boutures in vitro en terre dans des pots de 3 L (15,5 cm de diamètre, 20,5 cm de hauteur).
5. Remplissez le pot de 3 L avec 1850 g de substrat Klasmann 3:1 2 : Sable.
6. Placez les plantes dans la chambre de culture dans des conditions de lumière de 320 μmol·^m-2·^s-1 avec une combinaison de 55% de lumière blanche et de 81% d’infrarouge et réglez-la sur un régime de jours longs (photopériode de 16 h).
7. Réglez la température dans la chambre de culture à 22 °C/19 °C pour le jour et la nuit et l’humidité relative (HR) à 55 % pour toute l’expérience.
8. Maintenir le poids du pot à 60% de la teneur en eau relative du sol (SRWC) comme niveau de contrôle approprié pour maintenir la croissance et le rendement^30,31.
   REMARQUE : D’après des essais antérieurs, le maintien d’une teneur en eau volumétrique supérieure à 60 % a favorisé la croissance de mousse à la surface du sol et a augmenté le risque de maladies des plantes. De plus, la présence de mousse pourrait générer des signaux positifs trompeurs à partir de l’imagerie de fluorescence de la chlorophylle, qui est difficile à filtrer. Utilisez l’équation suivante : SRWC% = (FW-DW)/(TW-DW) × 100, où FW est le poids frais du sol, TW est le poids de turgescence et DW est le poids sec³².
   1. Prélever les échantillons de sol (100 g) dans trois sacs de mélange différents de substrat Klasmman 2 comme répétitions et peser le poids frais du sol.
   2. Saturez le sol avec de l’eau jusqu’à ce que les pots retiennent l’eau sans couler et pesez le poids de la turgescence du sol.
   3. Placez l’échantillon dans le four à 80 °C pendant 3 jours jusqu’à ce que les échantillons de sol soient complètement secs et pesez le poids sec du sol^{à 33}.
9. Placez les tapis bleus sur la surface du pot pour réduire l’évaporation.
   REMARQUE : La couleur bleue est nécessaire pour soustraire le fond de sol des pixels de la plante dans la segmentation de l’image.
10. Sélectionnez dix répétitions biologiques par traitement.
11. Randomisez les pots pendant l’irrigation (au total, 50 pots).
12. Ajoutez les supports bleus pour soutenir les plantes et éviter les dommages mécaniques lors de leur placement dans le système de phénotypage.

2. Application du stress

1. Au stade précoce de la tubération (28 jours après le repiquage des boutures in vitro ), divisez les plantes en cinq groupes de traitement et phénotypez dix plantes par traitement (figure 1).
2. Induire la contrainte simple et combinée à un niveau qui n’est pas préjudiciable comme suit :
   1. Dans la chambre de culture, gardez les plantes sous contrôle, les traitements de sécheresse et d’engorgement à 22 °C/19 °C jour/nuit (étape 1.7), avec différents pourcentages de SRWC :
      Contrôle (C) avec 60 % de SRWC pour l’ensemble de l’expérience.
      Sécheresse (D) avec 20% de SRWC progressivement pendant 7 jours, suivie d'1 jour de récupération.
      Engorgement (W) avec 160 % de SRWC pendant 5 jours, suivi de 10 jours de récupération.
   2. Pour maintenir le niveau d’eau au-dessus de la surface du sol lors du traitement d’engorgement, insérez un sac en plastique dans le pot vide, puis placez le pot principal avec de la terre dans le deuxième pot préparé.
   3. Placez les plantes dans une capsule de croissance à 30/28 °C jour/nuit et à 55 % HR pour les traitements thermiques. Imposer des contraintes thermiques simples et combinées comme suit :
      1. Pour la chaleur (H), maintenez la température de 30 à 28 °C avec 60 % de SRWC pendant 15 jours.
      2. Pour le triple stress Chaleur + Sécheresse + Engorgement (HDW), exposez les plantes à un traitement thermique à une température jour/nuit de 30 °C/28 °C pendant les 7 premiers jours (en gardant 60% de SRWC), suivi d’un traitement combiné sécheresse + chaleur pendant les 7 autres jours (20% de SRWC et 30 °C/28 °C) et enfin exposez les plantes à un stress d’engorgement pendant 1 jour. Pour ce dernier, remettez les plantes dans la chambre de culture (voir l’étape 1.7 pour les conditions) et induisez un engorgement à 160% de SRWC pendant 1 jour.
         REMARQUE : Les durées sélectionnées des stress induits ont été basées sur une expérience pilote qui a montré des effets de stress sans impacts néfastes avec un taux de survie de 100 % des plantes traitées. Dans l’environnement de la chambre de croissance, la variation des conditions environnementales était de l’ordre ± 0,2 °C pour la température et ± 3 % pour l’humidité.

3. Préparation de la plante pour le phénotypage

1. Après l’allumage des lumières à 6h00 dans les chambres de croissance, laisser les plantes s’acclimater dans des conditions de lumière de croissance constante (320 μmol·^m-2·^s-1) pendant au moins 2-3 h avant le début du protocole de phénotypage. Cela garantit que la photosynthèse et la régulation stomatique sont dans un état d’équilibre³⁴.
2. Avant la mesure, transférez les plantes de leur lieu de culture dans la zone tampon de croissance du système de phénotypage utilisé pour le chargement manuel des plantes dans le système pendant que l’évaluation automatisée est en mode veille et positionnée dans la serre (figure supplémentaire 1, figure supplémentaire 2 et figure supplémentaire 3).
   REMARQUE : Les plantes ont été maintenues dans la zone tampon de croissance pendant la période de phénotypage qui a duré 3,5 h. Dans la serre, la variation des conditions environnementales était de l’ordre de ± 2 °C pour la température, ± 5 % pour l’humidité et 20 % de fluctuation de l’intensité lumineuse. Ainsi, considérez que les mesures doivent commencer immédiatement et être courtes, en évitant l’influence des conditions de serre sur les plantes.
3. Dans la plate-forme de phénotypage, placez les pots dans les disques qui se déplacent automatiquement sur un tapis roulant à des intervalles donnés vers le capteur d’imagerie selon les protocoles de mesure spécifiés au paragraphe 4.
4. Étiquetez chaque plante/plateau avec un identifiant unique pour vous assurer que les données mesurées sont attribuées à la bonne plante tout au long de l’expérience.

4. Protocole de phénotypage

1. Optimiser le protocole de phénotypage à l’aide de plusieurs capteurs d’imagerie (fluorescence de la chlorophylle, IR thermique, RVB et imagerie hyperspectrale), permettant ainsi de mesurer simultanément les paramètres physiologiques et morphologiques des plantes (Figure 2).
   REMARQUE : Étant donné que les réponses des plantes reflètent les conditions environnementales et les effets diurnes, il est important d’envisager la randomisation des pots et d’effectuer le phénotypage au cours de la même période de la journée.
2. Dans la plateforme de phénotypage, assurez-vous que les plantes entrent dans le système par un tunnel d’adaptation (Figure 2A) où la hauteur de la plante est d’abord capturée, puis la hauteur de chaque capteur est ajustée en fonction de la distance de travail fixe.
3. Effectuer les mesures en deux tours comme justifié dans le protocole de mesure à l’aide du logiciel.
   1. Dans un premier temps, il s’agit de mesurer les réponses physiologiques quantifiées en tant que « réactions rapides » à l’aide de la fluorescence de la chlorophylle et de l’imagerie thermique.
   2. Commencez par mesurer les paramètres physiologiques sous traitement de stress thermique, puis le reste des traitements.
   3. Dans le deuxième tour, procédez à d’autres mesures pour évaluer les réponses plus lentes, y compris l’imagerie structurelle RVB et hyperspectrale, suivie d’une évaluation du poids et de l’arrosage.
4. Lors de l’étape de pesée et d’arrosage, définissez le poids de référence pour chaque plante afin de permettre l’arrosage et la pesée automatisés selon le traitement donné.
   1. Assurez-vous que le poids de référence total comprend le poids du disque, de l’insert situé sur la bande transporteuse, du support bleu, du tapis bleu, du pot, du sol et de la biomasse végétale dans le protocole défini.
   2. Pour une mesure précise de l’évapotranspiration lors de l’étape de pesée et d’arrosage, préparez des pots vides comme référence. De plus, préparez des pots supplémentaires pour corriger le poids de la biomasse végétale.
5. Pour mesurer 50 plantes, la durée totale du protocole de phénotypage est de 215 min (85 min au^1er tour et 130 min au^2ème tour).
6. Phénotypez quotidiennement toutes les plantes dans des conditions de contrôle (1 jour avant le traitement) puis induisez les traitements de stress pour surveiller les réponses dynamiques et évaluer les phases précoces et tardives du stress induit.

5. Ajuster les paramètres de chaque capteur d’image

1. Imagerie par fluorescence cinétique de chlorophylle
   REMARQUE : La fluorescence cinétique de la chlorophylle est utilisée pour étudier la capacité photosynthétique des plantes en réponse à différentes conditions environnementales, y compris les stress abiotiques, et pour fournir des informations précieuses sur l’efficacité quantique de la photochimie et de la dissipation de la chaleur (processus non photochimique).
   1. Effectuer une mesure de fluorescence de la chlorophylle sur des plantes adaptées à la lumière en utilisant un protocole de lumière courte pour discriminer les réponses des plantes sous différents traitements.
   2. Acclimater³⁵ plantes pendant 5 min sous la lumière dans le tunnel d’adaptation équipé de LED blanc froid (6500 K) à 500 μmol·^m-2·^s-1.
      REMARQUE : L’imagerie par fluorescence de la chlorophylle est la première mesure après l’adaptation à la lumière utilisée pour suivre l’évolution des capacités photosynthétiques des plantes.
   3. Sélectionner et optimiser le protocole prédéfini en fonction de la taille de l’installation et des intensités lumineuses requises.
   4. Optimisez les paramètres de mesure, y compris les paramètres de la caméra et de l’intensité lumineuse pour assurer l’acquisition d’un signal fort avec un rapport signal/bruit optimal.
      1. Ajustez les paramètres de l’appareil photo tels que l’obturateur (temps d’exposition, durée de mesure des flashs) et la sensibilité (gain électrique de l’appareil photo). Utilisez l’obturateur à 2 ms et la sensibilité à 12%.
         REMARQUE : Ces valeurs sont ajustées en fonction de la taille et de la forme des feuilles et de la distance définie entre le haut de la canopée et le capteur d’imagerie.
      2. Réglez l’intensité de la lumière actinique à 500 μmol·^m-2·^s-1 et réglez l’impulsion de saturation à 3200 μmol·^m-2·^s-1, ce qui est au moins 6 à 7 fois plus élevé que la lumière actinique.
   5. Pour mesurer les paramètres à l’état stationnaire de la lumière (Lss) (décrits ci-dessous), adaptez les plantes à la lumière pendant 5 minutes avant les mesures dans le tunnel d’adaptation à la lumière.
   6. Pour estimer le rendement quantique du photosystème à l’état stationnaire II (PSII) des plantes adaptées à la lumière, sélectionnez le protocole de lumière courte (Figure 3) et définissez le protocole comme suit.
      REMARQUE : La durée du protocole était de 10 s par plante.
      1. Commencez la mesure en allumant la lumière actinique blanc froid à 500 μmol·^m-2·^s-1 pendant 3 s pour mesurer la fluorescence à l’état d’équilibre dans la lumière (F_t_Lss aka. F_t')
      2. Appliquez une impulsion de saturation à 3200 μmol·^m-2·^s-1 pendant 800 ms pour mesurer la fluorescence maximale à l’état stationnaire dans la lumière (F_m_Lss aka. F_m')
      3. Éteignez la lumière actinique, puis allumez la lumière rouge lointain (735 nm) pour permettre au PSII de se détendre dans l’obscurité pendant 800 ms et de mesurer la fluorescence minimale à l’état stationnaire dans la lumière (F_o_Lss aka. F_o').
   7. Pour calculer les paramètres relatifs, utilisez un logiciel d’analyse de données qui soustrait l’arrière-plan et extrait les paramètres pertinents.
      REMARQUE : Les paramètres extraits du protocole utilisé sont : l’efficacité maximale de la photochimie PSII d’un échantillon adapté à la lumière à l’état stationnaire déterminée comme F_v/F_m_Lss aka. F_v'/F_m', rendement quantique du photosystème II ou efficacité de fonctionnement du photosystème II à l’état stationnaire léger défini comme QY_Lss aka. φ_PSII = F_q'/F_m', et la fraction de centres de réaction ouverts dans le PSII (QA oxydé) est déterminée comme qL_Lss = (F_q'/F_v') x (F₀'/F_t').
2. Imagerie thermique infrarouge (IR)
   REMARQUE : L’imagerie thermique IR est utilisée pour la mesure non invasive de la température réelle de la canopée, déterminant ainsi les différentes régulations stomatiques. Dans l’unité d’imagerie IR thermique, une caméra thermique est montée latéralement sur le bras robotique pour mesurer la température de la canopée depuis la vue latérale.
   1. Pour augmenter le contraste de la température de fond par rapport à la température de l’objet imagé pendant le traitement de l’image, utilisez une paroi chauffante contrôlée automatiquement sur le côté opposé de la caméra thermique pour augmenter le contraste. Réglez la température de la paroi à 8 °C au-dessus de la température de l’air dans l’unité d’imagerie.
      REMARQUE : Les images thermiques ont été acquises dans l’obscurité en utilisant le mode de balayage linéaire³⁵.
   2. Après l’acquisition de l’image, générez un masque végétal basé sur des données de vue latérale RVB et utilisez-le pour co-enregistrer les données thermiques dans l’analyse d’image. Cela garantit une identification précise de l’objet numérisé tout en éliminant les interférences de fond telles que le support de plante.
   3. Pour éviter l’influence des conditions environnementales fluctuantes tout au long de l’expérience, calculez le paramètre différence de température (delta T ou ΔT).
      REMARQUE : Le delta T (ΔT) est défini comme la différence entre la température mesurée de la surface de la feuille (la moyenne de tous les pixels de l’ensemble de la surface détectée de la plante) et la température de l’air ambiant à l’intérieur de la boîte d’imagerie.
3. Imagerie RVB
   REMARQUE : L’imagerie RVB est basée sur des caméras à système inertiel visuel (VIS) qui détectent la lumière dans la gamme visible de 400 à 700 nm, où elle est utilisée pour l’analyse approfondie de la morphologie et de l’architecture des plantes et l’extraction des caractéristiques de l’indice de couleur.
   1. L’unité d’imagerie contient une table tournante pour un positionnement précis du plateau et permet simultanément une imagerie multi-angle pour les vues latérales.
   2. Réglez l’imagerie RVB sur l’imagerie latérale pour capturer l’usine sous trois angles (0°, 120° et 240°), qui est prise en mode de balayage linéaire (RGB1) et l’imagerie de vue de dessus en mode instantané (RGB2).
   3. Les deux caméras disposent d’une source lumineuse à base de LED, assurant un éclairage homogène de la plante imagée et, ainsi, déterminant avec précision les caractéristiques morphologiques et de couleur.
   4. Extrayez les paramètres calculés à l’aide d’un logiciel d’analyse de données.
   5. Pour des paramètres supplémentaires basés sur les vues latérales et supérieures, calculez le volume de la plante (biomasse numérique)³⁶ :
      
   6. Calculez le taux de croissance relatif (RGR)³⁷ :
      
      Où T_n et T_n+1 indiquent l’intervalle de temps (jours).
4. Imagerie hyperspectrale
   REMARQUE : L’imagerie hyperspectrale est utilisée pour la visualisation de la réflectance spectrale des plantes. Les changements dans la réflectance des feuilles sont des indicateurs de l’état physiologique différent d’une plante donnée.
   1. Utilisez le capteur d’imagerie hyperspectrale pour quantifier la réflectance de la canopée dans la partie visible des spectres lumineux, avec une caméra hyperspectrale dans le visible-proche infrarouge (VNIR) dans la plage 380-900 nm et une caméra infrarouge à courte longueur d’onde (SWIR) dans la plage 900-1700 nm.
   2. Les caméras sont montées sur un bras robotique avec une source lumineuse à tube halogène (600 W) pour un éclairage homogène et spectralement approprié de l’échantillon lors de l’acquisition d’images se déplaçant sur la zone XZ.
   3. Les deux caméras fonctionnent en mode balayage linéaire et sont placées dans un boîtier d’imagerie isolé par la lumière.
   4. Avant chaque tour de mesure, effectuez deux mesures d’étalonnage (automatiquement) : l’étalonnage du courant d’obscurité et l’étalonnage radiométrique à l’aide de l’étalon de téflon à réflectance spectrale.
   5. L’image d’étalonnage sombre est soustraite de l’image d’étalonnage brute et blanche pour supprimer le bruit de courant d’obscurité. Ensuite, générez l’hyperspectral final en divisant l’image brute après soustraction par l’image d’étalonnage blanche.

6. Exportation de données et analyse d’images

1. Utilisez le logiciel d’analyse de données pour l’extraction automatique, la soustraction d’arrière-plan et la segmentation Plant Mask du pipeline de traitement d’images (Figure 2B).
2. Le logiciel effectue une analyse entièrement automatisée où l’application du masque, la soustraction de l’arrière-plan dans laquelle les plantes sont isolées de leur arrière-plan et le calcul des paramètres sont traités comme décrit pour l’imagerie RVB³⁸ et l’imagerie thermique²⁰.
3. Extrayez les paramètres mesurés et calculés à partir des pixels spécifiques à la plante, tels que définis par le masque de plante et le masque de plateau générés par l’image RVB.
4. Si les images n’ont pas été entièrement sélectionnées, ce qui peut se produire en raison de changements dans la verdure de la végétation au cours des stades ultérieurs du développement ou de l’effet du traitement du stress, ouvrez la partie analyse des données locales dans le logiciel et réajustez les paramètres du masque végétal dans le logiciel de l’analyseur de données en fonction de chaque capteur.
5. Dans le traitement d’images par fluorescence de chlorophylle, ajustez les paramètres d’analyse des paramètres de fluorescence de la chlorophylle (figure supplémentaire 4).
   1. Définissez le seuil sur True, ce qui signifie que la segmentation de l’usine est effectuée automatiquement.
   2. Définissez l’index de trame du masque sur False, ce qui signifie que pour la détection du masque végétal , la trame Time-Visual est utilisée comme défini dans le protocole de fluorescence de la chlorophylle.
6. Dans le traitement d’images thermiques, définissez les paramètres pour l’analyse du masque végétal (figure supplémentaire 5).
   1. Définissez le seuil automatique de génération du masque d’objet sur False.
   2. Définissez le masque de l’image latérale RVB sur True à utiliser pour l’analyse.
7. Dans le traitement d’image RVB, ajustez les paramètres d’analyse des paramètres du masque végétal (figure supplémentaire 6 et figure supplémentaire 7) en fonction de l’espèce et du stade de développement.
   1. Sélectionnez la formule 4*G-3*B-R, qui est la définition pour la génération du masque d’objet et définit le composant de couleur utilisé (composant Rouge, Vert, Bleu).
      REMARQUE : Cette formule standard et d’autres paramètres peuvent changer en fonction du type d’appareil photo utilisé (vue de dessus ou vue de la taille), du traitement appliqué et des différents cultivars.
   2. Ajustez le seuil utilisé pour la conversion de l’image en niveaux de gris avec un canal vert amélioré en une image binaire, déterminant la surface couverte par la plante.
   3. Ajustez la taille médiane du filtre utilisée pour réduire le bruit et les pixels invalides et remplissez ceux qui manquent.
   4. Ajustez la taille minimale de l’objet en pixels à inclure dans l’analyse.
   5. Ajustez la taille minimale des trous dans les objets de masque en pixels, généralement des dizaines de pixels. Les trous plus petits que cette valeur sont fermés et insérés dans les pixels de l’objet.
   6. Définissez l’option Utiliser la réduction de la réflexion sur True pour normaliser les valeurs RVB dans chaque pixel.
   7. Définissez l’option Ignorer les mauvais points exposés sur True pour recadrer les pixels sur/sous-exposés du masque végétal (par exemple, l’omission des reflets de surface ou des pixels sombres où le bruit est plus grand que le signal).
8. Colorez les paramètres de segmentation analysés à partir de RVB à l’aide du logiciel d’analyse de données pour fournir des informations sur les changements de couleur liés aux réponses au stress et à la sénescence des plantes.
   REMARQUE : La verdure est estimée à l’aide d’une gamme prédéfinie de couleurs représentant tous les stades de développement des plantes. L’intensité des canaux de couleur de tous les pixels correspondant à la surface de l’usine a été regroupée et regroupée pour être utilisée comme carte de couleurs source pour la segmentation des couleurs.
   1. Fournissez l’image RVB traitée (arrière-plan supprimé), la carte des couleurs - liste de teintes pour une analyse spécifique, en tant qu’entrée dans le logiciel.
   2. Pour obtenir des résultats impartiaux, effectuez la sélection des teintes à l’aide d’un ensemble de données d’entraînement et sélectionnez différents stades de développement et traitements.
      REMARQUE : L’analyse enregistre les valeurs R, V et B de chaque pixel de chaque image à partir de ce jeu de données d’entraînement.
   3. Définissez le nombre de teintes (sélectionnez 6 teintes) à l’aide du logiciel pour la sortie de définition des couleurs, allant de 0 à 255 pour chaque canal.
   4. Fournissez la liste des teintes générées dans le logiciel d’analyse de données (Couleurs).
9. Dans le traitement d’images hyperspectrales, traitez les données hyperspectrales acquises à l’aide de l’analyse pixel par pixel implémentée dans le logiciel de l’analyseur hyperspectral, avec étalonnage radiométrique et du bruit sombre, soustraction de l’arrière-plan et segmentation du masque végétal. Utilisez les spectres moyens et les indices de végétation pour une analyse plus approfondie.
10. Créez le masque pour l’extraction de données à partir de l’image hyperspectrale à partir de l’image VNIR prise par la caméra hyperspectrale VNIR. Pour le balayage hyperspectral SWIR, générez un masque végétal basé sur l’analyse SWIR.
    1. Dans VNIR Plant Mask, utilisez la formule 1.2*(2.5*(R740-R672)-1.3*(R740-R556)) pour visualiser les images où R représente la valeur de réflectance dans une longueur d’onde spécifique (Figure supplémentaire 8).
    2. Dans SWIR Plant Mask, utilisez la formule (R960-R1450)-(R960-R1200) dans le traitement d’images pour visualiser les images (Figure supplémentaire 9).

7. Pesée et arrosage

1. Stockez le poids (avant) l’arrosage pendant la procédure de pesée et d’arrosage. Appliquez ensuite un arrosage et gardez également le poids après l’arrosage.
2. Arrosez les plateaux en mode référence - chaque plateau avait une valeur de référence stockée dans la base de données, à laquelle il était toujours arrosé. Déterminer la référence en fonction du traitement.

8. Analyse des données

1. Analysez les données à l’aide de l’ANOVA et du test de Shapiro.
2. Effectuez des comparaisons par paires entre différents traitements par test Pairwise Wilcox dans R studio (version 4.2.3) à l’aide des packages (dplyr), (tidyverse), (rstatix) et (ggpubr).
   REMARQUE : L’analyse de l’image a été effectuée automatiquement à l’aide d’un logiciel d’analyse de données. Pour une analyse plus approfondie de l’acquisition d’images, utilisez un logiciel d’analyse de données spécifique au capteur.

Dans cette étude, le phénotypage automatisé basé sur l’image a été utilisé pour étudier les réponses morphologiques et physiologiques de la pomme de terre (cv. Lady Rosetta) en stress unique et combiné. L’approche appliquée a montré les réponses dynamiques des plantes en haute résolution spatio-temporelle lorsque le stress était induit au stade d’initiation des tubercules. Pour évaluer les phases précoces et tardives du stress, les résultats ont été présentés sur 3 périodes ([0-5 jours de phénotypage (DOP)], [6-10 DOP] et [11-15 DOP]) (Figure 1). Jusqu’à 0 DOP, toutes les plantes ont été cultivées dans des conditions de contrôle (C), puis de 1 à 5 DOP, où un stress d’engorgement (W) et un stress thermique (H) ont été appliqués. Ainsi, les réponses ont été observées comme suit : (i) dans 0-5 DOP, indique la chaleur et l’engorgement initiaux ; (ii) dans 6-10 DOP, reflétait la sécheresse précoce (D) et la chaleur et la sécheresse combinées (HD) ont été observées et (iii) dans 11-15 DOP, montrait les stress de chaleur tardive, de sécheresse et de chaleur combinée + sécheresse + engorgement (HDW). La récupération de l’engorgement a été observée dans 6-10 DOP et 11-15 DOP.

Traits morphologiques
L’imagerie RVB a été appliquée pour déterminer l’effet de différents stress et combinaisons sur la croissance des plantes au-dessus du sol. Les résultats de la figure 4 montrent que le traitement thermique et le stress d’engorgement (0-5 DOP) entraînent déjà une réduction du volume de la plante et du RGR par rapport au contrôle. Au cours de la période 6-10 DOP, le volume des plantes et le RGR des plantes témoins ont continuellement augmenté, tandis qu’en cas de chaleur, de chaleur combinée, de sécheresse et d’engorgement, cette augmentation du volume des plantes a été clairement réduite (Figure 4A). Comme les plantes sont très sensibles au stress d’engorgement, une diminution a été prononcée dans le RGR (figure 4B). Au cours du stress hydrique tardif (11-15 DOP), où le SRWC a été maintenu à 20%, une nette réduction du RGR a été observée par rapport au contrôle. Cependant, dans la phase tardive de l’HDW combiné, l’application du traitement de l’engorgement a provoqué une augmentation du RGR le dernier jour de stress.

Traits physiologiques
La combinaison de phénotypage structurel et physiologique a été appliquée pour révéler d’autres réponses au stress. L’utilisation de plusieurs capteurs d’imagerie permet de déterminer les réponses physiologiques dans la phase précoce du stress. Une analyse plus poussée des données de fluorescence de la chlorophylle a montré que l’engorgement affectait négativement l’efficacité photosynthétique où Fv'/Fm'(Fv/Fm_Lss) diminuait considérablement dans les degrés 0-5 et 6-10 DOP, mais une réponse récupératrice était observée dans les niveaux 11-15 DOP où Fv'/Fm’augmentait légèrement (figure 5A). Au cours de la phase de stress tardive (11-15 DOP), une réduction de Fv'/Fm’a été observée en cas de sécheresse et de chaleur et de sécheresse combinées. Dans les usines gorgées d’eau, l’efficacité opérationnelle des plantes (QY_Lss aliasφ PSII) était significativement plus faible par rapport aux autres traitements à 0-5 DOP et 6-10 DOP, mais une légère augmentation à 11-15 DOP, indiquant ainsi une récupération de la plante (figure 5B). De plus, les différents mécanismes de régulation de l’efficacité contribuant à la protection du PSII ont été déterminés en calculant la fraction de centres de réaction ouverts dans le PSII à l’état stationnaire léger (qL_Lss) (Figure 5C). Ce n’est qu’en période de sécheresse qu’une augmentation de qL a été observée, probablement en raison de la photoinhibition.

Ces résultats étaient conformes aux données IR qui reflétaient différents mécanismes sous-jacents sous stress (figure 6). Une augmentation du deltaT (ΔT) a été observée dans l’engorgement, réduisant le taux d’échange gazeux. En cas de sécheresse tardive et de stress combiné de chaleur et de sécheresse, une augmentation de ΔT était due à la fermeture des stomates, considérée comme l’une des principales réponses pour éviter les pertes d’eau excessives. D’autre part, une réduction de ΔT sous les traitements thermiques a été observée lorsque les stomates s’ouvrent pour améliorer l’efficacité de la transpiration et refroidir la surface des feuilles.

En examinant les données hyperspectrales, deux paramètres ont été sélectionnés à partir des données VNIR hyperspectrales pour évaluer les indices de réflectance foliaire, y compris le NDVI comme indicateur de la teneur en chlorophylle et le PRI comme indicateur de l’efficacité de la photosynthèse. Les résultats ont montré une diminution du NDVI et du PRI uniquement en cas d’engorgement, en lien avec la réduction observée des traits morphologiques (Figure 7A,B). De plus, d’après les données hyperspectrales SWIR utilisées pour évaluer la teneur en eau dans les plantes, une augmentation de l’indice d’eau dans l’engorgement a été observée pendant 0-5 DOP (Figure 7C). Cependant, lors des traitements thermiques, une réponse opposée a été observée lorsque l’indice d’eau était inférieur à celui du témoin. Ces résultats étaient en accord avec un examen de la végétation à partir de la segmentation des couleurs de la vue de dessus RVB. Les changements dans la proportion de teintes indiquent les réponses au stress au fil du temps (figure 8). L’indice de verdissement a montré une réduction de la teneur en pigments en cas de sécheresse et combinait HDW à la fin de la phase de stress et récupération progressive après le traitement par engorgement. Ainsi, l’utilisation des multiples capteurs d’imagerie a reflété la corrélation des traits morpho-physiologiques et a permis d’évaluer la performance globale de la plante sous stress abiotique.

Figure 1 : Chronologie de l’application des différents traitements, y compris l’âge des plantes en jours après le repiquage des boutures in vitro. Le jour 0 du phénotypage (DOP) a été mesuré dans des conditions de contrôle (C), puis les différents stress ont été induits avec des durées différentes. De 1 à 5 DOP, une contrainte d’engorgement (W) a été appliquée et la réponse initiale du traitement thermique (H). Les jours suivants, du 6 au 10 DOP, où la phase initiale du stress hydrique (D) et du stress combiné thermique et sécheresse (HD) a été présentée. Au cours de la période 11-15 DOP, la réponse des plantes à la phase tardive de sécheresse et aux traitements thermiques et à l’application de l’engorgement à l’HD (HDW) pendant 1 jour a été reflétée. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 2 : Schéma résumant le protocole de phénotypage et l’analyse des données. (A) Vue d’ensemble du protocole de phénotypage. Les plantes sont transportées vers le système de phénotypage à partir des conditions contrôlées de la chambre de croissance FS-WI (PSI). Les plantes ont été acclimatées à la lumière dans la chambre d’adaptation à la lumière pendant 5 min à 500 μmol.m-2.s-1 avant les mesures. Plusieurs capteurs d’imagerie ont été utilisés pour déterminer les traits morphologiques et physiologiques, suivis de la station de pesée et d’abreuvement. En fonction du traitement, les plantes ont été remises dans des conditions contrôlées, soit à 22 °C/19 °C, soit à 30 °C/28 °C. (B) Extraction et segmentation automatiques du pipeline de traitement d’images de chaque capteur d’imagerie. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 3 : Aperçu du protocole de lumière brève, pour l’imagerie par fluorescence de la chlorophylle. Le protocole de mesure a commencé par l’allumage de la lumière actinique blanc froid pour mesurer la fluorescence stationnaire dans la lumière (Ft_Lss), puis l’application d’une impulsion de saturation pour mesurer la fluorescence maximale stationnaire dans la lumière (Fm_Lss). La lumière actinique a été éteinte et la lumière rouge lointaine a été allumée pour déterminer la fluorescence minimale à l’état d’équilibre dans la lumière (Fo_Lss). La durée du protocole était de 10 s par plante. Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 4 : Imagerie RVB utilisée pour l’évaluation morphologique. (A) Volume de la plante calculé à partir de la zone des vues supérieures et latérales RVB. (B) Taux de croissance relatif (RGR) au stade d’initiation des tubercules. Les données représentent des valeurs moyennes ±écart-type (n = 10). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 5 : Imagerie par fluorescence de la chlorophylle sur des plantes adaptées à la lumière. (A) Efficacité maximale de la photochimie PSII d’un échantillon adapté à la lumière à l’état stationnaire (Fv/Fm_Lss). (B) Rendement quantique du photosystème II ou efficacité de fonctionnement du photosystème II à l’état stationnaire léger (QY_Lss). (C) Fraction de centres de réaction ouverts dans le PSII à l’état stationnaire léger (QA oxydé) (qL_Lss). Les données représentent des valeurs moyennes ±écart-type (n = 10). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 6 : L’imagerie IR thermique a été utilisée pour calculer la différence entre la température moyenne de la canopée extraite des images IR thermiques et la température de l’air (ΔT). Les données représentent des valeurs moyennes ±écart-type (n = 10). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 7 : Imagerie hyperspectrale pour la détermination des indices de végétation et de la teneur en eau. (A) Indice de végétation par différence normalisée (NDVI). (B) Indice de réflectance photochimique (PRI) calculé à partir de l’imagerie VNIR. (C) Indice d’eau calculé à partir de l’imagerie SWIR. Les données représentent des valeurs moyennes ±écart-type (n = 10). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure 8 : Indice de verdissement des plantes sous différents traitements. Le traitement de l’image est basé sur la transformation de l’image RVB d’origine en une palette de couleurs composée de 6 teintes définies. Les données représentent des valeurs moyennes ±écart-type (n = 10). Veuillez cliquer ici pour voir une version agrandie de cette figure.

Figure supplémentaire 1 : Intensité lumineuse mesurée pendant les jours de phénotypage (DOP). La durée des mesures est de 9h00 à 12h35. LI_Buff correspond aux données médianes de 5 capteurs de lumière répartis dans la serre. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 2 : Humidité relative (HR) mesurée pendant les jours de phénotypage (DOP). La durée des mesures de 9h00 à 12h35 RH_Buff correspond aux données médianes de 5 capteurs d’humidité répartis dans la serre. RH2 fait référence à l’humidité relative dans la chambre d’adaptation. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 3 : Température mesurée pendant les jours de phénotypage (DOP). La durée des mesures de 9h00 à 12h35 T_Buff fait référence aux données médianes de 5 capteurs de température répartis dans la serre. T2 fait référence à la température dans la chambre d’adaptation. T3 fait référence à la température de la paroi chauffante. T4 fait référence à la température dans l’unité d’imagerie IR thermique. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 4 : Capture d’écran d’un logiciel d’analyse de données montrant les paramètres ajustés pour l’analyse du masque végétal dans les capteurs d’imagerie par fluorescence de chlorophylle. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 5 : Capture d’écran d’un logiciel d’analyse de données montrant les paramètres ajustés pour l’analyse des masques végétaux dans les capteurs d’imagerie infrarouge thermique. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 6 : Capture d’écran d’un logiciel d’analyse de données montrant les paramètres ajustés pour l’analyse du masque végétal dans les capteurs d’imagerie RVB à 1 vue latérale. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 7 : Capture d’écran d’un logiciel d’analyse de données montrant les paramètres ajustés pour l’analyse du masque végétal dans les capteurs d’imagerie RGB2 à vue de dessus. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 8 : Capture d’écran d’un logiciel d’analyse de données montrant les paramètres ajustés pour l’analyse des masques végétaux dans les capteurs d’imagerie VNIR. Veuillez cliquer ici pour télécharger ce fichier.

Figure supplémentaire 9 : Capture d’écran d’un logiciel d’analyse de données montrant les paramètres ajustés pour l’analyse des masques végétaux dans les capteurs d’imagerie SWIR. Veuillez cliquer ici pour télécharger ce fichier.

Les auteurs déclarent qu’ils n’ont pas d’intérêts financiers concurrents connus ou de relations personnelles qui auraient pu sembler influencer les travaux rapportés dans cet article.