Skip to main content

Exploring tomato Solanum pennellii introgression lines for residual biomass and enzymatic digestibility traits



Residual biomass production for fuel conversion represents a unique opportunity to avoid concerns about compromising food supply by using dedicated feedstock crops. Developing tomato varieties suitable for both food consumption and fuel conversion requires the establishment of new selection methods.


A tomato Solanum pennellii introgression population was assessed for fruit yield, biomass phenotypic diversity, and for saccharification potential. Introgression lines 2–5, 2–6, 6–3, 7–2, 10–2 and 12–4 showed the best combination of fruit and residual biomass production. Lignin, cellulose, hemicellulose content and saccharification rate showed a wide variation in the tested lines. Within hemicellulose, xylose value was high in IL 6–3, IL 7–2 and IL 6–2, whereas arabinose showed a low content in IL 10–2, IL 6–3 and IL 2–6. The latter line showed also the highest ethanol potential production. Alkali pre-treatment resulted in the highest values of saccharification in most of lines tested, suggesting that chemical pretreatment is an important factor for improving biomass processability. Interestingly, extreme genotypes for more than one single trait were found, allowing the identification of better genotypes. Cell wall related genes mapping in genomic regions involved into tomato biomass production and digestibility variation highlighted potential candidate genes. Molecular expression profile of few of them provided useful information about challenged pathways.


The screening of S. pennellii introgression population resulted very useful for delving into complex traits such as biomass production and digestibility. The extreme genotypes identified could be fruitfully employed for both genetic studies and breeding.


Over the last decades, rising concerns upon depleting fossil fuels has resulted in an increased interest in fuels derived from bio-renewable sources including sugars, starch and lignocellulosic materials [1]. Lignocellulosic biomass materials constitute the most abundant renewable substrate for ethanol [2]. Currently, cellulosic feedstocks derived from dedicated biomass crops in the U.S., South America, Asia and Europe are from corn stover, sugarcane bagasse, or perennial low input-high yield crops such as miscanthus or switchgrass [3]. More interestingly, lignocellulosic biomass can be obtained in large-scale from agricultural residues, making their conversion into fuel more advantageous from the economic, environmental and strategic points of view [4]. In particular, the production of ethanol from biomass residuals represents a unique opportunity to avoid concerns about compromising food supply by using starch or sucrose based feedstocks [57]. So far, a few attempts have been conducted to investigate the potential of producing fuel from residual biomass of tomato (Solanum lycopersicum L.), a major vegetable crop worldwide.

Biomass production depends on several traits related to morphological and physiological processes controlling the plant vegetative growth. Developing new varieties for both food and fuel production, will require the establishment of new selection methods. Moreover, in order to elucidate genetic interactions between traits, it will be important to understand the correlations between traits and the extent to which they can be uncoupled, since positive and negative correlations can have profound effects on each other or comprising other aspects of crop production [8]. Genomic resources like wild introgression populations can facilitate the identification of tomato genotypes characterized by both high fruit and residual biomass production.

The tomato Solanum pennellii introgression population is a permanent mapping source for Quantitative Trait Loci (QTLs) analysis composed of a series of introgression lines, in which defined genomic segments of the S. pennellii genome replaced homologous regions in S. lycopersicum (cultivar M82) background. Such population can be very effective for identifying QTL, because any phenotypic difference between an introgression line and the recurrent parental line is attributed solely to donor parent genes within the introgressed chromosomal segment [9]. The assessment of S. pennellii introgression lines phenotypic and chemical traits can provide useful genetic information about tomato biomass production and potential fuel conversion.

Numerous structural and compositional features can have effects on lignocellulosic biomass processability. Cellulose is a polymer of 1–4 β linked glucose and forms crystalline fibrils within the cell walls. Hemicelluloses are complex-polymers of hexoses (mannose, glucose, galactose) and pentoses (xylose and arabinose), arranged in a non-crystalline manner, and interact with cellulose fibers [10]. The cellulose cristallinity and the heterogeneity of the hemicellulosic fraction confer recalcitrance to the biomass and represent a barrier for the utilization of the sugars locked in the polymers [11]. Pretreatments can be useful for modifying the architecture of the cell walls that compose the biomass, making it more accessible to hydrolytic enzymes. This involves the modification of lignin, removal of matrix polysaccharides, and reduction of cellulose crystallinity [12].

A way to improve feedstock amenability to fuel conversion is to delve into genetic variability within one species for breeding towards the best traits. In this work, we approach the study of a complex trait such biomass production and fuel conversion making use of S. pennellii introgression lines. Firstly, we assessed the tomato introgression lines biomass phenotypic diversity for identifying both traits contributing to high biomass yield and genotypes to employ in selection program for increasing biomass production. Then we characterised the residual biomass for lignin components, correlating them with the saccharification potential. In addition, we performed, on selected genotypes, a comparison among different cell wall polysaccharide components in order to identify chemical traits associated with tomato residual biomass digestibility. Finally, as proof of concept we assessed genes involved in cell wall biosynthesis in extreme genotypes by in silico analysis and by molecular assay, in order to investigate their role in phenotypic diversity.


Characterization of morphological traits in S. pennellii introgression lines

Thirty-seven tomato introgression lines, able to cover the full S. pennellii genome, have been assessed for fruit yield and biomass production. The frequency distribution of fruit yield and biomass parameters (Fig. 1) showed that the fruit yield of 50 % of the introgression lines felt in the range 369–754 g per plant, with a harvest index range of 45–61 and a residual biomass of 295–701 g. Overall, IL 12–4, 6–3 and 2–5 provided the highest fruit yields (1043, 978 and 958 g per plant respectively), whereas 4–4 and 6–2 gave the lowest (117 and 25 g per plant respectively), whilst the recurrent parent M82 showed a fruit yield of 665 g per plant and a harvest index of 61.5 (Additional file 1: Tables S1). Interestingly, IL 6-2 and IL 6–3 have both an indeterminate growth habit, but showed extremely different yield performance, in agreement with data reported by Eshed and Zamir [9].

Fig. 1

Frequency distribution of yield parameters, plant residual biomass and leaf area in 37 tomato introgression lines

