Literature DB >> 33497519

Red light imaging for programmed cell death visualization and quantification in plant-pathogen interactions.

Sergio Landeo Villanueva1, Michele C Malvestiti1, Wim van Ieperen2, Matthieu H A J Joosten1, Jan A L van Kan1.   

Abstract

Studies on plant-pathogen interactions often involve monitoring disease symptoms or responses of the host plant to pathogen-derived immunogenic patterns, either visually or by staining the plant tissue. Both these methods have limitations with respect to resolution, reproducibility, and the ability to quantify the results. In this study we show that red light detection by the red fluorescent protein (RFP) channel of a multipurpose fluorescence imaging system that is probably available in many laboratories can be used to visualize plant tissue undergoing cell death. Red light emission is the result of chlorophyll fluorescence on thylakoid membrane disassembly during the development of a programmed cell death process. The activation of programmed cell death can occur during either a hypersensitive response to a biotrophic pathogen or an apoptotic cell death triggered by a necrotrophic pathogen. Quantifying the intensity of the red light signal enables the magnitude of programmed cell death to be evaluated and provides a readout of the plant immune response in a faster, safer, and nondestructive manner when compared to previously developed chemical staining methodologies. This application can be implemented to screen for differences in symptom severity in plant-pathogen interactions, and to visualize and quantify in a more sensitive and objective manner the intensity of the plant response on perception of a given immunological pattern. We illustrate the utility and versatility of the method using diverse immunogenic patterns and pathogens.
© 2021 The Authors. Molecular Plant Pathology published by British Society for Plant Pathology and John Wiley & Sons Ltd.

Entities:  

Keywords:  RFP; chlorophyll; hypersensitive response; programmed cell death; thylakoid disassembly

Mesh:

Year:  2021        PMID: 33497519      PMCID: PMC7865082          DOI: 10.1111/mpp.13027

Source DB:  PubMed          Journal:  Mol Plant Pathol        ISSN: 1364-3703            Impact factor:   5.663


INTRODUCTION

In natural environments as well as in agricultural systems, plants are constantly threatened by pathogens, including nematodes, oomycetes, fungi, bacteria, and viruses. In order to invade, colonize, and obtain nutrients from their host, these pathogens must deal with plant defence mechanisms, and therefore they have evolved different lifestyles and strategies to overcome the plant immune system and cause disease. Given their sessile nature, plants rely in the first line of their defence on anatomical barriers and preformed phytochemical compounds to prevent pathogen invasion. However, when pathogens can overcome these first passive defence layers, recognition of the invaders by the plant can take place within the host tissue. To perceive danger, plants have evolved a wide range of receptors and their function consists of the recognition of so‐called immunogenic patterns (IPs), in both the apoplast and the cytoplasm (Gust et al., 2017; van der Burgh & Joosten, 2019). Triggers of plant immunity include damage‐associated molecular patterns (DAMPs) (Boller & Felix, 2009), structural components of the invaders (microbe‐associated molecular patterns, MAMPs), and pathogen‐secreted compounds known as effectors that serve to perturb plant defence responses and promote virulence (Kanyuka & Rudd, 2019). Apoplastic triggers of immunity are referred to as extracellular immunogenic patterns (ExIPs) (van der Burgh & Joosten, 2019). These are perceived by trans‐membrane pattern recognition receptors (PRRs) located on the cell surface, which have been grouped into receptor‐like kinases (RLKs, containing a cytoplasmic kinase domain) and receptor‐like proteins (RLPs, which lack such a domain) (Gust et al., 2017). By contrast, cytoplasmic danger signals are referred to as intracellular immunogenic patterns (InIPs) and are recognized by the activity of cytoplasmic receptors, usually from the category of nucleotide‐binding leucine‐rich repeat receptors (NB‐LRRs) (van der Burgh & Joosten, 2019). Extra‐ and intracellular receptors possess highly variable LRR domains that confer specific binding capacity to a wide array of ExIPs and InIPs (van der Burgh & Joosten, 2019). On perception of immunogenic patterns, PRRs and NB‐LRRs activate downstream signalling pathways, leading to the induction of immune responses to block the invader. Plant immune responses comprise multiple cellular processes, among others alkalization of the extracellular space (Granado et al., 1995; Jourdan et al., 2009), cell wall reinforcement via callose deposition (Luna Diez et al., 2010; Voigt, 2014), accumulation of reactive oxygen and nitrogen species (Turkan, 2018), an increase in intracellular calcium (Grant et al., 2000), biosynthesis of pathogenesis‐related (PR) proteins (Jain & Khurana, 2018) and phytoalexins (Jeandet, 2015), and the release of small RNAs to interfere with the transcriptional machinery of the invading pathogen (Cai et al., 2018). In an attempt to stop pathogen invasion and disease development, the affected plant cells can commit suicide by undergoing programmed cell death (PCD). Even though it represents a decisive and seemingly altruistic response of an individual plant cell to pathogen invasion, PCD can lead to either resistance or to disease development, depending on the lifestyle of the pathogen. In incompatible interactions with biotrophic pathogens, PCD is referred to as the hypersensitive response (HR) and manifests itself as an effective resistance response of the plant upon pathogen recognition. The HR represents a form of autophagic PCD, preventing the spread of biotrophic microbes (Kabbage et al., 2017; Michaeli et al., 2016). For example, race 5 of the biotrophic fungus Cladosporium fulvum, which causes leaf mould of tomato following a typical gene‐for‐gene interaction, cannot cause disease on tomato carrying the gene encoding the RLP Cf‐4 (Joosten et al., 1994; Thomas et al., 1997). This PRR recognizes the Avr4 effector protein that is secreted by C. fulvum race 5 and mediates the activation of the HR, leading to leaf mould resistance (Thomas et al., 1997). On the other hand, in interactions with necrotrophs, PCD occurs after induction of apoptosis of the affected host cells and becomes visible as a necrotic lesion (Curtis & Wolpert, 2004; van Baarlen et al., 2004). In this situation, PCD shows specific apoptosis hallmarks, which include cell shrinkage and membrane blebbing with retention of plasma membrane integrity, loss of mitochondrial membrane potential and release of cytochrome c into the cytoplasm, protease activation, chromatin condensation and nucleosome cleavage, decrease in ATP, and increase in reactive oxygen species (ROS; Dickman et al., 2017). Necrotrophic pathogens, such as fungi of the genus Botrytis, exploit the plant PCD pathway by actively triggering apoptosis in host cells through the activity of an array of effector compounds, allowing the pathogen to acquire nutrients from the dead plant cells and to colonize the plant tissue (Veloso & van Kan, 2018). To determine the level of resistance or susceptibility of a given plant genotype, assessment of disease symptoms in an infection assay is a key step to evaluate the outcome of an interaction between host and pathogen. Thus, visualization and quantification of PCD over a given period of time are essential to understand how the plant responds upon contact with the pathogen. Furthermore, a visual assessment is often performed to study the effect of individual pathogen IPs (including effectors) on the host, which is done by applying such molecules to a plant tissue by either agroinfiltration or IP injection. This approach allows plant genotypes to be compared for their responses to pathogens or specific IPs and for disease development or the establishment of resistance responses to be monitored over time. At present, the most widely used and commonly accepted method to visualize plant cell death is the trypan blue staining (TBS) technique, which is a cell viability test based on the principle that live cells, possessing intact and functional cell membranes, exclude trypan blue dye whereas dead cells do not. This staining method was first described by Ehrlich (1885) in a study in which trypan blue dye was injected into animal tissue and he observed that the brain tissue showed less staining. Subsequent research showed that when trypan blue was injected directly into the central nervous system instead, the brain tissue would stain equally well, but the stain would not travel to the rest of the body, showing a compartmentalization between the cerebrospinal fluid and the vasculature in the rest of the body (Ehrlich, 1885; Mott, 1913). In the course of the 20th century and until current times, TBS has been successfully implemented in plant biology and phytopathological studies because it was found that TBS allows plant vascular tissue (Escamez et al., 2017), dead Arabidopsis thaliana cells as a consequence of HR induced by Pseudomonas syringae effectors (Johansson et al., 2015), or cell death responses in Brassica rapa cotyledons to oomycete elicitins (Takemoto et al., 2005) to be visualized. TBS consequently became useful to assess the extent of plant tissue colonization by microbial pathogens and, for example, for detecting the presence of (micro‐)lesions in plant–microbe interactions (Vogel et al., 2005). Nevertheless, TBS has four major drawbacks. First, it is not understood how trypan blue dye can pass membranes of dead cells as well as of hyphae of live fungi and oomycetes, whereas the dye is specifically excluded from live plant cells with an intact plasma membrane. Second, the method requires the use of toxic chemicals, such as the dye itself, phenol, and chloral hydrate. Third, TBS represents a destructive assay that does not allow disease progression and/or the occurrence of plant cell death responses to be followed in the same sample over time, nor does it allow the tissue to be sampled for RNA, protein, or metabolite analyses after it is stained. Finally, the method provides a qualitative test of cell viability, but it does not easily allow PCD to be quantified in a standardized way. We here report a nondestructive, nontoxic, and reproducible method to visualize and quantify cell death in green plant tissues. This method can be implemented to study different pathosystems and allows the development of PCD to be followed over time. The protocol is based on the detection of light signals emitted by plant cells undergoing cell death. With the aid of a fluorescence imaging system, we observed red light emission by leaf tissue upon infiltration with known PCD‐inducing compounds, agroinfiltration, or infection by a pathogen. We provide evidence that the red light signal originates from the increase in chlorophyll fluorescence due to dismantling of chloroplast thylakoid membranes in cells undergoing PCD. We provide several illustrations that this method can be used as a versatile tool for quantitative and nondestructive evaluation of plant responses to microbial pathogens or their IPs.

