Literature DB >> 32265949

Photosynthetic Cost Associated With Induced Defense to Plasmopara viticola in Grapevine.

Antonio F Nogueira Júnior1, Merle Tränkner2, Rafael V Ribeiro3, Andreas von Tiedemann4, Lilian Amorim1.   

Abstract

Downy mildew caused by Plasmopara viticola is one of the most destructive diseases of Vitis vinifera worldwide. Grapevine breeding programs have introgressed P. viticola-resistant traits into cultivated V. vinifera genotypes and launched interspecific hybrids with resistance against downy mildew. In general, pathogen infection affects primary metabolism, reduces plant growth and development and modifies the secondary metabolism toward defense responses, which are costly in terms of carbon production and utilization. The objective of this work was to evaluate the photosynthesis impairment by inducible defenses at the leaf level in V. vinifera cultivars resistant to P. viticola. Photosynthetic limitations imposed by P. viticola in susceptible and resistant grapevine cultivars were evaluated. Histochemical localization of hydrogen peroxide and superoxide and the activity of ascorbate peroxidase were assessed. Measurements of leaf gas exchange, chlorophyll fluorescence and the response of leaf CO2 assimilation to increasing air CO2 concentrations were taken, and photosynthetic limitations determined in cultivars Solaris (resistant) and Riesling (susceptible). The net photosynthetic rates were reduced (-25%) in inoculated Solaris plants even before the appearance of cell death-like hypersensitive reactions ("HR"). One day after "HR" visualization, the net photosynthetic rate of Solaris was reduced by 57% compared with healthy plants. A similar pattern was noticed in resistant Cabernet Blanc and Phoenix plants. While the susceptible cultivars did not show any variation in leaf gas exchange before the appearance of visual symptoms, drastic reductions in net photosynthetic rate and stomatal conductance were found in diseased plants 12 days after inoculation. Decreases in the maximum Rubisco carboxylation rate and photochemical impairment were noticed in Riesling after inoculation with P. viticola, which were not found in Solaris. Damage to the photochemical reactions of photosynthesis was likely associated with the oxidative burst found in resistant cultivars within the first 24 h after inoculation. Both chlorophyll degradation and stomatal closure were also noticed in the incompatible interaction. Taken together, our data clearly revealed that the defense response against P. viticola causes a photosynthetic cost to grapevines, which is not reversible even 12 days after the pathogen infection.
Copyright © 2020 Nogueira Júnior, Tränkner, Ribeiro, von Tiedemann and Amorim.

Entities:  

Keywords:  Vitis vinifera; biotic stress; downy mildew; gas exchange; photosynthesis limitations

Year:  2020        PMID: 32265949      PMCID: PMC7098430          DOI: 10.3389/fpls.2020.00235

Source DB:  PubMed          Journal:  Front Plant Sci        ISSN: 1664-462X            Impact factor:   5.753


Introduction

