Expression of Exogenous GFP-CesA6 in Tobacco Enhances Cell Wall Biosynthesis and Biomass Production

Simple Summary Cellulose is synthesized at the plasma membrane by an enzymatic complex constituted by different cellulose synthase (CesA) proteins. The overexpression of CesA genes has been assessed for increasing cellulose biosynthesis and plant biomass. In this study, we analyzed transgenic tobacco plants (F31 line), stably expressing the Arabidopsis CesA6 fused to GFP, for possible variations in the cellulose biosynthesis. We found that F31 plants were bigger than the wild-type (wt), showing significant increases of stem height, root length, and leaf area. They bloomed about 3 weeks earlier and yielded more flowers and seeds than wt. In the F31 leaves, the expression of the exogenous GFP-CesA6 prompted the overexpression of all CesAs involved in the synthesis of primary cell wall cellulose and of other proteins responsible for plant cell wall building and remodeling. Instead, secondary cell wall CesAs were not affected. In the F31 stem, showing a 3.3-fold increase of the secondary xylem thickness, both primary and secondary CesAs expression was differentially modulated. Significantly, the amounts of cellulose and matrix polysaccharides increased in the transformed seedlings. The results evidence the potentiality to overexpress primary CesAs in tobacco for biomass production increase. Abstract Improved cellulose biosynthesis and plant biomass represent important economic targets for several biotechnological applications including bioenergy and biofuel production. The attempts to increase the biosynthesis of cellulose by overexpressing CesAs proteins, components of the cellulose synthase complex, has not always produced consistent results. Analyses of morphological and molecular data and of the chemical composition of cell walls showed that tobacco plants (F31 line), stably expressing the Arabidopsis CesA6 fused to GFP, exhibits a “giant” phenotype with no apparent other morphological aberrations. In the F31 line, all evaluated growth parameters, such as stem and root length, leaf size, and lignified secondary xylem, were significantly higher than in wt. Furthermore, F31 line exhibited increased flower and seed number, and an advance of about 20 days in the anthesis. In the leaves of F31 seedlings, the expression of primary CesAs (NtCesA1, NtCesA3, and NtCesA6) was enhanced, as well as of proteins involved in the biosynthesis of non-cellulosic polysaccharides (xyloglucans and galacturonans, NtXyl4, NtGal10), cell wall remodeling (NtExp11 and XTHs), and cell expansion (NtPIP1.1 and NtPIP2.7). While in leaves the expression level of all secondary cell wall CesAs (NtCesA4, NtCesA7, and NtCesA8) did not change significantly, both primary and secondary CesAs were differentially expressed in the stem. The amount of cellulose and matrix polysaccharides significantly increased in the F31 seedlings with no differences in pectin and hemicellulose glycosyl composition. Our results highlight the potentiality to overexpress primary CesAs in tobacco plants to enhance cellulose synthesis and biomass production.


Introduction
The cell wall is a dynamic polymeric network surrounding all plant cells and organized in layered structures each characterized by specific properties. Plant cell growth occurs simultaneously to the synthesis and deposition of a primary cell wall consisting of cellulose microfibrils embedded in a hydrated matrix of pectins, hemicelluloses, and proteins. This structural organization ensures tensile strength to oppose turgor pressure and the proper plasticity to allow water uptake and cell expansion [1]. Only some specialized cell types of the mechanical (fibers, sclereids) and vascular (xylem vessels) tissues appose a thickened secondary wall to the primary wall when expansion growth ceases. Strength, rigidity, and hydrophobicity are characteristics of the secondary wall, mainly consisting of cellulose, hemicelluloses, and lignin.
Cellulose is the most abundant polysaccharide of the biosphere with more than 10 11 tons estimated to be synthesized each year [2]. Cellulose consists of 1,4 linked β-glucan linear chains synthesized at the plasma membrane by the cellulose synthase complex (CSC), a multiprotein complex showing a symmetrical organization in six subunit clusters, each constituted by different cellulose synthase (CesA) proteins. The CSC architecture is not clearly defined, yet. Probably, three different CesA isoforms take part to the synthesis of cellulose in the primary or secondary cell wall. In Arabidopsis, CesA1, CesA3, and CesA6like groups (CesA2, CesA5, CesA6, and CesA9) are involved in primary cell wall growth, whereas the deposition of the secondary cell wall cellulose seems to be coordinated by the CesA4, CesA7, and CesA8 isoforms [3][4][5]. For a long time, cellulose microfibril was assumed to contain 36 β-1,4 glucan chains, supporting the CSC model in which each of the six subunits of the rosette contains six CesAs [6][7][8][9]. Accurate analytical techniques, including atomic force microscopy (AFM), nuclear magnetic resonance (NMR), single molecule imaging, living cell confocal microscopy, and the most advanced structural illumination microscopy (SIM), lean towards a structure consisting of 18-24 β-1,4 glucan chains [10][11][12][13]. Consequently, the "hexamer of hexamers" CSC model has been replaced with the "hexamer of trimers" in which each three CesA isoforms involved in the synthesis of primary or secondary wall are likely present in an equal stoichiometric ratio of 1:1:1 [14,15]. Moreover, new insights on the catalytic activity of single CesAs were recently obtained. The heterologous expression in yeast of Populus tremula × tremuloides CesA8 (PttCesA8), involved in secondary cell wall formation, was demonstrated to be sufficient to produce cellulose microfibrils in vitro [16]. Similarly, Physcomitrella patens CesA5, involved in primary cell wall formation, reconstituted into proteoliposomes, fostered the synthesis and self-assembly of cellulose microfibrils [17]. In addition, alternative 1:1:1 CSC models have been proposed as in the case of aspen (Populus tremula L.) secondary cell wall deposition during xylem vessel differentiation, where CesA8, CesA4, and CesA7 stoichiometry was found to shift from 3:2:1 in normal wood to 8:3:1 in tension wood driven by an increase of CesA8 expression [18]. Therefore, the organization of CesA proteins yet remains an open question and the possibility that it is species-and/or tissue-specific, as suggested by the studies on aspen, open up new perspectives for biotechnological applications, mainly related to increasing cellulose production and engineering of bioenergy crops.
Several studies have been performed to improve cellulose synthesis by increasing the expression levels of a single CesA obtaining unsatisfactory results, especially for secondary wall associated CesAs. The overexpression of PtCesA8 in aspen resulted in silencing of both the exogenous insert and its native counterpart, with a massive cellulose reduction in the secondary xylem and consequent alteration of plant growth [19]. A similar silencing response was reported in barley plants individually transformed with primary (CesA1, CesA2, and CesA6) and secondary (CesA4 and CesA8) cell wall CesA cDNAs. No increase of the total cellulose content in any of the transgenic barley lines was obtained and, in some lines, a significant decrease was evidenced other than some aberrant phenotypes [20]. Furthermore, in switchgrass (Panicum virgatum L.), CesA4 and CesA6 overexpression resulted in a decrease of cellulose content and crystallinity accompanied by increased xylan content [21]. On the contrary, in Arabidopsis, the overexpression of primary cell wall CesA2, CesA5, and CesA6 resulted in an improvement of cellulose biosynthesis and plant growth that seem to be related to changes in expression genes mainly related to cell division and expansion [22].
In a previous study [23], five independently transformed tobacco lines expressing the AtCesA6 gene fused at the C-terminus to the GFP fluorescent tag were obtained aimed at studying the internalization pattern of the resulting GFP-CesA6 chimeric protein [23]. The secretion pattern of the chimera GFP-CesA6 was characterized, highlighting that it recycled between plasma membrane and various subcellular compartments [23] following the CesA recycling pattern that underlies the dynamic remodeling of the cellulose synthase complex [3,24,25]. In the present study, we analyzed one of the stably transformed tobacco lines for possible variations in the synthesis of cellulose. Through different experimental approaches, we show that the expression of the exogenous GFP-CesA6 affects the biosynthetic machinery involved in the synthesis and remodeling of cell wall, causing an increase of cellulose and matrix polysaccharides, secondary xylem, and plant biomass.