For crop residual biomass and leaf expansion, IL 2–6 showed the highest values of 1136 g and 6.6 m2 per plant, respectively. It was followed by IL 4–4 (950 g) in the biomass ranking and by IL 12–4 (4.55 m2) with regard to leaf area. The lowest values of residual biomass and leaf area were showed by IL 3–2 and IL 1–1 (128 g and 0.47 m2 per plant respectively). M82 showed a biomass yield of 420 g per plant and a leaf area of 1.62 m2 per plant (Additional file 1: Tables S2 and Additional file 2: Table S3). Extreme genotypes were identified for each measured trait. The most skewed traits were fruit weight per plant, residual biomass per plant and leaf area per plant. The introgression lines 2–5, 2–6, 6–3, 7–2, 10–2 and 12–4 showed the highest combination of fruit and residual biomass production and the introgression lines 2–1, 3–2, 4–3, 5–1 and 7–1 the lowest.

Lignin content and saccharification rate

The highest lignin content was in IL 2-6 (21.5 g 100 g−1 of dry weight), the lowest in IL 1–1 and IL 8–3 (10.5 g 100 g−1), M82 display a lignin content of 17.6 g 100 g−1 of dry weight (Table 1; Additional file 2: Table S3). Saccharification rates higher than 10 mg glucose g−1 biomass h−1 (Table 1) were recorded in nine genotypes (IL 2–6, IL 3–3, IL 3–4, IL 4–3, IL 8–1, IL 8–2, IL 9–1, IL 10–2 IL 11-1 and M82), whereas a very low value was displayed by IL 8–3 (1.3 mg glucose g−1 biomass h−1). Fourty-seven per cent of genotypes evidenced a saccharification value ranging between 8 and 12 mg glucose g−1 biomass h−1. Only the 8 % of the introgression lines showed saccharification rates lower than 2 mg glucose g−1 biomass h−1 (IL 1–1, IL 3–2 and IL 8–3) and the 3 % higher than 13 mg glucose g−1 biomass h−1 (IL 3–4). Out of total lines analyzed, eight line, producing high amount of dry biomass (IL 1–3, IL 1–4, IL 2–2, IL 2–6, IL 6–1, IL 6–3, IL 7–2, IL 12–4), showed a correlation (R2 = 0.86) between saccharification rate and lignin content (Fig. 2).

Table 1 Lignin content and saccharification rate in 37 tomato introgression lines
Fig. 2

Linear regression between lignin content and saccharification rate in a subset of ILs with high dry biomass production

Cell wall composition and theoretical ethanol yield

Table 2 shows the analysis of cell wall polysaccharide components performed on selected ILs. Total cellulose ranged from 40.7 g 100 g−1 of dry weight (IL 4–1) to 55.9 g 100 g−1 of dry weight (IL 4–3). Within the total cellulose, the crystalline fraction accounted for 3.5 (IL 3–1) to 17.6 (IL 4–3) g 100 g−1 of dry weight. The highest hemicellulose content recorded in IL 4–3 (6.0 g 100 g−1 d.w.) was about three times greater than that showed by IL 2–5. Finally, the pectin content displayed a wide range comprised between 2.0 and 10.1 g 100 g−1 d.w. (in IL 9–3 and IL 10–3 respectively). M82 showed a low level of total cellulose (38.7 g 100 g−1 of dry weight) but the crystalline fraction and hemicellulose accounted for 15.2 and 5.2 g 100 g−1 of dry weight respectively. Among the introgression lines tested, IL 2–6 and IL 6–3 showed high content of cellulose (51.7 and 55.8 g 100 g−1 d.w. respectively) and hemicellulose (5.6 and 4.4 g respectively) but moderate crystalline cellulose (11.5 and 12.5 g). Conversely, IL 3–1 and IL 4–1 showed low content of cellulose (41.0 and 40.7 g 100 g−1 d.w. respectively), crystalline cellulose (3.5 and 5.4 g respectively) and hemicellulose (2.9 and 2.5 g respectively). Interestingly, lignin was positively correlated with crystalline cellulose (r = 0.52 at p < 0.05) and with hemicellulose (r = 0.51 at p < 0.05). Crystalline cellulose, in turn, was also correlated with hemicellulose (r = 0.87 at p < 0.01).

Table 2 Chemical composition of 13 tomato introgression lines residual biomass

The monosaccharide composition in the hemicellulosic fraction showed a large degree of variability between genotypes tested (Table 3). Glucose was the highest represented sugar, ranging from 26.9 g 100 g−1 d.w (IL 4–1) to 18.4 g 100 g−1 d.w (IL 4–3). Xylose followed glucose in the rank, showing a variation comprised between 13.4 g 100 g−1 d.w. in IL 3–1 and 26.5 g 100 g−1 d.w in IL 6–3. Fucose was the least represented monosaccharide, with values ranging from 0.27 (IL 2–6) and 0.92 g 100 g−1 d.w. (IL 10–3). Based on cell wall content of cellulose and hemicellulose monosaccharides (on dry weight) a theoretical ethanol yield was calculated (Fig. 3). IL 2–6 showed the highest ethanol potential production (2681 L ha−1) as it provided with the highest biomass yield (6 t ha−1) and a good conversion of biomass into ethanol (451 L t−1 d.w.). IL 4–3 resulted in a very low theoretical ethanol yield (539 L ha−1), in spite of the best biomass quality for ethanol production (486 L t−1 d.w.). Finally, IL 4–1 showed the lowest potential conversion of biomass into ethanol (340 L t−1 d.w).

Table 3 Monosaccharide composition of hemicellulose in 13 tomato introgression lines
Fig. 3

Theoretical ethanol yield in 13 tomato introgression lines plus M82

Comparison among different saccharification pretreatments

Three biomass pretreatments were assessed in order to understand if saccharification rate is affected by treatment in genotype specific manner. Alkali pre-treatment caused the highest values of saccharification in most lines, whereas water pre-treatment was generally the least effective (Fig. 4). IL 4–3, IL 6–3, and IL 2–6 resulted in the highest values of saccharification rate (41, 38 and 35 mg glucose g−1 biomass h−1 respectively) using the alkali pretreatment. Interestingly, the saccharification rate was positively correlated with cellulose content, upon alkali pre-treatment, and with the crystalline cellulose content and hemicellulose content upon all the three pretreatments tested (Additional file 3). Finally, IL biomass saccharification rate was not significantly correlated with the hemicellulose arabinose/xylose ratio, regardless of the pre-treatment type. Our results display that biomass enzymatic digestibility is affected both by genotype cell wall structures and pretreatment type employed.