RESULTS

This study was initiated after the observation that on exposure to a green light source, leaf tissue undergoing a typical HR emits light that could be detected using the red light filters of a fluorescence imaging system. Our initial experiments indicated that the intensity of the emitted light positively correlated with the severity of visual signs of cell death. We set out to explore whether this method would be suitable for the development of an alternative, fast, simple, and safe nondestructive tool to visualize cell death in plant tissue. This tool should allow users to quantify PCD in a standardized and reproducible manner and to follow the development of PCD over a given time course. Therefore, we first tested whether the areas yielding a light signal that was detected with red light imaging originated from the affected tissue that in PCD tests is normally stained with trypan blue. Nicotiana benthamiana:Cf‐4 leaves were infiltrated with pure Avr4 protein from C. fulvum and we compared the visible HR symptoms with the classic TBS and the red light signal on the same leaves (Figure 1). At 3 days after infiltration (dai), HR symptoms were observed on leaves infiltrated with 5 µM of Avr4 but not 1 µM (Figure 1, left column). Comparing TBS (Figure 1, middle column) to the red light signal (Figure 1, right column), it can be noticed that the areas stained with trypan blue exactly match with the areas from which the red light signal is detected by the imaging system. The signal intensity was clearly stronger at higher Avr4 concentration and showed a punctate pattern, with the most intense red light emission from areas that displayed the most obvious signs of cell death.
FIGURE 1

Cell death symptoms caused by Avr4 infiltration (a: 1 µM; b: 5 µM) in leaves of Nicotiana benthamiana:Cf‐4 plants, visualized at 3 days after infiltration with visible light (left column), trypan blue staining (middle column), and red light imaging system (right column)

Cell death symptoms caused by Avr4 infiltration (a: 1 µM; b: 5 µM) in leaves of Nicotiana benthamiana:Cf‐4 plants, visualized at 3 days after infiltration with visible light (left column), trypan blue staining (middle column), and red light imaging system (right column) Recently, several studies have used UV transillumination to monitor cell death responses (Guo et al., 2020; Zdrzalek et al., 2020). We compared the use of UV light as a source of excitation for imaging of tissue undergoing PCD to the imaging based on green light excitation described above (Figure S1). When the leaves were excited with UV light, a much stronger background signal was detected in the healthy tissue and in the leaf veins. Such a background signal interferes with the red light signal that derives from tissue undergoing PCD, thereby rendering the UV imaging system less sensitive for areas yielding lower signal intensities. Having shown that the necrotic areas exactly match with the areas from which the red light signal originates, we infiltrated different IPs into leaves of N. benthamiana to test whether the intensity of the red light signal correlates with the severity of the visible PCD symptoms (Figures 2, 3, 4). On infiltration of flg22, at 3 dai no cell death symptoms were observed on the leaves, and also no red light signal was detected (Figure 2a). On the other hand, a clear red light signal was detected on Avr4 infiltration in leaves of N. benthamiana stably expressing Cf‐4, but, apart from some chlorosis, no cell death symptoms were observed (Figure 2b). In contrast, BcNEP1 infiltration caused the swift development of a clearly visible area of dry necrotic tissue that matched with the area yielding a strong red light signal (Figure 2c). In a ROS burst assay, the incubation of leaf disks of N. benthamiana in a solution containing either flg22, Avr4, or BcNEP1 triggered a strong ROS burst (Figure 2, lower panel).
FIGURE 2

