Fenotipado basado en imágenes de alto rendimiento para determinar las respuestas morfológicas y fisiológicas a estreses simples y combinados en patata

Los efectos potenciales del cambio climático, incluido el aumento de la intensidad y la frecuencia de las olas de calor, las inundaciones y las sequías, tienen impactos negativos en los cultivos^{en crecimiento 1}. Es importante comprender la influencia del cambio climático en la variabilidad de los cultivos y las consecuentes fluctuaciones en la producción anual de cultivos². Con el aumento de la población y la demanda de alimentos, mantener el rendimiento de las plantas de cultivo es un desafío, por lo tanto, se requiere urgentemente encontrar cultivos resistentes al clima para el mejoramiento ^3,4. La papa (Solanum tuberosum L.) es uno de los cultivos alimentarios esenciales que contribuye a la seguridad alimentaria mundial debido a su alto valor nutricional y mayor eficiencia en el uso del agua. Sin embargo, la reducción del crecimiento y del rendimiento en condiciones desfavorables es un problema principal, particularmente en las variedades susceptibles ^5,6. Muchos estudios destacaron la importancia de investigar enfoques alternativos para mantener la productividad del cultivo de papa, incluyendo prácticas agrícolas, encontrando genotipos tolerantes y entendiendo el impacto del estrés en el desarrollo y el rendimiento ^7,8,9, que también es muy demandado por los productores (o agricultores) europeos de ^papa10.

Las plataformas de fenotipado automatizado, incluido el fenotipado basado en imágenes, permiten los análisis cuantitativos de la estructura y función de la planta que son esenciales para seleccionar rasgos relevantes de interés^11,12. El fenotipado de alto rendimiento es una técnica avanzada y no invasiva para determinar diversos rasgos morfológicos y fisiológicos de interés de forma reproducible y rápida ¹³. A pesar de que el fenotipo refleja diferencias genotípicas en relación con los efectos ambientales, la comparación de plantas en condiciones controladas con estrés permite vincular la extensa información fenotípica con una condición específica (estrés)¹⁴. El fenotipado basado en imágenes es esencial para describir la variabilidad fenotípica, y también es capaz de detectar un conjunto de rasgos a lo largo del desarrollo de la planta, independientemente del tamaño de la población¹⁵. Por ejemplo, la medición de los rasgos morfológicos, incluida la forma, el tamaño y el índice de color de las hojas mediante sensores de imágenes Rojo-Verde-Azul (RGB), se utiliza para determinar el crecimiento y el desarrollo de las plantas. Además, las mediciones de los rasgos fisiológicos, incluido el rendimiento fotosintético, la temperatura del dosel y la reflectancia de las hojas, se cuantifican utilizando múltiples tipos de sensores, como la fluorescencia de clorofila, el infrarrojo térmico (IR) y las imágenes hiperespectrales¹⁶. Estudios recientes en ambientes controlados mostraron el potencial del uso del fenotipado basado en imágenes para evaluar diferentes mecanismos y respuestas fisiológicas de plantas bajo estreses abióticos como el calor en la papa¹⁷, la sequía en la cebada¹⁸, el arroz¹⁹ y los tratamientos combinados de sequía y calor en el trigo²⁰. A pesar de que el estudio de las respuestas de las plantas a múltiples interacciones de estrés es complejo, los hallazgos revelan nuevos conocimientos en la comprensión de los mecanismos de las plantas para hacer frente a los rápidos cambios en las condiciones climáticas²¹.

Las respuestas fisiológicas y morfológicas de las plantas están directamente influenciadas por las condiciones de estrés abiótico (altas temperaturas, déficit hídrico e inundaciones), lo que resulta en una reducción del rendimiento²². A pesar de que las papas tienen una alta eficiencia en el uso del agua en comparación con otros cultivos, el déficit de agua afecta negativamente la cantidad y calidad del rendimiento debido a la arquitectura de las raíces poco profundas⁵. Dependiendo de la intensidad y duración del nivel de sequía, el índice de área foliar se reduce, y el retraso en el crecimiento del dosel con inhibición de la formación de nuevas hojas es pronunciado durante las etapas posteriores del estrés, lo que lleva a una reducción en la tasa fotosintética²³. El nivel umbral de agua es crítico con exceso de agua o períodos de sequía prolongados, lo que resulta en un efecto negativo en el crecimiento de las plantas y el desarrollo de los tubérculos debido a la limitación de oxígeno, la disminución de la conductividad hidráulica de las raíces y la restricción del intercambio de gases^24,25. Además, las papas son sensibles a las altas temperaturas, donde las temperaturas por encima de los niveles óptimos dan como resultado un retraso en las tasas de iniciación, crecimiento y asimilación de los tubérculos²⁶. Cuando los estreses aparecen en combinación, las regulaciones bioquímicas y las respuestas fisiológicas difieren de las respuestas individuales al estrés, lo que pone de relieve la necesidad de investigar las respuestas de las plantas a las combinaciones de estrés²⁷. Las tensiones combinadas pueden dar lugar a reducciones (aún más) graves en el crecimiento de las plantas y efectos determinantes en los rasgos relacionados con la reproducción²⁸. El impacto de la combinación de estrés depende de la dominancia de cada estrés sobre los demás, lo que lleva a una respuesta mejorada o suprimida de la planta (por ejemplo, la sequía generalmente conduce al cierre de los estomas mientras los estomas están abiertos para permitir el enfriamiento de la superficie de la hoja bajo estrés térmico). Sin embargo, la investigación combinada sobre el estrés aún está emergiendo, y se requieren más investigaciones para comprender mejor la compleja regulación que media las respuestas de las plantas en estas condiciones²⁹. Por lo tanto, este estudio tiene como objetivo resaltar y recomendar un protocolo de fenotipado utilizando múltiples sensores de imagen que pueden ser adecuados para evaluar las respuestas morfofisiológicas y comprender los mecanismos subyacentes del rendimiento general de la papa bajo tratamientos de estrés únicos y combinados. Según la hipótesis, la combinación de múltiples sensores de imagen demostró ser una herramienta valiosa para caracterizar las estrategias tempranas y posteriores durante la respuesta al estrés de las plantas. La optimización del protocolo de fenotipado basado en imágenes será una herramienta interactiva para que los investigadores y mejoradores de plantas encuentren rasgos de interés para la tolerancia al estrés abiótico.