Fig. 4

Saccharification rate in 13 tomato introgression lines plus M82, as a function of the pretreatment type

In silico identification of candidate QTL involved in biomass production

Based on information deposited in cell-wall literature studies and in database resources [13], we selected 23 protein families implicated in the construction and modification processes of plant cell walls. Using this dataset, the tomato genome was explored to identify the enzymes involved in the biosynthesis of key polysaccharides such as pectin, lignin, cellulose, hemicellulose and starch. In order to identify candidate QTL for biomass production by in silico analysis, we looked at selected ILs for cell wall related genes (see Additional file 2). Several genes located in the introgression regions of contrasting genotypes 2–5, 2–6, 3–1, 4–1 4–3, 6–2, 6–3, 7–2, 9–3, 10–2, 10–3, 11–3, and 12–4 were identified. Interestingly, introgression region 2–6 included a high number of Peroxidases, Laccases, Glycosyltransferases and MYB transcription factors in a region of 3.7 Mb. Moreover, two Cellulose synthase, two UDP-D-glucose dehydrogenase, a specific GDP-mannose dehydratase, a Sucrose synthase, a Galactosyltransferase, a GDP-mannose dehydratase and a GHMP kinase have been identified in this chromosome area (Table 4). Conversely, IL 4–3 in a much larger region includes only 30 genes coding for enzymes involved in cell wall biosynthesis.

Table 4 Proteins involved in the construction and modification processes of cell wall polysaccharides identified in tomato 13 introgression line chromosome regions

Molecular expression profile

RT- qPCR analysis was used to analyze the expression pattern of key genes in two contrasting ILs (12–4 and 4–3). Aldehyde dehydrogenase (ALDH), Xyloglucan endotransglucosylase (BRU1), Cellulose synthase-like glycosyltransferase (CslC1-1,CslC2, CslC6), Alpha-L-fucosidase1 (FUCA1) Glucose transporter 8 (GLUT8), GDP-mannose transporter (GMT) and Cell-wall invertase (INV2) gene expression were assessed in the two selected lines and in the M82 genotype used as control. The invertase INV2 gene showed an opposite expression pattern in the two ILs tested, in fact it was up-regulated in IL 12-4 and down-regulated in IL 4-2 (Fig. 5a). GMT and GLUT8, respectively mannose and glucose transporters showed an increased transcription in IL 4–3 compared to M82 control, while were both down regulated in IL 12–4 (Fig. 5b). The expression of the cellulose synthesis like genes was up-regulated in IL 12–4 and down-regulated in IL 4–3, except for the Cslc1.1 gene (Fig. 5c). Transcription levels of the BRU gene, involved in xyloglucan synthesis, was high in both IL, while transcription of FUCA, involved in hydrolysis L-fucose, was lower than in M82 in both ILs (Fig. 5d). Finally, the expression of ALDH was slightly up-regulated in IL 12–4 and down-regulated in IL 4-3 (Fig. 5e).

Fig. 5

Expression fold change of different genes among IL 12-4 (dark gray) and IL 4-3 (light gray) respect to M82. Panel (a) shows the expression of cell wall invertase (INV2); (b), GDP-mannose transporter (GMT) and Glucose transporter 8 (GLUT8); (c), cellulose synthase-like glycosyltransferase (CslC1-1,CslC2, CslC6); (d),  xyloglucan endotransglucosylase (BRU1) and α-L-fucosidase1 (FUCA1); and (e), Aldehyde dehydrogenase (ALDH). Error bars relative to RT–qPCR experiments represent the estimate of standard error of the mean (SEM) for the replicates


Large variations in fruit yield and biomass production among lines was evidenced in a two years field experiment conducted with S. pennellii IL population. Our results support previous findings [14] reporting that most of tomato introgression lines showed lower yield and fruit mean weight but higher vegetative biomass than the parent M82. Interestingly, in our study a number of genotypes with extreme phenotypes were identified. Analysis of extreme genotypes that exceed +/− 2.5 mean value of quantitative trait loci (QTL), provide nearly equivalent power to complete genotyping at a reduced cost [15, 16].

The introgression lines tested exhibited dramatic variations also in cell wall lignin content, with the highest value being more than two times greater than the lowest. A similar survey on a set of Arabidopsis thaliana mutants for lignin biosynthetic genes revealed that lignin content ranged from 6.0 g 100 g−1 dry weight to 15.9 [17]. Lignin composition could affect the saccharification efficiency and usually plants with increased biomass and reduced lignin have an improved fuel production [17]. However in our work a subset of ILs with high dry biomass and lignin content, also showing a high saccharification rate was identified. Different Mischanthus genotypes also display different patterns of correlation between lignin content and saccharification efficiency [18]. Transgenic plants with modified lignin composition enhanced biomass saccharification [19]. In tobacco, up-regulation of sucrose metabolism genes appears to directly impact primary growth and therefore biomass production, also with slight decrease of lignin content [20]. Determining sample composition, especially structural carbohydrate, and using this information to predict relative performance for fuel conversion could be very useful [21]. In our study a number of lines with high potential ethanol production were identified, and as outlined, ethanol production can depend on different cell wall structure variations. Indeed, the cell wall architecture can be assembled from many different types of polysaccharides, phenylpropanoids and structural proteins. In our research, crystalline cellulose was highly correlated with hemicelluloses; with this respect, the amount and composition of branches attached to the hemicellulose backbone can affect the cell wall plasticity and crystal structure [22]. Within the hemicelluloses, the amount of xylan was high in IL 2–6, IL 4–3, IL 6–2, IL 6–3, IL 7–2, IL 10–3 and IL 11–3 and low in IL 3-1 and IL 12–4. Cell wall recalcitrance varies among plant species and even within different genotypes of the same species [23]. The close relationship between recalcitrance and the chemical composition of the non cellulosic matrix suggests that cell wall strength could be tuned by carefully controlling the matrix composition [17, 24].