Visible cell death symptoms (left) and detected red light signal (right) on infiltration of 0.1 µM flg22 (a), 0.01 µM Avr4 (b), and 0.004 µM BcNEP1 (c) at 3 days after ifiltration in leaves of Nicotiana benthamiana:Cf‐4. Reactive oxygen species burst triggered by 0.1 µM flg22 and the Milli‐Q water (MQ) control in N. benthamiana wild‐type (WT) (d), 0.01 µM Avr4 in N. benthamiana:Cf‐4 and the N. benthamiana WT control plant (e), and 0.004 µM BcNEP1 and the MQ control in N. benthamiana WT (f). n = 8, error bars show the SE

FIGURE 3

Cell death symptoms on agroinfiltration of Avr4/Cf‐4 (a), AvrBS4/BS4 (b), NRC1K191R and NRC1D481V (c) in Nicotiana benthamiana leaves at 3 days after infiltration (dai), visualized under visible light and with the red light imaging system. (d) Red light signal intensity measured upon agroinfiltration of Avr4 after virus‐induced gene silencing of Cf‐4 in N. benthamiana:Cf‐4 at 3 dai, n = 8. Bars represent SD. ***Significant difference using a t test, p < .001

FIGURE 4

Dose‐dependent cell death symptoms caused by BcNEP1 (a, b, c) and Avr4 (d) protein infiltration in Nicotiana benthamiana:Cf‐4 at 3 days after infiltration, as visualized under visible light and with the red light imaging system. Quantification of red light signal intensities emitted by N. benthamiana:Cf‐4 leaves infiltrated with 0.01, 0.1, and 1 µM Avr4 protein and Milli‐Q water (MQ) as a control (e), and with 0.004, 0.04, and 0.4 µM BcNEP1 protein (f). Bars represent averages ± SD (n = 8). Significant difference using analysis of variance followed by Tukey's test, *p < .05, **p < .01, ****p < .0001

Visible cell death symptoms (left) and detected red light signal (right) on infiltration of 0.1 µM flg22 (a), 0.01 µM Avr4 (b), and 0.004 µM BcNEP1 (c) at 3 days after ifiltration in leaves of Nicotiana benthamiana:Cf‐4. Reactive oxygen species burst triggered by 0.1 µM flg22 and the Milli‐Q water (MQ) control in N. benthamiana wild‐type (WT) (d), 0.01 µM Avr4 in N. benthamiana:Cf‐4 and the N. benthamiana WT control plant (e), and 0.004 µM BcNEP1 and the MQ control in N. benthamiana WT (f). n = 8, error bars show the SE Cell death symptoms on agroinfiltration of Avr4/Cf‐4 (a), AvrBS4/BS4 (b), NRC1K191R and NRC1D481V (c) in Nicotiana benthamiana leaves at 3 days after infiltration (dai), visualized under visible light and with the red light imaging system. (d) Red light signal intensity measured upon agroinfiltration of Avr4 after virus‐induced gene silencing of Cf‐4 in N. benthamiana:Cf‐4 at 3 dai, n = 8. Bars represent SD. ***Significant difference using a t test, p < .001 Dose‐dependent cell death symptoms caused by BcNEP1 (a, b, c) and Avr4 (d) protein infiltration in Nicotiana benthamiana:Cf‐4 at 3 days after infiltration, as visualized under visible light and with the red light imaging system. Quantification of red light signal intensities emitted by N. benthamiana:Cf‐4 leaves infiltrated with 0.01, 0.1, and 1 µM Avr4 protein and Milli‐Q water (MQ) as a control (e), and with 0.004, 0.04, and 0.4 µM BcNEP1 protein (f). Bars represent averages ± SD (n = 8). Significant difference using analysis of variance followed by Tukey's test, *p < .05, **p < .01, ****p < .0001 To test whether this method can also be used in combination with Agrobacterium transient transformation assays (ATTA; agroinfiltration), genes encoding IPs and their corresponding immune receptors were transiently expressed in N. benthamiana leaves and the development of HR was monitored. Expression of Avr4 or Cf‐4 alone caused no symptoms of HR and no red light signal was detected outside the area of wounding caused by the syringe during the infiltration. Transient expression of both Avr4 and Cf‐4 in the same leaf area caused symptoms of HR and a clear red light signal was detected (Figure 3a). Similarly, transient expression of the protein Avr4BS4 from Xanthomonas campestris pv. vesicatoria and its matching tomato immune receptor, the NB‐LRR BS4, triggered symptoms of HR and a red light signal only when infiltrated in combination (Figure 3b) (Schornack et al., 2004). Moreover, the elicitor‐independent HR triggered by NRC1D481V was also tested, which is a constitutively active mutant of the tomato NB‐LRR NRC1. In the agroinfiltrated areas, HR symptoms were observed together with a strong emission of a red light signal. In contrast, expression of NRC1K191R, which is an inactive mutant of NRC1, caused no HR and no red light signal was detected (Figure 3c) (Gabriëls et al., 2007). Furthermore, quantification of HR upon agroinfiltration was assessed by transiently expressing Avr4 in N. benthamiana:Cf‐4 after virus‐induced silencing of the gene encoding the immune receptor Cf‐4. A significant decrease in the signal intensity was observed when Cf‐4 was silenced (Figure 3d). To explore whether the imaging system also functions for plants other than N. benthamiana, we infiltrated a different IP, namely the Botrytis cinerea polygalacturonase 3 (BcPG3) protein, in detached leaves of wild‐type (WT) Arabidopsis Col‐0 plants and in an Atsobir1 Arabidopsis mutant (Zhang et al., 2014). As depicted in Figure S2, necrotic tissue in combination with a clear red light signal was observed only in leaves of the WT plants, whereas neither necrotic spots nor a red light signal were detected on BcPG3 infiltration in the Atsobir1 mutant. In addition, a specific red light signal was detected in leaves of transgenic tomato seedlings systemically expressing both the C. fulvum Avr4 and the Cf‐4 resistance protein (Joosten et al., 1994), even before symptoms of necrosis were observed (Figure S3). Besides the qualitative visualization of PCD in different plant species, this methodology should also allow the quantification of the magnitude of PCD by measuring the intensity of the red light signal. To test this application, we infiltrated leaves of N. benthamiana:Cf‐4 with different concentrations of Avr4 protein from C. fulvum and BcNEP1 protein from B. cinerea (Figure 4). While even the highest concentration (1 µM) of Avr4 protein barely induced any sign of PCD as observed with the naked eye (Figure 4d), the imaging system detected increasing red light signal intensities emitted by the infiltrated tissue on infiltration of higher Avr4 concentrations (Figure 4d,e). By contrast, no red light signal was detected in the mock‐infiltrated areas, except at sites where the syringe was applied and pressed on the leaf to infiltrate the Milli‐Q water (MQ), thereby causing some wounding (Figure 4d, MQ). On the other hand, BcNEP1 infiltration caused visible necrotic spots even at the lowest concentrations used (Figure 4a), whereby the red light signal intensity measured on BcNEP1 infiltration was up to 25‐fold higher than in Avr4 infiltrated leaves (Figure 4e,f). In addition, it can be noticed that when BcNEP1 was infiltrated at higher concentrations no clear differences in the necrotic tissues could be observed (Figure 4b,c). However, when quantified with the imaging system, increasing red light signal intensities were recorded in leaves infiltrated with increasing BcNEP1 concentrations (Figure 4f). Furthermore, to monitor the development of PCD over a given time course of infection, leaves of N. benthamiana and Arabidopsis Col‐0 were inoculated with a conidial suspension of B. cinerea. In both cases, development of necrotic lesions was observed from 2 days postinoculation (dpi) and the imaging system detected a clear red light signal at the site of the lesion (Figure 5a,b). For N. benthamiana, pictures of the same detached leaves were taken at 2, 3, and 4 dpi. This was not possible for Arabidopsis leaves because they were too fragile after detachment, and therefore symptom development was followed on detaching different leaves at 2 and 3 dpi. Following incubation, we observed that the lesions expanded and that the red light signal was sometimes detected in irregular shapes, outside of the circular necrotic spot caused by the fungus. The latter light‐emitting zones were not visible to the naked eye (Figure 5a,b, yellow arrows). Moreover, the imaging system detected a red light signal for both Arabidopsis and N. benthamiana, which was visible as concentric rings of different intensities at the location where the necrotic lesions developed (Figure 5a,b, red arrows). In addition to leaves, the red light imaging system was also tested to detect disease symptoms on stem tissue. When B. cinerea conidia were inoculated onto N. benthamiana intact and cut stems, red light emission was observed in stem tissue showing visible disease symptoms (Figure S4).
FIGURE 5