1. Preparación del material vegetal y condiciones de crecimiento

1. Trasplante in vitro de esquejes de papa (Solanum tuberosum L., cv. Lady Rosetta) de cultivo de tejidos a macetas de 250 mL.
2. Llena las macetas con Klasmman Substrate 2 completamente saturado y mantenlas en la cámara de crecimiento en condiciones de poca luz durante 1 semana.
3. Ajuste las condiciones de luz en el nivel del dosel a 160 μmol·^m-2·^s-1 con una combinación de 25% de luz blanca y 35% de infrarrojos utilizando un fotómetro.
4. Trasplante las plantas después de 10 días de crecimiento de los esquejes in vitro en tierra a macetas de 3 L (15,5 cm de diámetro, 20,5 cm de altura).
5. Rellena la maceta de 3 L con 1850 g de Klasmann Sustrato 2: Arena 3:1.
6. Coloque las plantas en la cámara de crecimiento bajo condiciones de luz de 320 μmol·^m-2·^s-1 con una combinación de 55% de luz blanca y 81% de infrarrojos y configúrelo en un régimen de día largo (fotoperíodo de 16 h).
7. Ajuste la temperatura en la cámara de crecimiento a 22 °C/19 °C para el día/noche y la humedad relativa (HR) al 55% durante todo el experimento.
8. Mantener el peso de la maceta al 60% de contenido relativo de agua del suelo (SRWC) como el nivel de control adecuado para mantener el crecimiento y rendimiento^30,31.
   NOTA: Según ensayos anteriores, mantener el contenido volumétrico de agua por encima del 60% promovió el crecimiento de musgo en la superficie del suelo y elevó el riesgo de enfermedades de las plantas. Además, la presencia de musgo podría generar señales positivas engañosas a partir de imágenes de fluorescencia de clorofila, que es difícil de filtrar. Utilice la siguiente ecuación: SRWC% = (FW-DW)/(TW-DW) × 100, donde FW es el peso fresco del suelo, TW es el peso de turgencia y DW es el peso seco³².
   1. Seleccione las muestras de suelo (100 g) de tres bolsas de mezcla Klasmman Substrate 2 diferentes como réplicas y pese el peso fresco del suelo.
   2. Sature la tierra con agua hasta que las macetas retengan el agua sin gotear y pesen el peso de la turgencia del suelo.
   3. Coloque la muestra en el horno a 80 °C durante 3 días hasta que las muestras de suelo estén completamente secas y pesen el peso seco del suelo³³.
9. Coloque las alfombrillas azules en la superficie de la maceta para reducir la evaporación.
   NOTA: El color azul es necesario para restar el fondo del suelo de los píxeles de la planta en la segmentación de la imagen.
10. Seleccionar diez réplicas biológicas por tratamiento.
11. Aleatorizar las macetas durante el riego (en total, 50 macetas).
12. Agregue los soportes azules para sostener las plantas y evitar daños mecánicos al colocarlas en el sistema de fenotipado.

2. Aplicación de estrés

1. En la etapa temprana de tuberización (28 días después del trasplante de los esquejes in vitro ), dividir las plantas en cinco grupos de tratamiento y fenotipar diez plantas por tratamiento (Figura 1).
2. Induzca la tensión simple y combinada a un nivel que no sea perjudicial de la siguiente manera:
   1. En la cámara de crecimiento, mantener las plantas bajo control, los tratamientos de sequía y encharcamiento a 22 °C/19 °C día/noche (paso 1.7), con diferentes porcentajes de SRWC:
      Control (C) con 60% de SRWC durante todo el experimento.
      Sequía (D) con 20% de SRWC gradual durante 7 días, seguido de 1 día de recuperación.
      Encharcamiento (W) con 160% de SRWC durante 5 días, seguido de 10 días de recuperación.
   2. Para mantener el nivel de agua por encima de la superficie del suelo en el tratamiento de encharcamiento, inserte una bolsa de plástico en la maceta vacía y luego coloque la maceta principal con tierra en la segunda maceta preparada.
   3. Coloque las plantas en una cápsula de crecimiento a 30/28 °C día/noche y 55% HR para tratamientos térmicos. Imponga tensiones térmicas simples y combinadas de la siguiente manera:
      1. Para calor (H), mantener la temperatura entre 30 y 28 °C con un 60% de CSRW durante 15 días.
      2. Para el triple estrés de Calor + Sequía + Anegamiento (HDW), exponga las plantas a un tratamiento térmico a 30 °C/28 °C de temperatura día/noche durante los primeros 7 días (manteniendo un 60% de SRWC), seguido de un tratamiento combinado de sequía + calor durante los otros 7 días (20% de SRWC y 30 °C/28 °C) y finalmente exponga las plantas a un estrés de encharcamiento durante 1 día. Para esto último, vuelva a colocar las plantas en la cámara de crecimiento (consulte el paso 1.7 para conocer las condiciones) e induzca el anegamiento al 160% de SRWC durante 1 día.
         NOTA: Las duraciones seleccionadas de los estreses inducidos se basaron en un experimento piloto que mostró efectos de estrés sin impactos perjudiciales con una supervivencia del 100% de las plantas tratadas. En el ambiente de la cámara de crecimiento, la variación de las condiciones ambientales estuvo en el rango ± 0,2 °C para la temperatura y ± 3% para la humedad.

