Gas Exchange Rates Decrease and Leaf Temperature Increases in Nicotiana benthamiana Leaves Transiently Overexpressing Hemagglutinin in an Agrobacterium-Assisted Viral Vector System

In this study, gas exchange characteristics and temperature of Nicotiana benthamiana leaves transiently overexpressing hemagglutinin (HA), an influenza vaccine antigen, with an Agrobacterium tumefaciens-assisted viral vector were investigated. Inoculation of leaves with an empty viral vector not containing the HA gene decreased the net photosynthetic rate (Pn) and transpiration rate (T) from 2 to 3 days post-infiltration (DPI) in the A. tumefaciens suspension. Expression of HA with the vector decreased Pn and T to much lower levels until 4 DPI. Such significant decreases were not observed in leaves infiltrated with suspension of A. tumefaciens not carrying the viral vector or in uninfiltrated leaves. Thus, viral vector inoculation itself decreased Pn and T to a certain extent and the HA expression further decreased them. The decreases in Pn and T in empty vector-inoculated and HA expression vector-inoculated leaves were associated with decreases in stomatal conductance, suggesting that the reduction of gas exchange rates was caused at least in part by stomatal closure. More detailed gas exchange and chlorophyll fluorescence analyses revealed that in HA vector-inoculated leaves, the capacity of ribulose-1,5-bisphosphate carboxylase/oxygenase to assimilate CO2 and the capacity of photosynthetic electron transport in planta were downregulated, which contributed also to the decrease in Pn. Leaf temperature (LT) increased in viral vector-inoculated leaves, which was associated with the decrease in T. When HA vector-inoculated leaves were grown at air temperatures (ATs) of 21, 23, and 26°C post-infiltration, HA accumulated earlier in leaves and the days required for HA content to attain its peak became shorter, as AT was higher. The highest LT was found 1–2 days earlier than the highest leaf HA content under all post-infiltration AT conditions. This phenomenon could be applicable in a non-destructive technique to detect the optimum harvesting date for individual plants to determine the day when leaf HA content reaches its maximum level, irrespective of spatiotemporal variation of AT, in a plant growth facility.


Gas Exchange Rates Decrease and Leaf Temperature Increases in Nicotiana benthamiana Leaves Transiently Overexpressing Hemagglutinin in an Agrobacterium-Assisted Viral Vector System
In this study, gas exchange characteristics and temperature of Nicotiana benthamiana leaves transiently overexpressing hemagglutinin (HA), an influenza vaccine antigen, with an Agrobacterium tumefaciens-assisted viral vector were investigated. Inoculation of leaves with an empty viral vector not containing the HA gene decreased the net photosynthetic rate (P n ) and transpiration rate (T) from 2 to 3 days post-infiltration (DPI) in the A. tumefaciens suspension. Expression of HA with the vector decreased P n and T to much lower levels until 4 DPI. Such significant decreases were not observed in leaves infiltrated with suspension of A. tumefaciens not carrying the viral vector or in uninfiltrated leaves. Thus, viral vector inoculation itself decreased P n and T to a certain extent and the HA expression further decreased them. The decreases in P n and T in empty vectorinoculated and HA expression vector-inoculated leaves were associated with decreases in stomatal conductance, suggesting that the reduction of gas exchange rates was caused at least in part by stomatal closure. More detailed gas exchange and chlorophyll fluorescence analyses revealed that in HA vector-inoculated leaves, the capacity of ribulose-1,5-bisphosphate carboxylase/oxygenase to assimilate CO 2 and the capacity of photosynthetic electron transport in planta were downregulated, which contributed also to the decrease in P n . Leaf temperature (LT) increased in viral vector-inoculated leaves, which was associated with the decrease in T. When HA vector-inoculated leaves were grown at air temperatures (ATs) of 21, 23, and 26 • C post-infiltration, HA accumulated earlier in leaves and the days required for HA content to attain its peak became shorter, as AT was higher. The highest LT was found 1-2 days earlier than the highest leaf HA content under all post-infiltration AT conditions. This phenomenon could be applicable in a non-destructive technique to detect the optimum harvesting date for individual plants to determine the day when leaf HA content reaches its maximum level, irrespective of spatiotemporal variation of AT, in a plant growth facility.