Downy mildew, caused by the biotrophic oomycete Plasmopara viticola, is one of the most devastating diseases of grapevine. P. viticola is native to North America and was considered to be of minor importance on Vitis aestivalis, in which it was first described in 1834 (Wilcox et al., 2015). However, this pathogen was inadvertently introduced into Europe in 1878 and now it is endemic in the main grapevine growing regions worldwide (Gessler et al., 2011). The European grapevine (Vitis vinifera) evolved in the absence of this disease and it is extremely susceptible to P. viticola. After P. viticola introduction the pathogen quickly spread throughout the European vineyards (Rossi et al., 2013). P. viticola affects grapevine leaves, inflorescences, berries and shoots with high disease severity causing defoliation and producing low-quality, unsightly or entirely damaged grapes. The pathogen can also cause weakening, dwarfing and death of young shoots (Agrios, 2005). When the weather is conducive to the disease and in the absence of effective disease control, downy mildew can easily reduce the grapevine yield up to 100% (Agrios, 2005; Wilcox et al., 2015). P. viticola zoospores penetrate the host through stomata and develop intercellular mycelium in the mesophyll of grapevine leaves, generating globose haustoria responsible for the oomycete nutrition. After infection and colonization periods, sporangiophores and sporangia emerge from the stomata (Ash, 2000; Agrios, 2005; Gessler et al., 2011; Kamoun et al., 2015; Armijo et al., 2016). Currently, downy mildew is of paramount importance in all humid parts of the world where grapevine is grown. In subtropical regions of South America, the occurrence of rainy periods at the time of budbreak and flowering can lead to severe epidemics of downy mildew. Consequently, high frequency of fungicide sprays is necessary to control downy mildew during the rainy season, i.e., from summer to autumn (Cappello et al., 2017). Although European V. vinifera cultivars are highly susceptible to downy mildew, Muscadinia species and several American and Asian Vitis species exhibit varying levels of resistance to P. viticola. By the conventional breeding, V. vinifera has been crossed with grapevines showing resistance against P. viticola, and resistant interspecific hybrids have been subsequently found. In resistant accessions obtained by crosses with North American Vitis species, P. viticola completes its life cycle and releases less sporangia than on susceptible cultivars (Bellin et al., 2009). The genomic regions conferring resistance to downy mildew referred to as resistance to P. viticola (Rpv) are quantitatively inherited. To date, 27 quantitative trait loci (QTLs) with major effects on downy mildew resistance are known and described (Fischer et al., 2004; Welter et al., 2007; Bellin et al., 2009; Marguerit et al., 2009; Vezzulli et al., 2019). Rpv1 and Rpv2 are responsible for the resistance derived from Muscadinia rotundifolia (Peressotti et al., 2010) and map to chromosomes 12 and 18, respectively (Merdinoglu et al., 2003; Welter et al., 2007). The locus Rpv3, which was first described in cultivars Regent and Bianca, came from the American wild Vitis species and is located on chromosome 18 (Fischer et al., 2004; Welter et al., 2007; Bellin et al., 2009). The Asian resistance loci Rpv8, Rpv10, and Rpv12 originated from Vitis amurensis. Rpv8 and Rpv12 map to chromosome 14 (Venuti et al., 2013) and Rpv10 maps to chromosome 9 and is present in cultivar Solaris (Blasi et al., 2011; Schwander et al., 2012; Venuti et al., 2013). The table of loci traits in grapevine relevant for breeding and genetics with a complete description of associated markers, their chromosomal localization, and the donor genotype/species of P. viticola resistance is available online (www.vivc.de/loci). The mechanisms by which grapevine cells reduce downy mildew intensity are complex and not fully elucidated (Kortekamp, 2006; Liu et al., 2015). These mechanisms comprise pre-existing chemical defenses and inducible structural and biochemical defenses. High constitutive levels of antimicrobial compounds, such as inositol and caffeic acid, are observed in healthy leaves of resistant cultivars (Figueiredo et al., 2008). The induced structural defenses include formation of callose (Kortekamp et al., 1997; Gindro et al., 2003) and lignification (Dai et al., 1995). The accumulation of stilbenes and pathogenesis-related proteins (Kortekamp, 2006; Richter et al., 2006; Slaughter et al., 2008), the generation of reactive oxygen species, the induction of peroxidases (Kortekamp and Zyprian, 2003), and finally, cell death-like hypersensitive reactions (“HR”) represent induced biochemical defenses found in grapevine resistant cultivars after inoculation with P. viticola (Díez-Navajas et al., 2008). Pathogen infection in incompatible interactions may affect the plant primary and also the secondary metabolism due to the induction of defenses (Berger et al., 2007). In the incompatible interaction of Phytophthora nicotianae and Nicotiana tabacum, photosynthesis is switched off few hours after pathogen inoculation, and defense mechanisms are initiated. This metabolic shift is related to an early blockage of intercellular sugar transport by callose deposition and increases in apoplastic invertase activity (Scharte et al., 2005). Both constitutive and inducible defenses are costly in terms of carbohydrate production and utilization. However, inducible systems seem to be less costly than constitutive systems in the absence of pathogens because inducible defenses require allocation of both energy and resources away from growth and reproduction (Cipollini, 1998; DeWitt et al., 1998; Purrington, 2000; Heil and Baldwin, 2002). The cost of inducible defenses is a widely discussed topic among breeders and pathologists and there are few papers conveying a significant piece of novel information on the drawbacks of plant disease resistance. The development of new resistant varieties to P. viticola is a promising way to control downy mildew and the better understanding of the effect of P. viticola in the leaf carbon assimilation of resistant plants is an important issue for the scientific community interested in grapevine-pathogen interaction. Plasmopara viticola affects the photosynthesis of susceptible grapevine cultivars by reducing the diffusion of CO2 in leaf mesophyll (Jermini et al., 2010; Nogueira Júnior et al., 2019). In diseased plants, the concentration of photosynthetic pigments in areas nearby lesions and the primary photochemistry are reduced (Moriondo et al., 2005). Reductions and abundance in Rubisco activity have also been reported in susceptible cultivars of V. vinifera and Vitis labrusca infected by P. viticola (Gamm et al., 2011; Nogueira Júnior et al., 2019). However, the effect of downy mildew on photosynthetic capacity and leaf gas-exchange for resistant cultivars has not been characterized. Here, alterations in gas exchange, chlorophyll fluorescence and production of oxygen reactive species of leaves from cultivars with different resistance levels to downy mildew were examined. The objective of this work was to estimate the impairment of photosynthesis caused by inducible responses to the infection of P. viticola on downy mildew resistant cultivars and to elucidate the underlying limitations imposed by P. viticola to the photosynthesis of resistant and susceptible cultivars. The oxidative burst induced by P. viticola was also evaluated and related to the photosynthetic limitations in resistant cultivars.

Materials and Methods

The characterization of photosynthetic cost associated with inducible defenses to P. viticola in grapevine was performed in three steps. The first step, which is described in section “Gas Exchange Variables in Cultivars Resistant and Susceptible to Downy Mildew,” aimed to evaluate gas exchange in grapevine leaves of different cultivars, resistant and susceptible to downy mildew. The resistant cultivars Solaris, Cabernet Blanc and Phoenix and the susceptible cultivars Riesling, Niagara Rosada, Merlot and Moscato were used in this step. For the second step two representative cultivars, one resistant (Solaris) and one susceptible (Riesling) to P. viticola were selected. This step aimed to evaluate the diffusive, photochemical and biochemical limitations imposed by P. viticola on grapevine photosynthesis. The experiment at the last step aimed to characterize some induced responses in the cultivar Solaris and to correlate these responses with the early decrease on photosynthesis of leaves challenged with P. viticola.

Plant Material and Inoculum Maintenance

All experiments were performed with potted grapevine plants using the downy mildew resistant cultivars of V. vinifera Solaris, Cabernet Blanc and Phoenix; the susceptible cultivars Riesling, Merlot, Moscato; and the susceptible cultivar Niagara Rosada (V. labrusca). The cuttings were grown in pots containing sterilized substrate (manure, clay soil and sand at a ratio of 3:3:1) under greenhouse conditions where the average air temperature was 25 ± 5°C (error SEM), the RH was 60% and the photosynthetic active radiation (PAR) was 500 μmol m–2 s–1 with a 12-h photoperiod. The substrate was fertilized with 3 g of Osmocote® per gram of soil. After budbreak, the plants were conducted in a single stem and top pruned when they presented 6–7 fully expanded leaves. The plants received 200 mL of water daily and were fertilized monthly with 10 g of Hakaphos® Blau. Leaves with symptoms of downy mildew were collected in vineyards, and sporangia of P. viticola were harvested, dried at room temperature and maintained at −25°C in Eppendorf tubes. Sporangia were rehydrated in 50 mL of bidistilled water and suspensions of the inoculum at a concentration of 104 sporangia mL–1 were obtained using a Neubauer chamber. One month after budbreak, potted plants of cv. Merlot were inoculated by spraying inoculum suspension in the abaxial side of leaves up to runoff. Immediately after the inoculation, plants were kept in an incubation chamber at 25°C and covered with a black plastic bag to keep leaves in darkness and high air relative humidity (>90%) for 24 h. Three drops of 50 μL each of the sporangia suspension used for inoculation were placed in a polystyrene dish and kept in the same incubation chamber of inoculated plants to evaluate the viability of sporangia. Viability of sporangia was determined by the observation at light microscope of the opened sporangium 24 h after incubation. Afterward, the plastic bags were removed and plants were kept in the same incubation chamber with a 12-h photoperiod and PAR of 400 μmol m–2 s–1. The inoculum production was carried out by successive inoculations of P. viticola suspensions into new healthy potted plants every month.