Plants Material and Growth Conditions
Stable transgenic tobacco plants expressing the GFP-CesA6 construct of the fourth generation, F 3 1 line [23], have been used in this study. Nicotiana tabacum cv SR1 (wt) and F 3 1 tobacco seeds were germinated in a growth chamber (temperature 22 • C, humidity 60%, light intensity 25 µE/m 2 /s, photoperiod 16/8 h day/night) on solid MS (Murashige & Skoog, Sigma-Aldrich, St. Louis, MO, USA) medium and transplanted to the soil after two months.

Microscopic Observations and Chlorophyll and Carotenoid Content
Morphometric parameters were measured using ImageJ 1.53e from seed images acquired on a stereomicroscope (Stemi 508, Zeiss, Germany) and on leaves scanned with an EPSON XP-255 scanner. The epidermal pavement cell area was analyzed following the method of Li et al. [26] with some modifications. Epidermal peels were taken from the abaxial lamina of the most expanded leaves of 2-month-old wt and F 3 1 tobacco seedlings and gently pressed onto the slides. Observations were performed under a LSM710 microscope (Zeiss, Germany) at 40× magnification. The area value was measured using the Graphic tool of the Zen2012 Black Edition program.
For stem anatomy observations, the first internode of 2-month-old wt and F 3 1 seedlings was fixed overnight in cold (4 • C) FAA solution (Ethanol 50%, Glacial acetic acid 0.5%, Formaldehyde [37-40%] 10%, water up to the final volume) and then embedded in paraffin. Thin cross-sections (10 µm) were cut with a microtome (Leica RM2155) and picked up on a slide. Clearing procedure, fuchsin A (0.2%) and calcofluor white (0.1%) staining were performed as described by Ursache et al. [27]. Fuchsin A was detected at 561 nm excitation and acquired at 600-650 nm; calcofluor white was detected at 405 nm excitation and acquired at 425-475 nm.
The spectrophotometric determination of leaf content of chlorophylls (a, b, and total) and total carotenoids was carried out according to the Lichtenthaler and Buschmann (2001) method [28] and reported as nmol/g fresh weight (fw).

RNA Extraction and RT-qPCR Analysis
Total RNA was isolated from 2-month-old wt and F 3 1 leaves and stems as previously reported [29]. RT-qPCR analysis was performed using SYBR Green fluorescent detection in a CFX96 Real-Time System Cycler (Bio-Rad, Hercules, CA, USA) with three biological and three technical replicates per sample. The primer sequences used are reported in Table S1. The PCR program was as follows: 10 min at 95 • C; 50 cycles of 15 s at 95 • C; 20 s at 60 • C; and an increment of 0.5 • C every 0.5 s from 65 • C to 95 • C. The specificity of PCR products was checked in a melting-curve test. The Nicotiana tabacum Elongation factor EF-1α (EF-1α, AF120093), L25 ribosomal protein (L25, L18908), and Ubiquitin-conjugating enzyme E2 (Ntubc2, AB026056) genes were tested as reference genes [30]. In wt and F 3 1 samples, all genes had a variation coefficient below 0.1, according to Czechowski et al. [31]. The EF-1α gene was chosen as the reference gene because it showed the lowest variation coefficient among the genes (Table S2). Differences in gene expression between transformed and wt samples were considered significant when the expression was at least doubled (greater than or equal to two-fold upregulation) or halved (less than or equal to two-fold downregulation), according to Chen et al. [32].

Leaf Proteins Fractionations and Immunoblotting Analyses
Three wt and F 3 1 leaves were finely powdered in a pre-cooled mortar with liquid nitrogen and 100 mg of the powder were suspended in the homogenization buffer (40 mM Hepes-NaOH buffer, pH 7.5, 10 mM imidazole, 1 mM benzamidine, 5 mM 6-aminohexanoic acid, 10 mM dithiothreitol and 1 mM phenyl methyl sulfonyl fluoride) and centrifuged at 800× g (Allegra 2IR, Beckman, Brea, CA, USA) for 10 min at 4 • C. The resulting pellet, consisting of cell walls and cell debris, was sequentially washed and extracted to isolate the cell wall protein fraction (CW) according to De Caroli et al. [29]. The 800 g supernatant, containing the protein intracellular fraction, including soluble cargo and intrinsic membrane proteins, was layered on a Ficoll (40%) pad and ultracentrifuged at 100,000× g for 2 h at 4 • C (Optima XPN-1000, Beckman). The proteins in the formed interface were collected and delipidated to obtain the intrinsic membrane proteins fraction (MEM); the supernatant was concentrated by filtration on Centricon Plus 20 (Amicon, Merck KGaA, Darmstadt, Germany) to obtain the soluble cargo proteins (SOL) [29]. Protein concentration was determined as described by Bradford [33]. An identical amount of SOL, MEM, and CW protein fractions (10 µg per line) was subjected to SDS-PAGE and Western blotting analysis. Anti-GFP (1:5000 v/v) (Molecular Probes, now ThermoFisher Scientific, Waltham, MA, USA) and anti-XTH (Agrisera, Vännäs, Sweden) were used in TBS + 1% skimmed-milk powder. All primary antibodies were combined with anti-rabbit secondary antibodies coupled to peroxidase (1:10,000 v/v) (Sigma-Aldrich, St. Louis, MO, USA). Chemiluminescence reaction occurred using Clarity or Clarity Max (Biorad, Hercules, CA, USA) and the signal was captured by ChemiDoc MP Imaging Systems (Biorad). Densitometric analysis was performed using Image Lab software 6.0.1 (Biorad).

