# Physiological performance of transplastomic tobacco plants overexpressing aquaporin AQP1 in chloroplast membranes

Alicia Fernández-San Millán, Iker Aranjuelo, Cyril Douthe, Miquel Nadal, María Ancín, Luis Larraya, Inmaculada Farran, Jaume Flexas, Jon Veramendi

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Alicia Fernández-San Millán, Iker Aranjuelo, Cyril Douthe, Miquel Nadal, María Ancín, et al.. Physiological performance of transplastomic tobacco plants overexpressing aquaporin AQP1 in chloroplast membranes. Journal of Experimental Botany, 2018, 69 (15), pp.3661-3673. 10.1093/jxb/ery148 .
hal-05240234

HAL Id: hal-05240234

https://hal.inrae.fr/hal-05240234v1

Submitted on 4 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>
RESEARCH PAPER

# Physiological performance of transplastomic tobacco plants overexpressing aquaporin AQP1 in chloroplast membranes

Alicia Fernández-San Millán<fnref n="1" />, Iker Aranjuelo<fnref n="1" />, Cyril Douthe<fnref n="2" />, Miquel Nadal<fnref n="2" />, María Ancín<fnref n="1" />, Luis Larraya<fnref n="1" />, Inmaculada Farran<fnref n="1" />, Jaume Flexas<fnref n="2" />, and Jon Veramendi<fnref n="1" />,*

Received 19 December 2017; Editorial decision 12 April 2018; Accepted 12 April 2018

Editor: Susanne von Caemmerer, Australian National University, Australia

## Abstract

The leaf mesophyll CO₂ conductance and the concentration of CO₂ within the chloroplast are major factors affecting photosynthetic performance. Previous studies have shown that the aquaporin NtAQP1 (which localizes to the plasma membrane and chloroplast inner envelope membrane) is involved in CO₂ permeability in the chloroplast. Levels of NtAQP1 in plants genetically engineered to overexpress the protein correlated positively with leaf mesophyll CO₂ conductance and photosynthetic rate. In these studies, the nuclear transformation method used led to changes in NtAQP1 levels in the plasma membrane and the chloroplast inner envelope membrane. In the present work, NtAQP1 levels were increased up to 16-fold in the chloroplast membranes alone by the overexpression of *NtAQP1* from the plastid genome. Despite the high NtAQP1 levels achieved, transplastomic plants showed lower photosynthetic rates than wild-type plants. This result was associated with lower Rubisco maximum carboxylation rate and ribulose 1,5-bisphosphate regeneration. Transplastomic plants showed reduced mesophyll CO₂ conductance but no changes in chloroplast CO₂ concentration. The absence of differences in chloroplast CO₂ concentration was associated with the lower CO₂ fixation activity of the transplastomic plants. These findings suggest that non-functional pores of recombinant NtAQP1 may be produced in the chloroplast inner envelope membrane.

**Keywords:** Aquaporin, chloroplast envelope, CO₂ permeability, plastid transformation, protein targeting, tobacco.

## Introduction

It is predicted that future increases in the human population will require a 30% increase in crop yield rates (Edgerton, 2009). Improving the photosynthetic performance of crops is one way in which plant production might be increased (Parry et al., 2011; Reynolds et al., 2011; Parry et al., 2013; Flexas et al., 2013), and a number of strategies have been identified that, either individually or in combination, might achieve this (Long et al., 2006; Flexas et al., 2006; Parry

et al., 2013; Flexas et al., 2012). Photosynthetic performance is affected by two major factors: the concentration of CO₂ within the chloroplast and the efficiency of the carboxylation biochemistry. Availability of CO₂ at the carboxylation site in the chloroplast can be limited by its diffusion into the substomatal cavities, referred to as stomatal conductance (g_s), and by the conductance of CO₂ from the substomatal cavity to the chloroplast, referred to as mesophyll conductance (g_m).

<footnote n="1">¹ Instituto de Agrobiotecnología (Universidad Pública de Navarra-CSIC), Departamento de Producción Agraria, Campus Arrosadía, 31006 Pamplona, Spain</footnote>

<footnote n="2">² Research Group on Plant Biology under Mediterranean Conditions, Departament de Biologia, Universitat de les Illes Balears, Carretera de Valldemossa Km 7.5, 07122 Palma de Mallorca, Illes Balears, Spain</footnote>

<footnote>*Correspondence: jon@unavarra.es</footnote>
Classically, $g_m$ has been described not to limit photosynthesis, and the CO₂ concentration was thought to be similar in the substomatal cavity ($C_i$) and in the chloroplast stroma ($C_c$). However, over the past decade, a number of studies (Flexas et al., 2006; Scafaro et al., 2011; Kaldenhoff, 2012; Evans and von Caemmerer, 2013; Flexas et al., 2013) have shown that $g_m$ has a major influence on CO₂ diffusion into the chloroplast, with a consequent impact on the photosynthetic rate. At the cellular level, atmospheric CO₂ has to pass through the cell wall and three membranes (the plasma membrane and the two membranes of the chloroplast envelope) to reach the chloroplast stroma. The CO₂ permeability of the chloroplast envelope is low, probably due to its relatively large protein content (Priestley and Woolhouse, 1980); indeed, it was estimated that it may account for almost half of the internal leaf resistance to CO₂ (Uehlein et al., 2008). As a result, under light-saturated conditions, photosynthesis is limited by the availability of CO₂ within the chloroplast. Other studies have shown that the $g_m$ can change quickly in response to varying environmental conditions, such as leaf temperature (Bernacchi et al., 2002), water stress (Galmés et al., 2007), blue light (Loreto et al., 2009), and the external CO₂ concentration (Flexas et al., 2007a). This rapid modification of $g_m$ points to the existence of additional components, some of them probably proteins, controlling the conductance of the mesophyll to CO₂ diffusion. Proteins forming pore-like structures, such as aquaporins (AQPs), might help explain how these rapid variations in $g_m$ occur.

AQPs are small proteins that increase the permeability of cell membranes to water and certain small, neutral molecules, including CO₂ (Maurel et al., 2008; Gomes et al., 2009; Chaumont and Tyerman, 2014; Kaldenhoff et al., 2014; Groszmann et al., 2017). AQPs were discovered for the first time in plants in the vacuolar tonoplast of Arabidopsis (Maurel et al., 1993), and are present in the whole plant kingdom. AQPs are located in the plasma membrane and also in most of the intracellular membranes. Many isoforms of AQPs exist, which can be classified according to their sequence homologies and subcellular localization. The plasma membrane intrinsic protein (PIP) class includes isoforms that are most abundant in the plasma membrane. This class can be subdivided into subclasses PIP1 and PIP2 according to sequence similarity. Investigations on the mesophyll cells of tobacco leaves have shown that the plasma membrane protein NtAQP1 (a PIP1 member) facilitates CO₂ transport, and that it has important functions in photosynthesis and stomatal opening (Uehlein et al., 2003). Further studies revealed a dual localization of NtAQP1 in the plasma membrane and the inner envelope membrane (IEM) of the chloroplast (Uehlein et al., 2008). A mutation in the *Arabidopsis thaliana AtPIP1;2* gene was found to be associated with reduced $g_m$ and a reduction in the rate of photosynthesis (Heckwolf et al., 2011). Genetic engineering to modify NtAQP1 expression levels confirmed these results, revealing a function for NtAQP1 in CO₂ conductance. Antisense or RNA interference-mediated downregulation of *NtAQP1* resulted in a reduction of IEM CO₂ permeability, $C_c$ values, and photosynthetic performance (Uehlein et al., 2003; Flexas et al., 2006; Uehlein et al., 2008). Overexpression of NtAQP1

in tobacco and Arabidopsis, however, increased chloroplast membrane CO₂ permeability, the rate of photosynthesis, and plant growth (Aharon et al., 2003; Uehlein et al., 2003; Flexas et al., 2006). A similar positive effect on CO₂ permeability, plus an increase in leaf net photosynthesis, was observed in tomato and rice plants overexpressing AQP (Hanba et al., 2004; Sade et al., 2010). It was eventually suggested that the function of NtAQP1 might depend on its localization in the cell, and that it might provide a water channel in the plasma membrane and a CO₂ channel in the chloroplast envelope (Uehlein et al., 2008).

In the present study, the hypothesis that higher levels of AQPs in the chloroplast would increase CO₂ transport and the rate of photosynthesis was tested by overexpressing NtAQP1 from the chloroplast genome of tobacco. Compared with nuclear transformation, plastid transformation provides the advantage of high transgene expression levels (Bock, 2015). In addition, the recombinant protein is confined to the chloroplast, eliminating the effect of AQP1 modification in the plasma membrane. Therefore, the main objective of the present study was to evaluate the role of NtAQP1 overexpression specifically in the chloroplast membranes on CO₂ permeability and photosynthetic performance.

## Materials and methods

### Production of plants overexpressing NtAQP1 in the chloroplast