Necrotic lesions developing after Botrytis cinerea inoculation on leaves of (a) Nicotiana benthamiana and (b) Arabidopsis, visualized with visible light (left) and with the red light imaging system (right). Yellow arrows indicate the red light signal detected in tissue outside the necrotic lesions; red arrows indicate distinct concentric rings with different red light signal intensities within the lesion. dpi, days postinoculation

Necrotic lesions developing after Botrytis cinerea inoculation on leaves of (a) Nicotiana benthamiana and (b) Arabidopsis, visualized with visible light (left) and with the red light imaging system (right). Yellow arrows indicate the red light signal detected in tissue outside the necrotic lesions; red arrows indicate distinct concentric rings with different red light signal intensities within the lesion. dpi, days postinoculation We also followed the development of necrotic lesions on two lily cultivars inoculated with two different Botrytis elliptica isolates. Figure 6 shows the necrotic lesions that were observed at 3 and 4 dpi. In all cases the lesions expanded during the incubation. With both isolates, distinct areas were observed emitting a light signal of different intensities, with a fairly abrupt transition between them (Figure 6b). The outer area showed a lower red light intensity (Figure 6b) and this was not visible as necrotic tissue when the picture was taken with visible light. Moreover, for the most aggressive isolate (Be9401) in both cultivars the imaging system detected a much weaker red light signal in the centre of the lesion compared to the concentric rings around the inoculation site and the red light signal was also detected in leaf veins. In contrast, for the less aggressive isolate (Be9732), the signal remained intense at the inoculation site, while the lesions were still expanding and no red light signal was detected in the leaf veins. To illustrate the differences in symptom severity and to compare isolate virulence and cultivar susceptibility, the signal intensities over the lesion areas were quantified and the values were plotted on a bar chart (Figure S5).
FIGURE 6

(a) Necrotic lesions developed after inoculation of Botrytis elliptica isolates Be9401 and Be9732 on leaves of lily genotypes OT and LO, visualized with visible light (upper) and with the red light imaging system (bottom), at 3 and 4 days postinoculation (dpi). For each time point and cultivar, three representative leaves are shown. (b) Magnification of a necrotic lesion observed at 3 dpi on genotype LO, highlighted by a red rectangle in (a), showing an area with different light signal intensities. Yellow dashed lines indicate the borders between these areas

(a) Necrotic lesions developed after inoculation of Botrytis elliptica isolates Be9401 and Be9732 on leaves of lily genotypes OT and LO, visualized with visible light (upper) and with the red light imaging system (bottom), at 3 and 4 days postinoculation (dpi). For each time point and cultivar, three representative leaves are shown. (b) Magnification of a necrotic lesion observed at 3 dpi on genotype LO, highlighted by a red rectangle in (a), showing an area with different light signal intensities. Yellow dashed lines indicate the borders between these areas To explain the origin and nature of the observed red light signal in tissue undergoing PCD, we hypothesized that the signal might derive from the disassembly of thylakoid membranes of the chloroplasts. To test this hypothesis, chloroplasts were isolated from N. benthamiana leaves and treated with different concentrations of the membrane‐disrupting surfactants Nonidet P‐40 and sodium dodecyl sulphate (SDS), and the surfactant solutions were also directly infiltrated in detached N. benthamiana leaves (Figures S6 and S7). As shown in Figure 7a,c, red light emission was detected in the chloroplast suspension after treatment with both surfactants, whereas no red light was emitted by the untreated chloroplast suspension and the pure surfactant solutions. The graphs in Figure 7b,d show that the red light intensities increased when the chloroplasts were treated with a higher concentration of the surfactant, as well as with a longer incubation time with 10% Nonidet P‐40, and with 10% and 20% SDS. On direct leaf infiltration with the surfactant solutions, the tissue collapsed in a dose‐dependent manner and eventually showed a brownish colouration that resembled the necrotic tissue observed upon infiltration of various IPs and in the inoculation assays (Figures [Link], [Link], [Link]).
FIGURE 7