3. Preparación de la planta para el fenotipado

1. Después de encender las luces a las 6:00 a.m. en las cámaras de crecimiento, permita que las plantas se aclimaten bajo las condiciones de luz de crecimiento constante (320 μmol·^m-2·^s-1) durante al menos 2-3 h antes del inicio del protocolo de fenotipado. Esto asegura que la fotosíntesis y la regulación estomática estén en un estado estacionario³⁴.
2. Antes de la medición, transfiera las plantas de su lugar de cultivo al área de amortiguación de crecimiento del sistema de fenotipado utilizado para la carga manual de plantas en el sistema mientras el marcado automatizado está en modo de espera y posicionado dentro del invernadero (Figura complementaria 1, Figura complementaria 2 y Figura complementaria 3).
   NOTA: Las plantas se mantuvieron en el área de amortiguamiento de crecimiento durante el período de fenotipado que duró 3.5 h. En el invernadero, la variación de las condiciones ambientales estuvo en el rango de ± 2 °C para la temperatura, ± 5% para la humedad y 20% de fluctuación en la intensidad de la luz. Por lo tanto, considere que las mediciones deben comenzar de inmediato y ser cortas, evitando la influencia de las condiciones del invernadero en las plantas.
3. En la plataforma de fenotipado, coloque los recipientes en los discos que se mueven automáticamente en una cinta transportadora en intervalos determinados hasta el sensor de imagen de acuerdo con los protocolos de medición especificados en la sección 4.
4. Etiquete cada planta/bandeja con una identificación única para asegurarse de que los datos medidos se asignen a la planta correcta durante todo el experimento.

4. Protocolo de fenotipado

1. Optimizar el protocolo de fenotipado utilizando múltiples sensores de imagen (fluorescencia de clorofila, IR térmico, RGB e imágenes hiperespectrales), permitiendo así la medición simultánea de los parámetros fisiológicos y morfológicos de las plantas (Figura 2).
   NOTA: Dado que las respuestas de las plantas reflejan las condiciones ambientales y los efectos diurnos, es importante considerar la aleatorización de las macetas y realizar el fenotipado dentro del mismo período del día.
2. En la plataforma de fenotipado, asegúrese de que las plantas ingresen al sistema a través de un túnel de adaptación (Figura 2A) donde primero se captura la altura de la planta y luego se ajusta la altura de cada sensor en función de la distancia de trabajo fija.
3. Realice las mediciones en dos rondas según lo justificado en el protocolo de medición utilizando el software.
   1. En la primera ronda, comprenden las mediciones de respuestas fisiológicas cuantificadas como "reacciones rápidas" utilizando fluorescencia de clorofila e imágenes térmicas.
   2. Comienza midiendo los parámetros fisiológicos bajo tratamientos de estrés térmico y luego el resto de tratamientos.
   3. En la segunda ronda, proceda con otras mediciones para evaluar las respuestas más lentas, incluyendo RGB estructural e imágenes hiperespectrales, seguidas de la evaluación del peso y el riego.
4. Durante la etapa de pesaje y riego, defina el peso de referencia para cada planta para permitir el riego y el pesaje automatizados para el tratamiento dado.
   1. Asegúrese de que el peso total de referencia incluya el peso del disco, el inserto ubicado en la cinta transportadora, el soporte azul del soporte, la estera azul, la maceta, el suelo y la biomasa vegetal en el protocolo definido.
   2. Para una medición precisa de la evapotranspiración durante la etapa de pesaje y riego, prepare macetas vacías como referencia. Además, prepare macetas adicionales para corregir el peso de la biomasa de la planta.
5. Para medir 50 plantas, la duración total del protocolo de fenotipado es de 215 min (85 min en la^1ª ronda y 130 min en la^2ª ronda).
6. Fenotipar diariamente todas las plantas en condiciones de control (1 día antes del tratamiento) y luego inducir los tratamientos de estrés para monitorear las respuestas dinámicas y evaluar las fases tempranas y tardías del estrés inducido.

5. Ajuste de la configuración de cada sensor de imagen