Cell Wall Isolation, Polysaccharide Fractionation and Analyses
Cell walls were isolated, purified (delipidated, destarched, and deproteinated), and essiccated as described by Leucci et al. [34] starting from 10 g fw aliquots of wt and F 3 1 leaves from 2-month-old wt and F 3 1 seedlings. The chemical solubilization of cell wall polysaccharides was carried out according to Lenucci et al. [35]. Aliquots (100 mg) of the purified cell walls were sequentially extracted with (a) 50 mM cyclohexane-trans-1,2-diamine-N,N,N ,N tetraacetate (CDTA, Na salt, 10 mL), pH 6.5, at 20 • C for 6 h with shaking (CDTA-1); (b) 50 mM CDTA (10 mL), pH 6.5, at 20 • C for 2 h with shaking (CDTA-2); (c) 50 mM Na 2 CO 3 + 20 mM NaBH 4 (10 mL) at 1.0 • C for 16 h with shaking (Na 2 CO 3 -1); (d) 50 mM Na 2 CO 3 + 20 mM NaBH 4 (10 mL) at 20 • C for 4 h (Na 2 CO 3 -2); (e) 0.5 M KOH + 10 mM NaBH 4 (10 mL) at 1.0 • C for 3 h, under nitrogen; (f) 1 M KOH + 10 mM NaBH 4 (10 mL) at 1 • C for 6 h, under nitrogen; and (g) 4 M KOH + 10 mM NaBH 4 (10 mL) at 20 • C for 24 h, under nitrogen. At the end of the sequential extractions, an insoluble residue, usually considered α-cellulose, remained. After each extraction, the soluble polymers were separated from the insoluble residue by centrifugation (5000× g, 10 min in a Biofuge 15 Heraeus Sepatech centrifuge). The combined alkali extracts and the insoluble residue were acidified to pH 5.0 with glacial acetic acid. The CDTA-1, CDTA-2, Na 2 CO 3 -1, and Na 2 CO 3 -2 extracts were combined and referred to as CDTA + Na 2 CO 3 extract. Similarly, the 0.5, 1.0, and 4.0 M KOH extracts were combined and referred to as KOH extract. The CDTA + Na 2 CO 3 and KOH extracts as well as the insoluble residue were dialyzed exhaustively in the presence of 0.05% chlorbutol (with SnakeSkin Pleated Dialysis Tubing, Pierce; Mr cutoff 10,000) against several changes of distilled water for 7 days at 5.0 • C in a cold room, freeze-dried, and accurately weighted. Aliquots (2 mg) of the freeze-dried CDTA + Na 2 CO 3 and KOH extracts were subjected to acid hydrolysis with 2N trifluoroacetic acid (TFA), at 121 • C for 60 min. For CDTA + Na 2 CO 3 , a pretreatment with 90% formic acid at 100 • C for 6 h was required to improve the release of sugar residues from pectins. The hydrolyzates were centrifuged at 10,000× g and the supernatants were taken to dryness in a SpeedVac concentrator (Savant). The amount of cellulose present in the insoluble residue was quantified by the Updegraff [36] method.
To perform 1 H NMR analysis and spectra acquisition and processing, appropriate aliquots of the dry hydrolyzates (four independent replicas) of wt and F 3 1 transgenic line for both pectin and hemicellulose extracts were processed for 1 H NMR analysis. For each sample, 600 µL of deuterium oxide (D 2 O) containing 0.05% w/v TSP-d4 (sodium salt of trimethylsilyl propionic acid) as a chemical shift reference, δ = 0 ppm, was added to the lyophilized plant material. The content was mixed thoroughly with a vortex mixer at room temperature for 1 min and then filled in a 5 mm NMR tube. All measurements were performed at 300 K on a Bruker Avance III 600 Ascend NMR spectrometer (Bruker, Ettlingen, Germany), operating at 600.13 MHz for 1 H observation, equipped with a TCI cryoprobe incorporating a z axis gradient coil and automatic tuning matching (ATM). For each sample, a 1D sequence with pre-saturation and composite pulse for selection (zgcppr Bruker standard pulse sequence) was acquired, with 128 scans, 16 dummy scans, 5 s relaxation delay, size of fid of 64 K data points, a spectral width of 12,019.230 Hz (20.0276 ppm), and an acquisition time of 2.73 s. The resulting FIDs were multiplied by an exponential weighting function corresponding to a line broadening of 0.3 Hz before Fourier transformation, automated phasing, and baseline correction. Metabolite identifications were based on 1 H and 13 C assignment by 1D and 2D omo and heteronuclear experiments (2D 1 H Jres, 1 H COSY, 1 H-13 C HSQC, and HMBC) and by comparison with the literature data [37][38][39][40][41][42]. NMR data processing was performed by using TopSpin 3.6.1 (Bruker, Biospin, Milano, Italy). All spectra were referenced to the TSP signal (0.00 ppm). The change in the composition of selected components, identified by NMR, between wt and F 3 1 transgenic line hemicelluloses and pectins was determined by analyzing the integrals of selected distinctive unbiased NMR signals after spectra processing, using TSP for chemical shift calibration and quantification [43], after normalization to the amount of hydrolyzate used in each analysis. Differences in the metabolite content were represented as log2foldchange of the calculated average intensities of the corresponding selected signals [44,45].
A TLC analysis was also carried out to visualize the glycosidic residues from the pectin and hemicellulose extracts of wt and transformed leaves. The dried hydrolyzates were solubilized in pyridine and appropriate aliquots (typically 10 µL) deposited on a 20 cm × 20 cm TLC plate (Merck KGaA, Darmstadt, Germany) together with a standard solution containing 5µg/mL of rhamnose, arabinose, galactose, glucose, xylose, mannose, and galacturonic acid. The chromatographic run and the sugar staining were performed according to the method reported by O'Rourke et al. [46]. The ride lasted three hours using the following mobile phase: ethyl acetate:pyridine:acetic acid:water (6:3:1:1, v/v). The sugars were visualized by immersing the plate in a 0.5% thymol mixture in ethanol:sulfuric acid (19: 1 v/v) followed by incubation in a pre-heated oven at 105 • C for 5 min.

Protoplast Preparation, Radioactive Labeling and Assay of De Novo Synthesized Cell Wall Polysaccharides
The de novo synthesis of cell wall polysaccharides was assessed on leaf tobacco protoplasts prepared from 2-month-old wt and F 3 1 leaves as described by De Caroli et al. [47]. Freshly prepared protoplasts were incubated in 1 mL of K3M (K3 solution in which sucrose was replaced by mannitol) medium (26 h) in the presence of 430 kBq of D-[U-14 C]glucose (Amersham Biosciences, Amersham, UK) for the last 8 h of incubation. At the end of the incubation period, the protoplasts were treated as reported in Leucci et al. [48]. Furthermore, to remove the possible presence of callose, the cell walls were incubated with 1,3-β-endoglucanase (barley, Megazyme, Wicklow, Ireland) at 40 • C for 18 h and the residual material was centrifuged at 800× g for 10 min at 4 • C to collect purified cell walls. The quantification of newly synthesized matrix polysaccharides and cellulose was performed after nitric acid/acetic acid hydrolysis (Updegraff,[36]) of the protoplasts' regenerating cell walls and counting for radioactivity of both the hydrolyzates (matrix polysaccharides) and the acid resistant material (cellulose).

Statistical Analyses
Statistical analyses were based on the t-Student test. Statistical comparisons were performed using SigmaStat software, version 11.0 (Systat Software Inc., Chicago, IL, USA).

Expression of the Exogenous GFP-CesA6 Construct in Tobacco Plants
At the same growth stage, the five transformed lines of tobacco expressing the exogenous GFP-CesA6 gene were bigger than wt, a character clearly evident already a few weeks after germination. When observed at the confocal microscope, all lines showed similar fluorescence intensity indicating no differences in the expression of GFP-CesA6 [23].
Confocal observations were here confirmed through RT-qPCR analyses revealing comparable expression levels of AtCesA6 gene in all transformed lines ( Figure 1a). Western blot of the membrane protein fractions purified from the leaves of 2-month-old transformed seedlings showed the presence of bands with similar protein level in the five lines (Figures 1b and S1). For subsequent investigations aimed at studying the potential effects of the construct expression on plant growth and cellulose synthesis, among the five lines, we chose the transgenic line F 3 1 as representative of the fourth generation.
1,3-β-endoglucanase (barley, Megazyme, Wicklow, Ireland) at 40 °C for 18 h and th sidual material was centrifuged at 800× g for 10 min at 4 °C to collect purified cell w The quantification of newly synthesized matrix polysaccharides and cellulose was formed after nitric acid/acetic acid hydrolysis (Updegraff,[36]) of the protoplasts' re erating cell walls and counting for radioactivity of both the hydrolyzates (matrix pol charides) and the acid resistant material (cellulose).

