High Photosynthetic Rates in a Solanum pennellii Chromosome 2 QTL Is Explained by Biochemical and Photochemical Changes

Enhanced photosynthesis is strictly associated with to productivity and it can be accomplished by genetic approaches through identification of genetic variation. By using a Solanum pennellii introgression lines (ILs) population, it was previously verified that, under normal (CO2), IL 2–5 and 2–6 display increased photosynthetic rates by up to 20% in comparison with their parental background (M82). However, the physiological mechanisms involved in the enhanced CO2 assimilation exhibited by these lines remained unknown, precluding their use for further biotechnological applications. Thereby, here we attempted to uncover the physiological factors involved in the upregulation of photosynthesis in ILs 2–5 and 2–6 under normal (CO2) as well as under elevated (CO2). The results provide evidence for increased biochemical capacity (higher maximum carboxylation velocity and maximum electron transport rate) in plants from IL 2–5 and 2–6, whereas the diffusive components (stomatal and mesophyll conductances) were unaltered in these ILs in comparison to M82. Our analyses revealed that the higher photosynthetic rate observed in these ILs was associated with higher levels of starch as well as total protein levels, specially increased RuBisCO content. Further analyses performed in plants under high (CO2) confirmed that biochemical properties are involved in genetic variation on chromosome 2 related to enhanced photosynthesis.