Gas Exchange Variables in Cultivars Resistant and Susceptible to Downy Mildew

Seven experiments were performed, one per cultivar, with the three resistant and the four susceptible cultivars. Five plants of each cultivar were inoculated as described above and kept at humid chamber for 12 h. As a control, five plants were sprayed with distilled water. The experiments were performed in the full randomized design. Net photosynthetic rate (A), stomatal conductance (g), intercellular CO2 concentration (C) and transpiration (E) were measured on the inoculation day (day 0) and 1, 4, and 12 days after the inoculation (dai) of P. viticola in all resistant cultivars, except for cv. Solaris in which the gas exchange was evaluated at 0, 1, 2, 5, and 12 dai. In the susceptible cultivars the gas exchange variables were measured at 0, 1, 2, 6, and 12 dai. The instantaneous carboxylation efficiency (k) was estimated as A/Ci (Machado et al., 2005) in each measurement of leaf gas exchange. The evaluations were performed with portable infrared gas analyzers (GFS3000, Heinz Walz GmbH, Germany; and Li-6400-XT, LI-COR Inc., United States) equipped with fluorometers (PAM 3055 and 6400-40, respectively). Cuvette temperature was set to 25°C and relative humidity was 60% at the inlet of the cuvette. The cuvette sizes were 2 and 4 cm2, respectively for Li-6400-XT and GFS3000. CO2 concentration was kept constant at 400 μmol mol−1. A LED array provided PAR of 1000 μmol m–2 s–1. A, gs, Ci and E were recorded after fluxes had stabilized. The leaf areas evaluated for gas exchange were photographed at 4 and 12 dai and digital images were processed with Quant software (Vale et al., 2001) to estimate downy mildew severity. All experiments, except for cv. Phoenix, were repeated once. The average values of photosynthetic variables measured in healthy and diseased plants were compared by Student’s t-test (p ≤ 0.05) using STATISTICA 6.0 software (StatSoft Inc., Tulsa, OK, United States). The normal distribution of data was verified by the Shapiro-Wilk’s test and transformations were performed when necessary (Supplementary Table S1).

Photosynthetic and Photochemical Evaluations in Solaris and Riesling Cultivars

The photosynthetic limitations caused by P. viticola were evaluated in the grapevine cultivars Solaris and Riesling, which are resistant and susceptible to downy mildew, respectively. Three plants of each cultivar were inoculated, and three plants were sprayed with water as described previously. Six evaluations were performed in each experiment from the 5th day after inoculation onward (twice a day) using the GFS3000 (Heinz Walz, GmbH, Germany) equipped with a fluorometer (PAM 3055). Leaf CO2 assimilation response to increasing air CO2 concentration (Ca) was initially evaluated at 400 μmol mol–1, which was gradually changed to 250, 200, 150, and 50 μmol mol–1 and then gradually increased to 400, 600, 800, 1400, and 2000 μmol mol–1. Day respiration (Rd) was obtained through a linear regression between A (response variable) and Ci (predictor variable) under C <400 μmol mol–1 and corresponds to A when Ci is zero. Mesophyll conductance (gm) was estimated as follows (Flexas et al., 2007): where Γ∗ is the photosynthetic compensation point, i.e., the CO2 concentration at which the photorespiratory efflux of CO2 is equal to the CO2 photosynthetic assimilation rate; and J is the transport of electrons from chlorophyll fluorescence assessments. The calculated values of gm were used to convert A-Ci curves into A-Cc (Cc is CO2 concentration at carboxylation sites in chloroplasts) curves using the following equation (Flexas et al., 2007): From A-Cc curves, the maximum Rubisco carboxylation rate (Vcmax) and maximum rate of electron transport driving regeneration of ribulose-1,5-bisphosphate (Jmax) were estimated as follows (Farquhar et al., 1980; Flexas et al., 2007): where Kc and Ko are the Michaelis-Menten constants of Rubisco for carboxylation and oxygenation (Bernacchi et al., 2003), respectively, and O is the internal O2 concentration, which is considered equal to the external O2 concentration. Vcmax and Jmax were estimated by non-linear regressions using STATISTICA 6.0 software (StatSoft Inc., Tulsa, OK, United States). The maximum leaf CO2 assimilation (Amax) and the maximum stomatal conductance (gsmax) were also obtained in each A-Cc curve. The experiment was performed twice and the average values of photosynthetic variables measured in healthy and diseased plants were compared using Student’s t-test (p ≤ 0.05).

Detection and Quantification of Chlorophyll a Fluorescence, H2O2, O2– and APX Activity