1. Imágenes de fluorescencia cinética de clorofila
   NOTA: La fluorescencia cinética de clorofila se utiliza para investigar la capacidad fotosintética de las plantas en respuesta a diferentes condiciones ambientales, incluidos los estreses abióticos, y para proporcionar información valiosa sobre la eficiencia cuántica de la fotoquímica y la disipación de calor (proceso no fotoquímico).
   1. Realizar mediciones de fluorescencia de clorofila en plantas adaptadas a la luz utilizando un protocolo de luz corto para discriminar las respuestas de las plantas bajo diferentes tratamientos.
   2. Aclimatar³⁵ plantas durante 5 min bajo la luz en el túnel de adaptación equipado con LEDs de color blanco frío (6500 K) a 500 μmol·^m-2·^s-1.
      NOTA: Las imágenes de fluorescencia de clorofila son la primera medición después de la adaptación a la luz utilizada para monitorear los cambios en las capacidades fotosintéticas de las plantas.
   3. Seleccione y optimice el protocolo predefinido de acuerdo con el tamaño de la planta y las intensidades de luz requeridas.
   4. Optimice los ajustes de medición, incluidos los ajustes de la cámara y la intensidad de la luz para garantizar la adquisición de una señal fuerte con una relación señal-ruido óptima.
      1. Ajuste la configuración de la cámara, como el obturador (tiempo de exposición, duración de los flashes de medición) y la sensibilidad (ganancia eléctrica de la cámara). Utilice el obturador a 2 ms y la sensibilidad al 12%.
         NOTA: Estos valores se ajustan en función del tamaño y la forma de la hoja y la distancia definida entre la parte superior del dosel y el sensor de imagen.
      2. Ajuste la intensidad de la luz actínica a 500 μmol·^m-2·^s-1 y ajuste el pulso de saturación a 3200 μmol·^m-2·^s-1, que es al menos 6-7 veces mayor que la luz actínica.
   5. Para medir los parámetros en estado estacionario ligero (Lss) (descrito a continuación), las plantas de adaptación a la luz duran 5 minutos antes de las mediciones en el túnel de adaptación a la luz.
   6. Para estimar el rendimiento cuántico del fotosistema II en estado estacionario (PSII) de plantas adaptadas a la luz, seleccione el protocolo de luz corta (Figura 3) y establezca el protocolo de la siguiente manera.
      NOTA: La duración del protocolo fue de 10 s por planta.
      1. Comience la medición encendiendo la luz actínica de color blanco frío a 500 μmol·^m-2·^s-1 durante 3 s para medir la fluorescencia de estado estacionario en la luz (F_t_Lss aka. F_t')
      2. Aplique un pulso de saturación a 3200 μmol·^m-2·^s-1 durante 800 ms para medir la fluorescencia máxima en estado estacionario en luz (F_m_Lss aka. F_m')
      3. Apague la luz actínica, luego encienda la luz roja lejana (735 nm) para permitir que PSII se relaje en la oscuridad durante 800 ms y mida la fluorescencia mínima en luz en estado estacionario (F_o_Lss aka. F_o').
   7. Para calcular los parámetros relativos, utilice un software de análisis de datos que reste el fondo y extraiga los parámetros relevantes.
      NOTA: Los parámetros extraídos del protocolo utilizado son: máxima eficiencia de la fotoquímica PSII de la muestra adaptada a la luz en estado estacionario a la luz determinada como F_v/F_m_Lss aka. F_v'/F_m', rendimiento cuántico del fotosistema II o eficiencia operativa del fotosistema II en estado estacionario de luz definido como QY_Lss aka. φ_PSII = F_q'/F_m', y la fracción de centros de reacción abiertos en PSII (QA oxidado) se determina como qL_Lss = (F_q'/F_v') x (F₀'/F_t').
2. Imágenes térmicas infrarrojas (IR)
   NOTA: Las imágenes térmicas IR se utilizan para la medición no invasiva de la temperatura real del dosel, determinando así las diferentes regulaciones estomáticas. En la unidad de imágenes térmicas IR, una cámara térmica está montada lateralmente en el brazo robótico para medir la temperatura de la cubierta desde la vista lateral.
   1. Para aumentar el contraste de la temperatura de fondo sobre la temperatura del objeto fotografiado durante el procesamiento de la imagen, utilice una pared calefactada controlada automáticamente en el lado opuesto de la cámara térmica para aumentar el contraste. Regule la temperatura de la pared a 8 °C por encima de la temperatura del aire en la unidad de imagen.
      NOTA: Las imágenes térmicas se adquirieron en la oscuridad utilizando el modo de escaneo lineal³⁵.
   2. Después de la adquisición de la imagen, genere una máscara de planta basada en datos de vista lateral RGB y utilícela para corregistrar con los datos térmicos en el análisis de la imagen. Esto garantiza una identificación precisa del objeto escaneado al tiempo que elimina las interferencias de fondo, como el soporte de la planta.
   3. Para evitar la influencia de las condiciones ambientales fluctuantes a lo largo de todo el experimento, calcule el parámetro diferencia de temperatura (delta T o ΔT).
      NOTA: Delta T (ΔT) se define como la diferencia entre la temperatura medida de la superficie de la hoja (el promedio de todos los píxeles de toda la superficie detectada de la planta) y la temperatura del aire ambiente dentro de la caja de imágenes.
3. Imágenes RGB
   NOTA: Las imágenes RGB se basan en cámaras de sistema visual inercial (VIS) que detectan luz en el rango visible de 400 a 700 nm, donde se utilizan para el análisis en profundidad de la morfología de la planta, la arquitectura y la extracción de las características del índice de color.
   1. La unidad de imagen contiene una mesa giratoria para un posicionamiento preciso de la bandeja y, al mismo tiempo, permite la obtención de imágenes desde varios ángulos para vistas laterales.
   2. Configure las imágenes RGB basadas en imágenes de vista lateral para capturar la planta desde tres ángulos (0°, 120° y 240°), que se toman en el modo de escaneo en línea (RGB1) y las imágenes de vista superior en el modo de instantánea (RGB2).
   3. Ambas cámaras tienen una fuente de luz basada en LED, lo que garantiza una iluminación homogénea de la planta fotografiada y, por lo tanto, determina con precisión las características morfológicas y de color.
   4. Extraiga los parámetros calculados mediante el uso de software de análisis de datos.
   5. Para obtener parámetros adicionales basados en las vistas lateral y superior, calcule el volumen de la planta (biomasa digital)³⁶:
      
   6. Calcular la tasa de crecimiento relativo (RGR)³⁷:
      
      Donde T_n y T_n+1 indican el intervalo de tiempo (días).
4. Imágenes hiperespectrales
   NOTA: Las imágenes hiperespectrales se utilizan para la visualización de la reflectancia espectral de las plantas. Los cambios en la reflectancia de las hojas son indicadores de los diferentes estados fisiológicos de la planta en cuestión.
   1. Utilice el sensor de imágenes hiperespectrales para cuantificar la reflectancia del dosel en la parte visible de los espectros de luz, con una cámara hiperespectral de infrarrojo cercano (VNIR) en un rango de 380-900 nm y una cámara infrarroja de longitud de onda corta (SWIR) en un rango de 900-1700 nm.
   2. Las cámaras están montadas en un brazo robótico con una fuente de luz de tubo halógeno implementada (600 W) para una iluminación de muestra homogénea y espectralmente apropiada durante la adquisición de imágenes que se mueven a través del área XZ.
   3. Ambas cámaras funcionan en modo de escaneo lineal y están colocadas en una caja de imágenes aislada de la luz.
   4. Antes de cada ronda de medición, realice dos mediciones de calibración (automáticamente): calibración de corriente oscura y calibración radiométrica utilizando el estándar de teflón de reflectancia espectral.
   5. La imagen de calibración oscura se resta de la imagen de calibración en blanco y sin procesar para eliminar el ruido de corriente oscuro. A continuación, genere la hiperespectral final dividiendo la imagen en bruto después de la sustracción por la imagen de calibración blanca.

6. Exportación de datos y análisis de imágenes

1. Utilice el software del analizador de datos para la extracción automática, la sustracción de fondo y la segmentación de la máscara de planta de la canalización de procesamiento de imágenes (Figura 2B).
2. El software realiza un análisis totalmente automatizado en el que la aplicación de mascarillas, la sustracción de fondo en la que las plantas se aíslan de su fondo y el cálculo de parámetros se procesan como se describe para la imagen RGB³⁸ y la imagen térmica²⁰.
3. Extraiga los parámetros medidos y calculados de los píxeles específicos de la planta definidos por la máscara de planta y la máscara de bandeja generadas por imagen RGB.
4. Si las imágenes no se seleccionaron completamente, lo que puede ocurrir debido a cambios en el verdor de la vegetación durante etapas posteriores de desarrollo o el efecto del tratamiento de estrés, abra la parte de análisis de datos locales en el software y vuelva a ajustar la configuración de la máscara de planta en el software del analizador de datos dependiendo de cada sensor.
5. En el procesamiento de imágenes de fluorescencia de clorofila, ajuste los parámetros de análisis de la configuración de la máscara vegetal de fluorescencia de clorofila (Figura complementaria 4).
   1. Establezca el umbral como Verdadero, lo que significa que la segmentación de la planta se realiza automáticamente.
   2. Establezca el Índice de marco de máscara como Falso, lo que significa que para la detección de máscara de planta , se utiliza el marco Tiempo-Visual tal como se define en el protocolo de fluorescencia de clorofila.
6. En el procesamiento de imágenes térmicas, establezca los parámetros para el análisis de la máscara de plantas (Figura complementaria 5).
   1. Establezca el umbral automático para la generación de máscara de objeto como Falso.
   2. Establezca la máscara de la imagen lateral RGB como True para utilizarla en el análisis.
7. En el procesamiento de imágenes RGB, ajuste los parámetros de análisis de la configuración de la máscara de planta (Figura complementaria 6 y Figura complementaria 7) en función de la especie y la etapa de desarrollo.
   1. Seleccione la fórmula 4 * G-3 * B-R, que es la definición para la generación de máscaras de objetos y define el componente de color utilizado (componente rojo, verde, azul).
      NOTA: Esta fórmula estándar y otras configuraciones pueden cambiar según el tipo de cámara utilizada (vista superior o de tamaño), el tratamiento aplicado y los diferentes cultivares.
   2. Ajuste el umbral utilizado para la conversión de la imagen en escala de grises con un canal verde mejorado a una imagen binaria, determinando la superficie cubierta por la planta.
   3. Ajuste el tamaño medio del filtro utilizado para reducir el ruido y los píxeles no válidos y rellene los que faltan.
   4. Ajuste el tamaño mínimo del objeto en píxeles que se incluirá en el análisis.
   5. Ajuste el tamaño mínimo de los agujeros de los objetos de máscara en píxeles, normalmente decenas de píxeles. Los agujeros más pequeños que este valor se cierran y se toman en los píxeles del objeto.
   6. Establezca Usar reducción de reflexión como verdadero para normalizar los valores RGB en cada píxel.
   7. Establezca Omitir puntos expuestos incorrectos como Verdadero para recortar los píxeles expuestos por encima o por debajo de la máscara de la planta (por ejemplo, omisión de los reflejos de la superficie o los píxeles oscuros donde el ruido es mayor que la señal).
8. Colorea la configuración de segmentación analizada desde RGB utilizando el software analizador de datos para proporcionar información sobre los cambios de color relacionados con las respuestas al estrés y la senescencia de la planta.
   NOTA: El verdor se estima utilizando una gama predefinida de colores que representan todas las etapas del desarrollo de la planta. La intensidad en los canales de color de todos los píxeles correspondientes a la superficie de la planta se agrupó y agrupó para ser utilizada como mapa de color de origen para la segmentación de color.
   1. Proporcione la imagen RGB procesada (fondo eliminado), mapa de colores: lista de tonos para un análisis específico, como entrada en el software.
   2. Para obtener resultados imparciales, realice la selección de tonos utilizando un conjunto de datos de "entrenamiento" y seleccione diferentes etapas de desarrollo y tratamientos.
      NOTA: El análisis guarda los valores R, G y B de cada píxel de cada imagen de este conjunto de datos de entrenamiento.
   3. Defina el número de tonos (seleccione 6 tonos) utilizando el software para la salida de definición de color, que va de 0 a 255 para cada canal.
   4. Proporcione la lista de tonos generados en el software del analizador de datos (Colores).
9. En el procesamiento de imágenes hiperespectrales, procese los datos hiperespectrales adquiridos utilizando el análisis píxel por píxel implementado en el software del analizador hiperespectral, con calibración radiométrica y de ruido oscuro, sustracción de fondo y segmentación de máscaras de plantas. Utilice los espectros promedio y los índices de vegetación para un análisis más detallado.
10. Cree la máscara para la extracción de datos de la imagen hiperespectral de la imagen VNIR tomada por la cámara hiperespectral VNIR. Para el escaneo hiperespectral SWIR, genere una máscara de planta basada en el análisis SWIR.
    1. En VNIR Plant Mask, use la fórmula 1.2 * (2.5 * (R740-R672) -1.3 * (R740-R556)) para visualizar las imágenes donde R representa el valor de reflectancia en una longitud de onda específica (Figura complementaria 8).
    2. En SWIR Plant Mask, utilice la fórmula (R960-R1450)-(R960-R1200) en el procesamiento de imágenes para visualizar las imágenes (Figura complementaria 9).

7. Pesaje y riego

1. Guarde el peso (antes) de regar durante el procedimiento de pesaje y riego. Luego aplique riego y mantenga el peso después de regar también.
2. Riegue las bandejas en el modo de referencia: cada bandeja tenía un valor de referencia almacenado en la base de datos, al que siempre se regaba. Determine la referencia en función del tratamiento.

8. Análisis de datos

1. Analice los datos utilizando ANOVA y prueba de Shapiro.
2. Realice comparaciones por pares entre diferentes tratamientos mediante la prueba de Wilcox por pares en R Studio (versión 4.2.3) utilizando los paquetes (dplyr), (tidyverse), (rstatix) y (ggpubr).
   NOTA: El análisis de la imagen se realizó automáticamente mediante el uso de un software analizador de datos. Para un análisis más detallado de la adquisición de imágenes, utilice el software de análisis de datos específico del sensor.

En este estudio, se utilizó el fenotipado automatizado basado en imágenes para investigar las respuestas morfológicas y fisiológicas de la papa (cv. Lady Rosetta) bajo estrés simple y combinado. El enfoque aplicado mostró las respuestas dinámicas de las plantas en alta resolución espacio-temporal cuando se indujo estrés en la etapa de iniciación del tubérculo. Para evaluar las fases temprana y tardía del estrés, los resultados se presentaron en 3 períodos de tiempo ([0-5 días de fenotipado (DOP)], [6-10 DOP] y [11-15 DOP]) (Figura 1). Hasta 0 DOP, todas las plantas se cultivaron bajo condiciones de control (C), luego de 1-5 DOP, donde se aplicó estrés por encharcamiento (W) y estrés por calor (H). Así, las respuestas se observaron de la siguiente manera: (i) en 0-5 DOP, indicaron el calor inicial y el encharcamiento; (ii) en 6-10 DOP, se observó la sequía temprana (D) y se observó calor combinado y sequía (HD) y (iii) en 11-15 DOP, se observaron los estreses de calor tardío, sequía y calor combinado + sequía + anegamiento (HDW). La recuperación del anegamiento se observó en 6-10 DOP y 11-15 DOP.

Rasgos morfológicos
Se aplicaron imágenes RGB para determinar el efecto de diferentes tensiones y combinaciones en el crecimiento de las plantas aéreas. Los resultados de la Figura 4 muestran que el tratamiento térmico y el estrés de encharcamiento (0-5 DOP) ya causan una reducción del volumen de la planta y de la RGR en comparación con el control. Durante 6-10 DOP, el volumen de plantas y la RGR de las plantas de control aumentaron continuamente, mientras que bajo calor, calor combinado, sequía y anegamiento, este aumento en el volumen de plantas se redujo claramente (Figura 4A). Debido a que las plantas son muy susceptibles al estrés por encharcamiento, se produjo una disminución pronunciada de la RGR (Figura 4B). Durante el estrés por sequía tardía (11-15 DOP), donde el SRWC se mantuvo en un 20%, se observó una clara reducción de la RGR en comparación con el control. Sin embargo, en la fase tardía de la HDW combinada, la aplicación del tratamiento de encharcamiento provocó un aumento de la RGR en el último día de estrés.

Rasgos fisiológicos
Se aplicó la combinación de fenotipado estructural y fisiológico para revelar nuevas respuestas al estrés. El uso de múltiples sensores de imagen permite la determinación de las respuestas fisiológicas en la fase temprana del estrés. Un análisis más detallado de los datos de fluorescencia de la clorofila mostró que el anegamiento estaba afectando negativamente a la eficiencia fotosintética, donde Fv'/Fm' (Fv/Fm_Lss) disminuyó drásticamente en 0-5 DOP y 6-10 DOP, pero se observó una respuesta de recuperación en 11-15 DOP donde Fv'/Fm' aumentó ligeramente (Figura 5A). Durante la fase de estrés tardío (11-15 DOP), se observó una reducción de Fv'/Fm' en sequía y calor y sequía combinados. En las plantas anegadas, la eficiencia operativa de las plantas (QY_Lss también conocida comoφ PSII) fue significativamente menor en comparación con otros tratamientos en 0-5 DOP y 6-10 DOP, pero un ligero aumento a 11-15 DOP, lo que indica la recuperación de la planta (Figura 5B). Además, se determinaron los diferentes mecanismos en la regulación de la eficiencia que contribuye a la protección de PSII mediante el cálculo de la fracción de centros de reacción abiertos en PSII en un estado estacionario ligero (qL_Lss) (Figura 5C). Solo bajo sequía se observó un aumento en la CV, probablemente debido a la fotoinhibición.

Estos hallazgos estuvieron de acuerdo con los datos de RI que reflejaron diferentes mecanismos subyacentes bajo tensiones (Figura 6). Se observó un aumento de deltaT (ΔT) en el anegamiento, reduciendo la tasa de intercambio de gases. Bajo sequía tardía y estrés combinado de calor y sequía, un aumento en el ΔT se debió al cierre de los estomas, considerado una de las principales respuestas para evitar la pérdida excesiva de agua. Por otro lado, se observó una reducción de ΔT bajo tratamientos térmicos a medida que los estomas se abren para mejorar la eficiencia de la transpiración y enfriar la superficie foliar.

Al investigar los datos hiperespectrales, se seleccionaron dos parámetros de los datos VNIR hiperespectrales para evaluar los índices de reflectancia de las hojas, incluido el NDVI como indicador del contenido de clorofila y el PRI como indicador de la eficiencia de la fotosíntesis. Los resultados mostraron una disminución del NDVI y del PRI solo bajo encharcamiento, en relación con la reducción observada en los rasgos morfológicos (Figura 7A,B). Además, a partir de los datos hiperespectrales SWIR utilizados para evaluar el contenido de agua en las plantas, se observó un aumento en el índice de agua en el anegamiento durante 0-5 DOP (Figura 7C). Sin embargo, bajo tratamientos térmicos, se observó una respuesta opuesta donde el índice de agua fue menor que el control. Estos hallazgos estuvieron de acuerdo con un examen de la vegetación a partir de la segmentación de color de la vista superior RGB. Los cambios en la proporción de tonos indican las respuestas al estrés a lo largo del tiempo (Figura 8). El índice de reverdecimiento mostró una reducción en el contenido de pigmentos bajo sequía y HDW combinado en la fase de estrés tardía y una recuperación gradual del tratamiento de encharcamiento. Por lo tanto, el uso de los sensores de imagen múltiple reflejó la correlación de los rasgos morfofisiológicos y permitió la evaluación del rendimiento general de la planta bajo estreses abióticos.

Figura 1: Cronología de la aplicación de los diferentes tratamientos, incluyendo la edad de las plantas en días después del trasplante de los esquejes in vitro . El día 0 de fenotipado (DOP) se midió en condiciones de control (C), y luego se indujeron las diferentes tensiones con diferentes duraciones. De 1 a 5 DOP se aplicó tensión de encharcamiento (W) y la respuesta inicial del tratamiento térmico (H). Los días siguientes 6-10 DOP, donde se presentó la fase inicial de estrés por sequía (D) y estrés combinado por calor y sequía (HD). Durante 11-15 DOP se reflejó la respuesta de las plantas a la fase tardía de sequía y tratamientos térmicos y la aplicación de encharcamiento a HD (HDW) durante 1 día. Haga clic aquí para ver una versión más grande de esta figura.

Figura 2: Esquema que resume el protocolo de fenotipado y análisis de datos. (A) Descripción general del protocolo de fenotipado. Las plantas se transportan al sistema de fenotipado desde las condiciones controladas en la cámara de crecimiento FS-WI (PSI). Las plantas se aclimataron a la luz en la cámara de adaptación a la luz durante 5 min a 500 μmol.m-2.s-1 antes de las mediciones. Se utilizaron múltiples sensores de imagen para determinar los rasgos morfológicos y fisiológicos, seguidos de la estación de pesaje y riego. Dependiendo del tratamiento, las plantas se volvieron a colocar en condiciones controladas, ya sea a 22 °C/19 °C o a 30 °C/28 °C. (B) Extracción y segmentación automáticas de la tubería de procesamiento de imágenes de cada sensor de imágenes. Haga clic aquí para ver una versión más grande de esta figura.

Figura 3: Breve descripción del protocolo de luz para imágenes de fluorescencia de clorofila. El protocolo de medición comenzó encendiendo la luz actínica de color blanco frío para medir la fluorescencia de estado estacionario en la luz (Ft_Lss) y luego aplicando un pulso de saturación para medir la fluorescencia máxima en estado estacionario en la luz (Fm_Lss). La luz actínica se apagó y la luz roja lejana se encendió para determinar la fluorescencia mínima en estado estacionario en la luz (Fo_Lss). La duración del protocolo fue de 10 s por planta. Haga clic aquí para ver una versión más grande de esta figura.

Figura 4: Imágenes RGB utilizadas para la evaluación morfológica. (A) Volumen de planta calculado a partir del área de vistas superior y lateral RGB. (B) Tasa relativa de crecimiento (RGR) durante la etapa de iniciación del tubérculo. Los datos representan valores medios ± desviación estándar (n = 10). Haga clic aquí para ver una versión más grande de esta figura.

Figura 5: Imágenes de fluorescencia de clorofila en plantas adaptadas a la luz. (A) Máxima eficiencia de la fotoquímica PSII de una muestra adaptada a la luz en estado estacionario a la luz (Fv/Fm_Lss). (B) Rendimiento cuántico del fotosistema II o eficiencia operativa del fotosistema II en estado estacionario de luz (QY_Lss). (C) Fracción de centros de reacción abiertos en PSII en estado estacionario ligero (QA oxidado) (qL_Lss). Los datos representan valores medios ± desviación estándar (n = 10). Haga clic aquí para ver una versión más grande de esta figura.

Figura 6: Se utilizaron imágenes térmicas IR para calcular la diferencia entre la temperatura promedio del dosel extraída de las imágenes térmicas IR y la temperatura del aire (ΔT). Los datos representan valores medios ± desviación estándar (n = 10). Haga clic aquí para ver una versión más grande de esta figura.

Figura 7: Imágenes hiperespectrales para determinar los índices de vegetación y el contenido de agua. (A) Índice de Vegetación de Diferencia Normalizada (NDVI). (B) Índice de reflectancia fotoquímica (PRI) calculado a partir de imágenes VNIR. (C) Índice de agua calculado a partir de imágenes SWIR. Los datos representan valores medios ± desviación estándar (n = 10). Haga clic aquí para ver una versión más grande de esta figura.

Figura 8: Índice de enverdecimiento para plantas bajo diferentes tratamientos. El procesamiento de imágenes se basa en la transformación de la imagen RGB original en un mapa de color que consta de 6 tonos definidos. Los datos representan valores medios ± desviación estándar (n = 10). Haga clic aquí para ver una versión más grande de esta figura.

Figura complementaria 1: Intensidad lumínica medida durante los días de fenotipado (DOP). La duración de las mediciones desde las 9:00 am hasta las 12:35 pm. LI_Buff refiere a los datos medianos de 5 sensores de luz distribuidos en el invernadero. Haga clic aquí para descargar este archivo.

Figura suplementaria 2: Humedad relativa (HR) medida durante los días de fenotipado (DOP). La duración de las mediciones desde las 9:00 am hasta las 12:35 pm. RH_Buff refiere a los datos medianos de 5 sensores de humedad distribuidos en el invernadero. RH2 se refiere a la humedad relativa en la cámara de adaptación. Haga clic aquí para descargar este archivo.

Figura complementaria 3: Temperatura medida durante los días de fenotipado (DOP). La duración de las mediciones desde las 9:00 am hasta las 12:35 pm. T_Buff refiere a los datos medianos de 5 sensores de temperatura distribuidos en el invernadero. T2 se refiere a la temperatura en la cámara de adaptación. T3 se refiere a la temperatura de la pared calefactora. T4 se refiere a la temperatura en la unidad de imágenes térmicas IR. Haga clic aquí para descargar este archivo.

Figura complementaria 4: Captura de pantalla del software del analizador de datos que muestra los parámetros ajustados para el análisis de mascarillas vegetales en sensores de imágenes de fluorescencia de clorofila. Haga clic aquí para descargar este archivo.

Figura complementaria 5: Captura de pantalla del software del analizador de datos que muestra los parámetros ajustados para el análisis de máscaras de plantas en sensores de imágenes térmicas infrarrojas. Haga clic aquí para descargar este archivo.

Figura complementaria 6: Captura de pantalla del software del analizador de datos que muestra los parámetros ajustados para el análisis de máscaras de plantas en sensores de imágenes RGB de 1 vista lateral. Haga clic aquí para descargar este archivo.

Figura complementaria 7: Captura de pantalla del software del analizador de datos que muestra los parámetros ajustados para el análisis de máscaras vegetales en sensores de imagen de vista superior RGB2. Haga clic aquí para descargar este archivo.

Figura complementaria 8: Captura de pantalla del software del analizador de datos que muestra los parámetros ajustados para el análisis de la mascarilla vegetal en los sensores de imagen VNIR. Haga clic aquí para descargar este archivo.

Figura complementaria 9: Captura de pantalla del software del analizador de datos que muestra los parámetros ajustados para el análisis de máscaras de plantas en sensores de imágenes SWIR. Haga clic aquí para descargar este archivo.

Los autores declaran que no tienen intereses financieros o relaciones personales contradictorias que puedan haber influido en el trabajo reportado en este artículo.