INTRODUCTION
Plants have become a viable platform for producing recombinant proteins such as biopharmaceuticals. Compared with transgenic animal-based or mammalian cell culture-based systems, plantbased systems offer advantages for biopharmaceutical production that include low production cost, high production scalability, and minimum risk of contamination with human pathogens, whereas animal and plant cells have a similar eukaryotic protein synthesis pathway (Fischer and Emans, 2000;Ma et al., 2003;Twyman et al., 2003;Desai et al., 2010). Indeed, two plantmade biopharmaceuticals, taliglucerase alfa (ELELYSO R ) and interferon-α (Interberry α R ), have been approved for human and animal use, respectively, and are already available in the market (Drake et al., 2017). In addition, a seasonal quadrivalent influenza vaccine was recently subjected to a Phase 3 efficacy study in multiple countries 1 and several other candidates have shown their initial safety and efficacy in early-phase clinical trials (Fischer et al., 2012). There are two main biotechnologies for transgene expression in plants: stable transformation of the nuclear or chloroplast genome and transient gene expression using plant viral or non-viral vectors. Among them, transient gene overexpression in Nicotiana species such as Nicotiana benthamiana with viral vectors enables rapid mass production of recombinant proteins (Pogue et al., 2002(Pogue et al., , 2010Gleba et al., 2004Gleba et al., , 2007Lico et al., 2008;Matoba et al., 2011;Chen et al., 2013;Peyret and Lomonossoff, 2015). The accumulation level of a target recombinant protein using this technology is reported to reach as high as 40% of total soluble protein (Marillonnet et al., 2005) or 3-4 g kg −1 on a leaf fresh mass (FM) basis (Marillonnet et al., 2005;Hamorsky et al., 2015) within days post-vector inoculation.
It is expected that massive synthesis and accumulation of a foreign protein in leaves would affect the physiology of the plant. Recent studies have shown that transient overexpression of a recombinant protein using a viral vector could cause severe stress, such as the endoplasmic reticulum (ER) stress, to N. benthamiana plants, leading to reduction in protein yield (Hamorsky et al., 2015;Roychowdhury et al., 2018). Viral vector inoculation would also induce defensive responses of plants, including the hypersensitive response (Lam et al., 2001;Mandadi and Scholthof, 2013) and systemic acquired resistance (Fu and Dong, 2013), which affect plant physiology. Understanding the physiological changes caused by viral vector inoculation and recombinant protein expression in host plants can contribute to improving plant conditions that are favorable for recombinant protein production. However, compared to biotechnology for target protein expression or efficacy of target proteins as biopharmaceuticals, less attention has been paid to physiological aspects of plants in plant-based biopharmaceutical production.
Gas exchange, i.e., CO 2 and water vapor transport driven by photosynthesis, respiration, and transpiration, is a fundamental physiological process in leaves, and is involved in supplying the energy and compounds required for various metabolisms. Given that virus infection often impairs photosynthetic activity 1 http://medicago.com and stomatal opening (Chaerle and Van Der Straeten, 2001;Rolfe and Scholes, 2010;Grimmer et al., 2012), such negative effects may be observed in plants inoculated with a virusderived vector for transient gene expression. Massive synthesis and accumulation of a foreign protein in plant cells can also alter the composition of endogenous proteins, such as significant reduction of ribulose-1,5-bisphosphate carboxylase/oxygenase (RuBisCO; Marillonnet et al., 2004;Bally et al., 2009), which is the key enzyme for photosynthetic CO 2 assimilation. In addition, stomatal malfunction is expected to influence thermal conditions and temperature of leaves via alteration of latent heat dissipation by transpiration. However, to the best of our knowledge, no study has been undertaken to examine the effects of transient overexpression of biopharmaceutical proteins using viral vectors on gas exchange characteristics of leaves.
The objective of the present study was to reveal gas exchange characteristics and temperature of N. benthamiana leaves transiently overexpressing hemagglutinin (HA), an influenza vaccine antigen, in an Agrobacterium tumefaciens-assisted viral vector system. We performed a series of experiments. In experiments (Exp.) 1 and 2, time course of net photosynthetic rate (P n ), transpiration rate (T), stomatal conductance (g s ), and leaf temperature (LT) in leaves post-inoculation of the viral vector were examined. Post-inoculation plant growth was also evaluated. In Exp. 3, photosynthetic gas exchange characteristics and electron transport capacity of photosystem II (PSII) were analyzed in detail. In Exp. 4, changes in leaf HA content and LT were investigated under different air temperature (AT) conditions.

Plant Material and Growth Conditions
Seeds of N. benthamiana were sown into rockwool cubes (AO36/40, ROCKWOOL B.V., Roermond, Netherlands) and seedlings were grown in a growth room. Light was provided by white fluorescent lamps for 16 h d −1 at a photosynthetic photon flux density (PPFD) of 200 µmol m −2 s −1 at the surface of the rockwool cubes. ATs were 25/20 • C (day/night) and mean relative humidity was 30-70%. CO 2 concentration was not lower than an atmospheric level of 400 µmol mol −1 even during the day. The rockwool cubes were subirrigated with tap water for the first week and subsequently with a nutrient solution (prescription A, OAT Agrio Co., Ltd., Tokyo, Japan) at an electrical conductivity of 0.18 S m −1 and pH of 6. At 14 days after sowing, seedlings with the second true leaf approximately 3 mm in length were transplanted onto rockwool blocks (Delta 6.5G, ROCKWOOL B.V.) and grown further in the same room until 35, 36, or 38 days after sowing (depending on experiment) at which time they were subjected to the treatments described below. The rockwool blocks were subirrigated with the above nutrient solution. Lateral shoots and flower buds were removed once a week.

Vector Construction and Agroinfiltration
A "deconstructed" tobamoviral replicon system (Marillonnet et al., , 2005 [magnICON R , ICON Genetics GmbH, Halle (Saale), Germany] was used to overexpress HA in N. benthamiana. Vector construction and A. tumefaciens (GV3101::pMP90, Koncz and Schell, 1986) transformation have been described previously (Matsuda et al., 2012(Matsuda et al., , 2017a. Briefly, the HA construct was designed to target the ectodomain of HA derived from influenza A virus (subtype H1N1, strain A/California/07/2009) to the ER, using an N-terminal secretory signal peptide and a C-terminal HDEL ER-retention signal peptide. The construct was subcloned into the plasmid vector (pICH26212), which contained genes encoding tobamovirusderived RNA-dependent RNA polymerase and movement protein, to form the HA expression vector (pNM216). Suspension of A. tumefaciens was infiltrated into leaves for inoculation, as described previously (Matsuda et al., 2012(Matsuda et al., , 2018. The optical density of the bacterial suspension at 600 nm was 0.03.

Treatments
In Exp. 1, plants were inoculated with the HA expression vector (pNM216) or an empty vector without the HA gene (pICH26212) by vacuum infiltration of A. tumefaciens suspension at a gage pressure of −85 kPa for 2 min or were uninfiltrated. In Exp. 2, plants were inoculated with either of the two vectors or with a "mock" inoculum, i.e., infiltrated with a suspension of A. tumefaciens not carrying any plasmid vectors. In Exp. 3, the 8th or 9th oldest leaf from each plant was used. One half of the leaf lamina bordering the primary vein was inoculated with the HA vector by syringe infiltration. The other half was uninfiltrated. In Exp. 1-3, plants were thereafter grown in growth chambers. Light was provided by white light-emitting diodes (LEDs) for 16 h d −1 at 200 µmol m −2 s −1 PPFD, which was measured and adjusted at the tops of plants. AT was maintained at 22 ± 1 • C throughout the day. Mean relative humidity was 55-80%. The chambers were ventilated with an air pump at an air exchange rate of approximately 0.5 h −1 to keep CO 2 concentration close to the atmospheric level. In Exp. 4, plants were inoculated with the HA vector or were uninfiltrated. Plants were thereafter grown in the growth chambers under the same light environment at an AT of 21, 23, or 26 • C throughout the day. In Exp. 4, experiments were replicated twice for each AT treatment.

Gas Exchange Measurements
In Exp. 1-3, gas exchange rates of the 8th or 9th oldest leaves were measured using a portable gas-exchange measurement system (LI-6400XT, LI-COR Inc., Lincoln, NE, United States). Light was provided by white LEDs (NSPW310DS, Nichia Corp., Tokushima, Japan). In Exp. 1 and 2, measurements were taken at a PPFD of 200 µmol m −2 s −1 and an atmospheric CO 2 concentration (C a ) of 400 µmol mol −1 . In Exp. 3, measurements were taken at PPFDs of 250 and 1,200 µmol m −2 s −1 and at mean C a of 98,194,292,390,589,787,987,1,186, and 1,784 µmol mol −1 . The coefficient of variation of C a was within 1.2%, irrespective of the C a level. In all measurements, the leaf chamber was maintained at an AT equal to that for plant growth and a leaf-to-air vapor pressure deficit lower than 1.1 kPa. The air change rate of the leaf chamber was regulated at 0.2 s −1 . Gas exchange parameters including intercellular CO 2 concentration (C i ) were calculated based on methods described by von Caemmerer and Farquhar (1981).

Chlorophyll Fluorescence Measurements
In Exp. 3, chlorophyll (Chl) fluorescence levels of the 9th oldest leaves were measured using a Chl fluorescence measuring system (DUAL-PAM/F, Heinz Walz GmbH, Effeltrich, Germany). The fiber optics were connected to the leaf chamber of the gas exchange measurement system via an adaptor (6400-06, LI-COR Inc.). Chl fluorescence levels and gas exchange rates were measured simultaneously under the conditions described above. Before measurements were taken, leaves were maintained in the dark for at least 30 min. The steady-state and the maximum Chl fluorescence levels in the dark (F o and F m , respectively) were first measured, and subsequently those in the light (F and F m , respectively) were measured at actinic-light PPFDs of 250 and 1,200 µmol m −2 s −1 . The maximum photochemical quantum yield (F v /F m ) and effective photochemical quantum yield (Y II ) of PSII were calculated as (F m -F o )/F m (Kitajima and Butler, 1975) and (F m -F)/F m (Genty et al., 1989), respectively.

Temperature Measurements
In Exp. 1, 2, and 4, LT of the 8th or 9th oldest leaf of each plant was measured by continuously attaching the measuring junction of a calibrated type-K thermocouple (0.1 mm in diameter) to the leaf surface. AT in each growth chamber was measured using a calibrated type-T thermocouple. Data were collected and recorded every 30 s by a data logger (GL220, Graphtec Corp., Kanagawa, Japan). LT was evaluated as daily mean difference between LT and AT (LT -AT), which was calculated as the averaged difference between LT and AT over a 16-h photoperiod.

Microscopic Observation of Leaf Epidermis
For uninfiltrated plants and plants inoculated with the HA expression vector, replicas of the adaxial and abaxial epidermis of the 8th or 9th oldest leaves were created at 6 DPI using a liquid adhesive. The replicas were observed under a digital microscope (VHX-1000, Keyence Corp., Osaka, Japan) at 400times magnification.

Growth Analysis
Uninfiltrated plants and plants inoculated with the empty vector and the HA expression vector were destructively harvested at 6 days post-infiltration (DPI). Plants at 0 DPI were also harvested. Total leaves of each plant were oven-dried at 100 • C for 1 h followed by 80 • C for 3 days to determine dry mass (DM).

Hemagglutinin Quantification
In Exp. 4, the 8th or 9th oldest leaves of plants inoculated with the HA vector were destructively harvested at different DPI. Each leaf lamina was divided into two portions along the primary vein. One half of the leaf lamina was weighed to determine FM and subsequently subjected to HA quantification by sandwich enzyme-linked immunosorbent assay using a commercial kit (Sino Biological Inc., Beijing, China) as previously described by Fujiuchi et al. (2016b). The other half was weighed to determine FM and subsequently oven-dried as described above to determine DM. The ratio of DM to FM was used to convert HA content per unit FM into that per unit DM. Because the absolute HA levels tended to be different depending on the lot number of the kit, HA content was normalized against the maximum level in each AT treatment and replication for evaluation.

Statistical Analyses
In each experiment, individual plants were considered independent biological replicates. Numbers of replicates per treatment were as follows: 3 in Exp. 1 and 2, 3-4 in Exp. 3, and 4 for HA content and 6 for LT in Exp. 4, respectively. Significant differences between two means and among three means were tested by Welch's t test and Tukey-Kramer's HSD test, respectively, at P < 0.05 using statistical software (R 3.4.3, R Core Team, 2017).

Changes in Gas Exchange Rates and Leaf Temperature
In Exp. 1, P n decreased in leaves inoculated with the empty vector and HA vector from 2 DPI, whereas P n in uninfiltrated leaves was almost unchanged ( Figure 1A). P n in HA vectorinoculated leaves was significantly lower than that in empty vector-inoculated leaves at 4 and 5 DPI. Similar trends among treatments were observed in T and g s (Figures 1B,C). LT -AT in HA vector-inoculated leaves increased from 2 to 4 DPI and that in empty vector-inoculated leaves slightly increased from 2 to 3 DPI ( Figure 1D). LT -AT in those leaves thereafter gradually decreased. LT -AT in uninfiltrated leaves was almost constant between 2 and 8 DPI. LT -AT at 4-6 DPI in HA vector-and empty vector-inoculated leaves was approximately 1.2 and 0.4 • C higher, respectively, than that in uninfiltrated leaves.
To investigate the factor that generated the reduction of P n , T, and g s in leaves inoculated with the viral vectors, i.e., the empty and HA vectors, in Exp. 2, we compared mock-inoculated leaves, which were infiltrated with the suspension of A. tumefaciens not carrying any plasmid vectors, with empty vector-and HA vector-inoculated leaves. P n , T, and g s in mock-inoculated leaves did not markedly decrease throughout the experimental period, as did empty vector-and HA vector-inoculated leaves (Figures 2A-C). This indicates that the reduction of P n , T, and g s in leaves inoculated with the viral vectors did not result from the vacuum infiltration treatment with an A. tumefaciens suspension. In leaves inoculated with the empty and HA vectors, P n , T, and g s decreased from 2 to 3 or 4 DPI, being lower in HA vector-inoculated leaves than in empty vector-inoculated leaves   (Figures 2A-C), as was observed in Exp. 1. LT -AT in empty vector-and HA vector-inoculated leaves tended to be higher than that in mock-inoculated leaves at 3-5 DPI, although the difference was not significant in this experiment partly due to a relatively large variation among biological replicates in empty vector-inoculated leaves ( Figure 2D). g s is an index of the diffusivity of CO 2 and water vapor through the stomata and is largely influenced by stomatal density and stomatal opening. We made replicas of the adaxial and abaxial epidermis of uninfiltrated leaves and HA vectorinoculated leaves at 6 DPI and observed them with a digital microscope. There was no notable difference in stomatal density or the shape of guard cells between leaves (data not shown), suggesting that the lower g s in HA vector-inoculated leaves was due to impairment of stomatal opening.

Photosynthetic Characteristics
In Exp. 3, P n in pre-infiltrated leaves, uninfiltrated leaves, and HA vector-inoculated leaves was measured at −1, 7, and 6 DPI, respectively, under low (250 µmol m −2 s −1 ) and high (1,200 µmol m −2 s −1 ) PPFD at various C a levels. We previously found that leaf HA content reached the maximum level at approximately 6 DPI under the growth conditions similar to those in the present experiment (Matsuda et al., 2017a). Preinfiltrated and uninfiltrated leaves exhibited comparable C a -P n curves under both PPFD levels (Figures 3A,B). At any C a level, P n of HA vector-inoculated leaves was significantly lower than that of pre-infiltrated and uninfiltrated leaves. To evaluate P n excluding the effect of lower g s in HA vector-inoculated leaves (Figures 1C, 2C), P n was plotted against C i instead of C a (Figures 3C,D). P n at a given C i was lower in HA vectorinoculated leaves than in pre-infiltrated and uninfiltrated leaves under both PPFD conditions. Therefore, the lower P n in HA vector-inoculated leaves was not solely accounted for by the lower g s , and the downregulation of photosynthetic reactions in chloroplasts should occur. We also measured F v /F m and Y II at low and high PPFD, and they were all significantly lower in HA vector-inoculated leaves than in pre-infiltrated and uninfiltrated leaves ( Table 1). This indicates that photosynthetic electron transport capacity of PSII was downregulated in HA vectorinoculated leaves.

Plant Growth
We next confirmed the effect of reduction of leaf P n in leaves inoculated with the empty and HA vectors on plant growth. Plants that had been inoculated with either of the vectors or uninfiltrated were grown and harvested at 6 DPI, and leaf DM was compared ( Table 2). There was no significant difference in leaf DM among them, although their leaf DM was significantly greater than that pre-infiltration. Thus, the effect of leaf P n reduction due to the vector inoculation on plant growth was not clearly observed under the conditions of the present study.

Changes in Hemagglutinin Content and Leaf Temperature at Different Air Temperatures
In Exp. 4, changes in HA content and LT -AT in leaves were examined at different AT levels post-infiltration. At 21 • C AT, HA content increased from 5 to 6 DPI and thereafter decreased slightly ( Figure 4A) or was almost unchanged (Figure 4B).
At 23 • C AT, HA content increased from 4 to 5 DPI, attained high levels at 5 and 6 DPI, and decreased from 6 to 7 DPI (Figures 5A,B). At 26 • C AT, HA content increased from 3 to 4 DPI and decreased sharply from 4 to 6 DPI (Figures 6A,B). Thus, the greater the post-infiltration AT, the earlier the accumulation of HA in leaves and the timing at which HA content peaked. The LT -AT during the overall experimental period was higher as the AT became higher, both in uninfiltrated and HA vector-inoculated leaves (Figures 4-6C,D). In uninfiltrated . An asterisk ( * ) represents a significant difference between uninfiltrated leaves and HA vector-inoculated leaves at each DPI by Welch's t-test at P < 0.05. leaves, LT -AT was almost unchanged throughout the period irrespective of AT. In contrast, in HA vector-inoculated leaves, LT -AT first increased, attained its peak, and then gradually decreased for all AT levels. Furthermore, the DPI that was required for LT -AT to reach the maximum level appeared to decrease with increasing AT, i.e., it was approximately 5, 4, and 3-4 DPI at 21, 23, and 26 • C, respectively.

Gas Exchange Rates Are Reduced and Photosynthetic Capacity Is Downregulated in Leaves Transiently Overexpressing Hemagglutinin With the Viral Vector
Inoculation of N. benthamiana leaves with the HA expression vector significantly decreased P n and T from 2 to 5 DPI compared to uninfiltrated leaves (Figures 1A,B) and leaves that were vacuum infiltrated with the mock A. tumefaciens suspension and did not receive a T-DNA vector (Figures 2A,B). Inoculation of leaves with the empty vector not containing the sequence of HA gene decreased P n (Figures 1, 2A) and T (Figures 1, 2B) yet not as much as did inoculation with the HA expression vector. Thus, the viral vector inoculation itself decreased P n and T in leaves, and the HA expression further decreased them. Moreover, the difference in g s among treatments showed a similar trend to those observed in P n and T (Figures 1, 2), suggesting that the decreases in P n and T in empty vector-and HA vector-inoculated leaves were at least partly accounted for by the decrease in g s . Infection of pathogens including viruses generally leads to reduction of g s (Grimmer et al., 2012). Salicylic acid accumulation during the hypersensitive response is thought to be involved in the reduction of g s via stomatal closure (Chaerle et al., 1999). Such response to virus infection might at least partly explain what was observed in the viral vector-inoculated leaves during the present study. Hamorsky et al. (2015) reported that inoculating N. benthamiana leaves with the same empty vector as used in the present study and with vectors to express cholera toxin B subunit resulted in a greater than 100-fold increase in the expression of PR1a, a marker gene for the hypersensitive response, compared to uninfiltrated leaves.
According to a biochemical model of C 3 photosynthesis, the so-called FvCB model (Farquhar et al., 1980), under high PPFD conditions, P n measured at low and high C i is limited by the capacity of RuBisCO to assimilate CO 2 and the capacity of photosynthetic electron transport to regenerate ribulose-1,5-bisphosphate in the Calvin cycle, respectively. Under low PPFD conditions, P n is limited by the electron transport capacity including light harvesting. P n in HA vectorinoculated leaves was lower than that in uninfiltrated leaves irrespective of PPFD or C i (Figures 3C,D), suggesting that both capacities of CO 2 assimilation by RuBisCO and electron transport were downregulated in planta. In previous studies, the apparent reduction of RuBisCO content in leaves transiently overexpressing a foreign protein with viral vectors has been observed Matsuda et al., 2012). This might be related to the remarkable reduction of P n in HA vector-inoculated leaves. In addition, results of our fluorometric analysis (Table 1) suggest that components involved in the photosynthetic electron transport chain were also downregulated in HA vector-inoculated leaves. It is likely that transient HA overexpression with the viral vector downregulated overall photosynthetic machinery in HA vector-inoculated leaves. Thus, both the decrease in g s and the downregulation of photosynthesis in mesophyll cells should contribute to the decreased P n by the HA expression in N. benthamiana leaves.
Despite the significant decrease in P n in empty vector-and HA vector-inoculated plants, their leaf DM at 6 DPI was not significantly different from that of the uninfiltrated control ( Table 2). This was probably because the post-infiltration period was only for 6 days and therefore the integral of the decreased P n in the viral vector-inoculated leaves during that period was too small to appear as a reduction in leaf DM. However, it is possible that downregulation of photosynthesis negatively affects metabolisms in plant cells via a supply shortage of energy and substrates. Specification of such negative effects, if any, will be a future research subject.

Leaf Temperature Increases Prior to Hemagglutinin Accumulation
Leaf temperature increased post-infiltration in empty vector-and HA vector-inoculated leaves, whereas LT was almost constant in uninfiltrated leaves and leaves infiltrated with the mock suspension (Figures 1, 2D). It is likely that the decrease in T (Figures 1, 2B) accompanied by the decrease in g s (Figures 1,  2C) in vector-inoculated leaves primarily contributed to the increase in LT via the suppression of latent heat dissipation. Chaerle et al. (1999) observed an increase in LT due to stomatal Temperature is known as an environmental factor that significantly influences the accumulation of various recombinant proteins in plants in transient gene expression systems (for a review, see Fujiuchi et al., 2016a). When AT post-infiltration was maintained at 26, 23, and 21 • C, the highest levels of leaf HA content were observed at 4 DPI (Figures 6C,D), 5-6 DPI (Figures 5C,D), and 6 DPI and later (Figures 4C,D), respectively. We previously reported a similar temperaturedependent difference in the time course of leaf HA content (Matsuda et al., 2017a). The activity of RNA replicase of tobacco mosaic virus almost linearly increases with increasing temperature between 12 and 35 • C (White and Dawson, 1978), which might be a factor contributing to the early accumulation of HA at high AT. Peak LT was also attained earlier as postinfiltration AT increases (Figures 4-6C,D). We recently reported that mean LT rather than mean AT post-infiltration was better at accounting for temperature dependency of the leaf HA level at 6 DPI under various thermal conditions and identified the importance of LT post-infiltration (Matsuda et al., 2017b). Based on our results in the present study, the time course of LT as well as mean LT post-infiltration should be considered when evaluating the effect of temperature on recombinant protein accumulation.
The peak LT was attained 1-2 days earlier than the peak HA content irrespective of AT. This finding leads us to propose the possibility of using thermal information of leaves for the nondestructive estimation of an appropriate DPI of harvest at which time the leaf HA content reaches almost its maximum level at each AT. For example, the day at which time comparable LT to or lower LT than that on the day before was detected can be appropriate for harvesting to obtain a high HA yield. AT can vary spatiotemporally in a greenhouse or even in a fully contained, environmentally controlled facility used for post-infiltration plant cultivation. If LT is monitored with infrared thermography, for example, the optimum harvest date, which might vary depending on individual plants due to the AT variation, could be detected. Indeed, Buyel et al. (2016) proposed usefulness of a noninvasive contact-free method to determine thermal properties of leaves in plant-based biopharmaceutical protein production to meet a prerequisite for good manufacturing practice, i.e., monitoring plant growth, maturation, and batch-to-batch consistency. Although application of this concept to practical production systems remains to be demonstrated, the present study presents the possibility of introducing a non-invasive sensing technique into the management of the upstream plant production process in viral vector-based biopharmaceutical protein production systems.

CONCLUSION
Viral vector inoculation significantly reduced P n and T in N. benthamiana leaves, which was associated with the reduction of g s . Inoculation of leaves with the empty vector itself decreased P n , T, and g s to some extent and the HA expression further decreased them. In HA vector-inoculated leaves, downregulation of photosynthetic capacity including electron transport activity occurred, in addition to the impairment of stomatal opening. LT increased in viral vector-inoculated leaves, which was associated with a decrease in T. HA accumulated earlier in leaves as AT post-infiltration was increased, and the peak LT appeared to be attained 1-2 days earlier than the peak HA content irrespective of AT. This suggests a possibility to use LT information to determine harvest date of leaves, at which time leaf HA content was nearly maximum, without destructive measurements.

AUTHOR CONTRIBUTIONS
RM conceived and designed the study. AU and HN performed the experiments. All authors analyzed the data and reviewed the manuscript. RM wrote the manuscript.

FUNDING
This work was financially supported in part by JSPS KAKENHI (Grant Nos. 26712021 and 17H03893 to RM and 24658217 to KF).