Chlorophyll a fluorescence was measured using an Imaging-PAM Maxi (Heinz Walz GmbH, Germany) in cv. Solaris. The measurements were performed 6, 12, 30, 60, 120, and 288 h after P. viticola inoculation (hai). After dark-adapting leaves for 20 min, a saturation light pulse of 2700 μmol m–2s–1 was applied and the maximum PSII quantum yield [Fv/Fm = (Fm−Fo)/Fm] was measured (Maxwell and Johnson, 2000). Fm and Fo denote the maximum and minimum fluorescence of dark-adapted samples, respectively. The detection of H2O2 and O2– was performed in cvs. Solaris and Riesling and its quantification was only performed where H2O2 was detected. The samples were collected at 6, 12, 30, 60, 120, and 288 hai. For quantification of H2O2, five leaf discs of 0.46 cm2 were removed with a cork borer approximately 2 cm away from the central vein of healthy and inoculated leaves. All five discs were immediately transferred to 1 mL of acetone acidified with 25 mM H2SO4 and frozen in liquid nitrogen, where they remained until measurements. For measurements, 50 μL of thawed samples were added to 1 mL of ferrous ammonium sulfate xylenol orange (FOX) solution that contained 250 μM ferrous ammonium sulfate, 100 mM sorbitol, 100 μM xylenol orange and 25 mM H2SO4. The samples were then incubated at room temperature for 30 to 45 min. H2O2 concentration was determined at 550 nm (EPOCH, BioTec, United States) and 850 nm (8453 UV–vis Spectroscopy System, Agilent, United States) using a standard curve ranging from 0 to 100 μM H2O2 (Wolff, 1994). Leaf samples (0.5 g) were harvested and immediately frozen in liquid nitrogen. To assess ascorbate peroxidase (APX) activity, samples were homogenized in 5 mL phosphate buffer (pH 7.6) including 1% polyvinylpyrrolidone (PVP) and 0.1 mM EDTA and centrifuged for 20 min at 16,000 g at 4°C. The supernatant was collected and used as crude extract in the reaction mixtures of the enzyme activity assays. For the APX assay, the 0.3-mL reaction mixture contained 0.5 mM ascorbic acid, 50 mM phosphate buffer, 1 mM EDTA, 0.5 mM H2O and 10–15 μL of the supernatant. The reaction was started by adding 10 μL of 15 mM of H2O2, and APX was assayed spectrometrically (EPOCH, BioTec, United States/8453 UV-VIS Spectroscopy System, Agilent, United States) following the decrease in absorbance at 290 nm (Nakano and Asada, 1981). Histochemical detection of H2O2 and O2– was performed on leaf discs (3 cm2) removed from the same leaves used for the estimation of H2O2 concentration. The samples were collected at 6, 12, 30, 60, 120, and 288 hai from each plant of each treatment. The samples were placed in glass vials containing 3,3′−diaminobenzidine tetrahydrochloride (DAB), pH 5.5, 1%; Sigma−Aldrich, or nitroblue tetrazolium (NBT), pH 6, 0.1%; Sigma−Aldrich, and kept at 25°C for 2 h in darkness (Thordal-Christensen et al., 1997). After this period, the samples were cleared in 80% ethanol for 24 h, stored in glycerol solution (70%), and then observed with an optical microscope. The average values of H2O2 concentrations measured in healthy and diseased plants were compared by Student’s t-test (p ≤ 0.05).

Results

Differential Gas Exchange Responses of Resistant and Susceptible Grapevine Cultivars to Plasmopara viticola Infection

The death-like hypersensitive reactions “HR” were observed 4 dai with P. viticola on the resistant cultivars Solaris, Cabernet Blanc and Phoenix (Supplementary Figure S1). In the susceptible cultivars, the symptoms of downy mildew and P. viticola sporulation were observed 6 to 7 dai. The visual area affected by “HR” did not change for cvs. Cabernet Blanc and Phoenix after 4 dai, and a slow progression (from 5.9% at 4 dai to 10.2% at 12 dai) was exclusively noticed in cv. Solaris. After 12 dai, necrotic areas accounted for 2.3 and 3.3% (on average) of Cabernet Blanc and Phoenix, respectively (Supplementary Figures S1, S2). Typical symptoms of downy mildew, i.e., oil spots on the adaxial leaf surface and abundant sporulation of P. viticola on the abaxial leaf surface (Supplementary Figure S2), were observed on all susceptible cultivars. Disease severity reached 47, 53, 46, and 34% in Riesling, Niagara Rosada, Merlot and Moscato, respectively. On the day of inoculation with P. viticola, photosynthetic rates were similar for inoculated and non-inoculated leaves of resistant cultivars Solaris, Cabernet Blanc and Phoenix (Figure 1). However, significant reductions of leaf photosynthesis (A) were found at 2 dai in Solaris (Figure 1B) and at 1 dai for Cabernet Blanc and Phoenix (Figures 1D,F). At the appearance of “HR” and compared with healthy leaves, A values were 54, 35, and 77% lower on diseased leaves of Solaris, Cabernet Blanc, and Phoenix, respectively (Figures 1B,D,F). Then, the A values for inoculated leaves did not change significantly until 12 dai.
FIGURE 1

Net photosynthetic rate (A) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and * indicates significant differences (p < 0.05).

Net photosynthetic rate (A) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and * indicates significant differences (p < 0.05). Photosynthetic rates ranged from 8.9 to 13.3 μmol m–2 s–1 on susceptible cultivars at the day of inoculation when no differences were found between inoculated and non-inoculated leaves (Figures 1A,C,E,G). When visual symptoms appeared, significant reduction of A was only observed in Niagara Rosada (Figure 1C). After 12 days of inoculation, A values on inoculated leaves were reduced by 56, 56, 51, and 41% compared with healthy leaves of Riesling, Niagara Rosada, Merlot, and Moscato, respectively. No differences in Ci of inoculated and non-inoculated leaves were found before the appearance of “HR” in resistant cultivars (Figure 2), and this situation persisted until 12 dai for Cabernet Blanc and Phoenix. A slight increase (10.5%) was detected on C of inoculated leaves of Solaris 12 dai compared with non-inoculated leaves (Figure 2B). Similar C values were measured on inoculated and non-inoculated leaves of susceptible cultivars until symptoms appearance (Figures 2A,C,E,G). From this time, diseased leaves of Riesling and Niagara Rosada showed increases in C values of 14% compared with healthy leaves (Figures 2A,C). Similar to Cabernet Blanc and Phoenix (Figures 2D,F), C of Merlot and Moscato remained stable from 6 to 12 dai regardless of the treatment (Figures 2E,G).
FIGURE 2

Intercellular CO2 concentration (Ci) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and * indicates significant differences (p < 0.05).