Knowledge on cell wall composition can be useful for better direct genetic approaches and pretreatment design to render biomass more amenable to processing. In this respect, the effect of lignin as well as of cellulose on biomass digestibility has been described in previous research [25]. As for the effect of pretreatment on biomass saccharification rate, in our research alkali pretreatment showed the best performances. Consistently with our results, Lima et al. [26] reported that one-step alkali pretreatment improves the enzymatic digestibility of Eucalyptus bark compared to two-step pretreatment with HCl 1 % followed by 4 % NaOH. Acid and alkali pretreatments have distinct mechanisms for biomass modification [27]. Acid pretreatment involves the hydrolysis of hemicelluloses by breaking the glycosidic linkages of polysaccharides [28]. Alkali pretreatment, in turn, breaks down the intermolecular ester bonds that cross-link lignin with hemicelluloses, thereby solubilizing lignin and hemicelluloses [29]. The removal of hemicelluloses increases the mean pore size of the substrate, which facilitates the hydrolysis of cellulose [30]. Correlation among parameters can be affected by composition. Significant positive correlation between hexose sugar production and hemicellulose arabinose content or arabinose/xylose ratio resulted in enhanced biomass digestibility both with acid and alkaline pre-treatments [31]. The composition study conducted here allowed us to detect genotypes with wide differences in biomass production, as well as in cell wall composition. Such large differences in the composition of biomass have been previously observed by van Acker et al. [17] in both, cellulose and hemicellulose glucose content among Arabidopsis thaliana mutants. A similarly large variation was also observed by Marriott et al. (2014) in Brachypodium mutants. The variation in biomass composition between lines observed in the present work, although surprisingly wide, can be explained by the fact that these ILs originated from a genetically divergent biparental population. The identification of genes that affect such traits is still challenging. To this purpose a core panel of genes involved in metabolic pathways of cell wall biosynthesis and degradation was mapped in selected ILs with extreme phenotypes. These candidate genes encode proteins predicted to play a role in the synthesis, modification, assembly and disassembly of lignin, hemicellulose, cellulose and pectin [32]. A previous comparative approach showed that the differences in wall architecture between Arabidopsis and rice actually mirror the diversity of the individual gene families involved in the cell wall dynamics of the respective plant species [33]. Changes in wall composition or architecture could be due to mutations either in genes directly related to cell wall metabolism or in genes involved in regulation [34]. Morgan et al. [35] performed a biochemical analysis of a tomato introgression line with increased levels of fruit citrate identifying a target gene that was successfully tested in transgenic plants. Mutations with major effect more frequently occur in domesticated or artificially disturbed populations [36]. This supports the use of inbreds obtained between wild and domesticated species for identifying traits with strong and repeatable phenotypes [37]. In our work several specific candidate genes co-localizations could be underlined. For example, the genomic region delimitated by IL 2–6 contains 20 Peroxidase, 19 MYB factors, 6 Glycosyltransferase, 7 Laccase, 2 UDP-D-glucose dehydrogenase and 2 Cellulose synthase as well as other important polysaccharides biosynthesis enzymes. Different isoforms of UDP-glucose-consuming pathway have a regulatory role in carbon partitioning between cell wall formation and sucrose synthesis [38]. A Glycosyltransferase mutant produced plants deficient in ferulic and coumaric acid, aromatic compounds known to be attached to arabinosyl residues in xylan substituted with xylosyl residues. The mutant plants exhibit an increased extractability of xylan and increased saccharification, probably reflecting a lower degree of diferulic cross-links [39]. Genes and regulatory elements present in S. pennelli and S. lycopersicum could be involved in alterations of cell wall composition and biomass production.

The molecular expression profile suggests that the partioning between source and sink organ could be challenged in different lines. Indeed, IL 12–4 (a high biomass producer line) showed a major cleavage of sucrose compared to IL 4-3, due to the up-regulation of a cell wall invertase (INV2) and a lower activity of a GDP-mannose transporter and a of glucose transporter 8. Cleavage of sucrose by invertase is generally correlated with growth and cell expansion, associated with sucrose partioning [40]. Excess of sucrose is broken down into fructose and UDP-glucose, which is employed in the synthesis of cell wall polymers [20]. Cellulose synthase-like (CslC2 and CslC6) proteins, involved in the synthesis of various β-glycan polymers [41] by using GDP-mannose as substrate, result both up-regulated in IL 12-4 and down-regulated in IL 4-3 (low biomass producer line) whilst CslC1 is up-regulated in both. A decrease of mannose substrate to synthesize 1,4-β-glucan backbone of mannose could improve IL12-4 digestibility. The expression of ALDH, an enzyme involved in the synthesis of lignin components such us ferulate and sinapate [42] is also up-regulated in IL 12-4. Both ILs have a similar amount of lignin but different saccharification rate; this difference may be related either to their lignin composition or to the fact that saccharification is a multigenic trait which is affected by a large number of factors in the cell wall [43]. Genes involved in carbon metabolism mapped in IL extreme genotypes are part of a “network QTL”, where several elements of a metabolic network are affected by expression QTLs, enzyme activity QTLs, or metabolite QTLs [44, 45] fine tuned in any genotype.


The present work shows the potential for exploitation of tomato ILs residual biomass for fuel conversion. The tomato introgression lines showed high variability in biomass production, cell wall composition and saccharification rate and, consequently, potential ethanol yield also resulted in a wide range of values among the genotypes. The trait enhancement found in extreme genotypes could be compared with the behavior of recurrent parent through a combined genomic and chemical profiling. It is evident, even from literature, that an intricate network of relations among components govern biomass production and digestibility. For this reason it is difficult to identify genotypes that combine both characteristics, on the other hand expression pattern of genes related to biomass production provided interesting clues which should be further investigated. Precise gene mapping is needed in order to predict biomass quality based on genomic information, while the interrogation of contrasting lines could permit the identification of the most eligible alleles linked to saccharification efficiency. Indeed, selected lines will be further explored to identify the candidate genes involved.


Phenotypic analysis