Statistical Analyses
Statistical analyses were based on the t-Student test. Statistical comparisons were formed using SigmaStat software, version 11.0 (Systat Software Inc., Chicago, IL, US

Expression of the Exogenous GFP-CesA6 Construct in Tobacco Plants
At the same growth stage, the five transformed lines of tobacco expressing the e enous GFP-CesA6 gene were bigger than wt, a character clearly evident already a weeks after germination. When observed at the confocal microscope, all lines showed ilar fluorescence intensity indicating no differences in the expression of GFP-CesA6 Confocal observations were here confirmed through RT-qPCR analyses revealing co rable expression levels of AtCesA6 gene in all transformed lines ( Figure 1a). Western of the membrane protein fractions purified from the leaves of 2-month-old transfor seedlings showed the presence of bands with similar protein level in the five lines (Fig  1b and S1). For subsequent investigations aimed at studying the potential effects o construct expression on plant growth and cellulose synthesis, among the five lines chose the transgenic line F31 as representative of the fourth generation.

Morphometric and Growth Parameters of F 3 1 Seeds and Seedlings
No morphological or morphometric differences were observed between wt and F 3 1 mature dry seeds by stereomicroscope and confocal microscopy. In both cases, seeds appeared oblong-obovate with an eccentric hilum and an average area, perimeter, major and minor axis length of about 0.31 mm 2 , 2.20 mm, 0.72 mm, and 0.55 mm, respectively (Figure 2a,b). The outer integument of the seed coat was brown pigmented with a striate texture (Figure 2a,b) and constituted of normally expanded jigsaw-puzzle-shaped cells, as clearly visualized by the autofluorescence of their thickened cell walls (Figure 2c,d).
Biology 2022, 11, x 7 of 23 target according to Schmittgen and Livak [49]. (b) Western blot analyses of membrane protein fractions in leaves of 2-month-old wt and transformed tobacco seedlings stably expressing GFP-CesA6; bands were detected by anti-GFP. The results of three independent biological and three technical replicates are presented with box plots.

Morphometric and Growth Parameters of F31 Seeds and Seedlings
No morphological or morphometric differences were observed between wt and F31 mature dry seeds by stereomicroscope and confocal microscopy. In both cases, seeds appeared oblong-obovate with an eccentric hilum and an average area, perimeter, major and minor axis length of about 0.31 mm 2 , 2.20 mm, 0.72 mm, and 0.55 mm, respectively ( Figure  2a,b). The outer integument of the seed coat was brown pigmented with a striate texture (Figure 2a,b) and constituted of normally expanded jigsaw-puzzle-shaped cells, as clearly visualized by the autofluorescence of their thickened cell walls (Figure 2c,d). Although transformation did not affect neither the rate (>98%) nor the seed germination time (4 days after sowing), a slight increase in size was macroscopically noticeable in Although transformation did not affect neither the rate (>98%) nor the seed germination time (4 days after sowing), a slight increase in size was macroscopically noticeable in 30-day-old F 3 1 seedlings compared to wt, a difference which became clearly apparent 2 months after sowing and was remarkable in the 5-month-old plants (Figure 3a-c).
logy 2022, 11, x 8 of 30-day-old F31 seedlings compared to wt, a difference which became clearly apparent months after sowing and was remarkable in the 5-month-old plants (Figure 3a-c).  Table 1). A significant (p < 0.05) increase in th average of both stem height (from the root collar to the shoot apex) and root length about 61% and 135%, respectively, was evidenced in the 2-month-old F31 seedlings com pared to wt, as well as in the number of true leaves (67%) and roots (70%), and in lea root, and stem fresh weights (fw, 72%, 307%, and 113%, respectively) ( Figure 4). Table 1. Average values of parameters of wt and F31 tobacco flowers and seeds; the data are t mean ± SD of 10 measurements on each line. Asterisk indicates significant differences (p < 0.05; according to the t-Student test, between the two samples. wt F31 Anthesis 28 weeks ± 5 days 25 weeks ± 3 days * Flower numbers 9 ± 1 13 ± 2 * Seed numbers 9659 ± 147 11.516 ± 250 * Seed total weight (g) 0.749 ± 0.047 0.979 ± 0.093 *  Table 1). A significant (p < 0.05) increase in the average of both stem height (from the root collar to the shoot apex) and root length of about 61% and 135%, respectively, was evidenced in the 2-month-old F 3 1 seedlings compared to wt, as well as in the number of true leaves (67%) and roots (70%), and in leaf, root, and stem fresh weights (fw, 72%, 307%, and 113%, respectively) ( Figure 4). Table 1. Average values of parameters of wt and F 3 1 tobacco flowers and seeds; the data are the mean ± SD of 10 measurements on each line. Asterisk indicates significant differences (p < 0.05; *), according to the t-Student test, between the two samples.