INTRODUCTION
Photosynthesis (A) is the main driving force for plant growth and biomass production having a central position for breeders seeking to increase crop yield (Evans, 2013;Nunes-Nesi et al., 2016;Nuccio et al., 2017;Batista-Silva et al., 2020;Flexas and Carriquí, 2020). The improvement of photosynthetic performance under optimal and suboptimal conditions can be accomplished by genetic approaches through identification of genetic variation, as well as selection of accessions exhibiting higher photosynthetic rates. Moreover, understanding the mechanisms involved in the photosynthesis associated traits are required for improving the photosynthetic efficiency in crop species.
The photosynthetic efficiency can be enhanced, among others, by improving light capture through increasing pigment contents and by ameliorating light energy conversion (Ort et al., 2015). Furthermore, A can be improved by changes in stomatal properties (Rebetzke et al., 2013), ribulose 1,5 bisphosphate carboxilase oxigenase (RuBisCO) kinetic attributes (Parry et al., 2013), and also by modifying expression of others enzymes involved in carbon reduction reactions (Miyagawa et al., 2001;Lefebvre et al., 2005;Driever et al., 2017).
Alterations in A may be associated with diffusional limitations as those derived from changes in stomatal (g s ) and mesophyll (g m ) conductances which are dependent on leaf anatomical and physiological properties and influenced by environmental cues (Flexas et al., 2007(Flexas et al., , 2012Vrabl et al., 2009;Galmés et al., 2011). For instance, high levels of CO 2 at the substomatal cavities (C i ) tend to reduce g s in leaves by an intrinsic property of guard cells (Mott, 1988;Ainsworth and Rogers, 2007). Remarkably, at lower C i , the maximum carboxylation velocity (V cmax ) limits photosynthesis. On the other hand, when C i is higher, CO 2 assimilation is limited by maximum electron transport rate (J max ), or the regeneration of triose phosphate (Farquhar and Sharkey, 1982). Moreover, enrichment in the availability of atmospheric (CO 2 ) augments RuBisCO carboxylase activity leading to reduced photorespiratory process, improving carbon gains (Long et al., 2004;Peterhansel et al., 2010).
Natural genetic variation in diffusional and biochemical traits related to photosynthetic capacity have been verified in different plant species (Geber and Dawson, 1997;Galmés et al., 2014;Singh et al., 2014). For example, Mediterranean accessions of tomato (S. lycopersicum) cultivated under drought stress demonstrated higher transpiration efficiency provided by enchanced g m /g s ratio (Galmés et al., 2011). Furthermore, altered photosynthetic capacity can also be explained by components related to CO 2 assimilation capacity namely V cmax and J max (Dreyer et al., 2001;Driever et al., 2014;Walker et al., 2014). Thus, selection of accessions with higher J max contributing to ribulose 1,5 bisphosphate (RuBP) regeneration has been proposed as a trait for plant breeding (Driever et al., 2014;Simkin et al., 2017a).
The use of introgression lines (ILs) populations, an important genomic tool for Quantitative Trait Locus (QTL) mapping, has been proposed to identify candidate genes associated with improved A (Flood et al., 2011;Gu et al., 2012;Nunes-Nesi et al., 2016). QTL analysis allow the identification of genomic regions where detailed studies should be performed to identify genes and nucleotide changes responsible for a certain function (Asín, 2002;Gonzalez-Martinez et al., 2006;Jamil et al., 2016;de Oliveira Silva et al., 2018). In a carbon isotope composition (δ 13 C) study performed with a Solanum pennellii population demonstrated that variation in stomatal arrangement may improve internal CO 2 diffusion allowing greater photosynthetic efficiency in dry environments (Xu et al., 2008;Muir et al., 2014). In addition, most of the genetic variation in tomato for g s might be explained by differences in both width and anatomical traits of the stomatal pore (Fanourakis et al., 2015).
Recently, genomic regions involved in the regulation of photosynthesis and respiration were identified in tomato (de Oliveira Silva et al., 2018). In this study, a QTL for high A was identified in a unique overlapping genomic region, delimited by two ILs (IL 2-5 and IL 2-6) and named as BIN 2K. The IL 2-5 exhibited up to 20% increased A, 50% higher nitrate content and 21% higher starch accumulation at the end of the light period in comparison with the parental line M82. Similarly, the IL 2-6 exhibited 19% increased A, 22% higher g s , 42% higher levels of starch at the end of the light period, and 58% higher starch turnover in comparison with M82. In agreement, the IL 2-6 also displayed higher g s , when evaluating the variation in stomatal responsiveness to desiccation (Fanourakis et al., 2015). These two ILs were also evaluated for δ 13 C, in which IL 2-5 exhibited higher values of δ 13 C compared to M82 and IL 2-6 (−28.1, −28.4, and −28.8, respectively) indicating increased water use efficiency (WUE) in this line (Xu et al., 2008). Notably, the IL 2-5 was also described as a drought-tolerant line mainly through alterations in organ morphogenesis and biochemical pathways (Gong et al., 2010). In addition, the IL 2-5, and IL 2-6 have also been characterized by exihibiting higher total plant weight, number of flower, fruits and seeds than M82 (Eshed and Zamir, 1995;Semel et al., 2006). However, the fruits from these ILs have reduced size and fresh weight when compared to M82 (Eshed and Zamir, 1995;Causse et al., 2004).
In order to investigate the physiological basis for the increased A associated with the genomic region BIN 2K, a detailed physiological characterization of the IL 2-5, IL 2-6, and M82 was performed under atmospheric and elevated (CO 2 ). Given that A is highly affected by environmental (CO 2 ) and RuBisCO limitation is reduced by high (CO 2 ; von Caemmerer and Farquhar, 1981), these ILs might exhibit higher carbon assimilation capacity under this condition. Our results demonstrate that the higher photosynthetic efficiency in the IL 2-5 and IL 2-6 was not associated with diffusive component (g s and g m ), but rather to increased photobiochemical capacity. Our findings open important avenues for enhancing our understanding of genetic variation on chromosome 2 that is involved in the regulation of photosynthetic capacity in tomato plants.

Plant Material and Experimental Conditions
Seeds of the ILs 2-5 and 2-6, as well as the cultivar M82 were sown in trays containing the commercial substrate Tropstrato HT, Vida Verde R . The plants were allocated in a greenhouse under naturally fluctuating conditions of light intensity (±600 µmol photons m −2 s −1 daily) and ambient CO 2 (±400 µmol mol −1 ) as well as semicontrolled conditions of temperature (±28 • C) and relative air humidity (±60%) during spring. After germination, seedlings were transplanted in 1.16 L pots containing the same commercial substrate, supplemented with chemical fertilizer (N:4; P 2 O 5 :14; K 2 O: 8) in the proportion of 0.5 kg of NPK to 10.0 kg of substrate. The plants were irrigated daily as needed and no restriction of root development, and consequently photosynthetic limitation, was observed at the end of the experiment when the physiological analyses and sample harvesting took place (4-week-old plants). Four plants for each genotype were grown side by side and the position of each plant was changed daily.
In a second experiment, seeds from the same genotypes were germinated and cultivated as described previoulsy with the following modifications. Pots containing plants with 10 days after transplanting (DAT) were transferred to open top chambers (OTCs) with the following specifications 1.6 m diameter and 1.8 m height. The chambers were placed in a greenhouse under naturally fluctuating conditions of light intensity (±600 µmol photons m −2 s −1 daily) and ambient CO 2 (±400 µmol mol −1 ) as well as semicontrolled conditions of temperature (±28 • C) and relative air humidity (±60%). The experiment was carried out in two OTCs. One was continuously maintained at the current ambient level of CO 2 (∼400 µmol mol −1 ) and another continuously maintained at an elevated level of CO 2 (∼800 µmol mol −1 ). In the chamber supplemented with elevated (CO 2 ), 800 µmol mol −1 was the setpoint and the (CO 2 ) was monitored and adjusted every day with a portable CO 2 meter AZ-77535. CO 2 was only injected following the photoperiod, i.e., from 06:00 a.m. to 06:00 p.m. The obtained (CO 2 ) along the experiments are presented in Supplementary Figure S1. Air temperature averages (±30 • C) and relative air humidity (±77%) inside the OTCs were measured throughout the experiment and did not differ significantly between the two sets of OTCs during the experiment. The two OTCs were installed at 1 m of distance between each other. The plants remained under these conditions for 3 weeks until further analysis. Four plants for each genotype were grown inside each chamber and the positions of the plants inside the chamber were daily changed.

Growth and Fruit Analyses
Growth parameters were determined in 4-week-old plants. Total leaf area was measured using a LI-3100C Area Meter (LiCor, Lincoln, NE, United States) while specific leaf area (SLA) was measured as described previously (Hunt et al., 2002). The dry matter accumulated in leaves, stems and roots were determined as described previously (de Oliveira Silva et al., 2018). Total fruit production was determined at the end of the cultivation in plants with 60 DAT. The fruits were harvested when they presented the reddish color, characteristic of mature fruits. We determined the total number of fruits per plant, the mean fruit fresh weight per plant as well as the fruit size, by measuring the fruit length, and diameter using a digital caliper.

Measurements of Gas Exchange Parameters and Chlorophyll a Fluorescence
Gas exchange and chlorophyll a fluorescence analyses were performed simultaneously using an open-flow infrared gas exchange analyzer system equipped with an integrated fluorescence chamber (IRGA, LI-COR Inc. LI-6400XT; NE). These analyses were realized during the light period, from 8 h to 12 h (solar time) using the 2 cm 2 leaf chamber at 25 • C, flow rate of 300 mol s −1 , 0.5 stomatal ratio (amphistomatic leaves), and saturating light intensity of 1,000 µmol photons m 2 s −1 . The leaf-to-air vapor pressure deficit ranged from 1.2 to 2.0 kPa and the amount of blue light was set to 10% of photosynthetic photon flux density (PPFD) to optimize stomatal aperture. The initial fluorescence (F 0 ) was measured in illuminating dark-adapted leaves (for at least 1 h) with weak modulated measuring beams (0.03 µmol m −2 s −1 ). A saturating white light pulse (8,000 µmol m −2 s −1 ) was applied for 0.8 s to obtain the maximum fluorescence (F m ), from which the variable-to-maximum Chl fluorescence ratio was then calculated: In light-adapted leaves, the steady-state fluorescence yield was measured with the application of a saturating white light pulse (8,000 µmol m −2 s −1 ) to achieve the light adapted maximum fluorescence (F m ). A far-red illumination (2 µmol m −2 s −1 ) was applied after turn off the actinic light to measure the light-adapted initial fluorescence (F 0 ). The capture efficiency of excitation energy by open PSII reaction centers (F v /F m ) was estimated following a previously described procedure (Logan et al., 2007) and the actual PSII photochemical efficiency (ϕPSII) was estimated as ϕPSII = (F m -F s )/F m (Genty et al., 1989). Electron transport rate (ETR) was calculated as follows: φPSII × PARi × β × α, where PARi is the incident radiation, β is the absorbed quanta partition factor between PSII/I (assumed to be 0.5), and α is the leaf absorbance. Leaf absorbance (α) was estimated by the chlorophyll content per unit area as follow: α = χ/(χ + 76), where χ is the chlorophyll content per unit area (Evans and Poorter, 2001).
From combined measurements of fluorescence and gas exchange, we estimated the photorespiration rate (P r ) according to Valentini et al. (1995). Dark respiration rate (R d ) was measured by the same gas exchange system after plants being at least 2 h into the dark period.
Photosynthetic light-response curves (A/PPFD) were performed using atmospheric (CO 2 ; C a ) of 400 µmol CO 2 mol −1 and the plants were exposed to a range of PPFD in the sequence: 1,000, 1,500, 1,300, 1,100, 1,000, 800, 600, 400, 300, 200, 100, 50, and 0 µmol photons m −2 s −1 . Variables derived from the A/PPFD curves were estimated from adjustments of light response curve by the non-rectangular hyperbolic model (von Caemmerer, 2000). Measurements for CO 2 response curves were taken at light saturation point of 1,000 µmol m 2 s −1 , 25 • C at atmospheric CO 2 concentration (C a ) of 400 µmol mol −1 and once the steady state was reached, the C a was decreased to 300, 200, 100, and 50 µmol mol −1 . Then, C a was increased to 400, 500, 600, 800, 1,000, 1,200, 1,400, and 1,600 µmol mol −1 . Curves were obtained using the second terminal leaflet of the third fully expanded leaf from the apex of 4−week−old plants.
For the high (CO 2 ) experiments, instantaneous gas exchange and chlorophyll a fluorescence analyses were performed as previously described, with minor modifications. All measurements were performed during the light period between 8 h and 12 h (solar time) under 1,000 µmol photons m −2 s −1 FIGURE 1 | Growth phenotypes of 4-week-old plants from two ILs of Solanum pennellii into a genetic background of Solanum lycopersicum (M82) grown at 400 µmol CO 2 mol -1 . Total leaf area (A); Specific leaf area (B); Leaf dry weight (C); Branch dry weight (D); Root dry weight (E); Shoot dry weight (F); Root: shoot dry weight ratio (G); and Total dry weight (H). Values are presented as means ± SE (n = 4). Different letters accompany means that differ between the genotypes (P < 0.05) by the Tukey test.
Frontiers in Plant Science | www.frontiersin.org Values are presented as means ± SE (n = 4). Asterisks indicate differences of the ILs in relation to M82 (P < 0.05) according to the Tukey test.
at the leaf level (light saturation) of PPFD, determined by A/PPFD curves and at 25 • C. The reference (CO 2 ) was set at 400 µmol mol −1 for plants at atmospheric (CO 2 ) and 800 µmol mol −1 for plants under elevated (CO 2 ). Instantaneous measurements were obtained using the second terminal leaflet of the third fully expanded leaf from the apex of 4−week−old plants. R d was determined as described previously.

Estimation of Mesophyll Conductance (g m ), Maximum Rate of Carboxylation (V cmax ), Maximum Rate of Electron Transport (J max ), and Photosynthetic Limitations
The concentration of CO 2 in the carboxylation sites (C c ) was calculated following Harley et al. (1992) as gm = A / Ci − ( * (ETR + 8(A + Rd)))/(ETR − 4(A + Rd)) , and using * for tomato measured previously (Hermida-Carrera et al., 2016). The g m data were employed to convert the A-C i curves into A-C c curves using the following equation C c = C i −(A/g m ). ETR calculation using the proxy of α estimation, by the chlorophyll content empirical equation (Evans and Poorter, 2001), was additionally checked through a sensitivity analysis to ensure that the possible variability in α do not affect g m estimates. Values for g m were also estimated by the Ethier and Livingston (2004) method.
The parameters from A/C c in plants under ambient (CO 2 ), as V cmax and J max were calculated by fitting the mechanistic model of CO 2 assimilation (Farquhar et al., 1980), using the C c based on temperature of kinetic parameters of RuBisCO (K c and K o ). V cmax , J max , and g m were normalized to 25 • C using the temperature response and plug-in equations described by Sharkey et al. (2007). For the second experiment, the V cmax single point was only calculated for plants growing under atmospheric (CO 2 ; 400 µmol mol −1 ; De Kauwe et al., 2016).
Values of A, g s , g m , and ETR were employed to calculate the quantitative photosynthesis limitation analysis as described by Grassi and Magnani (2005) that permits the partitioning into the functional components of photosynthetic constraints related to stomatal (l s ), mesophyll (l m ), and leaf biochemical (l b ) limitations. According to A-C c curves obtained in the firt experiment, photosynthesis at the growing (CO 2 ) was found at the curvi-linear part of the response, associated to the energy limited RuBP regeneration (J max ). Since ETR is tightly coupled to J max , and should reflect gross photosynthesis (Genty et al., 1989;Valentini et al., 1995), calculations of l b were confirmed using ETR as a surrogate for leaf biochemistry (Galle et al., 2009;Pons et al., 2009).

Metabolite Analyses
Samples from fully expanded leaves were harvested in five different time points (beginning, middle and end of the light period, and as well as middle and end of the dark period) for the first experiment and two different time points (middle of the light period and end of the dark period) for the second experiment from 4−week−old plants. The material was frozen in liquid nitrogen and stored at -80 • C for further analyses. Samples were homogenized and aliquots of approximately 20 mg were subjected to hot ethanol extraction as described by Cross et al. (2006). The levels of starch, glucose, fructose, and sucrose were measured as described by Fernie et al. (2001). The levels of chlorophyll (a and b) were determined exactly as previously described (Porra et al., 1989). Nitrate contents were determined as detailed by Sulpice et al. (2009) and malate and fumarate as Nunes-Nesi et al. (2007). The total protein was measured according by Bradford (1976) and the total amino acids content was determined using a colorimetric assay according to Yemm et al. (1955).

Protein Extraction and Western Blot Analysis
Total soluble protein was extracted from 200 mg (fresh weight) of leaf tissue according to the protocol of Vasconcelos et al. (2005) These data are derived from A/C c response curve shown in Figure 2B at ambient CO 2 concentration (C a ) of 400 µmol mol −1 setpoint. Measurements of photosynthetic carbon fixation rates were determined at saturating-light levels of 1,000 µmol m −2 s −1 in the second terminal leaflet of the third fully expanded leaf from ILs and M82 plants. Values are presented as means ± SE (n = 4). Different letters accompany means that differ between the genotypes (P < 0.05) by the Tukey test. A, net CO 2 assimilation rate; A gross , A + P r + R d ; R d , dark respiration; P r , photorespiration; ETR, electron transport rate; g s , stomatal conductance; g m , mesophyll conductance to CO 2 ; C c , chloroplastic CO 2 concentration; V cmax _C c , maximum carboxylation velocity based on C c ; and J max _C c , maximum capacity for electron transport rate based on C c ; l s , stomatal limitation; l m , mesophyll limitation; and l b , biochemical limitation.
and quantified using the Bradford method (Bradford, 1976). Forty micrograms of total protein was subjected to a solution containing sample buffer (125 mM Tris-HCl, pH 6.8, containing SDS 4%, 20% glycerol, 4% 2-mercaptoethanol, and 0.1% bromophenol blue). The samples were separated using 12% glycine sodium dodecyl sulfate polyacrylamide gel electrophoresis (Glycine-SDS-PAGE). The proteins were electrophoretically transferred to nitrocellulose membrane (0.45 µm; Bio-Rad) by the method of Towbin et al. (1979). The membranes were blocked for 2 h with 3% (w/v) Bovine serum albumin in Tris-buffered saline (2 M Tris-HCl, pH 7.4, containing 5 M NaCl, and 10% tween). The blots were then incubated for 2 h with specific antibodies raised in rabbit (serum diluted 1:12000 large subunit -AS03 037 code). The membranes were treated for 2 h with bridging antibodies (Anti-rabbit HRP-conjugated secondary antibodies GE Healthcare, RPN2108 code) diluted 1:10000. Large RuBisCO subunit was visualized with substrate solution (10 mg of 3,3' -Diaminobenzidine mixed with 10 mL of Tris-buffered saline and 10 µL of H 2 O 2 ). All steps were performed at room temperature. Tween-20 in Tris-buffered saline was used for washing (five times for 5 min between all treatments mentioned above) as well as for the dilution of antibodies. The bands on immunoblots were quantified by densitometry using the semi-quantitative ImageJ software (Version 1.1.4) as previously described (Tiranti et al., 2009;Brandt et al., 2018). Data were normalized by standard marker (∼50 kDa) and expressed in fold change.

Statistical Analyses
The experiments were performed in a randomized complete block design. The data of the two experiments was submitted to analysis of variance (ANOVA). In the second experiment, two levels of the same categorical variable (CO 2 environment) were considered. The variability of the results is expressed as the mean ± standard error (SE) of n independent biological replicates. Four biological replicates (n = 4) were utilized for each analysis, except for Western blot, where three biological replicate (n = 3) were analyzed. In order to verify the statistical difference between the genotypes, the means were compared using the Tukey test (P < 0.05). Pearson correlation matrices were performed to determine the correlated structure of the different parameters which were considered to be significant when P < 0.05. All analyses were performed using the GENES program (Cruz, 2013).

Growth-Related Parameters
In order to evaluate growth related phenotypes, ILs 2-5, 2-6, and M82 plants were grown side by side in a greenhouse. Detailed examination revealed that both ILs tended to perform better than M82, whereas IL 2-5 exhibited increased total leaf area ( Figure 1A), IL 2-6 displayed higher total dry weight ( Figure 1H). Despite certain differences in leaf, branch, and root dry weights (Figures 1C-E), no alteration in biomass allocation was observed, as indicated by similar root:shoot ratios in all genotypes ( Figure 1G).

Photosynthesis Related Parameters
Since the ILs 2-5 and 2-6 previously showed higher photosynthetic rate (de Oliveira Silva et al., 2018), we decided to investigate the photosynthetic capacity of these ILs by varying Values for Western blot are expressed as fold change over the control M82 and presented as means ± SE (n = 3). Signal intensities were measured using ImageJ software and normalized to the amount of loaded marker (∼50 kDa). Different letters accompany means that differ between the genotypes (P < 0.05) by the Tukey test. FW: Fresh weight.
light intensity. We found a similar A response to different light intensities in all genotypes (Figure 2A), therefore no differences for light saturation and compensation points were observed (Supplementary Table S1). We next characterized A in these lines under variable CO 2 concentrations by performing A/C c curve analysis. First, there were no differences in calculations on C i basis (data not shown) but only on C c . Therefore, it was necessary to carry out a complete photosynthetic study of the calculations considering the g m , including to evaluate the photosynthetic limitations for each of the lines. Interestingly, increased photosynthetic rates were observed under specific conditions for the two ILs ( Figure 2B). From A/C c curves, both ILs displayed higher CO 2 assimilation rate at the final point of the curve when compared to M82 (Figure 2B).
Concerning photosynthetic characterization of plants under atmospheric (CO 2 ; Table 1), further analyses indicated that plants from the two ILs displayed similar A, R d , g s , and g m . On the other hand, IL 2-6 displayed increased gross photosynthesis (A gross ) and both ILs higher P r , ETR and lower chloroplastic CO 2 concentration (C c ), probably due to the higher carboxylation demand, driven by the maximum RuBisCO carboxylation velocity based on C c (V cmax _C c ). In addition, the J max was higher in the ILs in comparison to M82.

Protein and RuBisCO Content in Leaves
The total protein levels were consistently increased by up to 21% in the IL 2-5 in comparison to M82 cultivar ( Figure 3A). Thus, we next analyzed the levels of RuBisCO via western blot with antibodies against its large subunit (RbcL; Figures 3B,C and Supplementary Figure S3). By quantifying the band intensity, we verified a significant increase (∼20%) in the levels of RbcL in both ILs compared to M82 (Figure 3C).

Diurnal Pattern of Metabolite Levels in Leaves
Increased starch content in leaves at the end of the light period in the ILs 2-5 and 2-6 was observed when compared to M82 ( Figure 4A). On the other hand, at the end of dark period, the starch levels were reduced in IL 2-6, indicating higher starch degradation during the night period ( Figure 4A). The glucose levels were increased at the end of the light period and sucrose showed lower levels in the middle of the night in the IL 2-5 compared to M82 (Figures 4B,D). Regarding the levels of fructose, we did not observe differences among the three genotypes ( Figure 4C). The levels of malate and fumarate were generally similar to those observed in leaves of M82, except for a reduction in malate levels at the beginning of the day (Figures 4E,F) and in fumarate at the end of the night. Total chlorophyll and amino acids as well as nitrate levels did not differ among genotypes (Supplementary Figures S2A,C,D).

Fruit Related Parameters
We also evaluated whether the greater photosynthetic efficiency would reflect on total fruit production (Supplementary Figure S4). We observed that plants from IL 2-5 produced greater number of fruits with 88% more fruits than M82 plants (Supplementary Figure S4A). However, both ILs produced smaller fruits than M82 plants, which resulted in no differences in total fruit fresh weight per plant between ILs and M82 (Supplementary Figures S4B-D).

Elevated (CO 2 ) Modifies Photosynthetic Parameters as Well as Carbon and Nitrogen Metabolism in Tomato ILs
Given the increased photosynthetic capacity observed in the ILs, indicated by the higher assimilation rate at saturating (CO 2 ), we comparatively analyzed plants from M82 and the two ILs growing for 3 weeks under atmospheric and elevated (CO 2 ). We noticed under the OTC conditions, which are slightly different from the greenhouse conditions, such as higher air temperature and relative air humidity, a minor decrease in the photosynthetic capacity of all genotypes as they exhibited lower A gross and ETR (Table 2). However, the ILs still displayed higher ETR, V cmax values ( Table 2), and protein content (IL 2-5; Figure 6B), confirming the superior photobiochemistry capacity of these ILs in comparison to M82 plants. Under elevated (CO 2 ), an increase of 22% in assimilation rate was observed in the ILs in comparison to M82 plants ( Table 2). In addition, A gross , P r , and ETR were significantly higher in the ILs at high (CO 2 ; Table 2). No changes were observed for R d , g s , and g m at high (CO 2 ) among ILs and M82 plants ( Table 2). Whereas stomatal limitation was increased in the IL 2-5, biochemical limitation decreased in both ILs at elevated (CO 2 ; Table 2).
To further understand the physiological advantage provided by increased photobiochemistry in the ILs, we explored the relationship between diffusive and biochemical limitations (l d /l b ) with A and ETR comparing ILs to M82 as well as the two growing conditions [400 and 800 ppm (CO 2 )] (Figure 5). At 400 ppm (CO 2 ), as expected diffusive limitations were double than biochemicals (l d /l b > 2) where ILs obtained higher A mostly due to lower biochemical limitation (thus promoted higher l d /l b ratio than M82; Figure 5A) in comparison to M82. However, under 800 ppm (CO 2 ) diffusive limitations became less important with all limitations highly balanced (l d /l b ∼1), and again ILs even at high (CO 2 ) showed biochemical limitations more reduced than diffusive ones (thus promoted higher l d /l b ratio than M82) when Measurements of photosynthetic carbon fixation rates were determined at ambient and elevated CO 2 concentration (C a ) of 400 µmol mol −1 and 800 µmol mol −1 , respectively, and at saturating-light levels of 1,000 µmol m −2 s −1 in the second terminal leaflet of the third fully expanded leaf from ILs and M82 plants. Values are presented as means ± SE (n = 4). Capital letters compare means that differ between the genotypes within a treatment (P < 0.05) by the Tukey test. Lowercase letters compare means that differ in a single genotype between the two treatments (P < 0.05) by the Tukey test. A, Assimilation rate; A gross , A + P r + R d ; R d , dark respiration; P r , photorespiration; ETR, electron transport rate; g s, stomatal conductance g m , mesophyll conductance to CO 2 estimated according to the Harley method; C c , chloroplastic CO 2 concentration; V cmax _ C c , maximum carboxylation velocity single point based on C c; l s , stomatal limitation; l m , mesophyll limitation; and l b , biochemical limitation.
compared to M82 (Figure 5A). The very same trend observed for the gas-exchange A parameter was observed for the fluorescence ETR value where ILs showed higher values than M82 at both growing conditions ( Figure 5B). We next analyzed the levels of the major leaf metabolites in the ILs under elevated (CO 2 ; Figure 6). In agreement with previous results suggesting higher photosynthetic efficiency of both ILs above 800 ppm (CO 2 ), significant increase in the total chlorophyll (IL 2-5), and protein contents (IL 2-5) at the middle of the day, as well as fructose (both ILs) at the end of the dark period was verified (Figures 6A,B,E). Whilst sucrose at the end of the dark period increased in IL 2-5, it was reduced in IL 2-6 in ambient and elevated (CO 2 ) compared to M82 (Figure 6F). On the other hand, glucose, fructose, and sucrose, total amino acids, and malate contents were unaltered at the middle of the light period in both ambients (Figures 6D-F and Supplementary  Figures S5B,C). Similar to the previous observations under atmospheric (CO 2 ), an increase in starch levels was observed in FIGURE 5 | Relationships between the ratio between diffusive and biochemical photosynthetic limitations (l d/ l b ) and photosynthesis (A; A) and electron transport rate (ETR; B) of 4-week-old plants from two ILs of Solanum pennellii into a genetic background of Solanum lycopersicum (M82) grown at 400 and at 800 µmol CO 2 mol −1 . The l d/ l b is the ratio between the diffusive limitations (g s + g m ) divided by the biochemical limitations. g s and g m, stomatal conductance and mesophylic conductance, respectively. Each point represents biological replicate. Circle symbol indicates genotypes cultivated under ambient (CO 2 ) in the legend and triangle symbol, elevated (CO 2 ). All correlations are significant at P < 0.005. R 2 : determination coefficient. Values are presented as means ± SE (n = 4). Capital letters compare means that differ between the genotypes within a treatment (P < 0.05) by the Tukey test. Lowercase letters compare means that differ in a single genotype between the two treatments (P < 0.05) by the Tukey test. FW: Fresh weight.
plants from both ILs compared to M82 at the middle of the light period ( Figure 6C). Noteworthy, there was an increase in starch levels at the middle of the light period in all genotypes under elevated (CO 2 ), nevertheless the ILs exhibited higher starch levels in relation to M82 (Figure 6C).
Pearson correlation analysis was applied to further comprehend how metabolic changes and structural components are driving A in the ILs (Supplementary Figure S6). No significant correlation was observed relating A and general leaf biochemical and structural parameters as total chlorophyll and protein content at both (CO 2 ) treatments (Supplementary Figure S6), indicating that the higher A observed in both ILs is mostly driven by specific changes at the photobiochemistry level ( Figure 5).

In silico Expression Analysis of Genes From BIN 2K
Based on leaf expression data from a recent study (Ranjan et al., 2016), we examined the BIN 2K differentially expressed genes in leaves from IL 2-5 and IL 2-6 plants (Supplementary Table S2). We observed that a subset of 30 genes were differentialy expressed in the ILs in comparison with M82. Out of 30 genes, nine were upregulated while 19 genes were downregulated only in IL 2-5. Regarding the IL 2-6, three genes were upregulated and 11 negatively regulated. We also verified that three genes were upregulated in plants from both ILs namely an AMPactivated protein kinase (Solyc02g091530.2.1), a polypeptide inferred from GFF3 feature (Solyc02g091770.1.1), and an adiponectin receptor 2 (Solyc02g092230.2.1). Noteworthy, four genes were strongly downregulated in plants from both IL 2-5 and 2-6 including a transcription factor described as CBF_NF-Y_archaeal histone (Solyc02g090930.2.1), a multidrug resistance protein mdtK from MATE family transporter related proteins (Solyc02g091070.2.1), an epidermal patterning factorlike protein 2 (Solyc02g091910.1.1), and a formamidase-like protein (Solyc02g092530.2.1). It is important to mention, however, that caution should be taken as not to over−interpret this comparison given that the RNAseq information was obtained from an entirely independent experiment and solely for leaves.

DISCUSSION
Given our previous demonstration that CO 2 assimilation rates of both ILs 2-5 and 2-6 plants were on average 20% higher compared to M82 (de Oliveira Silva et al., 2018), we decided to deeper investigate the physiological mechanisms responsible for their enhanced CO 2 assimilation. After evaluating photosynthetic parameters, we observed that both ILs exhibited substantially higher A (∼22%), and ETR under high (CO 2 ; Table 2), whereas no alterations were observed under ambient (CO 2 ) in the lines compared to their respective M82. Regarding photosynthetic results under ambient (CO 2 ), it is well established that modulations of photosynthesis mainly occur in response to environmental factors including light and temperature. In fact, the conditions of the current experiment were characterized by higher temperatures and light intensities in comparison with those previously described (de Oliveira Silva et al., 2018). Possibly, such conditions led to higher P r (∼12 µmol CO 2 m −2 s −1 ; Table 1) than those previously observed for ILs (∼8 µmol CO 2 m −2 s −1 ; de Oliveira Silva et al., 2018). As a result, even with an increased V cmax and J max (and, consequently, A gross ), increases in A could not be observed in the ILs at ambient (CO 2 ), as the lower C c led to higher P r under these conditions (Table 1). Nevertheless, the two ILs exhibited increased levels of RuBisCO, while only IL 2-5 displayed significantly higher protein content under ambient (CO 2 ; Figure 3). Increased protein levels can be overall associated with higher V cmax values providing an important source of variation in carbon uptake capacity (Hikosaka, 2004) and it is in agreement with the elevated biomass production and total leaf area of the ILs (Figures 1A,H).
The increased photosynthesis in the ILs under high (CO 2 ) could be explained, at least partially, by higher ETR values ( Table 2). Considering that the regeneration of RuBP is determinant of the limitation of A under high (CO 2 ; Mitchell et al., 2000), even small increments in photochemical components (e.g., ETR) may have resulted in higher photosynthetic activity in the ILs (Figure 5). The manipulation of the photosynthetic electron transport chain has been proposed as a valuable tool to enhance the carbon assimilation in plants (Simkin et al., 2017b(Simkin et al., , 2019. In agreement, transgenic Arabidopsis and tobacco plants expressing algal Cytochrome c 6 gene exhibited increase in ETR, A as well as chlorophyll and starch contents, and further increment in biomass production (Chida et al., 2007;Yadav et al., 2018).
Close association between carbon and nitrogen metabolisms has been demonstrated in plants (Foyer et al., 1998;Araújo et al., 2008;Gauthier et al., 2010). Thus, any alteration in carbon abundance is expected to reflect in nitrogen levels (Gutschick, 1999;Nunes-Nesi et al., 2010). Among the metabolites measured in plants under elevated (CO 2 ), chlorophyll and protein accumulated in the IL 2-5 as well as starch in leaves from both ILs (Figures 6A-C), which must be interpreted in light of photosynthetic data. Chlorophyll is an important part of the light-dependent reactions, and it is responsible for light harvesting during photosynthesis, therefore it has been considered an important predictor of photosynthetic potential (Croft et al., 2017;Liu et al., 2019). In its turn, the higher A was determining for the increased starch content, without, however, leading to the retroinhibition of the photosynthetic process ( Figure 6C). In good agreement, the high protein levels of ILs could mean a higher demand for proteins involved in the metabolism and export of carbohydrates resulting from the higher photosynthetic rate.
Based both on the available leaf expression data (Chitwood et al., 2013;Ranjan et al., 2016) and in our previously published work (de Oliveira Silva et al., 2018), we verified that several genes are differentially expressed in the IL 2-5 and IL 2-6. In addition to genes mentioned here (Supplementary Table S2), three other genes appears as potential candidates involved in the responses that lead to higher photosynthetic performance observed in the ILs 2-5 and 2-6, namely the ATP−dependent Clp protease (Solyc02g091280), NADH dehydrogenase ironsulfur protein 4 (Solyc02g092270), and beta subunit of an ATP synthase (Solyc02g091130; de Oliveira Silva et al., 2018). The proteins encoded by these genes are largely involved in energy metabolism and therefore might be involved with photosynthetic capacity. Moreover, genes enconding proteins related with photosynthesis, chlorophyll biosynthesis, and response to light stimulus were also upregulated mainly in the IL 2-5 (Ranjan et al., 2016). Noteworthy, the candidate genes for the higher photosynthetic capacity observed in plants from ILs 2-5 and 2-6 (Supplementary Table S2; and de Oliveira Silva et al., 2018) should be further investigated in more details. In this regard, introgression sublines with introgressed fragments smaller than the original set of ILs (Alseekh et al., 2013) are an alternative to reduce genetic variance and the amount of genes to be further analyzed; consequently, shortening the detection of the genes influencing photosynthetic capacity in tomato plants.
By analyzing the effects of the genetic variation on chromosome 2 of ILs 2-5 and 2-6 that modulates the upregulation of photosynthetic capacity, we have here provided novel insights into the mechanisms that govern photosynthetic capacity in tomato. Considering the higher values of J max in ILs under ambient (CO 2 ), followed by the enhanced A and ETR values in ILs exposed to elevated (CO 2 ), it seems reasonable to assume that both ILs optimize their investment in the components involved in RuBP regeneration. Thus, genetic manipulation in these components may be an opportunity to increase photosynthesis reflecting positively on growth in a future environment with high (CO 2 ; Ding et al., 2016;Driever et al., 2017;Simkin et al., 2017b), or currently in greenhouses that employ CO 2 fertilization to improve productivity (Mortensen, 1987). Noteworthy, considering that the control of the CO 2 assimilation in the Calvin-Benson cycle is shared among several enzymes, an increase in photosynthetic potential could be achieved combining multigene manipulation (Raines, 2003;Simkin et al., 2019). Recent studies demonstrated that increasing simultaneously the activity of different enzymes of the C3 cycle in the same plant resulted in an increased CO 2 fixation and biomass yield dramatically (Simkin et al., 2015(Simkin et al., , 2017a. Recent challenges associated with an engineering improved photosynthesis leading to the optimized utilization of CO 2 to maximize crop production have been discussed elsewhere . In this regard, the S. pennellii ILs population represents an opportunity to explore the large quantity of traits that affect plant survival as well as yield potential once ILs facilitate the identification of large number of genes affecting single or several plant phenotypes. Perhaps more importantly, our findings not only open up a novel fundamental research avenue, given that they identify a genomic region that lead to significant enhanced photosynthetic capacity, but also have important biotechnological implications, specifically with regard to the potential for increased yield, without increasing the need for arable land.

DATA AVAILABILITY STATEMENT
All datasets generated for this study are included in the article/Supplementary Material.