Research was carried out for two years in Naples, southern Italy (40°50' N, 14°15' E, 17 m a.s.l.), on the commercial tomato variety M82 (Solanum lycopersicum L.), Solanum pennellii LA716 and on 37 Solanum pennellii tomato introgression lines, kindly provided by Dr Dani Zamir (Hebrew University of Jerusalem) and reported in Table 1. Information on such lines can be found at website A randomized complete block design with three replicates was arranged and each plot had a 4.73 m2 (1.75 x 2.70 m) surface area. Transplanted plants were arranged in single rows spaced by 0.90 m from each other and the spacing was of 30 cm along the rows (3.2 pt · m−2). Experimental research on plants were conducted in accordance with local legislation. Plants were grown under standard tomato field procedures used for the area and fruits were harvested at full ripening. S. pennellii plants failed to grow properly in our climatic conditions and the few samples obtained were not included in the following analysis.

General analytical methods

Plant samples were randomly selected to assess the maximum leaf surface extension using a bench top LI-COR leaf area meter. At harvest, the following determinations were made: a) weight and number of ripe undamaged fruits, classified as marketable; fruit mean weight on 50 unit samples; b) residual biomass, including leaves, shoots, stems; and c) immature or damaged fruits. Harvest index was calculated as a ratio between marketable fruits and total plant weight. After harvest, residual biomass showed no fungal symptoms; therefore samples were randomly collected in each plot and immediately transferred to the laboratory, where they were dried in an oven at 70 °C under vacuum until they reached constant weight. After assessing the dry residue, samples were carefully milled, avoiding mixing of materials belonging to different plant organs. The final material, composed of particles ≤ 1 mm diameter, was stored in air-tight bags at −20 °C and further dried just before being processed.

Chemical analyses

After harvesting and weighing out, the residual plant biomass collected was oven-dried at 60 °C. Lignin determination and saccharification assay using water pre-treatment were performed in all the 38 genotypes tested for yield and biomass production. Other analyses (lignin, cellulose, hemicellulose, pectin, hemicellulose monosaccharides, crystalline cellulose and saccharification assay with acid or alkaline pretreatment) were performed on 13 selected genotypes (5 genotypes with biomass production < 300 g pt−1 of fresh weight; 3 genotypes with biomass production from 301 to 800 g pt−1 of fresh weight; 5 genotypes with biomass production > 800 g pt−1 of fresh weight).

Lignin determination: acetyl bromide method

Biomass powder was weighed out (4 mg) into 2 mL tubes. The biomass was heated at 50 °C for 3 h after adding 250 μL of acetyl bromide solution (250 μL of acetyl bromide and 750 μL of glacial acetic acid in volume) and vortexing every 15 min. After the samples were cooled to room temperature, the content was transferred into 5 mL volumetric flasks. A further 1 mL of NaOH (2 mol L−1) was used to rinse the tubes pouring the NaOH into the 5 mL flasks. 175 μL of hydroxylamine HCl (0.5 mol L−1) was added to the volumetric flasks and, after vortexing, the latter were filled up to 5 mL with glacial acetic acid and mixed several times. Finally, in order to measure the 280 nm UV adsorption by spectrophotometer, 100 μL of each sample was diluted in 900 μL of glacial acetic acid. The amount of lignin was calculated using the following formula: \( \begin{array}{l}\left[\mathrm{absorbance}/\left(\mathrm{coefficient}\ \mathrm{pathlength}\right)\right]\ \left[\left(\mathrm{total}\ \mathrm{volume}\ 100\%\right)/\mathrm{biomass}\ \mathrm{weight}\right],\ \mathrm{where}\ \mathrm{coefficient} = \\ {}15.69,\ \mathrm{pathlength} = 1,\kern0.2em \mathrm{total}\ \mathrm{volume} = 5,\ \mathrm{biomass}\ \mathrm{weight} = 4.\end{array} \)

Cellulose, hemicellulose and pectin determination


A mixture of 240 mL of water, 0.75 mL of glacial acetic acid and 2.25 g of sodium chlorite were added to 7.5 g of extracted and dried sample and kept at 75 °C for 3 h.

At hourly intervals, a volume equivalent to the initial amounts of glacial acetic acid and sodium chlorite was added to the biomass. The sample obtained was filtered and washed up first with cold water, then with warm water and finally with acetone. The residue was oven-dried at 105 °C for 24 h and then weighed to calculate the content of holocellulose.


1.3 g amount of the resulting holocellulose was treated with 26 mL of potassium acetate (0.6 mol L−1) and incubated at 75 °C for 3 h before adding 26 mL of ammonium oxalate (0.04 mol L−1). The suspension was kept at 75 °C for 3 h. Then, the samples were filtered and washed up with excess of water before the residue was oven-dried at 105 °C for 24 h. The pectin content was calculated as the difference between the holocellulose fraction and the above residue.

Cellulose and hemicellulose

A sample of holocellulose (3.8 g) was treated with 100 mL of sodium hydroxide (4.4 mol L−1) at room temperature for 30 min and filtered. Then, it was washed up sequentially with warm water (200 mL), 5 mL of acetic acid (2 mol L−1) and 500 mL of water. Next, the residue was oven-dried at 105 °C for 24 h and weighed, providing the cellulose fraction. The hemicellulose content was calculated by subtracting the cellulose and pectin amount from that of holocellulose. The filtration process as well as the subsequent residue drying reported for pectin, cellulose and hemicellulose were accurately performed and, in fact, neither material loss nor different water content in the dried cell walls compared to the initial samples were assessed.

Non cellulosic monosaccharide determination