wt F 3 1
Anthesis 28 weeks ± 5 days 25 weeks ± 3 days * Flower numbers 9 ± 1 13 ± 2 * Seed numbers 9659 ± 147 11.516 ± 250 * Seed total weight (g) 0.749 ± 0.047 0.979 ± 0.093 * Moreover, when the morphometric parameters (area, perimeter, major and minor axes, and circularity) of embryonal (cotyledons) and true leaves (progressively numbered according to their emergence order) were measured, all traits, except circularity, appeared significantly increased in the 2-month-old F 3 1 seedlings with respect to wt ( Figure 5). Moreover, when the morphometric parameters (area, perimeter, major and minor axes, and circularity) of embryonal (cotyledons) and true leaves (progressively numbered according to their emergence order) were measured, all traits, except circularity, appeared significantly increased in the 2-month-old F31 seedlings with respect to wt ( Figure 5). An amount of 30 leaves for each sample were examined for each line. The data represent the mean ± SD. Asterisk indicates significant differences (p < 0.05; *), according to the t-Student test, between the two samples.
Remarkably, the two cotyledonary leaves, emerging during seed germination, also appeared more than twice as large in the transgenic line than in wt. In the leaf numbered 2 ( Figure 5), the size and shape of wt and F31 abaxial epidermal cells were compared with the corresponding equal areas of proximal and distal zone. In the proximal zone, the cells  Moreover, when the morphometric parameters (area, perimeter, major and minor axes, and circularity) of embryonal (cotyledons) and true leaves (progressively numbered according to their emergence order) were measured, all traits, except circularity, appeared significantly increased in the 2-month-old F31 seedlings with respect to wt ( Figure 5). Remarkably, the two cotyledonary leaves, emerging during seed germination, also appeared more than twice as large in the transgenic line than in wt. In the leaf numbered 2 ( Figure 5), the size and shape of wt and F31 abaxial epidermal cells were compared with the corresponding equal areas of proximal and distal zone. In the proximal zone, the cells Remarkably, the two cotyledonary leaves, emerging during seed germination, also appeared more than twice as large in the transgenic line than in wt. In the leaf numbered 2 ( Figure 5), the size and shape of wt and F 3 1 abaxial epidermal cells were compared with the corresponding equal areas of proximal and distal zone. In the proximal zone, the cells showed the same pattern and approximately constant average cell area; in the distal zone, whereas the shape does not change, the area of the F 3 1 epidermal cells significantly increased by expansion growth with respect to the wt ( Figure 6).
showed the same pattern and approximately constant average cell area; in the dis whereas the shape does not change, the area of the F31 epidermal cells signific creased by expansion growth with respect to the wt ( Figure 6). On a fw basis, the leaves of 2-month-old F31 seedlings showed a slightly content of chlorophyll a (Chl a, 123.9 nmol/g fw) and total chlorophylls (Chl (a+ nmol/g fw) than wt (143.1 and 181.9 nmol/g fw, respectively). Instead, no signific iations were found in the levels of chlorophyll b (Chl b) and total carotenoids ( Furthermore, the ratio Chl a:Chl b slightly increased in F31 (3.9) with respect t (3.7), while the ratio between Chl (a+b) and total carotenoids (Car) revealed a c crease in the F31 (2.86) than wt (3.53) leaves. Thus, apart from the "giant" phenotype and larger fresh biomass productio errant morphological traits (i.e., curled or asymmetrical leaves, necrotic or differe mented areas, defective flower morphology) were observed in F31 plants. On a fw basis, the leaves of 2-month-old F 3 1 seedlings showed a slightly reduced content of chlorophyll a (Chl a, 123.9 nmol/g fw) and total chlorophylls (Chl (a+b), 155.6 nmol/g fw) than wt (143.1 and 181.9 nmol/g fw, respectively). Instead, no significant variations were found in the levels of chlorophyll b (Chl b) and total carotenoids (Table 2). Furthermore, the ratio Chl a:Chl b slightly increased in F 3 1 (3.9) with respect to the wt (3.7), while the ratio between Chl (a+b) and total carotenoids (Car) revealed a clear decrease in the F 3 1 (2.86) than wt (3.53) leaves. Table 2. Chlorophyll and carotenoid content in wt and F 3 1 tobacco leaves. The data represent the mean ± SD of four independent experiments. The asterisks indicate significant (p < 0.05; *) or highly significant (p < 0.001, **) differences, according to the t-Student test, between the two samples. Thus, apart from the "giant" phenotype and larger fresh biomass production, no aberrant morphological traits (i.e., curled or asymmetrical leaves, necrotic or differently pigmented areas, defective flower morphology) were observed in F 3 1 plants.

Exogenous GFP-CesA6 Increases the Expression of Endogenous Primary CesAs and of Other Genes Involved in the Biosynthesis and Remodeling of Cell Wall Polymers in Tobacco Leaves
To investigate whether the enhanced growth of F 3 1 seedling was supported by a general overexpression of cell wall synthesis-, deposition-, and/or remodeling-related genes, we performed a RT-qPCR on RNA samples isolated from the leaves of 2-monthold transformed and wt seedlings. In particular, we examined the expression levels of the endogenous CesAs involved in cellulose synthesis of primary (NtCesA1, NtCesA3, and NtCesA6) and secondary (NtCesA4, NtCesA7, and NtCesA8) cell walls, as well as of other genes pivotal for matrix polysaccharide (pectins and hemicelluloses) metabolism and cell wall expansive growth (namely NtXyl4, NtGal10, NtExp11, NtPIP1.1, and NtPIP2.7). In particular, NtXyl4 encodes for an UDP-glucose:xyloglucan 1,4-β-Dglucosyltransferase putatively involved in backbone elongation of xyloglucans; the product of NtGal10 is a protein with putative galacturonosyltransferase activity catalyzing pectin biosynthesis, while NtExp11 codes for an expansin-A11 causing loosening and extension of plant cell walls disrupting the non-covalent bonding between cellulose microfibrils and matrix glucans. NtPIP1.1 and NtPIP2.7 code for two plasma membrane aquaporins involved in roles typically reported for PIPs across plants species, including leaf cell expansion [50].
Compared to wt, F 3 1 seedlings showed an overexpression of NtCesA1 (log2Foldchange between 3 and 4), NtCesA3, and NtCesA6 (log2Foldchange between 2 and 3) but not of NtCesA4, NtCesA7, and NtCesA8. NtGal10 and NtExp11 were also both strongly upregulated (log2Foldchange between 3 and 4), while the expression level of NtXyl4 was only slightly above the typical significance threshold of 2 for log2Foldchange. NtPIP1.1 and NtPIP2.7 were also upregulated with a log2Foldchange between 2 and 3 ( Figure 7). and secondary (NtCesA4, NtCesA7, and NtCesA8) cell walls otal for matrix polysaccharide (pectins and hemicelluloses) pansive growth (namely NtXyl4, NtGal10, NtExp11, NtPIP1. NtXyl4 encodes for an UDP-glucose:xyloglucan 1,4-β-D-glu volved in backbone elongation of xyloglucans; the produc putative galacturonosyltransferase activity catalyzing pecti codes for an expansin-A11 causing loosening and extension the non-covalent bonding between cellulose microfibrils and NtPIP2.7 code for two plasma membrane aquaporins ported for PIPs across plants species, including leaf cell exp Compared to wt, F31 seedlings showed an overexpr change between 3 and 4), NtCesA3, and NtCesA6 (log2Fold not of NtCesA4, NtCesA7, and NtCesA8. NtGal10 and NtExp regulated (log2Foldchange between 3 and 4), while the ex only slightly above the typical significance threshold of 2 and NtPIP2.7 were also upregulated with a log2Foldchange  Additionally, Western blot analyses of protein extracts from purified soluble, membrane, and wall fractions isolated from the leaves of 2-month-old tobacco wt and F 3 1 seedlings were performed to get direct indications on the total abundance of the xyloglucan endotransglicosylase/hydrolases (XTHs). XTHs are a class of enzymes with a role in the assembly, growth, and remodeling of plant cell walls catalyzing the molecular grafting between xyloglucans of the hemicellulose-cellulose network [29,51]. Immunoblots with the anti-XTH serum revealed a single protein band present exclusively in the cell wall fraction with an approximate MW of 40 KDa compatible with that of the target proteins and with intensity about 4-fold higher in the transformed line than in wt (Figures 8 and S2). quantification of XTH abundance, reported as average of pixel numb sented with box plots. Immunolocalization was carried out with an anti-XTH serum. As dicate highly significant differences (p < 0.001; **), according to the t-Student test, betwee samples.

In the Stem of GFP-CesA6 Transformed Tobacco Seedlings the Secondary Xylem T Was Enhanced
From the transmitted light confocal observations of transverse sections (10 µ of first internode of 2-month-old seedlings, no anomalies in the histological eleme revealed (Figure 9a,b). The same sections, stained with fuchsin A and calcoflu (specific for lignin and cellulose, respectively), showed a 3.3-fold increase of the s xylem thickness in the F31 line (874 ± 45 µ m) compared to wt (262 ± 38 µ m) (Figu Gene expression analyses of the stem first internode revealed differentially prim secondary CesAs expression profiles (Figure 9f). While the expression of the NtCesA1 and NtCesA3 was significantly reduced with respect to the wt (log2Fo between −2 and −4), the NtCesA6 expression was increased (log2Foldchange a With regard to the secondary NtCesAs, the expression level of NtCesA8 and NtCes increased (log2Foldchange around 2 for both genes) while NtCesA4 expression duced (log2Foldchange between −2 and −3).

In the Stem of GFP-CesA6 Transformed Tobacco Seedlings the Secondary Xylem Thickness Was Enhanced
From the transmitted light confocal observations of transverse sections (10 µm thick) of first internode of 2-month-old seedlings, no anomalies in the histological elements were revealed (Figure 9a,b). The same sections, stained with fuchsin A and calcofluor white (specific for lignin and cellulose, respectively), showed a 3.3-fold increase of the secondary xylem thickness in the F 3 1 line (874 ± 45 µm) compared to wt (262 ± 38 µm) (Figure 9c-e). Gene expression analyses of the stem first internode revealed differentially primary and secondary CesAs expression profiles (Figure 9f). While the expression of the primary NtCesA1 and NtCesA3 was significantly reduced with respect to the wt (log2Foldchange between −2 and −4), the NtCesA6 expression was increased (log2Foldchange around 2). With regard to the secondary NtCesAs, the expression level of NtCesA8 and NtCesA7 genes increased (log2Foldchange around 2 for both genes) while NtCesA4 expression was reduced (log2Foldchange between −2 and −3).

The Expression of the GFP-CesA6 Construct Affects the Relative Amounts of Cell Wall Polysaccharides
To understand if the stable transformation with the GFP-CesA6 construct was associated with alterations in the levels of cellulose and/or matrix polysaccharides, the cell walls from the leaves of 2-month-old seedlings (wt and F 3 1) were isolated, purified, and sequentially extracted into three fractions enriched in pectins (CDTA + Na 2 CO 3 extracts), hemicelluloses (0.5-4.0 M KOH extracts), and cellulose (insoluble residue). The dry weight (dw) of the purified cell walls, expressed in mg/g fw, was significantly different (p < 0.05) between F 3 1 (17.6 ± 1.3) and wt (7.3 ± 0.6), with up to a 2.4-fold increase (Table 3).

The Expression of the GFP-CesA6 Construct Affects the Relative Amounts of Cell Wall Polysaccharides
To understand if the stable transformation with the GFP-CesA6 construct was associated with alterations in the levels of cellulose and/or matrix polysaccharides, the cell walls (e) Quantification of secondary xylem in the stem of wt and F 3 1 tobacco seedlings. The measurements are reported with box thickness plot graph. Asterisks indicate highly significant differences (p < 0.001; **), according to the t-Student test, between the two samples. (f) RT-qPCR analyses of NtCesA1, NtCesA3, NtCesA6, NtCesA4, NtCesA7, and NtCesA8 genes in the stem of 2-month-old wt and F 3 1 tobacco seedlings. The gene expression is reported as transcript inhibition level in F 3 1 leaves (log2foldchange) with respect to wt. The results of three independent biological and three technical replicates are presented with box plots. The horizontal red dashed lines indicate the significant threshold of log2foldchange. Table 3. Weights of the purified cell walls and extracts obtained from the leaves of wt and F 3 1 2-month-old tobacco seedlings. The data, expressed in mg/g fw, represent the mean ± standard deviation of four independent experiments. Percentage values refer to the sum of total weight of extracts. In parenthesis are reported the amounts of cellulose quantified in the insoluble residues by the Updegraff method [36]. Asterisks indicate significant (p < 0.05; *) or highly significant (p < 0.001 **) differences, according to the t-Student test, between wt and F 3 (3.29) Moreover, the weight of CDTA + Na 2 CO 3 and KOH extracts obtained from F31, expressed as mg dw/g tissue fw, were slightly higher (10% and 23%, respectively) than those obtained from wt leaves, while the weight of the insoluble residue was increased over 3 times in the transformed line (15.0 ± 2.28 mg/g fw) compared to wt (4.18 ± 0.56). In percentage terms, pectins and hemicelluloses of wt leaves contributed to about 23% and 11% of the total dw of the cell wall, while the remaining 66% consists of the insoluble residue, which turned out to contain about 18.2% of crystalline cellulose as estimated by the Updegraff method [36]. The weight of the insoluble residue from F31 was much higher (84.9%, containing 23.7% crystalline cellulose), at the expense of pectins and hemicelluloses, which contributed, respectively, to 9.9% and 5.2% of the purified cell wall dw.
Significant differences were also detected calculating the percentage of pectins, hemicelluloses, and cellulose relative to their total weight (Table 4). In the walls of wt seedling leaves, pectins represented the prevalent cell wall polysaccharide (44.2%), followed by cellulose (35.2%) and hemicelluloses (20.6%). Transformation, instead, determined a considerable increase in the percentage of cellulose, which, in the F 3 1 leaves, contributes to 61.1% of the total weight of cell wall polysaccharides, followed by pectins (25.6%) and hemicelluloses (13.3%). Table 4. Percentage distribution of pectins, hemicelluloses, and cellulose in the purified cell walls from the leaves of wt and F 3 1 2-month-old tobacco seedlings. The percentage values refer to the total weight of the polysaccharide material expressed in mg/g of dry weight (dw) of purified cell walls and represent the mean of four independent experiments. The asterisk (*) indicates significant (p < 0.05) difference, according to the t-Student test, between wt and F 3 1 samples. For a rough indication of the glycosidic composition of the structural polysaccharides, the CDTA + Na 2 CO 3 and KOH extracts were subjected to acid hydrolysis, and the resulting products were assayed both by 1 H-NMR and thin layer chromatography (TLC). Both analytical techniques showed qualitatively similar profiles of the wt and F 3 1 hydrolyzates. The relative concentrations of the various metabolites in the different 1 H-NMR spectra were calculated by integrating the area of selected distinctive unbiased NMR signal. In particular, signals at 7.14 ppm (ferulic acid), 5.32 ppm (α-galacturonic acid monomers), 5.22 ppm (galacturonic acid oligomers), 5.10 ppm (rhamnose), 4.64 ppm (β-galacturonic acid monomers), 3.30 ppm (methoxy groups), and 2.22 ppm (acetyl groups) were selected and integrated for pectin wt and F 3 1 samples (Figure S3a,b). Moreover, signals at 5.25 ppm (α-fucose), 5.22 ppm (α-glucose, 5.18 ppm (α-xylose), 5.16 ppm (α-mannose), 4.88 ppm (β-mannose), 4.63 ppm (β-glucose), 4.56 ppm (β-xylose), and 4.50 ppm (β-fucose) were selected and integrated for hemicellulose wt and F 3 1 samples ( Figure S4). CDTA + NaC 2 O 3 hydrolyzates revealed the presence of signals attributable to galacturonic acid and rhamnose, as well as of a mixture of oligomeric fragments, confirming the low susceptibility of pectins to acid hydrolysis [52]. The presence of methoxy and acetyl groups, probably esterified to the carboxyl and hydroxyl groups of the C2 and/or C3 of galacturonic acid residues [53], and ferulic groups, known to contribute to the formation of cross-links between the pectic polysaccharides and between them and structural wall proteins [54], were also evidenced. Mannose, glucose, xylose, and fucose were the main glycosyl residues identified in the KOH extracts ( Figure S4). The data, expressed as log2foldchange, showed a general trend towards a reduction of the areas of 1 H-NMR signals of galacturonic acid, rhamnose, oligogalacturonides, methoxy, acetyl, and ferulic groups between F 3 1 and wt lines, although most of these variations (with the exception of ferulic acid) were not significant at the chosen threshold (+2, −2 log2foldchange) (Figure 10a). On the contrary, an upward trend was found for glucose and xylose in the KOH extracts (Figure 10b).

Relative Percent (%)
lyzates. The relative concentrations of the various metabolites in the different 1 H-NMR spectra were calculated by integrating the area of selected distinctive unbiased NMR signal. In particular, signals at 7.14 ppm (ferulic acid), 5.32 ppm (α-galacturonic acid monomers), 5.22 ppm (galacturonic acid oligomers), 5.10 ppm (rhamnose), 4.64 ppm (β-galacturonic acid monomers), 3.30 ppm (methoxy groups), and 2.22 ppm (acetyl groups) were selected and integrated for pectin wt and F31 samples ( Figure S3a,b). Moreover, signals at 5.25 ppm (α-fucose), 5.22 ppm (α-glucose, 5.18 ppm (α-xylose), 5.16 ppm (α-mannose), 4.88 ppm (β-mannose), 4.63 ppm (β-glucose), 4.56 ppm (β-xylose), and 4.50 ppm (β-fucose) were selected and integrated for hemicellulose wt and F31 samples ( Figure S4). CDTA + NaC2O3 hydrolyzates revealed the presence of signals attributable to galacturonic acid and rhamnose, as well as of a mixture of oligomeric fragments, confirming the low susceptibility of pectins to acid hydrolysis [52]. The presence of methoxy and acetyl groups, probably esterified to the carboxyl and hydroxyl groups of the C2 and/or C3 of galacturonic acid residues [53], and ferulic groups, known to contribute to the formation of cross-links between the pectic polysaccharides and between them and structural wall proteins [54], were also evidenced. Mannose, glucose, xylose, and fucose were the main glycosyl residues identified in the KOH extracts ( Figure S4). The data, expressed as log2foldchange, showed a general trend towards a reduction of the areas of 1 H-NMR signals of galacturonic acid, rhamnose, oligogalacturonides, methoxy, acetyl, and ferulic groups between F31 and wt lines, although most of these variations (with the exception of ferulic acid) were not significant at the chosen threshold (+2, −2 log2foldchange) ( Figure  10a). On the contrary, an upward trend was found for glucose and xylose in the KOH extracts ( Figure 10b). The TLC ( Figure S5) of the glycosidic residues deriving from the CDTA + Na 2 CO 3 extracts obtained from wt and transformed leaves confirmed the presence, in both samples and with no substantial differences, of galacturonic acid and rhamnose, as well as of glucose, xylose, and mannose. Similarly, the hydrolyzates from the extracts in KOH does not show evident differences in the glycosidic composition of the hemicelluloses related to the transformation. In agreement with 1 H-NMR, the most represented sugars were glucose and xylose, main constituents of the xyloglucans, principal hemicellulose of dicotyledonous primary cell walls. The presence of spots corresponding to arabinose, mannose, and rhamnose residues was also evident.
The de novo synthesis of α-cellulose was also followed by radiolabeling of freshly prepared leaf protoplasts from 2-months-old seedlings in the presence of D-[U-14 C]glucose for 8 h. At the end of the incubation period, the regenerating cell walls were isolated and purified (destarched, deproteinated, and treated with 1,3-β-endoglucanase to remove eventual callose). The radioactivity incorporated into the insoluble residue resulting after acetic acid/nitric acid hydrolysis of the wt and F 3 1 purified cell walls, represented, respectively, 0.39% and 0.49% of the total radioactivity of the homogenate, corresponding to a 25.6% increase of cellulose biosynthesis in the F 3 1 protoplasts (Figure 11). not show evident differences in the glycosidic composition of the hemicelluloses rela to the transformation. In agreement with 1 H-NMR, the most represented sugars were g cose and xylose, main constituents of the xyloglucans, principal hemicellulose of dico ledonous primary cell walls. The presence of spots corresponding to arabinose, manno and rhamnose residues was also evident.
The de novo synthesis of α-cellulose was also followed by radiolabeling of fres prepared leaf protoplasts from 2-months-old seedlings in the presence of D-[U-14 C]g cose for 8 h. At the end of the incubation period, the regenerating cell walls were isola and purified (destarched, deproteinated, and treated with 1,3-β-endoglucanase to remo eventual callose). The radioactivity incorporated into the insoluble residue resulting af acetic acid/nitric acid hydrolysis of the wt and F31 purified cell walls, represented, resp tively, 0.39% and 0.49% of the total radioactivity of the homogenate, corresponding t 25.6% increase of cellulose biosynthesis in the F31 protoplasts ( Figure 11). Significantly, the radioactivity associated to the acetic acid/nitric acid hydrolyza (matrix polysaccharides) was also increased of approximately 21% in the F31 line, be their amount 2.8% and 3.4% in the wt and in the transformed seedling protoplasts, resp tively. These results once again showed the enhanced effect of the insertion of GFP-Ces gene on the whole cell wall synthesis and in general on the plant growth. Significantly, the radioactivity associated to the acetic acid/nitric acid hydrolyzates (matrix polysaccharides) was also increased of approximately 21% in the F 3 1 line, being their amount 2.8% and 3.4% in the wt and in the transformed seedling protoplasts, respectively. These results once again showed the enhanced effect of the insertion of GFP-CesA6 gene on the whole cell wall synthesis and in general on the plant growth.

Discussion
Cellulose, the main plant cell wall polysaccharide, is used in several industrial applications including paper, textile, food/feed, and bioethanol production [55]. In flowering plants, more than 80% of dry biomass consists of cell wall polymers, mainly polysaccharides, with cellulose accounting for approximately 35-45% of its composition [56]. Increasing cellulose biosynthesis, thus, is of considerable economic and social importance and represents an important goal to improve plant growth, biomass production, and wood quality and quantity. Among other experimental approaches, the overexpression of CesA genes has been appraised to increase cellulose production and, in general, to obtain biomass-producing crop plants. Economically important trees (poplar), as well as herbaceous (barley) and model plants (Arabidopsis), have been transformed for this purpose with contradicting results. Indeed, genetic manipulation of CesA members, particularly those involved in secondary wall biosynthesis, remains difficult to handle yet [19,20,57,58].
In this study, we demonstrate that the transgenic tobacco line F 3 1, expressing the exogenous GFP-CesA6 construct, shows an increased plant biomass, matrix and cellulose biosynthesis, and a "giant" phenotype, with no apparent other morphological aberrations. Compared to wt, F 3 1 seeds did not show variations neither in morphology nor in size ( Figure 2). Moreover, no significant differences were observed in seed germinability and sprout morphology. Changes in growth were visible already one month after sowing and, in a much more relevant way, at later growth stages (Figure 3). The 2-month-old F 3 1 seedlings exhibited higher values of all evaluated growth parameters (Figures 3 and 4), with an approximately 22% and 43% significant increase of stem and root length, respectively, as well as of leaves morphometric parameters. Significantly, aberrant forms were never detected during the development of stems, roots, and leaves, indicating that the F 3 1 line has no morphological or physiological anomalies. An anticipation of the anthesis and an increase of the average flower number per plants, total seed number, and weight ( Figure 3, Table 1) also characterized the F 3 1 line.
A relationship between photosynthetic activity and transgenic plant biomass increase was hypothesized in the poplar hybrid clone "Nanline895" (Populus deltoids × Populus euramericana) overexpressing CesA2 gene from Pinus massoniana under the control of constitutive 35S promoter. The transgenic plants showed increased chlorophyll content, secondary wall thickening, and xylem width, as well as improved cellulose synthesis, plant growth, and biomass accumulation [59]. Differently from transgenic poplar, the leaves of 2-month-old F 3 1 tobacco seedlings showed a slight but significant reduction in the content of chlorophyll a, the main pigment providing energy for oxygenic photosynthesis and emitting high energy electrons into the P680 and P700 photosystems. The amount of chlorophyll b and total carotenoids, both playing a role in light harvesting as accessory pigments, as well as the total content of chlorophylls remained instead unaffected.
The conflicting results are likely related to the different experimental systems investigated, i.e., expression of different genes involved in the biosynthesis of the primary cell wall (CesA6 and CesA2 in tobacco and poplar, respectively), use of diverse promoters (the endogenous CesA6 promoter in tobacco, the strong CaMV 35S constitutive promoter in poplar), and distinct growth habitus of the two species (herbaceous vs. arboreal). Indeed, the parameters related to chlorophylls concentrations significantly differed among plant functional groups (trees < shrubs < herbs) [60]. Moreover, several factors seem to regulate the expression of genes involved in chlorophyll biosynthesis. In N. benthamiana, for example, the expression of a DNA virus-encoded protein has been reported to cause structural and functional damages to chloroplasts and to inhibit photosynthesis [61].
The slight increase of Chl a: Chl b and Chl (a+b):Car ratios in the F 3 1 leaves suggests a process of light acclimatization possibly resulting from an indirect effect of transformation. Thus, taller F 3 1 plants, getting closer to the artificial light, activate adaptation and photoprotection mechanism reducing chlorophylls and simultaneously increasing carotenoid amounts [61]. Anyway, Chl a: Chl b ratio of both samples falls within the range of physiological variations, confirming that there is no indication of possible damage in F 3 1 plants [62,63].
Transformation with the exogenous GFP-CesA6 construct also enhanced the expression of endogenous primary NtCesAs (NtCesA1, NtCesA3, and NtCesA6) and of proteins not strictly related to cellulose biosynthesis but, more generally, responsible for plant cell wall building and modification (NtXyl4, NtGal10, NtExp11, NtPIP1.1, and NtPIP2.7). On the contrary, all cell wall secondary CesAs (NtCesA4, NtCesA7, and NtCesA8) were not significantly affected in the leaves of the F 3 1 line, indicating that the two sets of CesA genes are not related, at least in this organ. In Arabidopsis, the overexpression of CesA2, CesA5, and CesA6 was reported to enhance the expression of other primary CesAs and of the celluloserelated genes KORRIGAN2, SUS1, COBRA, CSI1, and CSI3 [22]. A simultaneous increase of cell division and expansive growth and of secondary cell wall deposition rate was also revealed, corroborating the idea that primary wall CesA engineering could affect plant growth in other ways than simply increasing cellulose biosynthesis. In close agreement with the report of Hu et al. [22], our results show that various components of the cell wall biosynthetic machinery are affected by the insertion of GFP-CesA6 construct, including those involved in the biosynthesis of non-cellulosic polysaccharides (xyloglucans and galacturonans) and in cell wall remodeling (i.e., NtExp11 and XTHs).
The expression of the GFP-CesA6 construct seems to also affect the relative amounts of cell wall polysaccharides. Indeed, a significant 2.4-fold increase of the purified cell wall material obtained per g of fw tobacco leaves was detected in the F 3 1 line. This was mainly attributable to an increase of the insoluble residue following the extraction of pectins and hemicelluloses and of the cellulose contained therein (Table 3). Interestingly, the contribution of total polysaccharides to the dw of the purified cell walls (Table 4) was slightly but significantly lower in the F 3 1 line than in the wt, suggesting that the biosynthesis of non-polysaccharide material may be stimulated by the transformation. Furthermore, the transformed plants had a considerable increase in the percentage of cellulose which, in the leaves, contributes to 61.1% of the total weight of the wall polysaccharides, followed by pectins (25.6%) and hemicelluloses (13.3%); in wt, pectins represent the prevalent polysaccharide (44.2%), followed by cellulose (35.2%) and hemicelluloses (20.2%). Similarly, a higher percentage of cellulose associated with a reduction of monosaccharides involved in the biosynthesis of hemicelluloses and pectins and a giant phenotype was previously reported by Hu et al. [22] in Arabidopsis plants overexpressing different CesA subunits. No significant difference emerged, instead, in the monosaccharide composition and in the degree of methylation and acetylation of the analyzed fractions which, in both wt and F 3 1 leaves, are compatible with the presence of pectic polysaccharides of the type of homogalacturonans and rhamnogalacturonans I and II, and of hemicelluloses of the type of xyloglucans characteristic of the type I primary wall of dicotyledonous Angiosperms. A significant decrease of polysaccharide feruloylation was also evidenced in the F 3 1 line, even though ferulic acid is reported to be a minor constituent of the tobacco leaf cell walls [64]. In addition, we also found that F 3 1 leaf protoplasts incubated in the presence of radioactive glucose showed an increased amount of cellulose biosynthesis and, significantly, also of matrix polysaccharides. Therefore, our data evidenced a correlation of the increased expression level of the cellulosic and non-cellulosic polysaccharides encoding genes to an actual increase of all the cell wall polysaccharides.
Although in leaves the exogenous expression of GFP-CesA6 had no effect on the expression level of wall secondary CesAs, in the stem of 2-month-old F 3 1 seedlings, variations of CesAs expression were revealed showing a downregulation of primary CesA1 and CesA3 and secondary CesA8, and an increased expression of primary CesA6 and secondary CesA4. Furthermore, confocal microscope observations of cross-sections of F 3 1 first internode showed a remarkable increase in the diameter of the secondary lignified xylem with respect to wt, suggesting an organ specificity of the transformation effects. In Arabidopsis, all primary CesAs are reported to physically interact with secondary CesAs, both in vitro and in planta. The possibility that mixed complexes of primary and secondary CesAs may occur at particular times has been suggested, mainly during the transition between primary and secondary cell wall [65]. Our differential expression of primary and secondary CesAs would seem to corroborate this hypothesis.
In Arabidopsis, the overexpression of CesA6-like genes has been suggested to influence both cell expansion and division [22]. Furthermore, it has been reported that several cellulose synthase-like (Csl) genes promote cell division and affect cell cycle, including CslD5 during Arabidopsis stomatal lineage development [66] and in meristems [67], CslD1 in maize [68], CslDs (CslD1, CslD2a, CslD2b, CslD3a, CslD5) in Dendrobium catenatum [69], and CslD4 in rice [70]. The increase of all parameters related to the dimension of F 3 1 stem, leaves, and roots is in line with the above-mentioned results, suggesting that both cell division and cell expansion are increased in the F 3 1 line. Indeed, a stimulus of the vascular cambium meristematic activity is strictly related to the sharp increase of secondary lignified xylem thickness observed in F 3 1 stem. Moreover, the average area of F 3 1 abaxial epidermal pavement cells from the distal lamina sector of the most expanded leaves was increased compared to wt congruently with an effect of GFP-CesA6 expression on cell expansion. Moreover, in F 3 1 leaves, we also observed the overexpression of PIP1 and PIP2; these genes code for two aquaporin isoforms assembling in heterotrimers at the plasma membrane and able to mediate significant transmembrane water fluxes [71]. PIP1 and PIP2 were reported to be involved in the expansion of rose petals [72] and pomegranate flowers [73], so that the overexpression of these two PIPs in F 3 1 line can be considered a further indication on the effect of GFP-CesA6 expression on the expansion growth.

Conclusions
Not all attempts to increase the biosynthesis of cellulose through the overexpression of the CesA genes have given positive results. In transgenic barley lines, overexpressing primary and secondary CesAs, no increase in the cellulose content was obtained and, in some lines, it decreased significantly. A decrease of cell wall thickness, dwarf plants, and a "brittle node" phenotype were also found and justified as the need for a coordinated overexpression of all three HvCesA genes for the formation of a functional complex [20]. Taken together, the results presented here highlight that the AtCesA6 GFP-tagged expression led to the formation of tobacco plants with a "giant" phenotype with an evident increase of plant biomass. AtCesA6 affected cell wall biosynthetic machinery involved in lignocellulosic biomass production: (i) enhancing the expression of other tobacco primary wall CesAs (NtCesA1, NtCesA3, and NtCesA6) to produce more cellulose and those of other genes involved in the cell wall construction and modification (NtXyl4, NtGal10, NtExp11, NtPIP1.1, and NtPIP2.7); (ii) causing the relative amounts of cell wall polysaccharides and possibly of non-polysaccharide material in transformed tobacco seedlings; and (iii) increasing secondary cell wall deposition clearly visible in tobacco secondary growth of the stem. We prospect the possibility of using tobacco plants expressing the exogenous gene GFP-CesA6 to increase plant biomass for sustainable development and the diversification of the energy sources.

Data Availability Statement:
The data presented in this study are available upon request from the corresponding authors.