Red light emission in a 24‐well plate containing a chloroplast suspension obtained from Nicotiana benthamiana leaves after treatment with Nonidet P‐40 (NP40) (a, b) or sodium dodecyl sulphate (SDS) (c, d) at different concentrations. Red light signal intensities, determined at three time points after the indicated treatments, were quantified and bars represent averages ± SD (n = 4). Significant difference using analysis of variance followed by Tukey's test, *p < .05, ***p < .001, ****p < .0001

Red light emission in a 24‐well plate containing a chloroplast suspension obtained from Nicotiana benthamiana leaves after treatment with Nonidet P‐40 (NP40) (a, b) or sodium dodecyl sulphate (SDS) (c, d) at different concentrations. Red light signal intensities, determined at three time points after the indicated treatments, were quantified and bars represent averages ± SD (n = 4). Significant difference using analysis of variance followed by Tukey's test, *p < .05, ***p < .001, ****p < .0001

DISCUSSION

Advantages and applications of red light imaging for studies in plant pathology

Various imaging systems have been developed to investigate plant defence responses. These include quantification methods for ion flux changes and ROS accumulation (Grant & Loake, 2000), nitric oxide generation (Delledonne et al., 1998), and biophoton generation (Bennett et al., 2005). Despite focusing on plant defence responses, none of these methods allows PCD to be investigated as a determinant of the outcome of the interaction in a given pathosystem, which usually requires a destructive staining method, TBS. In this study we describe a simple methodology to visualize and quantify PCD in green plant tissues using the red light filter of a fluorescence–luminescence imaging system, which was originally purchased for imaging DNA or protein gels, and western blots. We show that the system can also be applied for imaging plant cell death during interactions with pathogens or their IPs. Compared to classical TBS, the red light imaging method has several advantages. First of all, it is safer because it does not require the use of toxic chemicals for preparing the staining and destaining solutions, and it generates no chemical waste. The method is fast and simple because it is solely based on the detection of a light signal that is automatically processed by the imaging system. This allows the analysis of a large number of samples in a short time, which makes it useful for large‐scale screenings, with up to hundreds of leaves per day. Importantly, in contrast to TBS, this method is nondestructive and the same samples can be analysed at different time points to monitor the progression of PCD or the imaging may assist in selecting different regions in the sample for extraction of RNA, proteins, or metabolites. Both in the case of leaf infiltration with different IPs and in inoculation assays with pathogens, the method allows the plant response to be observed at different time points, thereby enabling the development of symptoms to be followed during the interaction over a given time span (Figures 5 and 6). The red light imaging system allows PCD to be visualized and quantified in a standardized manner. Different parameters, such as the exposure time and the focus on particular sites of the sample, can be set in advance and with the aid of software the red light intensities can be quantified for comparison among different samples. On the other hand, for lesion size measurements, it is possible to quantify the area of the tissue for which the red light signal exceeds background levels. In contrast, in TBS the staining intensity may depend on how long the tissue is treated with the staining and destaining solutions, thereby only permitting a qualitative visualization of PCD, but hampering a quantitative comparison across multiple samples. It is important to note that the imaging system requires placing the samples on a flat surface, therefore it is not possible to acquire images from three‐dimensional objects such as whole plants because in that case the device is unable to bring the complete area to be studied into focus. More specialized and expensive equipment would be needed for such purposes. This simple and nondestructive method allows new research aimed at studying the role of PCD in plant–pathogen interactions to be pursued. PCD symptoms can be detected before becoming visible to the naked eye and different stages of PCD might be able to be distinguished depending on light intensities. In necrotic lesions caused by Botrytis spp., distinct concentric areas were observed showing a red light signal of different intensity. The absence of a red light signal at the centre of the lesion (Figure 6b) indicates that at that location the light signal emitting compounds have been completely degraded. In some cases, when moving distally from the centre of the lesion, the first area surrounding the inoculation site showed a higher signal intensity as compared to the lesion periphery, resulting in the formation of a concentric ring (Figure 6b). This observation might indicate that PCD is more advanced in the brighter region and is only just induced in the peripheral area. It remains to be established whether higher signal intensities correlate with more advanced stages of PCD. Moreover, by inoculating plant tissue with transgenic pathogen strains expressing green fluorescent protein (GFP), it should be possible to investigate whether host tissue colonization has already occurred in the areas displaying low‐intensity red light signals, or whether PCD has been induced by effector compounds that are secreted by the pathogen and have moved into the surrounding host tissue. We tested the versatility of the methodology for investigating PCD development by infiltration of different IPs. Upon flg22 infiltration, we observed neither a red light signal nor development of a necrotic lesion, despite the fact that flg22 recognition took place in the leaf tissue and triggered a ROS burst as previously reported (Figure 2) (Gómez‐Gómez et al., 1999). This observation confirms that PCD activation requires specific signalling mechanisms to be induced and that not all IPs equally trigger PCD (van der Burgh & Joosten, 2019). There was a dose‐ and time‐dependent cell death response on infiltration of various IPs at different concentrations (Figures 1 and 4). Because the symptoms are not always clearly visible, the red light filter of the imaging system provides a more reliable and objective tool to quantitatively determine the contribution of given IPs to PCD induction in the host tissue. Finally, this method will help to investigate the function of pathogen‐secreted effectors. For example, by labelling pathogen‐secreted effectors with fluorescent tags, it will be possible to study their mobility in the tissue and localization over time during host infection, in conjunction with the host response to these effectors. Additionally, by identifying areas that are in an early stage of PCD, showing low‐intensity red light signals, it should be possible to specifically dissect and process tissue areas to isolate proteins or RNAs at the appropriate time and location. Possibly, in this way also specific effector targets in the host cells and PCD signalling components can be identified via coimmunoprecipitation assays, as the appropriate tissue can be collected using the imaging system.

PCD, chlorophyll fluorescence, and the red light signal