Biomass dry powder (4 mg) was partially hydrolyzed by adding 0.5 mL of trifluoroacetic acid (2 mol L−1). Then, the vials were flushed with dry argon, mixed and heated at 100 °C for 4 h, mixing periodically. The vials were then cooled to room temperature and dried in centrifugal evaporator with fume extraction overnight. The pellets were washed twice with 500 μL of 2-propanol and vacuum dried. Finally, the samples were resuspended in 200 μL of deionised water, filtered with 0.45 μm PTFE filters, and analyzed by HPAEC. Monosaccharides were separated and quantified by HPAEC using a Dionex ICS-3000 with integrated amperometry detection. Chromatographic separation was performed on a CarboPac PA20 (3 x 150 mm) column (Thermo) using a gradient elution. The mobile phase consisted of solution A: 100 % water, solution B: 200 mM NaOH, and solution C: 0.1 M Sodium Hydroxide, 0.5 M Sodium Acetate. A flow rate of 0.5 ml min-1 was used and the gradient was as follows: 0 min: 100 % A; 5 min: 99 % A, 1 % B; 15 min: 99 % A, 1 % B; 22 min: 47.5 % A, 22.5 % B, 30 % C; 30 min: 47.5 % A, 22.5 % B, 30 % C. The column was then washed as follows: 30.1 min: 100 % B; 37 min: 100 % B; 37.1 min: 99 % A, 1 % B; 50 min: 100 % A; 55 min: 100 % A. The separated monosaccharides were quantified by using external calibration with a mixture of nine monosaccharide standards at 100 μM (arabinose, fucose, galactose, galacturonic acid, glucose, glucuronic acid, mannose, rhamnose, and xylose) that were subjected to acid hydrolysis in parallel with the samples.

Crystalline cellulose

Biomass dry pellets after TFA hydrolysis were washed once with 1.5 ml of water, and twice using 1.5 ml of acetone. The dried pellets were left to air dry overnight before complete hydrolysis by adding 90 μl of 72 % (p/v) sulfuric acid, incubating at room temperature for 4 h. 1.89 ml of water was subsequently added and the sample was heated for 4 h at 120 °C. The glucose content of the supernatant was assessed using the colorimetric Anthrone assay, using a glucose standard curve.

Theoretical ethanol yield calculation

The theoretical ethanol yield was calculated considering the total cellulose conversion in the sample, according to the National Renewable Energy Laboratory standards [46, 47]. The theoretical ethanol yield was expressed also taking into account agronomical traits such as the biomass yield per surface area unit.

Saccharification assay

Formatting of plant materials

Loading of plant powder into 96-well plates, using a custom-made robotic platform (Labman Automation, Stokesley, North Yorkshire, UK), and saccharification assays were performed according to Gomez et al. (2010) [48] after water, acid or alkali pretreatment. Enzymatic hydrolysis was carried out using an enzyme cocktail with a 4:1 ratio of Celluclast and Novozyme 188.

In silico search for cell wall related genes

Marker sequence coordinates that delimitated introgression lines (2–5, 2–6, 3–1, 4–1, 4–3, 6–2, 6–3, 7–2, 9–3, 10–2, 10–3, 11–3, 12–4) were downloaded by SGN website [13]. ITAG annotated proteins retrieved in the genomic areas bounded by the markers chosen were recorded in an excel file. Proteins involved in the construction and modification processes of cell wall polysaccharides (Additional file 2) have been searched for each IL line genomic area and stored in a separated file.

Gene expression analysis

Total RNA of IL 12–4, IL 4–3 and M82 genotypes was extracted from leaf tissues, using a Kit Spectrum plant total RNA (Sigma), and treated with DNase I Digestion (Sigma). Total RNA was quality checked and cDNA synthesis was performed with oligo (dT) and SuperScript III Reverse Transcriptase (Invitrogen). Specific primers for candidate genes (Additional file 2: Table S5) were designed using Primer3 software. RT–qPCR was performed in a 12.5 μL reaction volume using the SensiFAST SYBR Hi-ROX Kit (Bioline) with 4.5 μL cDNA as a template. Each reaction was carried out in triplicate and run on the 7900HT Fast Real-Time PCR System (Applied Biosystems). Fold change of each transcript, normalized to EF (elongation factor), was calculated relative to expression in the M82 sample, using the 2–ΔΔCt method.

Statistical analysis

Data were processed by analysis of variance and mean separations were performed through the Duncan multiple range test, with reference to 0.05 and 0.01 probability levels, using SPSS software version 17. Data expressed as percentage were subjected to angular transformation before processing. Correlations were performed with all pairs of chemical parameters using the software mentioned above.

Availability of data and materials