Intercellular CO2 concentration (Ci) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and * indicates significant differences (p < 0.05). Considering downy mildew resistant cultivars, g and E were always lower in diseased leaves than in healthy leaves from 4 dai (Figures 3, 4). Differences in g between inoculated and non-inoculated leaves of susceptible cultivars were observed only at 12 dai in Riesling, Merlot, and Moscato (Figures 3A,E,G). Nevertheless, no reductions in E were found in susceptible cultivars with the exception of Riesling at 12 dai when the average E value of diseased leaves was reduced by 40% compared with healthy leaves (Figure 4A).
FIGURE 3

Stomatal conductance (g) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and ∗ indicates significant differences (p < 0.05).

FIGURE 4

Transpiration rate (E) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and * indicates significant differences (p < 0.05).

Stomatal conductance (g) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and ∗ indicates significant differences (p < 0.05). Transpiration rate (E) in Vitis spp. cultivars. Downy mildew symptoms were observed in Riesling (A), Niagara Rosada (C), Merlot (E) and Moscato (G) 6 days after inoculation. Cell death-like hypersensitive reactions “HR” were observed in Solaris (B), Cabernet Blanc (D) and Phoenix (F) 4 days after inoculation (Inoc.). The average values of healthy and diseased leaves were compared using the Student’s t test for each cultivar (n = 6, ± SE), and * indicates significant differences (p < 0.05). Significant reductions in the instantaneous carboxylation efficiency (k) due to downy mildew were measured from 1 dai for resistant cultivars and at 12 dai for susceptible cultivars (Supplementary Figure S3). In Niagara Rosada, reduction in k was also observed in diseased leaves at 6 dai (Supplementary Figure S3C). Healthy plants were more responsive to increasing chloroplastic CO2 concentrations when compared with P. viticola inoculated plants regardless of the cultivar (Figure 5). The maximum leaf CO2 assimilation (A) was 18.7 and 12.7 μmol m–2 s–1 in healthy and diseased leaves of Solaris, respectively, and a reduction of 23% was noticed due to P. viticola infection. V was similar in healthy and diseased leaves of Solaris; however, J was reduced by 34% in diseased leaves (Table 1). The maximum stomatal conductance (g) was reduced (−27%) in Solaris without any significant change in g (Table 1). In Riesling, A was 22% higher in healthy leaves than in diseased leaves. Unlike Solaris, biochemical limitations were higher in diseased leaves of Riesling, with plants showing low V and J (Table 1). In addition, no significant diffusive limitation was found in diseased leaves compared with healthy leaves of Riesling, as suggested by g and g (Table 1).
FIGURE 5

Response of leaf CO2 assimilation (A) to increasing CO2 concentration at carboxylation sites in chloroplasts (C) of Vitis vinifera cv. Solaris (A) and Riesling (B). Open symbols represent mean values (n = 3, ± SE) of healthy leaves, while filled symbols indicate diseased leaves.

TABLE 1

Maximum Rubisco carboxylation rate (V), maximum rate of electron transport driving regeneration of ribulose-1,5-bisphosphate (J), mesophyll conductance (g), maximum stomatal conductance (g), and maximum photosynthetic rate (A) in Solaris and Riesling cultivars as affected by downy mildew (Plasmopara viticola).

VariableHealthy plantsInfected plantst-test
Solaris
Vcmax (μmol m–2 s–1)83.69.670.511.1ns
Jmax (μmol m–2 s–1)131.29.286.47.7*
gm (mol m–2 s–1)0.110.030.050.02ns
gsmax (mmol m–2 s–1)143.56.5105.15.6*
Amax (μmol m–2 s–1)18.70.712.71.3*
Riesling
Vcmax (μmol m–2 s–1)58.45.933.00.8*
Jmax (μmol m–2 s–1)88.15.565.70.2*
gm (mol m–2 s–1)0.100.020.080.01ns
gsmax (mmol m–2 s–1)117.212.690.52.4ns
Amax (μmol m–2 s–1)13.10.610.10.2*
Response of leaf CO2 assimilation (A) to increasing CO2 concentration at carboxylation sites in chloroplasts (C) of Vitis vinifera cv. Solaris (A) and Riesling (B). Open symbols represent mean values (n = 3, ± SE) of healthy leaves, while filled symbols indicate diseased leaves. Maximum Rubisco carboxylation rate (V), maximum rate of electron transport driving regeneration of ribulose-1,5-bisphosphate (J), mesophyll conductance (g), maximum stomatal conductance (g), and maximum photosynthetic rate (A) in Solaris and Riesling cultivars as affected by downy mildew (Plasmopara viticola). Chlorophyll fluorescence from healthy and inoculated leaves of Solaris indicated similar maximum PSII quantum efficiency (F/F) at 12 hai (0.770 and 0.766). Reduced F/F in inoculated leaves was found at 5 dai (0.735) when “HR” symptoms were visible (Figure 6). In non-inoculated leaves, F/F was 0.752. Although the effects of P. viticola infection are clearly visible by imaging of chlorophyll a fluorescence, no significant difference was detected on any day (Supplementary Figure S4).
FIGURE 6

Digital and chlorophyll fluorescence images of grapevine leaves of cv. Solaris after 12 h and 5 days after infection (dai) with Plasmopara viticola. Chlorophyll fluorescence images show the maximum PSII quantum efficiency (F/F). Circles indicate the area used for calculation of F/F. Leaves were dark-adapted prior to measurements. The false color code depicted at the bottom of images represents the range of 0 (black) to 1 (white).

Digital and chlorophyll fluorescence images of grapevine leaves of cv. Solaris after 12 h and 5 days after infection (dai) with Plasmopara viticola. Chlorophyll fluorescence images show the maximum PSII quantum efficiency (F/F). Circles indicate the area used for calculation of F/F. Leaves were dark-adapted prior to measurements. The false color code depicted at the bottom of images represents the range of 0 (black) to 1 (white).