This study shows that tissue undergoing PCD emits red light that can be detected by means of the imaging system. We hypothesized that the red light signal derives from increased chlorophyll fluorescence, as a consequence of the disassembly of the thylakoid membranes in the chloroplasts. This hypothesis was tested by treating isolated chloroplasts with membrane‐disrupting surfactants (Figure 7) and by leaf infiltration with such surfactants (Figures S6 and S7). The membrane‐disrupting properties of the surfactants caused a concentration‐dependent emission of the red light signal from isolated chloroplasts and infiltrated leaves. From a photochemical point of view, this phenomenon can be explained by the abortion of photochemical activity due to the dismantling of the chloroplast thylakoid membranes in cells undergoing PCD, resulting in an increase in chlorophyll fluorescence. In photochemically active plant cells, chlorophyll molecules are tightly packed together in light harvesting complexes (LHCs) that are located at the thylakoid membranes of the chloroplasts. Chlorophyll fluorescence is quenched by photochemical activity when LHC‐photosystem II (PSII) complexes are assembled in thylakoids. Chlorophylls have two distinct absorption wavelength peaks in the visible light spectrum: the highest is found in the blue spectrum at c.450 nm and the lowest at 650 nm in the red spectrum. Moreover, chlorophylls possess a fluorescence emission peak in the red spectrum at c.680 nm. During PCD or on perturbation of thylakoid membrane stability, LHC‐PSII supercomplexes are disassembled from chloroplast thylakoid membranes. This causes cessation of the electron flow required for photochemical activity because the connection between PSII and photosystem I via plastoquinone, cytochrome‐b6f, and plastocyanin is lost. As consequence, light‐excited chlorophyll molecules will dissipate photon energy through heat and fluorescence (Baker, 2008).

Conclusion

In this study we show that red light fluorescence detection by a multipurpose imaging system can be used to visualize plant tissue undergoing cell death. Quantifying the intensity of the red light signal enables the magnitude of cell death to be evaluated and thereby provides a readout of the plant immune response in a faster, safer, and nondestructive manner when compared to previously developed chemical staining methodologies. This application can be implemented to screen for differences in symptom severity in plant–pathogen interactions, and to visualize and quantify in a more sensitive and objective manner the intensity of the plant response on perception of a given immunogenic pattern.

EXPERIMENTAL PROCEDURES

Plant material and growth conditions

Transgenic N. benthamiana plants stably expressing Cf‐4 (N. benthamiana:Cf‐4) (Gabriëls et al., 2006) and WT N. benthamiana plants were grown under 16 hr light at 24 °C and 8 hr darkness at 22 °C and RH of 75%. Arabidopsis ecotype Col‐0 and the Atsobir1 mutant line in the Col‐0 background (Gao et al., 2009), were grown under 12 hr light at 21 °C and 12 hr darkness at 19 °C at RH of 75%. Lilium genotypes from the hybrid types OT (Oriental × Trumpet) and LO (Longiflorum × Oriental) were planted as bulbs and grown in the greenhouse under natural day/night light regime and temperatures. Mature leaves of 2‐month‐old plants were used for the inoculation assays.

Red light imaging and signal quantification

Treated plant material was placed in a ChemiDoc MP imaging system, model Universal Hood III (Bio‐Rad), and images were acquired on excitation by a light source in the green visible spectrum (Green LED Module kit no. 1708284) or in the UV spectrum (UV‐B Transillumination no. 1001361), with filters capturing the light emitted in the red visible spectrum (filter 605/50). The exposure time was adjusted to avoid image saturation. For comparison purposes, the images shown were autoscaled according to the treatment that yielded the lowest signal intensity. For the IP infiltration assays, the mean signal intensity was calculated using Image Lab v. 6.0.1 software by manually selecting the treated areas and subtracting the background signal intensity value. For the disease assay in Lilium, the spots yielding a light signal were manually highlighted and the signal was calculated as the intensity over the selected area among 12 inoculations sites on three leaves per plant cultivar and per B. elliptica isolate at 2, 3, and 4 dpi.

Statistical analysis

GraphPad Prism v. 8.4.0 was used to perform an analysis of variance (ANOVA) with a post hoc Tukey's test on the IP infiltration experiments, chloroplast treatment with surfactants, and the Botrytis disease assays in Lilium.

Immunogenic pattern‐triggered cell death assays

NEP1 protein from B. cinerea (BcNEP1) was produced in Pichia pastoris (Schouten et al., 2008) and dissolved in 10 mM phosphate buffer, pH 5.5. The peptide flg22 (EZBiolab) and the protein Avr4 (produced in P. pastoris; van den Burg et al., 2001) from C. fulvum were dissolved in water. Protein infiltration was carried out in the first fully expanded leaf of 5‐ to 6‐week‐old N. benthamiana:Cf‐4 and WT plants, and the cell death response was scored at 3 dai. The infiltrated areas were photographed and the average red light signal intensities were calculated. B. cinerea polygalacturonase 3 (BcPG3) protein produced in P. pastoris (Zhang et al., 2014) was dissolved in 10 mM sodium acetate buffer pH 4.2 and infiltrated at a concentration of 1.5 µM in fully expanded rosette leaves of 6‐ to 7‐week‐old Arabidopsis plants. The necrotic areas were photographed and analysed by red light imaging at 8 dai. Tomato (Solanum lycopersicum) seedlings stably expressing both Avr4 and Cf‐4 were generated by crossing transgenic tomato cv. Moneymaker (MM) expressing C. fulvum effector Avr4 to MM‐Cf‐4 plants (Etalo et al., 2013; Stulemijer et al., 2007). After sowing, plants were grown for 4 weeks at 33 °C and 100% RH to suppress systemic PCD. Synchronized systemic PCD was eventually induced by transferring the seedlings to 20 °C and 70% RH. The first fully developed tomato leaves were harvested 4 hr after transfer and were photographed.

Agrobacterium‐mediated transient transformation and virus‐induced gene silencing in N. benthamiana

Agrobacterium‐mediated transient transformation (agroinfiltrations) were performed as previously described (van der Hoorn et al., 2000). The constructs pMOG800‐Cf‐4, pGWB20‐BS4, pVsF‐AvrBS4, pBINKS‐NRC1D481V, and pBINKS‐NRC1K191R‐GFP were infiltrated with Agrobacterium tumefaciens cultures at OD600 = 1. The construct pMOG800‐Avr4 was infiltrated with A. tumefaciens cultures at OD600 = 0.03. Symptoms were evaluated at 3 dai. Virus‐induced gene silencing (VIGS) using tobacco rattle virus (TRV)‐based vectors was performed in N. benthamiana:Cf‐4 as previously described (Liebrand et al., 2012).