Supporting data are included as additional files.


  1. 1.

    Gnansounou E. Production and use of lignocellulosic bioethanol in Europe: current situation and perspectives. Bioresour Technol. 2010;101:4842–50.

    CAS  Article  PubMed  Google Scholar 

  2. 2.

    Somerville C, Youngs H, Taylor C, Davis SC, Long SP. Feedstocks for lignocellulosic biofuels. Science. 2010;329:790–2.

    CAS  Article  PubMed  Google Scholar 

  3. 3.

    Cheng KK, Cai BY, Zhang JA, Ling HZ, Zhou YH, Ge JP, Xu JM. Sugarcane bagasse hemicelluloses hydrolysate for ethanol production by acid recovery process. Biochem Eng J. 2008;38:105–9.

  4. 4.

    Cardona CA, Sanchez OJ. Fuel ethanol production: process design trends and integration opportunities. Bioresour Technol. 2007;98:2415–57.

    CAS  Article  PubMed  Google Scholar 

  5. 5.

    Schievano A, D’Imporzano G, Adani F. Substituting energy crops with organic wastes and agro-industrial residues for biogas production. J Environ Manag. 2009;90:2537–41.

    CAS  Article  Google Scholar 

  6. 6.

    Sun Y, Cheng J. Hydrolysis of lignocellulosic materials for ethanol production: a review. Bioresour Technol. 2002;83:1–11.

    CAS  Article  PubMed  Google Scholar 

  7. 7.

    Wheals AE, Basso LC, Alves DMG, Amorim HV. Fuel ethanol after 25 years. Trends Biotechnol. 1999;17:482–7.

    CAS  Article  PubMed  Google Scholar 

  8. 8.

    Farrar K, Bryant DN, Turner L, Gallagher JA, Thomas A, Farrell M, et al. Breeding for bio-ethanol production in Lolium perenne L.: association of allelic variation with high water-soluble carbohydrate content. Bioenergy Res. 2012;5:149–57.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  9. 9.

    Eshed Y, Zamir D. An introgression line population of Lycopersicon pennellii in the cultivated tomato enables the identification and fine mapping of yield-associated QTL. Genetics. 1995;141:1147.

    CAS  PubMed  PubMed Central  Google Scholar 

  10. 10.

    Gírio FM, Fonseca C, Carvalheiro F, Duarte LC, Marques S, Bogel-Łukasik R. Hemicelluloses for fuel ethanol: a review. Bioresour Technol. 2010;101:4775–800.

    Article  PubMed  Google Scholar 

  11. 11.

    Liu H, Zhu JY, Fu S. Effects of lignin-metal complexation on enzymatic hydrolysis of cellulose. J Agric Food Chem. 2010;58:7233–8.

    CAS  Article  PubMed  Google Scholar 

  12. 12.

    Saha BD. Hemicellulose bioconversion. J Ind Microbiol Biotechnol. 2003;30:279–91.

    CAS  Article  PubMed  Google Scholar 

  13. 13.

    SGN website (

  14. 14.

    Gur A, Semel Y, Osorio S, Friedmann M, Seekh S, Ghareeb B, et al. Yield quantitative trait loci from wild tomato are predominately expressed by the shoot. Theor Appl Genet. 2011;122:405–20.

    Article  PubMed  PubMed Central  Google Scholar 

  15. 15.

    Lander ES, Botstein D. Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics. 1989;121:185–99.

    CAS  PubMed  PubMed Central  Google Scholar 

  16. 16.

    Darvasi A, Soller M. Selective genotyping for determination of linkage between a marker locus and a quantitative trait locus. Theor Appl Genetics. 1985;1992:353–9.

    Google Scholar 

  17. 17.

    Van Acker R, Vanholme R, Storme V, Mortimer JC, Dupree P, Boerjan W. Lignin biosynthesis perturbations affect secondary cell wall composition and saccharification yield in Arabidopsis thaliana. Biotechnol Biofuels. 2013;6:1–17.

    Article  Google Scholar 

  18. 18.

    De Souza P, Kamei CL, Torres F, Pattathil S, Hahn MG, Trindade LM, et al. How cell wall complexity influences saccharification efficiency in Miscanthus sinensis. J Exp Bot. 2015;66:4351–65.

    Article  PubMed  PubMed Central  Google Scholar 

  19. 19.

    Chen F, Dixon RA. Lignin modification improves fermentable sugar yields for biofuel production. Nature Biotechnol. 2007;25:759–61.

    CAS  Article  Google Scholar 

  20. 20.

    Coleman HD, Beamish L, Reid A, Park JY, Mansfield SD. Altered sucrose metabolism impacts plant biomass production and flower development. Transgenic Res. 2010;19:269–83.

    CAS  Article  PubMed  Google Scholar 

  21. 21.

    Lorenz AJ, Anex RP, Isci A, Coors JG, de Leon N, Weimer PJ. Forage quality and composition measurements as predictors of ethanol yield from maize (Zea mays L.) stover. Biotechnol Biofuels. 2009;2:1–8.

    Article  Google Scholar 

  22. 22.

    Slavov G, Allison G, Bosch M. Advances in the genetic dissection of plant cell walls: tools and resources available in Miscanthus. Front Plant Sci. 2013;4:217. doi:10.3389/fpls.2013.00217.

    Article  PubMed  PubMed Central  Google Scholar 

  23. 23.

    Scheller HV, Ulvskov P. Hemicelluloses. Ann Rev Plant Biol. 2010;61:263–89.

    CAS  Article  Google Scholar 

  24. 24.

    Silveira RL, Stoyanov SR, Gusarov S, Skaf MS, Kovalenko A. Plant Biomass Recalcitrance: Effect of Hemicellulose Composition on Nanoscale Forces that Control Cell Wall Strength. J Am Chem Soc. 2013;135:19048–51.

    CAS  Article  PubMed  Google Scholar 

  25. 25.

    Yoshida M, Liu Y, Uchida S, Kawarada K, Ukagami Y, Ichinose H, et al. Effects of cellulose crystallinity, hemicellulose, and lignin on the enzymatic hydrolysis of Miscanthus sinensis to monosaccharides. Biosci Biotechnol Biochem. 2008;72:805–10.

    CAS  Article  PubMed  Google Scholar 

  26. 26.

    Lima MA, Lavorente GB, da Silva HKP, Bragatto J, Rezende CA, Bernardinelli OD, et al. Effects of pretreatment on morphology, chemical composition and enzymatic digestibility of eucalyptus bark: a potentially valuable source of fermentable sugars for biofuel production – part 1. Biotechnol Biofuels. 2013;6:1–17.

    Article  Google Scholar 

  27. 27.

    Xu N, Zhang W, Ren SF, Liu F, Zhao CQ, Liao HF, et al. Hemicelluloses negatively affect lignocellulose crystallinity for high biomass digestibility under NaOH and H2SO4 pretreatments in Miscanthus. Biotechnol Biofuels. 2012;5:1–12.

    Article  Google Scholar 

  28. 28.

    Saha BD, Iten LB, Cotta MA, Wu YV. Dilute acid pretreatment, enzymatic saccharification, and fermentation of rice hulls to ethanol. Biotechnol Prog. 2005;21:816–22.

    CAS  Article  PubMed  Google Scholar 

  29. 29.

    Macdonald DG, Bakhshi NN, Mathews JF, Roychowdhury A, Bajpai P, Moo-Young M. Alkali treatment of corn stover to improve sugar production by enzymatic hydrolysis. Biotechnol Bioeng. 1983;25:2067–76.

    CAS  Article  PubMed  Google Scholar 

  30. 30.

    Adani F, Papa G, Schievano A, Cardinale G, D’Imporzano G, Tambone F. Nanoscale structure of the cell wall protecting cellulose from enzyme attack. Environ Sci Technol. 2011;45:1107–13.

    CAS  Article  PubMed  Google Scholar 

  31. 31.

    Li F, Ren S, Zhang W, Xu Z, Xie G, Chen Y, et al. Arabinose substitution degree in xylan positively affects lignocellulose enzymatic digestibility after various NaOH/H2SO4 pretreatments in Miscanthus. Bioresour Technol. 2013;130:629–37.

    CAS  Article  PubMed  Google Scholar 

  32. 32.

    The Tomato Genome Consortium. The tomato genome sequence provides insights into fleshy fruit evolution. Nature. 2012;485:635–41.

    Article  Google Scholar 

  33. 33.

    Yokoyama R, Nishitani K. Genomic basis for cell-wall diversity in plants. A comparative approach to gene families in rice and Arabidopsis. Plant Cell Physiol. 2004;45:1111–21.

    CAS  Article  PubMed  Google Scholar 

  34. 34.

    Gardner SL, Burrell MM, Fry SC. Screening of Arabidopsis thaliana stems for variation in cell wall polysaccharides. Phytochemistry. 2002;60:241–54.

    CAS  Article  PubMed  Google Scholar 

  35. 35.

    Morgan MJ, Osorio S, Gehl B, Baxter CJ, Kruger NJ, Ratcliffe RG, Fernie AR, Sweetlove LJ. Metabolic engineering of tomato fruit organic acid content guided by biochemical analysis of an introgression line. Plant Physiol. 2013;161(1):397–407. doi:10.1104/pp.112.209619157:998-1014e.

  36. 36.

    Lande R. The response to selection on major and minor mutations affecting a metrical trait. Heredity. 1983;50:47–65.

    Article  Google Scholar 

  37. 37.

    Eshed Y, Zamir D. A genomic library of Lycopersicon pennellii in L. esculentum: A tool for fine mapping of genes. Euphytica. 1994;79:175–9.

    CAS  Article  Google Scholar 

  38. 38.

    Klinghammer M, Tenhaken R. Genome-wide analysis of the UDP-glucose dehydrogenase gene family in Arabidopsis, a key enzyme for matrix polysaccharides in cell walls. J Exp Bot. 2007;58:3609–21.

    CAS  Article  PubMed  Google Scholar 

  39. 39.

    Chiniquy D, Sharma V, Schultink A, Baidoo EE, Rautengarten C, Cheng K, et al. XAX1 from glycosyltransferase family 61 mediates xylosyltransfer to rice xylan. Proc Natl Acad Sci USA. 2012;109:17117–22.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  40. 40.

    Tang GQ, Lüscher M, Sturm A. Antisense repression of vacuolar and cell wall invertase in transgenic carrot alters early plant development and sucrose partitioning. Plant Cell. 1999;11:177–89.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  41. 41.

    Cocuron JC, Lerouxel O, Drakakaki G, Alonso AP, Liepman AH, Keegstra K, et al. A gene from the cellulose synthase-like C family encodes a beta-1,4 glucan synthase. Proc Ntl Academy Sci USA. 2007;104:8550–5.

    CAS  Article  Google Scholar 

  42. 42.

    Boerjan W, Ralph J, Boucher M. Lignin biosynthesis. Annual Rev Plant Biol. 2003;54:519–46.

    CAS  Article  Google Scholar 

  43. 43.

    Marriott PE, Sibout R, Lapierre C, Fangel JU, Willats WG, Hofte H, et al. Range of cell-wall alterations enhance saccharification in Brachypodium distachyon mutants. Proc Natl Academy Sci USA. 2014;111:14601–6.

    CAS  Article  Google Scholar 

  44. 44.

    Carli P, Arima S, Fogliano V, Tardella L, Frusciante L, Ercolano MR. Use of network analysis as tool to capture key traits affecting tomato organoleptic quality. J Exp Bot. 2009;60:3379–86.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  45. 45.

    Steinhauser MC, Steinhauser D, Gibon Y, Bolger M, Arrivault S, Usadel B, et al. Identification of enzyme activity quantitative trait loci in a Solanum lycopersicum x Solanum pennellii introgression line population. Plant Physiol. 2011;157:998–1014.

    CAS  Article  PubMed  PubMed Central  Google Scholar 

  46. 46.

    Dowe N, Mcmillan J. SSF experimental protocols - lignocellulosic biomass hydrolysis and fermentation. National Renewable Energy Laboratory (NREL), Technical Report TP-510-42630). 2008;1:116.

  47. 47.

    Capareda S. Introduction to Biomass Energy Conversions. Science, Chapter 4. CRC Press, Boca Raton, FL, 2013.

  48. 48.

    Gomez L, Whitehead C, Barakate A, Halpin C, McQueen-Mason SJ. Automated saccharification assay for determination of digestibility in plant materials. Biotechnol Biofuels. 2010;3:1–12.

    Article  Google Scholar 