Oxidative Burst in the Resistant Cultivar Solaris

The presence of H2O2 as indicated by DAB staining was detected as brown spots close to the secondary veins of diseased leaves (Figure 7B). Brown spots were frequent and intense at 12 and 30 hai (Supplementary Figure S5). Necrotic tissue formed on diseased leaves from 120 hai, and no reaction to DAB staining was observed in diseased leaves at this time (data not shown). Brown spots were not detected in healthy leaves. H2O2 concentrations were lower than 0.5 μmol g–1 fresh mass (FM) in both healthy and diseased leaves of Solaris at 6 hai and from 120 hai. However, increases in leaf H2O2 concentration and data variability were noticed at 12 and 60 hai (Figure 7). O2 was detected as dark blue spots similar in shape to those spots observed in H2O2 detection and near secondary veins. Similar to H2O2, the frequency and intensity of blue spots were high at 12 hai (Figure 7C). No blue spots were observed at 30 hai in diseased leaves, and no O2 was detected in healthy leaves (data not shown). No significant difference was observed for APX activity between healthy and diseased leaves, although the highest APX activity (2.33 μmol H2O2 g–1 FM min–1) was noticed in diseased leaves at 288 hai (Supplementary Figure S6).
FIGURE 7

Time-course of leaf H2O2 concentration (A) in healthy and P. viticola infected leaves of cv. Solaris after inoculation and visual localization of H2O2 (B) and O2– (C) in P. viticola inoculated leaves, 12 h after inoculation. Scale bar = 50 μm.

Time-course of leaf H2O2 concentration (A) in healthy and P. viticola infected leaves of cv. Solaris after inoculation and visual localization of H2O2 (B) and O2– (C) in P. viticola inoculated leaves, 12 h after inoculation. Scale bar = 50 μm.

Discussion