Total RNA from *Nicotiana tabacum* L. (cv. Petite Havana SR1) leaves was extracted using the Ultraspec RNA kit (Biotecx Laboratories, Houston, TX, USA), and cDNA was synthesized using the SuperScript III system (Invitrogen, Carlsbad, CA, USA). The NtAQP1 gene (GenBank Accession AJ001416) was amplified by PCR with the primers NTAQP1for: AAGCTTTTGCAAGTATATT TTCCATGGCAGAAAACAAAGAA GAAGATGTTAAGCTCGG and NTAQP1rev: GCGGCCGCTTAA GACGACTTG TGGAATGGAATGGCTCTG. The full-length cDNA was then cloned into the pGEMTeasy vector (Promega, Madison, WI, USA) and sequenced. The tobacco NtAQP1 gene was subsequently cloned into the pAF chloroplast transformation vector (Fernández-San Millán et al., 2008) under the control of the *psbA* promoter and 5'-untranslated region, to obtain the expression vector pAF-AQP1. The Tic40 transit peptide sequence (240 bp) was amplified by PCR using cDNA from *A. thaliana* using the primers AtTic40TPfor: CCATGGAGAACCTTACCCTAGTTTC and AtTic40TPrev: GCGGCCGCAAGCTTTGCTTCTCTGTTTC. It was then fused with NtAQP1 at an NcoI restriction site to produce the expression vector pAF-TicAQP1.

*N. tabacum* L. (cv. Petite Havana SR1) was also used in plastid transformations. The PDS-1000/He biolistic system (Bio-Rad, Hercules, CA, USA) was used for the integration of transgenes as previously described (Daniell, 1997). The *aadA* gene, conferring resistance to spectinomycin, was used as a selectable marker gene. Two rounds of selection and shoot development on RMOP medium containing 500 mg/l spectinomycin were performed. The transplastomic plants produced were named AQP1 and TicAQP1.

### Southern and northern blotting

Southern and northern blotting experiments were performed as previously described (Sanz-Barrio et al., 2013). For Southern blotting, the flanking sequence P1 probe generated by PCR was used. After Southern blot confirmation of the T₀ generation, selected plants were transplanted into pots and grown in the greenhouse for seed production. T₁ plants were
used in further experiments. Northern blotting was performed using the AQP1-specific P2 probe (515 bp) obtained by NcoI digestion of AQP1.

### Protein extraction, separation, and western blotting

Leaf samples (100 mg) from transformed and untransformed 70-day-old plants were ground in liquid nitrogen, homogenized in 300 µl of 2× Laemmli buffer (0.5 M Tris-HCl, pH 6.5, 4% SDS, 20% glycerol, and 10% β-mercaptoethanol) and heated at 95 °C for 5 min. After 5 min of centrifugation at 20 000 g, the supernatant was deemed to represent the total protein (TP) content. TP was quantified using the RC-DC protein assay (Bio-Rad) with BSA as a standard. Proteins were separated by SDS-PAGE on 12% polyacrylamide gels and transferred to a polyvinylidene fluoride (PVDF) membrane for immunoblotting. The primary antibodies used were anti-PIP1 (Agrisera, Vännäs, Sweden) and anti-NtAQP1 (kindly provided by R. Kaldenhoff) (dilution 1:3000). Peroxidase-conjugated goat anti-rabbit or anti-chicken immunoglobulin G (Sigma-Aldrich, St Louis, MO, USA) (both at a dilution of 1:3000) were used as secondary antibodies with the anti-PIP1 and anti-NtAQP1 primary antibodies, respectively. Detection was performed using the chemiluminescence ECL western blotting system (GE Healthcare, Fairfield, CT, USA).

Relative quantification of NtAQP1 monomers and oligomers was performed by comparing dilution series of TP from wild-type (WT) plants and both types of transplastomic plant (three replicates were analysed). For each line, adequate amounts were loaded on to an SDS-PAGE gel, electrophoretically separated, and then analysed by western blotting. Immunoblots were quantified using GeneTools Analyzer software (SynGene, Cambridge, UK).

### Plasma membrane isolation

Plasma membranes from 50-day-old WT and transplastomic plants grown in a growth chamber [16 h light/8 h dark; 200 µmol m⁻² s⁻¹ photosynthetic photon flux density (PPFD) and a day/night temperature regime of 28 °C/25 °C] were obtained as previously described (Santoni, 2007).

### Chloroplast isolation, fractionation, and immunoblotting

For chloroplast isolation, leaves from 50-day-old tobacco plants were cut into 1–3 cm² pieces and homogenized in a blender. The isolation buffer (330 mM sorbitol, 20 mM MOPS, 13 mM Tris, 3 mM MgCl₂, 0.1% BSA) was six times (v/w) the fresh mass weight of the leaf samples. The homogenate was passed through a filter mesh and centrifuged for 5 min at 1000 g and 4 °C. The pellet fraction was resuspended in isolation buffer and chloroplasts were isolated by 80–40% Percoll gradient fractionation after centrifugation for 10 min at 7700 g and 4 °C. Isolated chloroplasts were washed in 3 volumes of washing buffer (330 mM sorbitol, 50 mM HEPES/KOH, pH 7.6, 3 mM MgCl₂).

For fractionation, the chloroplasts were lysed by freeze-thawing in hypotonic TE buffer [10 mM Tris, 2 mM EDTA, pH 7.5, including a cocktail of protease inhibitors from Roche (Mannheim, Germany)]. Stroma, envelopes, and thylakoids were separated by using discontinuous sucrose gradients (0.93/0.6/0.3 M) after 2 h of centrifugation at 20 000 rpm in a swing-out rotor. Stroma was collected from the upper fractions and one volume of extracts was combined with one volume of 2× Laemmli buffer. The thylakoid membranes sedimented out and were resuspended in 10 mM TE buffer, to which one volume of 2× Laemmli buffer was added. The chloroplast envelopes were collected at the interface between 0.9 and 0.6 M sucrose. The envelope proteins were concentrated by methanol/chloroform extraction (Ferro et al., 2002) and resuspended in 10 mM TE buffer with one volume of 2× Laemmli buffer added. All three fractions were heated for 1 h at 37 °C. Protein concentrations were determined by the Bradford method. All samples were resolved by 9% SDS-PAGE and transferred to PVDF membranes for western blotting. Antisera to ADP-glucose pyrophosphorylase (AGPase), LHC chlorophyll a/b binding protein 1 (Lhcb1), and Tic40 (Agrisera) proteins were used at dilutions of 1:1000.

### Plants used for gas exchange

Plants were grown in 2 litre pots (organic soil/perlite, 70/30 v/v) in two places, Pamplona and Mallorca (Spain). In Pamplona, plants were grown in a greenhouse at 24–28 °C and relative humidity ~40%. Plants were watered with water by drip irrigation and twice a week with 50% diluted Hoagland's solution. In Mallorca, plants were grown in a growth chamber at PPFD ~350 µmol m⁻² s⁻¹ at the top of the plants, daily temperature of 24–26 °C, relative humidity ~40%, and watered twice a week with water and once a week with 50% diluted Hoagland's solution.

### Gas exchange and chlorophyll fluorescence analyses

Gas exchange measurements were performed with a calibrated Li-6400 XT portable gas analyser (Licor, Lincoln, NE, USA) equipped with the 2 cm² Li-6400-40 Leaf Chamber Fluorometer. Determinations were conducted in apical fully developed leaves. Three independent experiments (two in Pamplona and one in Mallorca) were performed. For the measurements made in Pamplona, plants were transferred to a growth chamber with similar environmental conditions to those of the Mallorca site. For each plant, the same procedure was followed: first, stabilization until a steady state of stomatal conductance was reached (typically ~20–30 minutes) in ambient conditions (CO₂ concentration=400 µmol mol⁻¹, 1500 PPFD, and 25 °C). After stabilization, the Aₙ/Cᵢ curve was performed by changing the concentration of CO₂ entering the leaf chamber with, the following steps: 400, 300, 250, 200, 150, 100, 50, 400, 400, 500, 600, 700, 800, 1000, 1200 and 1500 µmol mol⁻¹, with typically 2–3 minutes between each step. Each Aₙ/Cᵢ curve was corrected for leaks by following the protocol described by Flexas et al., (2007b). In all three experiments, the results of net CO₂ assimilation and stomatal conductance were very consistent (data not shown). In the third experiment, performed in Mallorca, chlorophyll fluorescence was measured together with gas exchange to estimate the mesophyll conductance to CO₂. Therefore, all the results shown in the present paper correspond to this latter experiment.