Download references


This research was funded by the Ministry of Agriculture (BIOMASSVAL project) and a short-mobility grant was provided by the University of Naples Federico II for carrying out analyses at the Department of Biology of York University (United Kingdom).

Author information



Corresponding authors

Correspondence to L. D. Gomez or M. R. Ercolano.

Additional information

Competing interest

The authors declare that they have no competing interests.

Authors’ contributions

GC was involved in analysis, interpretation of data and in manuscript writing; LDG in interpretation of chemical data and paper drafting; FF in experimental field design and phenotypic data analysis; CB contributed to gene annotation and to molecular analysis; AA was involved in chemical monosaccaride analysis process and interpretation; AE was involved in chemical monosaccharide data interpretation; RS contributed to lignin and saccharification assays; SJMM critically revised the manuscript; CD was involved in drafting and revision of manuscript; LF was involved in data interpretation and discussion; MRE conceived the study and was mainly involved in interpretation of data and in manuscript writing. All authors read and approved the final manuscript.

Additional files

Additional file 1: Table S1.

Fruit yield and harvest index of 37 tomato introgression lines. Table S2. Residual biomass and leaf area of tomato introgression lines. (DOCX 113 kb)

Additional file 2: Table S3.

ANOVA analysis on phenotypic and chemical traits. Table S4. Protein families involved in the construction and modification processes of cell wall polysaccharides. Table S5. Tomato cell wall related proteins retrieved in S. pennelli genome region of selected introgression lines. Table S6. List of primers used in qPCR assay. (XLSX 55 kb)

Additional file 3: Figure S1.

Correlation between saccharification rate and crystalline cellulose and between saccharification rate and hemicellulose in 13 tomato introgression lines plus M82. (PDF 13 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Caruso, G., Gomez, L.D., Ferriello, F. et al. Exploring tomato Solanum pennellii introgression lines for residual biomass and enzymatic digestibility traits. BMC Genet 17, 56 (2016).

Download citation


  • Solanum pennellii population
  • Cell wall components
  • Saccharification
  • Biomass conversion
  • Extreme genotypes
  • Candidate genes