Trypan blue staining

N. benthamiana leaves infiltrated with BcNEP1 protein were harvested at 3 dai and boiled for 2–5 min in staining solution, consisting of 16.6% (wt/vol) of trypan blue in an ethanol/lactic acid/phenol/water (3:1:1:1) mixture. After this, the leaves were destained with chloral hydrate and photographed.

ROS burst measurements

For each treatment, we used eight leaf discs (5 mm) punched from the first fully expanded leaf of 5‐ to 6‐week‐old N. benthamiana:Cf‐4 plants. The discs were floated on 50 μl of Milli‐Q water (MQ) in a 96‐well plate and kept in the dark at room temperature for 5 hr. The MQ was then removed using tissue paper, replaced by 50 μl of MQ, and kept in the dark for 1 hr at room temperature. After this, 50 μl of solution containing the IPs, 50 μM luminol L‐012 (FUJIFILM), and 10 μg horseradish peroxidase were added to each well. The luminescence values were measured using a CLARIOstar plate reader (BMG LABTECH) over a period of 4 hr.

Disease assays with Botrytis spp.

Leaves, intact and cut stems of N. benthamiana plants, and leaves of WT Arabidopsis Col‐0 were inoculated with conidia of B. cinerea B05.10. The fungus was grown on malt extract agar (50 g/L; Oxoid) and sporulation was induced by illumination with UV‐A lamps. After harvesting, the concentration of the conidia was adjusted to 106/ml using potato dextrose broth (PDB; Oxoid) at 24 g/L PDB, and inoculated in 2‐μl droplets on the adaxial side of the leaves. Mock inoculation was carried out by applying 2‐μl droplets of PDB solution. Inoculated plants were then incubated in moist plastic boxes and were kept under the same light conditions under which the plants were grown before. Inoculated leaves were detached at 2 dpi, photographed, and analysed by red light imaging. The leaves were then stuck into wet floral foam and returned to the plastic boxes to allow disease progression to take place. This procedure was repeated at 3 and 4 dpi to visualize symptom development. Given the fragility of inoculated A. thaliana leaves, it was only possible to record symptom development for N. benthamiana leaves. Symptom development on inoculated N. benthamiana stems was evaluated at 4 dpi. Detached leaves from two different cultivars of Lilium spp. (OT and LO) were abaxially inoculated with 2‐μl droplets of B. elliptica conidia from isolates Be9401 and Be9732, respectively (105 conidia/ml in 12g/L PDB). Mock inoculation was carried out with 2‐μl droplets of PDB. Inoculated leaves were incubated in moist plastic boxes and pictures of red light imaging were taken at 3 and 4 dpi, with 2 s exposure. The signal intensities were analysed in a manually drawn area, encompassing the necrotic lesion for all the cultivar × pathogen combinations. An area of variable size was selected in each control leaf for subtraction of the background signal of healthy tissue.

Thylakoid disassembly test

Intact chloroplasts were isolated from N. benthamiana by stirring cut leaf pieces and following the protocol as described by Ho and Theg (2016). After ultracentrifugation, the Percoll layer containing the chloroplasts was transferred to a 24‐well plate with 150 µl aliquots per well. Two surfactants with membranolytic activity were added to each well to perturb thylakoid membrane stability. The surfactants were used at final concentrations of 0.01%, 0.1%, 1%, and 10% for Nonidet P‐40 (wt/vol = g/100 ml) and 0.1%, 1%, 10%, and 20% for SDS (wt/vol = g/100 ml). One hundred and fifty microlitres of the respective surfactant at the different concentrations was added to the wells containing the chloroplast suspension. Untreated chloroplast suspension and pure surfactants at the highest concentration were used as controls. After gently mixing and 20 min of incubation in the dark, the microtitre plates were placed in the imaging system to detect the light signal using red light imaging, at an exposure time of 10 s. FIGURE S1 Cell death symptoms caused by Avr4 infiltration for (a) 1 µM and (b) 5 µM) in leaves of N. benthamiana:Cf‐4 plants, visualized at 3 days postinocuation with visible light (left column), trypan blue staining (middle‐left column), and red light imaging using green (middle‐right column) and UV light (right column) as excitation wave lengths Click here for additional data file. FIGURE S2 Cell death symptoms monitored at 8 days postinoculation caused by BcPG3 protein infiltration in an Arabidopsis wild‐type Col‐0 plant (left panels) and in a Atsobir1 deletion mutant (right panels). Each leaf was visualized under visible light with visible light (left) and red light imaging (right) Click here for additional data file. FIGURE S3 Red light signal detected in leaves of Avr4/Cf‐4 tomato seedlings (a) and in control leaves of only Cf‐4‐expressing tomato seedlings (b) harvested 4 hr after changing growth conditions to a permissive temperature (see Experimental Procedures) Click here for additional data file. FIGURE S4 Necrotic lesions developing on Botrytis cinerea inoculation on cut stems (a) and intact stems (b) of Nicotiana benthamiana at 4 days postinoculation, visualized with visible light (left) and red light imaging (right). Red arrows indicate the inoculation sites on the stem Click here for additional data file. FIGURE S5 Red light signal intensity over the lesion areas measured at 2, 3, and 4 days postinoculation from necrotic lesions developed on lily OT and LO inoculated with Be9401 and Be9732. Bars represent averages ± SD (n = 12). Asterisks indicate significant differences using analysis of variance followed by Tukey’s test Click here for additional data file. FIGURE S6 (a) Leaf response of Nicotiana benthamiana infiltrated with NP40 at the indicated concentrations, visualized under visible light (upper panel) and with red light imaging (lower panel) at 20 min after infiltration of the surfactant. (b) Red light signal intensity values measured at 20 min after infiltration with NP40 at the indicated concentrations. Bars represent averages ± SD (n = 3). Asterisks indicate significant differences using analysis of variance followed by a Tukey test Click here for additional data file. FIGURE S7 (a) Leaf response of Nicotiana benthamiana infiltrated with sodium dodecyl sulphate (SDS) at the indicated concentrations, visualized under visible light (upper panel) and with red light imaging (lower panel) at 20 min after infiltration of the surfactant. (b) Red light signal intensity values measured at 20 min after infiltration with SDS at the indicated concentrations. Bars represent averages ± SD (n = 4). Asterisks indicate significant differences using analysis of variance followed by a Tukey test Click here for additional data file.
  44 in total