After performing the Aₙ/Cᵢ curve, the leaf was kept in the chamber and N₂ from a tank (Air Liquide) was piped into the Li-6400 inlet to remove O₂ from the entering air in the leaf chamber, to allow measurements to be made in non-photorespiratory conditions (Valentini et al., 1995). We then performed a light curve at ambient CO₂ concentration (400 µmol mol⁻¹) with the following PPFD steps: 1500, 2000, 1750, 1500, 1250, 1000, 800, 550, 300, 150, 100, 75, 50, 25 and 0 µmol m⁻² s⁻¹. These measurements were used to estimate the product of leaf absorption (α) and the partitioning of absorbed quanta between photosystems I and II (β) (see Valentini et al., 1995 and Pons et al., 2009 for details). We used only the first points of the curve, with PPFD >400 µmol m⁻² s⁻¹, to estimate αβ (Martins et al., 2013), avoiding non-linearity of ΦPSII versus ΦCO₂ due to changes/higher influence of leaf respiration at low PPFD. Values of (αβ) were 0.36 ± 0.04 (TicAQP1), 0.35 ± 0.03 (AQP1), and 0.38 ± 0.03 (WT), with no significant differences between genotypes (see Supplementary Fig. S1 at JXB online). Night respiration rate (Rₐ) was estimated by measuring leaf gas exchange in darkness, 1 h after the lights of the growing chamber were turned off (night).

Mesophyll conductance (gₘ) was estimated by the method developed by Harley et al. (1992), as follows:

$$g_m = A_N / (C_i - \{ \Gamma^* [ETR + 8(A_N + R_1)] / [ETR - 4(A_N + R_1)] \})$$

where Aₙ is the net CO₂ assimilation rate, Cᵢ is the CO₂ concentration in the subdomal cavity, Γ^* is the CO₂ compensation point in the absence of Rₐ (assumed to be 40 µmol mol⁻¹, from Walker et al., 2013), R₁ is the respiration rate in light (estimated as Rₐ/2), and ETR is the electron transport rate, estimated as follows:

$$ETR = \alpha \times \beta \times \Phi \text{PSII} \times \text{PPFD}$$

where ΦPSII is the yield of photosystem 2. ΦPSII was estimated using the 'Multiphase Flash' method described by Loriaux et al. (2013).
### Discrimination against ¹³CO₂

The ¹³C isotope discrimination (Δ, ‰) was calculated as:

$$ \Delta = \frac{\delta^{13}C_{\text{atm}} - \delta^{13}C_{\text{sample}}}{\delta^{13}C_{\text{sample}} + 1} $$

where δ¹³C_atm is the carbon isotope composition in atmospheric CO₂ in the greenhouse (−10.8‰) and δ¹³C_sample is the carbon isotope composition of leaf total organic matter (TOM).

The ¹³C/¹²C ratio (R) in plant material was expressed in δ notation (δ¹³C) with respect to Vienna Pee Dee Belemnite calcium carbonate (V-PDB), and measured with an analytical precision of 0.1‰:

$$ \delta^{13}C = \left( \frac{R_{\text{sample}}}{R_{\text{standard}}} \right) - 1 $$

δ¹³C accuracy was monitored using international secondary standards of known ¹³C/¹²C ratios (IAEA-CH7 polyethylene foil, IAEA-CH6 sucrose, and USGS-40 glutamic acid, IAEA, Austria).

TOM and gas δ¹³C determinations were conducted at the Serveis Científic-Tècnics of the University of Barcelona. For TOM analyses, 1 mg of dry ground leaf material was analysed using an elemental analyser (EA1108, Series 1, Carbo Erba Instrumentazione, Milan, Italy) coupled to an isotope ratio mass spectrometer (Delta C, Finnigan MAT, Bremen, Germany) operating in continuous flow mode. Air δ¹³C samples were analysed by gas chromatography (Agilent 6890 Gas Chromatograph, Agilent Technologies, Spain) coupled to an isotope ratio mass spectrometer Deltaplus via a GC-C Combustion III interphase (ThermoFinnigan, Thermo, Barcelona, Spain).

### Determination of Rubisco, starch, and chlorophyll contents

Samples from the same leaves used for gas exchange and chlorophyll fluorescence measurements were analysed for their Rubisco, starch, and chlorophyll contents. Three leaf discs (2.1 cm²) per plant from WT, AQP1, and TicAQP1 plants were frozen in liquid nitrogen and ground in a Mikrodismembrator (Braun, Melsungen, Germany). The volume of the extrac tion buffer (phosphate buffer, pH 7.0, 100 mM) was three times the fresh weight (v/w) of the powdered leaf sample obtained from the three leaf discs. Samples were mixed in a vortex and, after 15 min on ice, centrifuged for 5 min at 20000 g at 4 °C. Protein fractions recovered from the supernatants were quantified by the Bradford method. For the separation of these proteins, one volume of the protein fraction was combined with one volume of 2× Laemmli buffer, boiled for 5 min, and then centrifuged at 20000 g for 5 min. Samples (15 µg) of the proteins in these supernatants were separated by SDS-PAGE (10%) and the gels were stained with Coomassie brilliant blue G-250. The Rubisco levels of the transplastomic plants were compared with those of the WT plants using GeneTools Analyzer software (SynGene).

Starch was determined using an amyloglucosidase-based test kit (R-Biopharm AG, Darmstadt, Germany).

The leaf chlorophyll contents of the transplastomic and WT plants [measured with a SPAD 502 chlorophyll meter (Minolta Optics Inc, Tokyo, Japan)] were recorded in the same leaves used for photosynthetic rate measurements.

Leaf area of flowering plants grown in a growth chamber was determined after scanning with ImageJ.

### Electron microscopy

Samples from the same leaves used to measure the photosynthetic rate were fixed in Karnovsky fixative (4% formaldehyde and 5% glutaraldehyde in 0.025 M cacodylate buffer, pH 6.7) by vacuum infiltration and further prepared for examination by transmission electron microscopy at the Microscopy Service of the University of Navarre, Spain.

### Statistical analysis

One-way ANOVA was used to analyse differences in the measured variables between the control and transplastomic plants. Differences among

means were analysed using the Tukey test (P<0.05). All calculations were performed using SPSS 10.0 software.

## Results

### Generation of tobacco transplastomic plants and determination of homoplasmy

Tobacco plants expressing *NtAQP1* from the plastid genome were obtained by biolistic bombardment of the leaves with the engineered pAF vector (Fernández-San Millán et al., 2008), which inserted the transgenes between the *trnI* and *trnA* regions of the plastid genome (Fig. 1A). Two different transformation vectors were designed, both with the transgene controlled by the promoter and the 5'-untranslated region of the *psbA* gene. The pAF-AQP1 vector included the full coding sequence of *NtAQP1*. In the pAF-TicAQP1 vector, the 76 amino acid sequence transit peptide of *A. thaliana* Tic40 protein was fused to the N-terminus of NtAQP1. Two independent transplastomic lines for each construction, developed after two rounds of selection on spectinomycin, were analysed by Southern blotting (Fig. 1B). As predicted for the correct homologous recombination of the transgenes, the flanking region P1 probe hybridized to a 10.4 or 10.7 kb BamHI DNA fragment in the AQP1 and TicAQP1 transplastomic plants, respectively. A 7.1 kb band was detected only in the WT plants, indicating that all four lines were homoplasmic.

### Expression of aquaporin in the chloroplast

Analysis at the transcriptional level in the transplastomic plants revealed transcripts of the expected size in both cases (Fig. 2A). Monocistrons of 1.4 and 1.7 kb derived from the *psbA* promoter were detected in the AQP1 and TicAQP1 plants, respectively. Dicistrons transcribed from the upstream *rrn* promoter were also present. A greater abundance of transcripts was observed in the AQP1 plants than in the TicAQP1 plants, probably owing to the TicAQP1 transcripts being less stable. The endogenous AQP1 mRNA was below the detection limit.

The overexpression of AQP1 protein in the transformed chloroplasts was confirmed by immunoblotting (Fig. 2B). A faint 30 kDa band of the expected size for the AQP1 protein monomer was observed in the WT plants. A stronger signal of the same electrophoretic mobility was detected in the AQP1 and TicAQP1 plants. Thus, the TicAQP1 recombinant protein was correctly processed in the stroma of the chloroplast following cleavage of the Tic40 transit peptide. Higher molecular weight signals were mainly present in transplastomic plants, indicating the presence of abundant oligomeric structures despite the denaturing conditions used during sample preparation and electrophoresis. It is also possible that the increased AQP1 protein production or inadequate post-translational modifications caused misfolding or the formation of non-specific aggregates with other proteins that resulted in abnormal migration patterns in the gel. The putative non-specific protein aggregates of high molecular weight could indicate that only a proportion of the recombinant protein equates to functional AQP1 complexes. Relative AQP1 protein levels were estimated by densitometry of different protein extract dilutions in western
A

<figure><img src="image_3.png" /></figure>

<figure><img src="image_4.png" /></figure>

B