The results of this study indicate that P. viticola infection is associated with a cost in terms of CO2 assimilation per unit of leaf area, i.e., resistant cultivars reduced A earlier and more severely than susceptible cultivars. On average, 55% reduction in leaf CO2 assimilation was noticed few hours after pathogen contact with its host, and such a decrease remained up to 12 days after P. viticola inoculation in resistant cultivars Solaris, Cabernet Blanc and Phoenix. This behavior could result from the downregulation of proteins related to photosynthesis and carbohydrate metabolism as previously reported in resistant grapevine cultivars (Figueiredo et al., 2017; Nascimento-Gavioli et al., 2017). A fast and massive transcription of signaling genes affected by P. viticola has also been reported on Rpv-bearing in grapevine cultivars (Fröbel et al., 2019). This response presumably is initiated by a receptor gene that triggers the signal transduction cascade such as that noted in a pathogen-associated molecular pattern (PAMP, Fröbel et al., 2019). In contrast, no effect on leaf gas exchange was observed during the infection and colonization of P. viticola in susceptible cultivars. The net photosynthetic rate on susceptible cultivars was reduced only when symptoms appeared at approximately 6 days after pathogen inoculation. In the susceptible cultivar Trincadeira there was no reduction of photosynthetic pigments content 24 h after P. viticola inoculation (Nascimento et al., 2019). This suggests that the photosynthetic rate remains unaltered on the first stages of the interaction where there are no visible lesions or symptoms (Nascimento et al., 2019). Similarly to our results, in the biotroph interaction between barley and powdery mildew (Blumeria graminis), no reduction of photosynthesis was observed in compatible interactions up to 3 dai. However, the photosynthesis in those leaves was progressively reduced up to 60% compared with healthy leaves at 7 dai (Swarbrick et al., 2006). Plasmopara viticola infection increased stomatal closure and decreased transpiration rates of resistant cultivars from 4 dai. Stomatal closure is a mechanism related to grapevine resistance to downy mildew, impairing pathogen penetration and sporulation (Allègre et al., 2009). As soon as P. viticola zoospores are released from sporangia, they move toward stomata by chemotaxis, encyst and rapidly penetrate leaf tissue through the stomata in susceptible cultivars. Partial stomatal closure contributes to avoiding P. viticola penetration in resistant grapevine cultivars, although it is insufficient to limit infection (Allègre et al., 2009). Stomatal conductance was reduced only 12 days after inoculation in most susceptible cultivars and was less affected when compared with photosynthesis. In cv. Marselan, a susceptible cultivar to downy mildew, P. viticola affects guard cell function and interferes with stomatal opening, as observed for other oomycetes (Allègre et al., 2007). For instance, Phytophthora infestans can cause abnormal stomata functioning in potato leaves to facilitate the emergence of sporangiophores and sporulation in dark conditions (Farrell et al., 1969). Although drastic alterations were observed in A, g and E of inoculated resistant grapevines, no variation was noticed in C of Cabernet Blanc and Phoenix leaves while it increased in inoculated leaves of Solaris at 12 dai. This was somehow unexpected since decreases in C are expected in leaves due to decreases in stomatal conductance. However, the low variation in C of resistant cultivars and the reduction in the instantaneous carboxylation efficiency after the infection (Supplementary Figure S3) in parallel to stomatal closure suggest that another physiological process controls intercellular CO2 concentrations of P. viticola-infected leaves. A rapid increase in the respiration rate of barley leaves resistant to Erysiphe graminis was reported after pathogen contact, while such increase in respiration was found only 3 days after inoculation in compatible interactions (Smedegaard-Petersen and Tolstrup, 1985). A late increase in C after the development of disease symptoms was observed in downy mildew-susceptible cultivars, which is in agreement with the linear increase in C with increasing rust severity observed in V. labrusca infected with Phakopsora euvitis (Nogueira Júnior et al., 2017). The activation of signal transduction pathways following the recognition of P. viticola by resistant grapevines leads to the production of reactive oxygen species (Aziz et al., 2003), such as O2– and H2O2 (Kortekamp and Zyprian, 2003). H2O2 production is one of the earliest detectable cytological responses to downy mildew in resistant cultivars, as noticed in Solaris 12 h after inoculation (our data and Trouvelot et al., 2008). The production of reactive oxygen species is associated with localized cell death and is crucial for delimiting the growth of pathogen into leaf tissues (Aziz et al., 2003). Reactive oxygen species cause oxidative damage to a variety of macromolecular targets, such as DNA, lipids and proteins, including enzymes of the Calvin–Benson Cycle and proteins located in the thylakoid membrane, leading to irreversible damage and ultimately tissue necrosis (Halliwell, 2006). Oxidative stress blocks key steps in chlorophyll biosynthesis by directly or indirectly inhibiting the enzymes involved in chlorophyll synthesis (Aarti et al., 2006). Herein, necrotic areas were observed 4 days after inoculation and ranged from 1.8 to 6.0% of total leaf area in resistant cultivars Solaris, Cabernet Blanc and Phoenix. Solaris inoculated with P. viticola showed low maximum PSII quantum efficiency 5 days after inoculation and low regeneration of RuBP driven by electron transport (Figure 6 and Table 1). This photochemical limitation is presumably related to the production of H2O2 and the degradation of cell photosynthetic pigments. As shown in Figure 6, the occurrence of distinct spots indicative of low photosynthetic efficiency is clearly visible upon infection with P. viticola. These spots might contribute to overall reduction on photosynthetic performance. In grapevine leaves of cv. Marselan infected with P. viticola, chlorophyll concentrations were reduced by approximately 30% (Gamm et al., 2011). Measurements conducted in spots and adjacent areas without symptoms confirmed that F/F was reduced in the spots, but the asymptomatic adjacent areas had F/F values comparable to non-infected leaves (Gamm et al., 2011). In the present study, F/F tends to decrease at 5 dai (although not statistically significant). This might be due to the occurrence of necrotic spots that are present in the leaf area considered for calculation of F/F. Reduced F/F is an indicator for impaired PSII functionality due to damage to reaction center proteins by reactive oxygen species, and this finding is consistent with reduced assimilation rates. In susceptible cultivars, such as Riesling (this study), Sangiovese (Moriondo et al., 2005) and Niagara Rosada (Nogueira Júnior et al., 2019), the increase in photochemical limitations is also observed during downy mildew symptom expression. Low concentrations of photosynthetic pigments in the yellow halos surrounding P. viticola lesions are related with the reduction in photochemical reactions in these cultivars (Moriondo et al., 2005). Expression of plant defenses is often assumed to be costly, requiring diversion of resources away from plant growth and development (Walters and Boyle, 2005), and P. viticola greatly enhances carbohydrate hydrolysis and represses photosynthesis-associated proteins in resistant cultivars (Nascimento-Gavioli et al., 2017). In the incompatible interaction of barley and Erysiphe graminis f. sp. hordei, pathogen recognition occurs during haustoria development into host cells. After recognition, several histochemical and biochemical alterations are induced in hosts (Smedegaard-Petersen and Tolstrup, 1985). The cost to the net photosynthetic rate in barley leaves challenged by E. graminis f. sp. hordei is also similar to the cost observed in downy mildew-resistant grapevines. However, in the incompatible interaction of Eucalyptus grandis and Puccinia psidii, the host reacts immediately during pathogen appressoria formation. Consequently, only small flecks are observed in the host before pathogen penetration and photosynthetic activity is unaffected (Alves et al., 2011). Although resistant cultivars decrease the photosynthetic activity when infected by P. viticola, its use provides a substantial contribution to the sustainability of viticulture while reducing pesticide applications (Buonassisi et al., 2017). P. viticola can complete its life cycle in resistant cultivars. However, low oomycete sporulation is noted compared with the sporulation in susceptible cultivars, and sporangiophores are abnormal. The reduced sporulation affects the secondary spread of the pathogen and reduces the rate of disease progression in field (Parlevliet, 1979). In the last 20 years, newly bred grapevine cultivars with disease resistance, desirable agronomic attributes and good enological characteristics have been introduced to the market (Buonassisi et al., 2017). These cultivars bear different Rpv resistance loci to downy mildew and several studies have been developed to identify the resistance loci on these materials. Locus Rpv 3 is present in all resistant cultivars used in the present work. However, the cultivar Solaris has also the Rpv10 locus[1]. In general, the effect of P. viticola on photosynthesis was similar for cvs. Solaris, Cabernet Blanc and Phoenix.

Conclusion

In conclusion, P. viticola differently affects leaf gas exchange in resistant and susceptible cultivars, with resistant cultivars exhibiting faster responses. While net photosynthetic rate and stomatal conductance of resistant cultivars are rapidly and negatively affected by P. viticola, susceptible cultivars only exhibit reduced leaf gas exchange at late stages of disease development. These results would suggest that the recognition followed by induced resistance responses is costly for grapevine carbon assimilation. In the absence of recognition, the pathogen colonizes the leaf tissue without affecting photosynthesis until the appearance of visual symptoms. Although the induction of defenses in resistant grapevine cultivars is costly, the use of genetic resistance to control downy mildew as an alternative to chemical control is highly desirable and complementary for a sustainable viticulture.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

AN and LA devised the project and the main conceptual ideas. AN and MT performed the experiments. AN, LA, RR, and AT revised critically the project and the manuscript for important intellectual content. All authors discussed the results and contributed to the final manuscript.

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
  34 in total

Review 1.  Fitness costs of induced resistance: emerging experimental support for a slippery concept.

Authors:  Martin Heil; Ian T Baldwin
Journal:  Trends Plant Sci       Date:  2002-02       Impact factor: 18.313

2.  Genetic dissection of sex determinism, inflorescence morphology and downy mildew resistance in grapevine.

Authors:  Elisa Marguerit; Christophe Boury; Aurélie Manicki; Martine Donnart; Gisèle Butterlin; Alice Némorin; Sabine Wiedemann-Merdinoglu; Didier Merdinoglu; Nathalie Ollat; Stéphane Decroocq
Journal:  Theor Appl Genet       Date:  2009-02-24       Impact factor: 5.699