1.  The RPM1 plant disease resistance gene facilitates a rapid and sustained increase in cytosolic calcium that is necessary for the oxidative burst and hypersensitive cell death.

Authors:  M Grant; I Brown; S Adams; M Knight; A Ainslie; J Mansfield
Journal:  Plant J       Date:  2000-08       Impact factor: 6.417

2.  Characterization of the tomato Cf-4 gene for resistance to Cladosporium fulvum identifies sequences that determine recognitional specificity in Cf-4 and Cf-9.

Authors:  C M Thomas; D A Jones; M Parniske; K Harrison; P J Balint-Kurti; K Hatzixanthis; J D Jones
Journal:  Plant Cell       Date:  1997-12       Impact factor: 11.277

Review 3.  Plant Immunity: Thinking Outside and Inside the Box.

Authors:  Aranka M van der Burgh; Matthieu H A J Joosten
Journal:  Trends Plant Sci       Date:  2019-06-03       Impact factor: 18.313

Review 4.  The Life and Death of a Plant Cell.

Authors:  Mehdi Kabbage; Ryan Kessens; Lyric C Bartholomay; Brett Williams
Journal:  Annu Rev Plant Biol       Date:  2017-01-11       Impact factor: 26.379

5.  Fungal endopolygalacturonases are recognized as microbe-associated molecular patterns by the arabidopsis receptor-like protein RESPONSIVENESS TO BOTRYTIS POLYGALACTURONASES1.

Authors:  Lisha Zhang; Ilona Kars; Bert Essenstam; Thomas W H Liebrand; Lia Wagemakers; Joyce Elberse; Panagiota Tagkalaki; Devlin Tjoitang; Guido van den Ackerveken; Jan A L van Kan
Journal:  Plant Physiol       Date:  2013-11-20       Impact factor: 8.340

6.  Differences in cell death induction by Phytophthora Elicitins are determined by signal components downstream of MAP kinase kinase in different species of Nicotiana and cultivars of Brassica rapa and Raphanus sativus.

Authors:  Daigo Takemoto; Adrienne R Hardham; David A Jones
Journal:  Plant Physiol       Date:  2005-06-24       Impact factor: 8.340

Review 7.  Many Shades of Grey in Botrytis-Host Plant Interactions.

Authors:  Javier Veloso; Jan A L van Kan
Journal:  Trends Plant Sci       Date:  2018-04-30       Impact factor: 18.313

8.  Plants send small RNAs in extracellular vesicles to fungal pathogen to silence virulence genes.

Authors:  Qiang Cai; Lulu Qiao; Ming Wang; Baoye He; Feng-Mao Lin; Jared Palmquist; Sienna-Da Huang; Hailing Jin
Journal:  Science       Date:  2018-05-17       Impact factor: 47.728

9.  Functional analysis of RXLR effectors from the New Zealand kauri dieback pathogen Phytophthora agathidicida.

Authors:  Yanan Guo; Pierre-Yves Dupont; Carl H Mesarich; Bo Yang; Rebecca L McDougal; Preeti Panda; Paul Dijkwel; David J Studholme; Christine Sambles; Joe Win; Yuanchao Wang; Nari M Williams; Rosie E Bradshaw
Journal:  Mol Plant Pathol       Date:  2020-07-08       Impact factor: 5.663

10.  The rice NLR pair Pikp-1/Pikp-2 initiates cell death through receptor cooperation rather than negative regulation.

Authors:  Rafał Zdrzałek; Sophien Kamoun; Ryohei Terauchi; Hiromasa Saitoh; Mark J Banfield
Journal:  PLoS One       Date:  2020-09-15       Impact factor: 3.240

View more
  8 in total

1.  Metabolic pathway genes for editing to enhance multiple disease resistance in plants.

Authors:  Ajjamada C Kushalappa; Niranjan G Hegde; Kalenahalli N Yogendra
Journal:  J Plant Res       Date:  2022-08-29       Impact factor: 3.000

2.  The Intracellularly Acting Effector Foa3 Suppresses Defense Responses When Infiltrated Into the Apoplast.

Authors:  Nico Tintor; Gea A M Nieuweboer; Ilse A W Bakker; Frank L W Takken
Journal:  Front Plant Sci       Date:  2022-05-23       Impact factor: 6.627

3.  Red light imaging for programmed cell death visualization and quantification in plant-pathogen interactions.

Authors:  Sergio Landeo Villanueva; Michele C Malvestiti; Wim van Ieperen; Matthieu H A J Joosten; Jan A L van Kan
Journal:  Mol Plant Pathol       Date:  2021-01-26       Impact factor: 5.663

Review 4.  Sensor-based phenotyping of above-ground plant-pathogen interactions.

Authors:  Florian Tanner; Sebastian Tonn; Jos de Wit; Guido Van den Ackerveken; Bettina Berger; Darren Plett
Journal:  Plant Methods       Date:  2022-03-21       Impact factor: 5.827

5.  Phytophthora capsici sterol reductase PcDHCR7 has a role in mycelium development and pathogenicity.

Authors:  Weizhen Wang; Fan Zhang; Sicong Zhang; Zhaolin Xue; Linfang Xie; Francine Govers; Xili Liu
Journal:  Open Biol       Date:  2022-04-06       Impact factor: 6.411

6.  Cytotoxic activity of Nep1-like proteins on monocots.

Authors:  Maikel B F Steentjes; Andrea L Herrera Valderrama; Laetitia Fouillen; Delphine Bahammou; Thomas Leisen; Isabell Albert; Thorsten Nürnberger; Matthias Hahn; Sébastien Mongrand; Olga E Scholten; Jan A L van Kan
Journal:  New Phytol       Date:  2022-05-03       Impact factor: 10.323

7.  Histological Methods to Detect Early-stage Plant Defense Responses during Artificial Inoculation of Lolium perenne with Epichloë festucae.

Authors:  Mostafa Rahnama; Damien J Fleetwood; Richard D Johnson
Journal:  Bio Protoc       Date:  2021-05-05

8.  A novel robust and high-throughput method to measure cell death in Nicotiana benthamiana leaves by fluorescence imaging.

Authors:  Yuxuan Xi; Vincent Chochois; Thomas Kroj; Stella Cesari
Journal:  Mol Plant Pathol       Date:  2021-08-23       Impact factor: 5.663

  8 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.