<figure><img src="image_5.png" /><figcaption>**Fig. 1.** Integration of *Nicotiana tabacum* AQP1 into the tobacco chloroplast genome. (A) Map of the wild-type (WT) and AQP1-transformed plastid (pt) genomes. The transgenes were targeted to the intergenic region between *trnL* and *trnA*. The selectable marker gene *aadA* (encoding aminoglycoside 3'-adenylyltransferase) was driven by the 16S ribosomal RNA operon promoter (*Prrm*). AQP1 was driven by the *psbA* promoter and 5'-untranslated region (*PpsbA*). Arrows within boxes show the direction of transcription. Numbers below each ptDNA indicate the predicted size of hybridizing fragments when total DNA was digested with *BamH* I. A 0.8 kb fragment of the targeting region for homologous recombination was used as a probe (P1) for Southern blot analysis. *TpsbA*, 3'-untranslated region of the *psbA* gene. (B) Southern blot analysis of two independent lines (1 and 2) for each transformation cassette.</figcaption></figure>

blots, analysing both monomeric and multimeric signals. The expression level of AQP1 protein in the AQP1 and TicAQP1 plants was approximately 10-fold and 16-fold greater, respectively, than in the WT plants.

As expected, analysis of purified plasma membranes by immunoblotting indicated that the AQP1 protein levels in the WT and transplastomic plants were similar (Fig. 2C), confirming that the AQP1 protein synthesized in the chloroplast was not exported out of the chloroplast.

*NtAQP1 was mainly localized to the chloroplast envelope*

The distribution of AQP1 in the chloroplast of the AQP1 and TicAQP1 transplastomic plants was examined by chloroplast purification and suborganellar fractionation followed by immunoblotting (Fig. 3). The equal loading of gels for each fraction

and the purity of the fractions were assessed using specific anti-bodies: anti-ADP-glucose pyrophosphorylase (AGPase) for the stroma (ap Rees, 1995), anti-Tic40 for the envelope, and anti-LHC chlorophyll a/b binding protein 1 (Lhcb1) for the thylakoid membrane (Farmaki et al., 2007). AQP1 was detected in the envelope and thylakoid membrane fractions but no signal was seen for the stroma soluble fraction (Fig. 3). As expected, a stronger signal was detected in the AQP1 and TicAQP1 transplastomic plants than in the WT plants. Monomeric and oligomeric forms were detected in the envelope and thylakoid fractions of both types of transplastomic plants. The AQP1 signal was most intense in the envelope fraction; note that 10-fold more thylakoid fraction protein was loaded to enable the detection of AQP1. Monomers and dimers were present mainly in the envelope fraction, with faint signals detected in the thylakoid fraction. Trimers and tetramers were detected in the envelope fraction but not in the thylakoid fraction.
<figure><img src="image_6.png" /></figure>

<figure><img src="image_7.png" /></figure>

<figure><img src="image_8.png" /><figcaption>**Fig. 2.** Analysis of *AQP1* expression in wild-type (WT) and AQP1 and TicAQP1 transplastomic plants. (A) Northern blot analysis of leaf samples. The expected transcript sizes of the mono- and dicistrons originating from different promoters are indicated below the map of the transformed plastid genome. The 515 bp *AQP1* sequence (P2) was used as a probe. A 10 µg aliquot of total RNA was loaded per well. Ethidium bromide-stained rRNA was used to assess loading. (B) Western blot analysis of total protein from leaf samples (two independent lines for each construction). The lower panel was overexposed to show the 30 kDa AQP1 monomer, which was not detected in the upper panel. A 30 µg aliquot of protein was loaded per well. (C) Western blot analysis of proteins extracted from the plasma membrane. A 3 µg aliquot of protein was loaded per well. The positions and sizes of molecular weight protein standards are indicated. The blots in B and C were detected using anti-NtAQP1 as the primary antibody.</figcaption></figure>

A higher proportion of high molecular weight aggregates was observed in the thylakoid fraction of TicAQP1 plants relative to the envelope fraction. Immunoblotting could not provide an accurate estimate of the relative distribution of NtAQP1 in each membrane fraction owing to the disproportionate-ate method by which each membrane fraction was purified. However, the higher *NtAQP1* expression level in TicAQP1 transplastomic plants correlated with the higher AQP1 content

in the thylakoid and envelope membranes of these plants relative to AQP1 transplastomic plants (Fig. 3).

## Photosynthetic performance and protein and starch metabolism were impaired by chloroplast NtAQP1 overexpression

The transplastomic plants, particularly TicAQP1 plants, showed retarded growth in comparison to the WT plants during the first 3 weeks following transplantation into pots, but thereafter they caught up and reached a similar height under standard growth conditions (data not shown).

Net photosynthesis (A_N) analyses showed both the AQP1 and the TicAQP1 plants to have lower CO₂ fixation rates than the WT plants (Fig. 4A), with no differences between the transplastomic plants. As expected, no significant differences were observed between genotypes for leaf stomatal conductance (g_s) (Fig. 4D). In contrast, g_m was diminished by 50% in AQP1 and TicAQP1 plants relative to WT plants (Fig. 4E). A_N/C_i curve determinations highlighted that the Rubisco maximum carboxylation capacity (V_Cmax) and the maximum electron transport rate contributing to ribulose 1,5-bisphosphate regeneration (J_max) values measured in both transplastomic plants were lower than those in WT plants (Fig. 4B, C). Gas exchange analyses also showed that the transplastomic plants had a higher C_i value (Fig. 5A), while no significant differences in C_c were detected (Fig. 5B). The C_c/C_i ratio was not altered in the transplastomic plants relative to the WT plants (Fig. 5C). In contrast, compared with the WT plants, AQP1 and TicAQP1 plants showed higher ¹³CO₂ discrimination (Δ) values (Fig. 4F). Transplastomic plants had lower soluble protein levels than their WT counterparts (Table 1). Moreover, the Rubisco levels were strongly reduced in the AQP1 and TicAQP1 plants (reductions of 23% and 41%, respectively). A 4-fold reduction in the leaf starch content, and reduced chlorophyll levels, were detected in the transplastomic plants relative to the WT plants (Table 1). In addition, the leaf mass area was reduced in the transplastomic plants (12.5% in AQP1 and 15.9% in TicAQP1 plants).

The ultrastructure of the mesophyll cell chloroplasts was analysed by transmission electron microscopy. Major differences were observed between the WT and transplastomic plants (Fig. 6). The WT plants showed the normal architecture of the thylakoid network, arranged in grana and lamellae (Fig. 6A), while the AQP1 and TicAQP1 plants showed abnormal granal stacking with a reduced number of appressed thylakoids (Fig. 6B, C). In comparison to the normal granal structure of the WT plants (Fig. 6D), the transplastomic plants (especially the TicAQP1 plants) also showed defective grana and swelling of the thylakoid lumen (Fig. 6E, F). Large protein aggregates detected by western blot (Figs 2B and 3), particularly in TicAQP1 plants, could have affected the thylakoid membrane integrity, resulting in abnormal thylakoid architecture. No differences were observed in the envelope membranes (Fig. 6G–I).

The integrity of the thylakoid membranes was assessed by treating isolated chloroplasts with SDS, a product commonly used in cell permeation assays in bacteria (Griffith and Wolf,
<figure><img src="image_9.png" /><figcaption>**Fig. 3.** Localization of AQP1 in the thylakoid and envelope membranes. Envelope, thylakoid, and stroma fractions were isolated from wild-type (WT), AQP1, and TicAQP1 leaves, and separated by SDS-PAGE. Samples of 2, 20, and 30 µg of protein from the envelope, thylakoid, and stroma, respectively, were loaded per well. Representative western blots performed with antibodies to AQP1, the inner-membrane Tic40 protein, the thylakoid membrane-specific LHC chlorophyll a/b binding protein 1 (Lhcb1), and the stroma-specific ADP-glucose pyrophosphorylase (AGPase) are shown. Asterisks indicate the positions of monomer (*), dimer (**), trimer (***), and tetramer (****) AQP1.</figcaption></figure>

For the chloroplasts of the AQP1 and especially the TicAQP1 plants, total chlorophyll solubilization was obtained with lower SDS concentrations than those required to achieve the same effect with WT plant chloroplasts (Supplementary Fig. S2A, C). The same pattern was observed for the solubilization of AQP1 and the thylakoid membrane Lhcb1 proteins, but not for the stromal AGPase or thylakoid-lumen-associated TL29 proteins (Supplementary Fig. S2B).

## Discussion

### Targeting of recombinant AQP1 to chloroplast membranes

The present study shows that NtAQP1 can be overexpressed from the plastid genome and that it localizes to the chloroplast membranes. Native NtAQP1 localizes to both the plasma membrane and the chloroplast IEM. No transit peptide for its chloroplast targeting has been identified, and the sorting mechanisms responsible for this dual localization are unknown (Luu and Maurel, 2013). The expression of NtAQP1 from the plastid genome results in the AQP1 polypeptide being synthesized in the stroma and not in the cytosol as in WT plants, conditioning its subsequent import pathway to the IEM. Hence, in addition to the AQP1 coding sequence, and given the uncertainty that the recombinant AQP1 from the chloroplast stroma would reach its target location correctly, a second construct with the transit peptide of the IEM Tic40 protein fused to AQP1 was prepared for plastid transformation. Tic40 is an integral IEM protein involved in protein translocation across this membrane (Chou et al., 2003) that follows the post-import pathway for

IEM targeting (Li and Schnell, 2006). Both constructs resulted in the incorporation of AQP1 in the chloroplast envelope and the thylakoid membranes. This result shows that the topology-determining sequence information within NtAQP1 is sufficient for its integration into the envelope membranes from the stroma (Fig. 3). Chloroplast IEM proteins that follow the post-import pathway are first imported from the cytoplasm into the chloroplast stroma in the form of a soluble, processed, intermediate product, and subsequently reinserted from the stroma into the IEM. NtAQP1 expressed from the plastid genome could putatively use the second part of this pathway from the stroma to the IEM. It has been shown that conserved proline residues in the N-terminal region of IEM proteins such as Tic40 and Tic110 are required for stromal reinsertion (Chiu and Li, 2008). The six proline residues present at positions 21, 35, 36, 37, 39, and 43 of NtAQP1 suggest that it uses a common import mechanism from the stroma.

Native IEM proteins utilize two different pathways for their targeting. It is unknown whether the native AQP1, with six membrane-spanning alpha helices, uses the post-import pathway or the stop-transfer pathway. Very little is known about the insertion of polytopic IEM proteins. For instance, the Cor413im1 membrane protein, with five or six transmembrane domains, is incorporated in the IEM via the stop-transfer pathway (Okawa et al., 2014), while Tic110, which has six transmembrane domains, utilizes the post-import pathway (Lübeck et al., 1997; Li and Schnell, 2006). If native and plastidial AQP1 use different import pathways, this could result in an improper location within the IEM that negatively affects its functionality. Tic40, a component of the TIC complex, is involved in the reinsertion process for proteins that use the post-import pathway.
<figure><img src="image_10.png" /><figcaption>**Fig. 4.** (A) Net photosynthesis, (A~N~), (B) maximum carboxylation velocity of Rubisco (V~Cmax~), (C) maximum electron transport rate contributing to ribulose 1,5-bisphosphate regeneration (J~max~), (D) stomatal conductance (g~s~), (E) mesophyll conductance (g~m~), and (F) ¹³C isotope discrimination (Δ) of wild-type (WT) and AQP1 and TicAQP1 transplastomic plants. Representative data from two independent experiments are presented. Values are means ±SE (n=7). Different letters indicate significantly different values (ANOVA, P<0.05).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 5.** (A) Intercellular CO₂ concentration (C~i~), (B) chloroplast CO₂ concentration (C~c~), (C) and C~c~/C~i~ ratio of wild-type (WT) and AQP1 and TicAQP1 transplastomic plants. Representative data from two independent experiments are presented. Values are means ±SE (n=7). Different letters indicate statistically different values (ANOVA, P<0.05).</figcaption></figure>

pathway (Chiu and Li, 2008). It could be expected that the recombinant TicAQP1 protein, which includes the Tic40 transit peptide, is inserted close to and potentially interacts with the

TIC complex. It remains to be elucidated whether the IEM import pathways of recombinant TicAQP1 and AQP1 proteins are the same, but the physiological performance of both
transplastomic plants was similar. The putative drawback related to IEM import seems to be equivalent in both transplastomic plants, irrespective of the presence of the Tic40 transit peptide.

However, integration of NtAQP1 into the thylakoid membranes was unexpected, and the mechanism of protein sorting remains unknown. Plastid transformation has also allowed the successful integration of other foreign proteins, with or without signal peptides, to the thylakoid membranes (Henig et al., 2007; De Marchis et al., 2011; Ahmad et al., 2012; Shanmugabalaji et al., 2013; Scotti et al., 2015), indicating that different import mechanisms might be used. The expression of the *Synechococcus BicA* bicarbonate transporter in tobacco plastids unexpectedly resulted in dual targeting of the protein to the thylakoid membranes and, in a smaller proportion, to the chloroplast envelope (Pengelly et al., 2014). A model for contact zones between

<table><caption>**Table 1.** Biochemical variables measured in young leaves from wild-type (WT) and transplastomic plants grown in a growth chamber</caption><thead><tr><th></th><th>WT</th><th>AQ1</th><th>TicAQP1</th></tr></thead><tbody><tr><td>Starch (µmol Glu g FW<sup>-1</sup>)</td><td>21.7 ± 2.9<sup>a</sup></td><td>5.5 ± 0.2<sup>b</sup></td><td>5.7 ± 0.4<sup>b</sup></td></tr><tr><td>Soluble protein (mg g FW<sup>-1</sup>)</td><td>16.1 ± 0.7<sup>a</sup></td><td>13.3 ± 1.0<sup>ab</sup></td><td>12.3 ± 0.4<sup>b</sup></td></tr><tr><td>Insoluble protein (mg g FW<sup>-1</sup>)</td><td>13.0 ± 1.7<sup>a</sup></td><td>10.5 ± 1.3<sup>a</sup></td><td>10.2 ± 2.4<sup>a</sup></td></tr><tr><td>Rubisco (% relative to WT)</td><td>100 ± 3.4<sup>a</sup></td><td>76.7 ± 3.9<sup>b</sup></td><td>58.8 ± 5.4<sup>c</sup></td></tr><tr><td>Chlorophyll content (SPAD)</td><td>42.9 ± 1.0<sup>a</sup></td><td>34.2 ± 0.9<sup>b</sup></td><td>30.0 ± 0.6<sup>c</sup></td></tr></tbody></table>

Values are means±SE (n=5-7).

Different superscript letters denote significant differences (ANOVA, P<0.05).

The chlorophyll content was measured using a Minolta SPAD 502 chlorophyll meter.

plasma and thylakoid membranes, allowing protein trafficking in short-lived connection assemblies, has been proposed for the cyanobacterium *Synechocystis* (Pisareva et al., 2011). In addition, functional thylakoid membranes were developed in association with the chloroplast envelope in *Chlamydomonas* under certain conditions (Hoober et al., 1991). These and other investigations have suggested a role of the IEM for thylakoid biogenesis in vascular plants (Celedon and Cline, 2013). This mechanism might tentatively explain the dual targeting of recombinant AQP1 to the envelope and thylakoid membranes.

Physiological performance of transplastomic plants

The present study sought to determine whether CO₂ transport to the chloroplast could be boosted by increasing the amount of AQP1 in the chloroplast membranes. Much higher AQP1 protein levels (up to 16-fold higher than in the WT) were obtained in this study than in another study using nuclear transformation, in which double the levels in WT were obtained (Flexas et al., 2006); this difference was probably due to the plastid transformation method. Despite the integration of AQP1 into the chloroplast envelope membranes, the transplastomic plants overexpressing NtAQP1 showed lower photosynthetic rates than the WT plants. Associated with the low values of A<sub>N</sub>, transplastomic plants showed both a reduction in CO₂ diffusion capacity (associated with g<sub>m</sub> but not g<sub>s</sub>) and a lower photosynthetic capacity (V<sub>Cmax</sub> and J<sub>max</sub>). Because the different lines differed in their C<sub>i</sub>, and g<sub>m</sub> responds to C<sub>i</sub> (Flexas et al., 2007a), the observed differences in g<sub>m</sub> could be attributable to differences in C<sub>i</sub>. However, g<sub>m</sub>/C<sub>i</sub> ratios (obtained from

<figure><img src="image_12.png" /><figcaption>**Fig. 6.** Changes in chloroplast ultrastructure due to AQP1 overexpression from the plastid genome. Transmission electron microscopic images of chloroplasts from (A) wild-type (WT), (B) AQP1, and (C) TicAQP1 plants. (D–F) Detail of the thylakoids from (D) WT, (E) AQP1, and (F) TicAQP1 plants. (G–I) Detail of the chloroplast envelope (delineated by two arrowheads) in (G) WT, (H) AQP1, and (I) TicAQP1 plants. C, cytosol; s, starch granule; st, stroma; v, vacuole. Scale bar=1 µm.</figcaption></figure>
data shown in Figs 4E and 5A) were 0.3 × 10⁻³ for the two
transplastomic lines, compared with ~0.7 × 10⁻³ for the WT,
suggesting that $g_m$ is reduced in the transplastomic lines regard-
less of their $C_i$. Considering that between 200 and 400 µmol
CO₂ mol⁻¹ air $g_m$ decreases at an approximate rate of 0.1% per
µmol CO₂ mol⁻¹ air (shown for different species, including
tobacco, by Flexas et al., 2008), if WT plants had had the same
$C_i$ as the transplastomic plants, their $g_m$ would have decreased
by ~5%, that is, it would have still been significantly higher.

The realized photosynthesis achieved by a plant can be
manipulated by two means (Gago et al., 2014; Flexas et al.
2016): either by modifying the photosynthetic capacity of
the plant (i.e. changing the rate of photosynthesis for a given
substrate availability) or by changing the diffusion capacity of
the leaf (which, in turn, will modify the quantity of substrate
available for photosynthesis). In this study, increasing AQP1
expression dramatically decreased the photosynthesis of the
transplastomic plants. The reduced A_N was an unexpected
result, especially given the role of AQP1 in CO2 diffusion, as
described in tobacco (Uehlein et al., 2003; Flexas et al., 2006)
and in other species (Hanba et al., 2004; Sade et al., 2010).
In a previous study, Flexas et al. (2006) showed that, com-
pared with the corresponding WT plants, photosynthetic
rates were lower in NtAQP1-deficient plants and higher in
NtAQP1-overexpressing plants, suggesting that variations
in the photosynthetic rate were certainly linked to changes
in C_c (Flexas et al., 2006). In the present study, the higher
C_i detected in the transplastomic (compared with the WT)
plants indicates that stomatal opening was not involved in
the lower A_N of the AQP1 and TicAQP1 plants. Indeed, the
same stomatal conductance in the three genotypes, along
with a clearly lower CO2 fixation by the photosynthetic
machinery in the transplastomic plants, might be related to
the higher C_i in AQP1 and TicAQP1 plants. Following from
this observation, it could be tentatively concluded that the
lower A_N in the transplastomic plants was caused by their
lower $g_m$ (by ~50%) compared with the WT plants, since $g_m$
is now recognized to play a major role in CO2 diffusion into
the chloroplast (Flexas et al., 2006; Flexas et al., 2007b; Scafaro
et al., 2011; Kaldenhoff, 2012; Evans and von Caemmerer,
2013; Flexas et al., 2013). Nevertheless, the hypothesis on
limited diffusion (reduced $g_m$) in AQP1 and TicAQP1 plants
cannot explain alone their lower photosynthetic rates. In
fact, $C_c$ was not different between WT and transplastomic
plants, and nor was the $C_c/C_i$ ratio, despite very different
rates of photosynthesis (Figs 4 and 5). Although $g_s$ was not
affected by the overexpression of AQP1 in transplastomic
plants, the overall CO2 supply was reduced as a consequence
of the reduction of $g_m$. However, because the CO2 demand
was also reduced due to impaired $V_{Cmax}$ and $J_{max}$, the overall
result was a relatively higher $C_c$ in the transplastomic plants
compared with WT plants; the difference was more evident
according to the Δ measurements (Fig. 4F) than according
to chlorophyll fluorescence-based estimates (Fig. 5B). If $C_c$
(i.e. the substrate for photosynthesis) was the same between
the WT, AQP1, and TicAQP1 plants, then differences in A_N
were more likely related to differences in photosynthetic
capacity.

Evidence that reduced $V_{Cmax}$ and $J_{max}$ are the true factors responsible for decreased $A_N$ in the two transplastomic lines arises from reverse photosynthesis modelling and from carbon isotope discrimination. Using reverse modelling, it is possible to estimate how large $C_c$ and $g_m$ should be for the transplastomic lines to reach WT $A_N$ if their $V_{Cmax}$ and $J_{max}$ is reduced. It turns out from this simulation that $C_c$ should increase from the estimated ~230 µmol mol⁻¹ to ~330 µmol mol⁻¹, which would require $g_m$ values of 8.8 and 4.3 µmol CO₂ m⁻² s⁻¹ for AQP1 and TicAQP1, respectively; these are unrealistically large values, far out of the range of estimates for any plant species (Flexas et al., 2012). On the other hand, stable isotopes, such as δ¹³C, have been proposed as indicators of stomatal opening and CO₂ diffusion (Farquhar et al., 1989; Araus et al., 2003). Changes in Δ have been linked to changes in CO₂ availability and/or Rubisco carboxylation activity (Farquhar et al. 1989; Brugnoli and Farquhar, 2000; Ghashghaie et al., 2003). The fact that CO₂ supply was reduced in transplastomic plants due to the reduction of $g_m$, but Δ was still larger in these plants than in WT plants, suggests a reduction of carboxylase activity relative to CO₂ supply around Rubisco.

The localization of AQP1 to the thylakoid membranes
might negatively affect their functionality, perhaps via inter-
action with proteins of the photosynthetic apparatus that
affects the protein dynamics within the thylakoid, with
a consequent impact on photosynthetic capacity. Indeed,
the reduced membrane resistance to SDS (Supplementary
Fig. S2) suggests that the thylakoid membrane was damaged
in some way, possibly related to the presence of recombinant
AQP1. It has been reported that targeting foreign membrane
proteins to thylakoid membranes by chloroplast transforma-
tion causes mutant phenotypes with reduced growth and
photosynthetic capacity, altered thylakoid ultrastructure, and
impairment of the integrity of the photosystems (Henig et al.,
2007; Gnanasekaran et al., 2016). A similar response could
be produced by the recombinant AQP1 located in the thy-
lakoid membranes, particularly considering the presence of
large protein aggregates (Fig. 3). However, despite the higher
content of recombinant AQP1 in the thylakoid membranes
of TicAQP1 plants, the photosynthetic parameters were simi-
lar in both AQP1 and TicAQP1 plants (Fig. 4). Perhaps a
threshold of chloroplast damage was reached in the AQP1
plants and the greater content of the recombinant protein in
TicAQP1 plants was inconsequential. Another possibility is
that other processes, in addition to the functionality of the
thylakoids, might be involved in the photosynthetic impair-
ment of these plants.

Other factors affecting recombinant AQP1 functionality

The gas exchange data highlight the fact that, in contrast to
what was expected, *NtAQP1* overexpression from the chloro-
plast genome constrains CO₂ diffusion from the substomatal
cavity to the chloroplast. An explanation for this harmful effect
is not easy to provide, but other factors, in addition to the
above-mentioned factors related to the import pathway and
thylakoid targeting, could be involved.
The regulation of the function of AQP proteins depends on several processes, including post-translational modifications and protein interactions, that affect both their activity and their subcellular localization (Hove and Bhave, 2011; Chaumont and Tyerman, 2014; Verdoucq et al., 2014). It is possible that differences between the stromal and cytosolic environments prevent the required AQP post-translational modifications and/or interactions in the stroma.

Another important aspect to be considered is that AQP5 assemble as homo- and/or heterotetramers in the membranes. The AQP monomer is the functional unit for water transport, but the tetramer and its composition may be important for CO2-related transport. Further, in tobacco, CO2 diffusion is greater when tetramers consist only of NtAQP1 from the PIP1 family (Otto et al., 2010), and different proportions of PIP1 and PIP2 subunits in the tetramer may promote either water or CO2 transfer, or both (Flexas et al., 2012). If improper AQP1 monomers (due to incorrect post-translational modifications) synthesized in the chloroplast stroma interact with normal WT AQP1 synthesized in the cytosol, non-functional or partially functional homotetramers might find their way to the IEM. It must be noted that unusual oligomeric and possible non-specific protein aggregates were detected (Figs 2B and 3). Given that the transplastomic plants in this study had 10–16 times the AQP1 content of the WT plants, defective homotetramers might have prevailed over any homotetramers composed exclusively of WT AQP1 monomers, thus compromising AQP1-mediated CO2 transport in the IEM.

Finally, despite the engineered modification to AQP1 expression levels being able to alter $g_m$ (suggesting a role for AQP5 in CO2 transport), the molecular bases controlling $g_m$ remain unclear and are affected by components such as the cell walls, plasma and chloroplastic membranes, and carbonic anhydrases. Thus, other regulatory mechanisms that control $g_m$ might be unmasked in plants with altered levels of AQP5 (Verdoucq et al., 2014), such as those analysed in the present study.

In conclusion, in an attempt to increase CO2 diffusion across the chloroplast membranes we engineered tobacco plants with substantially increased levels of the CO2-permeable NtAQP1 protein within the chloroplast membranes. However, there was no improvement in intracellular CO2 transfer capacity in the transplastomic plants. In fact, we observed impairment in photosynthetic capacity, which could be partially attributed to changes in the thylakoid ultrastructure, and reductions in $V_{Cmax}$ and $J_{max}$, which were the main limiting factors. This study serves to highlight obstacles that need to be overcome in work towards engineering improved CO2 transfer capacity for enhanced photosynthesis.

## Supplementary data

Supplementary data are available at JXB online.

Fig. S1. Relationship between $\Phi$PSII versus $\Phi$CO2 used to correct ETR.

Fig. S2. Membrane disruption of chloroplasts from wild-type (WT), AQP1, and TicAQP1 plants with increasing concentrations of SDS.

## Acknowledgements

The authors thank Ralf Kaldenhoff (Darmstadt University of Technology, Germany) for providing the NtAQP1 antibody and Iván Jáuregui for technical assistance. MA is a recipient of a Formación de Profesorado Universitario fellowship from the Ministerio de Educación, Cultura y Deporte, Spain. MN was supported by a predoctoral fellowship BES-2015-072578 from the Ministerio de Economía y Competitividad (MINECO), Spain, co-financed by the European Science Foundation.

## References

Aharon R, Shahak Y, Wininger S, Bendov R, Kapulnik Y, Galili G. 2003. Overexpression of a plasma membrane aquaporin in transgenic tobacco improves plant vigor under favorable growth conditions but not under drought or salt stress. The Plant Cell **15**, 439–447.

Ahmad N, Michoux F, McCarthy J, Nixon PJ. 2012. Expression of the affinity tags, glutathione-S-transferase and maltose-binding protein, in tobacco chloroplasts. Planta **235**, 863–871.

ap Rees T. 1995. Where do plants make ADP-Glc? In: Pontis HG, Salerno G, Echeverría E, eds. Sucrose metabolism, biochemistry, physiology and molecular biology. Rockville: American Society of Plant Physiologists, 143–155.

Araus JL, Villegas D, Aparicio N, del Moral L, Hani El S, Rharrabti Y, Ferrio JP, Royo C. 2003. Environmental factors determining carbon isotope discrimination and yield in durum wheat under Mediterranean conditions. Crop Science **43**, 170–180.

Bernacchi CJ, Portis AR, Nakano H, von Caemmerer S, Long SP. 2002. Temperature response of mesophyll conductance. Implications for the determination of Rubisco enzyme kinetics and for limitations to photosynthesis in vivo. Plant Physiology **130**, 1992–1998.

Bock R. 2015. Engineering plastid genomes: methods, tools, and applications in basic research and biotechnology. Annual Review of Plant Biology **66**, 211–241.

Brugnoli E, Farquhar GD. 2000. Photosynthetic fractionation of carbon isotopes. In: Leegood RC, Sharkey TD, von Caemmerer S, eds. Photosynthesis: physiology and metabolism. Dordrecht: Kluwer Academic Publishers, 399–434.

Celedon JM, Cline K. 2013. Intra-plastid protein trafficking: how plant cells adapted prokaryotic mechanisms to the eukaryotic condition. Biochimica et Biophysica Acta **1833**, 341–351.

Chaumont F, Tyerman SD. 2014. Aquaporins: highly regulated channels controlling plant water relations. Plant Physiology **164**, 1600–1618.

Chiu CC, Li HM. 2008. Tic40 is important for reinsertion of proteins from the chloroplast stroma into the inner membrane. The Plant Journal **56**, 793–801.

Chou ML, Fitzpatrick LM, Tu SL, Budziszewski G, Potter-Lewis S, Akita M, Levin JZ, Keegstra K, Li HM. 2003. Tic40, a membrane-anchored co-chaperone homolog in the chloroplast protein translocon. The EMBO Journal **22**, 2970–2980.

Daniell H. 1997. Transformation and foreign gene expression in plants by microprojectile bombardment. Methods in Molecular Biology **62**, 463–489.

De Marchis F, Pompa A, Mannucci R, Morosinotto T, Bellucci M. 2011. A plant secretory signal peptide targets plastome-encoded recombinant proteins to the thylakoid membrane. Plant Molecular Biology **76**, 427–441.

Edgerton MD. 2009. Increasing crop productivity to meet global needs for feed, food, and fuel. Plant Physiology **149**, 7–13.

Evans JR, von Caemmerer S. 2013. Temperature response of carbon isotope discrimination and mesophyll conductance in tobacco. Plant, Cell & Environment **36**, 745–756.

Farmaki T, Sanmartín M, Jiménez P, Paneque M, Sanz C, Vancanneyt G, León J, Sánchez-Serrano JJ. 2007. Differential distribution of the lipoxygenase pathway enzymes within potato chloroplasts. Journal of Experimental Botany **58**, 555–568.

Farquhar GD, Ehleringer JR, Hubick KT. 1989. Carbon isotope discrimination and photosynthesis. Annual Review of Plant Physiology and Plant Molecular Biology **40**, 503–537.

Fernández-San Millán A, Ortigosa SM, Hervás-Stubbs S, Corral-Martínez P, Seguí-Simarro JM, Gaétan J, Coursaget P, Veramendi J. 2008. Human papillomavirus L1 protein expressed in tobacco chloroplasts
self-assembles into virus-like particles that are highly immunogenic. Plant Biotechnology Journal **6**, 427–441.

**Ferro M**, **Salvi D**, **Riviere-Rolland H**, **Vermat T**, **Seigneurin-Berny D**, **Grunwald D**, **Garin J**, **Joyard J**, **Rolland N.** 2002. Integral membrane proteins of the chloroplast envelope: identification and subcellular localization of new transporters. *Proceedings of the National Academy of Sciences, USA* **99**, 11487–11492.

**Flexas J**, **Barbour MM**, **Brendel O**, *et al.* 2012. Mesophyll diffusion conductance to CO₂: an unappreciated central player in photosynthesis. *Plant Science* **193-194**, 70–84.

**Flexas J**, **Díaz-Espejo A**, **Berry JA**, **Cifre J**, **Galmés J**, **Kaldenhoff R**, **Medrano H**, **Ribas-Carbó M.** 2007b. Analysis of leakage in IRGA's leaf chambers of open gas exchange systems: quantification and its effects in photosynthesis parameterization. *Journal of Experimental Botany* **58**, 1533–1543.

**Flexas J**, **Díaz-Espejo A**, **Conesa MA**, *et al.* 2016. Mesophyll conductance to CO₂ and Rubisco as targets for improving intrinsic water use efficiency in C3 plants. *Plant, Cell & Environment* **39**, 965–982.

**Flexas J**, **Díaz-Espejo A**, **Galmés J**, **Kaldenhoff R**, **Medrano H**, **Ribas-Carbo M.** 2007a. Rapid variations of mesophyll conductance in response to changes in CO₂ concentration around leaves. *Plant, Cell & Environment* **30**, 1284–1298.

**Flexas J**, **Niinemets U**, **Gallé A**, *et al.* 2013. Diffusional conductances to CO₂ as a target for increasing photosynthesis and photosynthetic water-use efficiency. *Photosynthesis Research* **117**, 45–59.

**Flexas J**, **Ribas-Carbó M**, **Díaz-Espejo A**, **Galmés J**, **Medrano H.** 2008. Mesophyll conductance to CO₂: current knowledge and future prospects. *Plant, Cell & Environment* **31**, 602–621.

**Flexas J**, **Ribas-Carbó M**, **Hanson DT**, **Bota J**, **Otto B**, **Cifre J**, **McDowell N**, **Medrano H**, **Kaldenhoff R.** 2006. Tobacco aquaporin NtAQP1 is involved in mesophyll conductance to CO₂ in vivo. *The Plant Journal* **48**, 427–439.

**Gago J**, **Douthe C**, **Florez-Sarasa I**, **Escalona JM**, **Galmes J**, **Fernie AR**, **Flexas J**, **Medrano H.** 2014. Opportunities for improving leaf water use efficiency under climate change conditions. *Plant Science* **226**, 108–119.

**Galmés J**, **Medrano H**, **Flexas J.** 2007. Photosynthetic limitations in response to water stress and recovery in Mediterranean plants with different growth forms. *New Phytologist* **175**, 81–93.

**Ghashgaheja J**, **Badeck FW**, **Lanigan G**, **Nogués S**, **Tcherkez G**, **Deléens E**, **Cornic G**, **Griffiths H.** 2003. Carbon isotope fractionation during dark respiration and photorespiration in C3 plants. *Phytochemistry Reviews* **2**, 145–161.

**Gnanasekaran T**, **Karcher D**, **Nielsen AZ**, *et al.* 2016. Transfer of the cytochrome P450-dependent dhurrin pathway from *Sorghum bicolor* into *Nicotiana tabacum* chloroplasts for light-driven synthesis. *Journal of Experimental Botany* **67**, 2495–2506.

**Gomes D**, **Agasse A**, **Thiébaud P**, **Delrot S**, **Gerós H**, **Chaumont F.** 2009. Aquaporins are multifunctional water and solute transporters highly divergent in living organisms. *Biochimica et Biophysica Acta* **1788**, 1213–1228.

**Griffith KL**, **Wolf RE Jr.** 2002. Measuring β-galactosidase activity in bacteria: cell growth, permeabilization, and enzyme assays in 96-well arrays. *Biochemical and Biophysical Research Communications* **290**, 397–402.

**Groszmann M**, **Osborn HL**, **Evans JR.** 2017. Carbon dioxide and water transport through plant aquaporins. *Plant, Cell & Environment* **40**, 938–961.

**Hanba YT**, **Shibasaka M**, **Hayashi Y**, **Hayakawa T**, **Kasamo K**, **Terashima I**, **Katsuhara M.** 2004. Overexpression of the barley aquaporin HvPIP2;1 increases internal CO₂ conductance and CO₂ assimilation in the leaves of transgenic rice plants. *Plant & Cell Physiology* **45**, 521–529.

**Harley PC**, **Loreto F**, **Di Marco G**, **Sharkey TD.** 1992. Theoretical considerations when estimating the mesophyll conductance to CO₂ flux by analysis of the response of photosynthesis to CO₂. *Plant Physiology* **98**, 1429–1436.

**Heckwolf M**, **Pater D**, **Hanson DT**, **Kaldenhoff R.** 2011. The *Arabidopsis thaliana* aquaporin AtPIP2;1 is a physiologically relevant CO₂ transport facilitator. *The Plant Journal* **67**, 795–804.

**Henig A**, **Bonfig K**, **Roitsch T**, **Warzecha H.** 2007. Expression of the recombinant bacterial outer surface protein A in tobacco chloroplasts leads to thylakoid localization and loss of photosynthesis. *The FEBS Journal* **274**, 5749–5758.

**Hoober JK**, **Boyd CO**, **Paavola LG.** 1991. Origin of thylakoid membranes in *Chlamydomonas reinhardtii* y-1 at 38°C. *Plant Physiology* **96**, 1321–1328.

**Hove RM**, **Bhave M.** 2011. Plant aquaporins with non-aqua functions: deciphering the signature sequences. *Plant Molecular Biology* **75**, 413–430.

**Kaldenhoff R.** 2012. Mechanisms underlying CO₂ diffusion in leaves. *Current Opinion in Plant Biology* **15**, 276–281.

**Kaldenhoff R**, **Kai L**, **Uehlein N.** 2014. Aquaporins and membrane diffusion of CO₂ in living organisms. *Biochimica et Biophysica Acta* **1840**, 1592–1595.

**Li M**, **Schnell DJ.** 2006. Reconstitution of protein targeting to the inner envelope membrane of chloroplasts. *The Journal of Cell Biology* **175**, 249–259.

**Long SP**, **Ainsworth EA**, **Leakey AD**, **Nösberger J**, **Ort DR.** 2006. Food for thought: lower-than-expected crop yield stimulation with rising CO₂ concentrations. *Science* **312**, 1918–1921.

**Loreto F**, **Tsonev T**, **Centritto M.** 2009. The impact of blue light on leaf mesophyll conductance. *Journal of Experimental Botany* **60**, 2283–2290.

**Loriaux SD**, **Avenson TJ**, **Welles JM**, **McDermitt DK**, **Eckles RD**, **Riensche B**, **Genty B.** 2013. Closing in on maximum yield of chlorophyll fluorescence using a single multiphase flash of sub-saturating intensity. *Plant, Cell & Environment* **36**, 1755–1770.

**Luu DT**, **Maurel C.** 2013. Aquaporin trafficking in plant cells: an emerging membrane-protein model. *Traffic* **14**, 629–635.

**Lübeck J**, **Heins L**, **Soll J.** 1997. A nuclear-coded chloroplastic inner envelope membrane protein uses a soluble sorting intermediate upon import into the organelle. *The Journal of Cell Biology* **137**, 1279–1286.

**Martins SC**, **Galmés J**, **Molins A**, **DaMatta FM.** 2013. Improving the estimation of mesophyll conductance to CO₂: on the role of electron transport rate correction and respiration. *Journal of Experimental Botany* **64**, 3285–3298.

**Maurel C**, **Reizer J**, **Schroeder JI**, **Chrispeels MJ.** 1993. The vacuolar membrane protein gamma-TIP creates water specific channels in *Xenopus oococytes*. *The EMBO Journal* **12**, 2241–2247.

**Maurel C**, **Verdoucq L**, **Luu DT**, **Santoni V.** 2008. Plant aquaporins: membrane channels with multiple integrated functions. *Annual Review of Plant Biology* **59**, 595–624.

**Okawa K**, **Inoue H**, **Adachi F**, **Nakayama K**, **Ito-Inaba Y**, **Schnell DJ**, **Uehara S**, **Inaba T.** 2014. Targeting of a polytopic membrane protein to the inner envelope membrane of chloroplasts in vivo involves multiple transmembrane segments. *Journal of Experimental Botany* **65**, 5257–5265.

**Otto B**, **Uehlein N**, **Sdorra S**, *et al.* 2010. Aquaporin tetramer composition modifies the function of tobacco aquaporins. *The Journal of Biological Chemistry* **285**, 31253–31260.

**Parry MA**, **Andralojc PJ**, **Scales JC**, **Salvucci ME**, **Carmo-Silva AE**, **Alonso H**, **Whitney SM.** 2013. Rubisco activity and regulation as targets for crop improvement. *Journal of Experimental Botany* **64**, 717–730.

**Parry MA**, **Reynolds M**, **Salvucci ME**, **Raines C**, **Andralojc PJ**, **Zhu XG**, **Price GD**, **Condon AG**, **Furbank RT.** 2011. Raising yield potential of wheat. II. Increasing photosynthetic capacity and efficiency. *Journal of Experimental Botany* **62**, 453–467.

**Pengelly JJ**, **Förster B**, **von Caemmerer S**, **Badger MR**, **Price GD**, **Whitney SM.** 2014. Transplastomic integration of a cyanobacterial bicarbonate transporter into tobacco chloroplasts. *Journal of Experimental Botany* **65**, 3071–3080.

**Pisareva T**, **Kwon J**, **Oh J**, **Kim S**, **Ge C**, **Wieslander A**, **Choi JS**, **Norling B.** 2011. Model for membrane organization and protein sorting in the cyanobacterium *Synechocystis* sp. PCC 6803 inferred from proteomics and multivariate sequence analyses. *Journal of Proteome Research* **10**, 3617–3631.

**Pons TL**, **Flexas J**, **von Caemmerer S**, **Evans JR**, **Genty B**, **Ribas-Carbo M**, **Brugnoli E.** 2009. Estimating mesophyll conductance to CO₂: methodology, potential errors, and recommendations. *Journal of Experimental Botany* **60**, 2217–2234.

**Priestley DA**, **Woolhouse HW.** 1980. The chloroplast envelope of *Phaseolus vulgaris* L..1. Isolation and compositional characteristics. *Journal of Experimental Botany* **31**, 437–447.

**Reynolds M**, **Bonnett D**, **Chapman SC**, **Furbank RT**, **Manès Y**, **Mather DE**, **Parry MA.** 2011. Raising yield potential of wheat. I. Overview of a consortium approach and breeding strategies. *Journal of Experimental Botany* **62**, 439–452.
Sade N, Gebretsadik M, Seligmann R, Schwartz A, Wallach R, Moshelion M. 2010. The role of tobacco Aquaporin1 in improving water use efficiency, hydraulic conductivity, and yield production under salt stress. Plant Physiology **152**, 245–254.

Santoni V. 2007. Plant plasma membrane protein extraction and solubilization for proteomic analysis. Methods in Molecular Biology **355**, 93–109.

Sanz-Barrio R, Corral-Martinez P, Ancin M, Segui-Simarro JM, Farran I. 2013. Overexpression of plastidial thioredoxin f leads to enhanced starch accumulation in tobacco leaves. Plant Biotechnology Journal **11**, 618–627.

Scafaro AP, von Caemmerer S, Evans JR, Atwell BJ. 2011. Temperature response of mesophyll conductance in cultivated and wild Oryza species with contrasting mesophyll cell wall thickness. Plant, Cell & Environment **34**, 1999–2008.

Scotti N, Sannino L, Idoine A, Hamman P, De Stradis A, Giorio P, Maréchal-Drouard L, Bock R, Cardi T. 2015. The HIV-1 Pr55 gag polyprotein binds to plastidial membranes and leads to severe impairment of chloroplast biogenesis and seedling lethality in transplastomic tobacco plants. Transgenic Research **24**, 319–331.

Shanmugabalaji V, Besagni C, Piller LE, Douet V, Ruf S, Bock R, Kessler F. 2013. Dual targeting of a mature plastoglobulin/fibrillin fusion

protein to chloroplast plastoglobules and thylakoids in transplastomic tobacco plants. Plant Molecular Biology **81**, 13–25.

Uehlein N, Lovisolo C, Siefritz F, Kaldenhoff R. 2003. The tobacco aquaporin NtAQP1 is a membrane CO₂ pore with physiological functions. Nature **425**, 734–737.

Uehlein N, Otto B, Hanson DT, Fischer M, McDowell N, Kaldenhoff R. 2008. Function of *Nicotiana tabacum* aquaporins as chloroplast gas pores challenges the concept of membrane CO₂ permeability. The Plant Cell **20**, 648–657.

Valentini R, Epron D, De Angelis P, Matteucci G, Dreyer E. 1995. In situ estimation of net CO₂ assimilation, photosynthetic electron flow and photorespiration in Turkey oak (*Quercus cerris* L.) leaves: diurnal cycles under different levels of water supply. Plant Cell and Environment **18**, 631–640.

Verdoucq L, Rodrigues O, Martinière A, Luu DT, Maurel C. 2014. Plant aquaporins on the move: reversible phosphorylation, lateral motion and cycling. Current Opinion in Plant Biology **22**, 101–107.

Walker B, Ariza LS, Kaines S, Badger MR, Cousins AB. 2013. Temperature response of in vivo Rubisco kinetics and mesophyll conductance in *Arabidopsis thaliana*: comparisons to *Nicotiana tabacum*. Plant, Cell & Environment **36**, 2108–2119.