3.  Histological responses to downy mildew in resistant and susceptible grapevines.

Authors:  Ruiqi Liu; Lan Wang; Jiali Zhu; Tingting Chen; Yuejin Wang; Yan Xu
Journal:  Protoplasma       Date:  2014-07-16       Impact factor: 3.356

4.  Induced defenses and phenotypic plasticity.

Authors:  D Cipollini
Journal:  Trends Ecol Evol       Date:  1998-05       Impact factor: 17.712

5.  Early stage metabolic events associated with the establishment of Vitis vinifera - Plasmopara viticola compatible interaction.

Authors:  Rui Nascimento; Marisa Maia; António E N Ferreira; Anabela B Silva; Ana Ponces Freire; Carlos Cordeiro; Marta Sousa Silva; Andreia Figueiredo
Journal:  Plant Physiol Biochem       Date:  2019-01-25       Impact factor: 4.270

6.  Proteome of Plasmopara viticola-infected Vitis vinifera provides insights into grapevine Rpv1/Rpv3 pyramided resistance to downy mildew.

Authors:  Maria Carolina Andrade Nascimento-Gavioli; Sarah Zanon Agapito-Tenfen; Rubens Onofre Nodari; Leocir José Welter; Fernando David Sanchez Mora; Luciano Saifert; Aparecido Lima da Silva; Miguel Pedro Guerra
Journal:  J Proteomics       Date:  2016-05-26       Impact factor: 4.044

7.  Construction of a reference linkage map of Vitis amurensis and genetic mapping of Rpv8, a locus conferring resistance to grapevine downy mildew.

Authors:  Paule Blasi; Sophie Blanc; Sabine Wiedemann-Merdinoglu; Emilce Prado; Ernst H Rühl; Pere Mestre; Didier Merdinoglu
Journal:  Theor Appl Genet       Date:  2011-03-15       Impact factor: 5.699

8.  Metabolic consequences of susceptibility and resistance (race-specific and broad-spectrum) in barley leaves challenged with powdery mildew.

Authors:  Philip J Swarbrick; Paul Schulze-Lefert; Julie D Scholes
Journal:  Plant Cell Environ       Date:  2006-06       Impact factor: 7.228

9.  Laminarin elicits defense responses in grapevine and induces protection against Botrytis cinerea and Plasmopara viticola.

Authors:  Aziz Aziz; Benoît Poinssot; Xavier Daire; Marielle Adrian; Annie Bézier; B Lambert; Jean-Marie Joubert; Alain Pugin
Journal:  Mol Plant Microbe Interact       Date:  2003-12       Impact factor: 4.171

10.  The Rpv3-3 Haplotype and Stilbenoid Induction Mediate Downy Mildew Resistance in a Grapevine Interspecific Population.

Authors:  Silvia Vezzulli; Giulia Malacarne; Domenico Masuero; Antonella Vecchione; Chiara Dolzani; Vadim Goremykin; Zeraye Haile Mehari; Elisa Banchi; Riccardo Velasco; Marco Stefanini; Urska Vrhovsek; Luca Zulini; Pietro Franceschi; Claudio Moser
Journal:  Front Plant Sci       Date:  2019-03-06       Impact factor: 5.753

View more
  6 in total

1.  Blocked at the Stomatal Gate, a Key Step of Wheat Stb16q-Mediated Resistance to Zymoseptoria tritici.

Authors:  Mélissa Battache; Marc-Henri Lebrun; Kaori Sakai; Olivier Soudière; Florence Cambon; Thierry Langin; Cyrille Saintenac
Journal:  Front Plant Sci       Date:  2022-06-27       Impact factor: 6.627

2.  Transcriptional Profiling of Resistant and Susceptible Cultivars of Grapevine (Vitis L.) Reveals Hypersensitive Responses to Plasmopara viticola.

Authors:  Peijie Gong; Jun Kang; Ehsan Sadeghnezhad; Ruoxuan Bao; Mengqing Ge; Yaxian Zhuge; Lingfei Shangguan; Jinggui Fang
Journal:  Front Microbiol       Date:  2022-04-25       Impact factor: 6.064

3.  The Molecular Priming of Defense Responses is Differently Regulated in Grapevine Genotypes Following Elicitor Application against Powdery Mildew.

Authors:  Chiara Pagliarani; Amedeo Moine; Walter Chitarra; Giovanna Roberta Meloni; Simona Abbà; Luca Nerva; Massimo Pugliese; Maria Lodovica Gullino; Giorgio Gambino
Journal:  Int J Mol Sci       Date:  2020-09-15       Impact factor: 5.923

4.  Proteomic analysis of early-stage incompatible and compatible interactions between grapevine and P. viticola.

Authors:  Guo-Tian Liu; Bian-Bian Wang; David Lecourieux; Mei-Jie Li; Ming-Bo Liu; Rui-Qi Liu; Bo-Xing Shang; Xiao Yin; Li-Jun Wang; Fatma Lecourieux; Yan Xu
Journal:  Hortic Res       Date:  2021-05-01       Impact factor: 6.793

Review 5.  Powdery Mildew Resistance Genes in Vines: An Opportunity to Achieve a More Sustainable Viticulture.

Authors:  Viviana Sosa-Zuniga; Álvaro Vidal Valenzuela; Paola Barba; Carmen Espinoza Cancino; Jesus L Romero-Romero; Patricio Arce-Johnson
Journal:  Pathogens       Date:  2022-06-18

6.  Mono-Locus and Pyramided Resistant Grapevine Cultivars Reveal Early Putative Biomarkers Upon Artificial Inoculation With Plasmopara viticola.

Authors:  Ramona Mihaela Ciubotaru; Pietro Franceschi; Luca Zulini; Marco Stefanini; Domen Škrab; Marcia Denise Rossarolla; Peter Robatscher; Michael Oberhuber; Urska Vrhovsek; Giulia Chitarrini
Journal:  Front Plant Sci       Date:  2021-07-01       Impact factor: 5.753

  6 in total

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