Integrating Gene Expression Analysis and Ecophysiological Responses to Water Deficit in Leaves of Tomato Plants

Integrating Gene Expression Analysis and Ecophysiological Responses to Water Deficit in Leaves of Tomato Plants | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Integrating Gene Expression Analysis and Ecophysiological Responses to Water Deficit in Leaves of Tomato Plants View ORCID Profile G Bortolami , View ORCID Profile T de Werk , View ORCID Profile M Larter , View ORCID Profile A Thonglim , View ORCID Profile B Mueller-Roeber , View ORCID Profile S. Balazadeh , View ORCID Profile F. Lens doi: https://doi.org/10.1101/2024.07.05.602262 G Bortolami 1 Naturalis Biodiversity Center, Research Group Functional Traits , PO Box 9517, 2300 RA Leiden, The Netherlands 2 Plant Ecology Research Laboratory, School of Architecture, Civil and Environmental Engineering , 1015, Lausanne, Switzerland Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for G Bortolami T de Werk 3 University of Potsdam, Institute of Biochemistry and Biology , Karl-Liebknecht-Straße 24-25, Haus 20, 14476 Potsdam, Germany 4 Max-Planck Institute of Molecular Plant Physiology , Am Mühlenberg 1, 14476 Potsdam, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for T de Werk M Larter 1 Naturalis Biodiversity Center, Research Group Functional Traits , PO Box 9517, 2300 RA Leiden, The Netherlands 5 BIOGECO, INRA, Université de Bordeaux , 33615 Pessac, France Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for M Larter A Thonglim 1 Naturalis Biodiversity Center, Research Group Functional Traits , PO Box 9517, 2300 RA Leiden, The Netherlands Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for A Thonglim B Mueller-Roeber 3 University of Potsdam, Institute of Biochemistry and Biology , Karl-Liebknecht-Straße 24-25, Haus 20, 14476 Potsdam, Germany 4 Max-Planck Institute of Molecular Plant Physiology , Am Mühlenberg 1, 14476 Potsdam, Germany Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for B Mueller-Roeber S. Balazadeh 4 Max-Planck Institute of Molecular Plant Physiology , Am Mühlenberg 1, 14476 Potsdam, Germany 6 Leiden University, Institute Biology Leiden, Sylvius Laboratory , Sylviusweg 72, 2333 BE Leiden, The Netherlands Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for S. Balazadeh For correspondence: s.balazadeh{at}biology.leidenuniv.nl frederic.lens{at}naturalis.nl F. Lens 1 Naturalis Biodiversity Center, Research Group Functional Traits , PO Box 9517, 2300 RA Leiden, The Netherlands 6 Leiden University, Institute Biology Leiden, Sylvius Laboratory , Sylviusweg 72, 2333 BE Leiden, The Netherlands Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for F. Lens For correspondence: s.balazadeh{at}biology.leidenuniv.nl frederic.lens{at}naturalis.nl Abstract Full Text Info/History Metrics Preview PDF Abstract Soil water deficit (WD) is one of the most important abiotic stresses affecting plant survival and crop yield. Despite its economic relevance, many gaps remain in our understanding of how crops respond to WD, especially concerning the synergistic coordination of molecular and ecophysiological adaptations delaying plant damage and mortality. In this study, we investigated the gene expression imposed by a progressive WD and combined it with measurements pointing to key ecophysiological thresholds in leaves of tomato plants. We uncovered the transcriptomic changes in mature leaves at four stages defined by physiological markers relating to different WD intensities: partial stomatal closure, complete stomatal closure, after leaf wilting, and beginning of embolism development in the veins. By identifying key transcription factors (TFs) across these progressively worsening WD stages, we investigated the timing and impact of ABA-(in)dependent gene regulatory pathways during WD. In addition, we compared the transcriptome in young developing versus mature leaves and explored the physiological mechanisms that may explain the higher tolerance to dehydration in younger leaves. By correlating the transcriptomic changes to precise ecophysiological measurements, the combined dataset will serve as a framework for future studies comparing leaf molecular and physiological responses to WD at specific intensities. Highlight Integrated ecophysiological and gene expression analyses identify key mechanisms underlying the different thresholds of tomato responses to water deficit Introduction The environmental conditions, encompassing all abiotic and biotic factors, define the niche boundaries where a plant species can actively adapt and grow ( Chase & Leibold, 2003 ; Treurnicht et al. , 2020 ). Any deviation in these conditions from the optimal niche range leads to plant stress and a reduction in growth and yield ( Bailey-Serres et al. , 2019 ; Zhang et al. , 2022 ). In the present age of anthropogenic climate change, it is evident that global hydrological patterns are swiftly evolving, leading to diminished precipitation in numerous geographical areas and making them less suitable for agriculture ( Greve et al. , 2014 , 2017 , 2018 ; Wartenburger et al. , 2017 ; Xu et al. , 2019). In a world where plant growth will become more demanding, there is an urgent need to focus more on mechanisms underlying drought responses in crops to safeguard food production ( Pingali, 2012 ; Ray et al. , 2013 ). Tomato ( Solanum lycopersicum L.) is one of the major crop species cultivated in the Mediterranean area where droughts are projected to become more intense and more frequent ( Guerreiro et al. , 2018 ; Raymond et al. , 2018 ; Lionello & Scarascia, 2020 ). Previous research investigating tomato response to water deficit (WD) is scattered and does not provide a complete picture of plants’ sensibility to stress. Available drought studies in tomato either focused on only one single WD intensity ( Gong et al. , 2010 ; Iovieno et al. , 2016 ; Mishra et al. , 2016 ; Zhou et al. , 2019 ; Diouf et al. , 2020 ; Bian et al. , 2021 ; Pirona et al. , 2023 ), used osmotic stress or dehydration as a substitute for diminishing soil water availability ( Dong et al. , 2023 ; Pirona et al. , 2023 ), examined the WD response only during the seedling stage ( Gong et al. , 2010 ; Zhou et al. , 2019 ; Pirona et al. , 2023 ), or investigated the physio- anatomical drivers of xylem embolism development ( Lamarque et al. , 2020 ; Harrison Day et al. , 2022 ; Haverroth et al. , 2024 ). Traditionally, when plant responses to WD are studied, two main approaches are applied to investigate the sequence of events leading to drought-induced plant mortality. The first set of studies focuses on the ecophysiological and/or anatomical approach, looking at modifications of the long-distance water transport pathway and/or changes in xylem anatomical traits imposed by different levels of soil water content ( Lens et al. , 2011 ; Martin- StPaul et al. , 2017; Corso et al. , 2020 ; Lamarque et al. , 2023 ; Lens et al. , 2023). A second line of publications is molecular-oriented and investigates the finer-scale changes in the gene, protein, and chemical chain of reactions during different stages of WD ( Hummel et al. , 2010 ; Taiz et al. , 2024 ). Surprisingly, these two scientific disciplines are typically separated from each other and focus on different plant groups: the hydraulic-anatomical approach deals almost exclusively with tree species, studying the sequence of events triggering the forest decline ( Choat et al. , 2018 ; Mantova et al. , 2022 ), while the molecular approach mainly focuses on non-woody model species (one above all, Arabidopsis thaliana ) subjected to induced stress under controlled growth chamber conditions ( Lamesch et al. , 2012 ; Singh et al. , 2022 ). Cross-pollination between the two disciplines is only sporadic, mostly limited to the quantification of the expression of key genes during the ecophysiological monitoring of species subjected to WD ( Zhang et al. , 2020 ; Thonglim et al. , 2023 ), or constrained to the measurements of leaf water potential (Ψ leaf ) and gas exchange when investigating gene regulatory pathways imposed by WD ( Iovieno et al. , 2016 ; da Silva et al. , 2022 ). Finally, only a few studies performed ecophysiological investigations in plants over- or underexpressing selected genes ( Lamarque et al. , 2020 , Thonglim et al. , 2021 ). Consequently, achieving a deeper comprehension of plant responses to WD demands a more integrated approach, wherein the comparison of gene expression at specific time intervals during WD - based on ecophysiological monitoring - is essential. Here, we performed a holistic drought study on the leaves of tomato plants. We outlined three key objectives: firstly, we coupled precise ecophysiological monitoring (water potential, gas exchange, stomatal and xylem anatomy, xylem embolism spread, turgor loss point) of fully developed leaves in mature (around 50 days old) tomato plants with quantification of gene expression in leaves at four WD intensities: 1) partial stomatal closure, 2) full stomatal closure, 3) leaf wilting and the start of embolism formation, and 4) the point where embolism starts to exponentially increase corresponding to 12% loss of hydraulic conductivity in the main veins. Secondly, we identified and compared the gene expression at these four physiological thresholds. Thirdly, we focused on a leaf developmental stage-specific response, providing a putative transcriptomic basis for the higher resilience of younger tomato leaves to WD. Materials and methods Plant material, sampling, and water deficit experimental design Tomato ( Solanum lycopersicum cv. MoneyMaker) seeds were germinated on full-strength Murashige-Skoog media containing 3% sucrose (MS 3%). After three weeks, plants with the most developed hypocotyls were transferred into 3L pots in a mixture of potting soil (basis biomix, Lensli® substrates, Bleiswijk, The Netherlands), vermiculite and sand in a ratio of 25:8:2, along with three tablespoons of osmocote fertilizer (ICL Growing Solutions, Geldermal, The Netherlands). At the Institute of Biology Leiden (Leiden University), all potted plants were placed in the same growth chamber at 24°C and 70% relative humidity, with a daily 16-h light and 8-h night cycle at 200 µmol/m 2 s -1 photosynthetic active radiation and watered every other day until the beginning of the water deficit experiment. For the RNA-Seq sampling, we divided 50-day-old plants into two groups: a well-watered (WW) control group and a water deficit (WD) group that was fully deprived of water from this point forward. To safeguard against potential plant harm that could impact gene expression analysis while monitoring WD intensity, leaf water potential (Ψ leaf ) was assessed every other day in a subgroup of plants (three WW and three WD individuals). This was done until critical thresholds were reached, based on already available leaf vulnerability curves linking Ψ leaf with drought-induced embolism formation in leaf veins using a previous batch of plants with equal age (see section “Non-invasive optical determination of embolism vulnerability in leaves”). The monitored plants were not utilized for RNA extraction. The Ψ leaf measurements were conducted using a pressure chamber (Model 1000, PMS Instrument Company, Albany, OR, USA) on a leaflet from a fully matured leaf, following the procedure recommended by Rodriguez-Dominguez et al. (2022) . When a specific WD intensity was reached (i.e., at 4, 5, 8, and 10 days after water deficit; Fig. 1 , Suppl. Fig. S1), a mature leaf from the 6-8 th node from the base was used for ecophysiological monitoring, and RNA extraction as detailed in Supplementary Fig. S2. Download figure Open in new tab Figure 1. Sampling strategy and physiological thresholds in 50-day-old S. lycopersicum cv MoneyMaker plants. A) Leaf water potential (Ψ leaf ) variation over time for well-watered (WW, green circles) and water-deficit (WD, yellow circles) plants during the 10d drought experiment. The four sampling points for RNA extraction (at 4, 5, 8, and 10 days after water deficit, d-WD) are indicated with black vertical lines. Symbols represent averages ± SE (n=3-6 per sampling point). B) Stomatal safety margin during Ψ leaf decline, defined by the difference between the best-fitted blue curve for stomata dynamics ( gs , n=33), and the red vulnerability curve (average ± SE) based on leaf percentage of embolized pixels (leaf PEP, n=8). The purple vertical line indicates the turgor loss point (Ψ TLP ), and the black vertical lines indicate the four sampling dates for RNA extraction. Single-leaflet gas exchange measurements and stomata anatomy Maximal leaf gas exchange measurements were registered between 9:00 a.m. and 2:00 p.m. on a well-exposed leaflet close to the one used for RNA extraction (Suppl. Fig. S2) using a TARGAS-1 portable photosynthesis system (PP Systems, Amesbury, MA, USA). Tomato plants (n=33, including the ones used for RNA-Seq analyses) were grown in multiple batches from January until November 2022 as described above, and measured between 45 and 60 days after sowing. Optimal conditions of photosynthetic active radiation were set in the cuvette (1500 μmol photons m -2 s -1 ). The following parameters were recorded on each selected leaflet: stomatal conductance ( gs , mmol H 2 O m -2 s -1 ), assimilation ( A , μmol CO 2 m -2 s -1 ), and water use efficiency ( WUE = A gs -1 , μmol CO 2 mmol H 2 O -1 ). The same leaflet from four well-watered plants was later stored in 70% ethanol for anatomical observations of stomata density and size. Before these stomatal measurements could take place, the samples were further dehydrated using a series of 70%, 96%, and 100% ethanol, submerged in acetone, and placed in the critical point dryer (Leica EM CPD300, Leica Microsystems, Wetzlar, Germany), mounted on scanning electron microscope (SEM) stubs (12.2x10 mm, Ted Pella Inc., Redding, CA, USA) using carbon adhesive tabs (12 mm diameter, Electron Microscopy Sciences, Hatfield, PA, USA), and coated with a Quorum Q150T S Platinum Palladium sputter coater (Quorum Technologies, Laughton, UK). The abaxial leaf side was photographed with an SEM (JSM-6480 LV, JEOL, Tokyo, Japan) at 5 kV. Stomata were counted in six photos per sample, covering a total area of 1.84 mm 2 , and their length was measured on 36 random stomata per sample. Intervessel pit membrane thickness The thickness of the intervessel pit membranes in tomato leaves was measured in the leaflets xylem of five different plants before the determination of vulnerability to embolism, see below. We detached a leaflet from a node close to the one used for embolism detection, we collected a 1 cm-long section and immediately fixed it in Karnovsky’s fixative. As described in Thonglim et al. (2021) , the samples were cleaned three times in 0.1 mM cacodylate buffer, then post-fixed with 1% buffered osmium tetroxide, rinsed again with buffer solution, stained with 1% uranyl acetate, and dehydrated in a series of ethanol: 1% uranyl acetate replacement, with increasing concentration of ethanol (30, 50, 70, 96%, and twice in ≥99%). The samples were then infiltrated with Epon 812n (Electron Microscopy Sciences, Hatfield, UK) and placed at 60°C for 48 h in an oven. The Epon blocks were trimmed to a thickness of 2 μm using a rotary microtome with a glass knife. Subsequently, the cross-sections with many vessel–vessel contact areas were cut into ultrathin sections of 90–95 nm using a Leica EM UC7 ultramicrotome with a diamond knife. The sections were dried and mounted on film-coated copper slot grids with Formvar coating (Agar Scientific, Stansted, UK), and post-stained with uranyl acetate and lead citrate. Ultrastructural observations of intervessel pits were performed and photographed using a JEM-1400 Plus TEM (JEOL, Tokyo, Japan) equipped with an 11-megapixel camera (Quemesa, Olympus). The photos were analyzed with ImageJ (Schneider et al. , 2012), and the thickness was measured in three different spots of the membrane. We measured the thickness of six to fourteen intervessel pit membranes from five different leaflets, for a total of fifty-two membranes. Non-invasive optical determination of embolism vulnerability in leaves Air embolism formation and propagation in the leaf xylem were measured by adapting an optical technique ( Brodribb et al. , 2016 ) monitoring light transmission through the leaf veins. Tomato plants (n=8) were grown in multiple batches from September 2020 until January 2022 as described above and measured between 45 and 60 days after sowing. Plants were moved to the nearby hydraulics laboratory at Naturalis Biodiversity Center, where the soil was gently washed out without damaging the root system. One leaflet (still attached to the plant and in the same position as the one used for RNA-Seq or water potential measurements, Suppl. Fig. S2) was fixed to a scanner (Perfection V850, Seiko Epson Corporation, Suwa, Japan) and scanned every five minutes during one week of dehydration using an automating IT script developed by platform Phenobois Bordeaux, while the stem water potential (Ψ stem ) was measured at least once a day with the pressure chamber using a bagged leaflet. The resulting images were analyzed with ImageJ software (Schneider et al. , 2012) following instructions from http://www.opensourceov.org (last consulted November 2023). The formation of embolisms over time was determined by subtracting the differences in pixels in the 1 st - to 3 rd -order veins among the subsequent image scans. Background noise, mainly caused by tissue shrinkage, was removed using mild filters and manual image inspection. Vulnerability curves corresponding to the percentage of embolized pixels (PEP) as a function of Ψ were fitted for every sample based on the following equation ( Pammenter and Van der Willigen, 1998 ): PEP = 100 / {1 + exp [S/25 x (Ψ-P 50 )]} where P 50 is the xylem pressure at which 50% of embolized pixels were observed and S is the slope of the curve at the inflection point. Pressure-volume curves and determination of turgor loss point We determined the leaf turgor loss point (Ψ TLP ) in four leaflets from two well-watered plants right after the last sampling point using the pressure-volume curve method ( Sack et al. , 2010 ). Briefly, plants were watered at pot capacity, and the next day a leaflet from a mature fully developed leaf from the 6-8 th node (comparable to the one used for RNA-Seq, Suppl. Fig. S2) was detached and equilibrated inside a sealed plastic bag. The leaf was weighed with a precision scale (Practum224-1S, Sartorius, Göttingen, Germany), its water potential was measured with the pressure chamber, and afterward air-dehydrated in the lab. Starting from the onset of dehydration, its weight and water potential were measured two to three times per day for three days until complete dehydration. Pressure-volume relationships were constructed for every sample by plotting the inverse of leaf water potential (–1/Ψ leaf ) against the relative water content (RWC), and Ψ TLP estimated as the point of transition between the linear and nonlinear portions of the pressure-volume relationship. RNA extraction and library preparation A leaflet from a mature and fully developed leaf (Suppl. Fig. S2) was taken from the 6-8 th node from the base and flash frozen in liquid nitrogen. In the latest sampling date (10 d-WD ), also young developing leaves were sampled from the topmost (closest to the shoot apical meristem) secondary shoots. All samples were collected between 2 p.m. and 3:30 p.m. In each case, the target leaflets from two plants were combined to form a single biological replicate, and each time point contained two to three biological replicates for the WW and WD batch. Tissues were ground into a fine powder using a Tissuelyzer (QIAGEN), ensuring the samples remained frozen. Total RNA was extracted using TRIzol® Reagent (Thermo Fisher) according to the manufacturer’s instructions. Total RNA purity was assessed using 260/280 and 260/230 nm wavelength absorption ratios, using a NanoDrop™ spectrophotometer (DeNovix DS-11 FX+). Library preparation and paired-end sequencing were performed by BGI Genomics (China). In total, 27 samples had RNA of sufficient quality and quantity to be sequenced (Suppl. Table S1). Differential expression and co-expression analyses Read trimming, alignment of reads to the genome, and mapping of reads to gene models were performed using CLC Genomics Workbench (Qiagen). Reads were aligned to the SL4.00 genome release, using ITAG4.1 gene annotation data (obtained from the solgenomics FTP site: https://solgenomics.net/ftp/tomato_genome/ , accessed in June 2023). Automated human readable gene descriptions (AHRD) corresponding to the ITAG4.1 genome were downloaded from the Solgenomics website and added manually (Solgenomics, 2020). Differential expression analyses were performed in R (version 4.3.1), using the DESeq2 package ( Love et al. , 2014 ) between WD and WW samples on every sampling date. Correction for the overestimation of effect sizes for lowly expressed genes was done with a log 2 fold change (LFC) shrink, as described by Zhu et al. (2019) . Differentially expressed genes were defined as having an LFC difference of at least 1 in the WD plants compared to the control, and a Benjamini-Hochberg FDR adjusted p-value (p adj ; as given by DESeq2) below 0.05 (Suppl. Table S2). Differential gene regulation between developing and mature leaves was determined using the DESeq2 package, calculating differential regulation on the interaction between leaf age and treatment while controlling separately for both leaf age and treatment. Genes were defined as differentially regulated using a stringent filter (LFC > |2|; p adj < 0.05, Suppl. Table S3) to account for the increased number of independent variables and reduce the chance for false positives. K-means clustering was performed on Variance Stabilizing Transformation (VST)-normalized counts, as obtained from the raw read counts using DESeq2 in R. Clustering accuracy was validated manually by plotting data grouped by clusters along its three principal components (Suppl. Fig. S3). Weighted gene correlation network analysis (WGCNA), co-expression analysis, and identification of water deficit markers Normalized counts were obtained via a VST from the unaltered read counts using the DESeq2 package. The weighted gene co-expression network analysis (WGCNA) was performed on differentially expressed genes using the WGCNA package in R (version 4.3.1), using a soft power threshold of 20, generating a signed network on the previously normalized counts (chosen as described in Langfelder and Horvath, 2008 ; Suppl. Fig. S4). Cluster-Trait relationships (using watering regime, sampling date, position of the leaf on the plant, Ψ leaf , gs , A , and leaf PEP as functional traits) were calculated from the cluster eigengenes as described in Langfelder and Horvath (2008) . Genes encoding for transcription factors (TFs) were identified from the GO term “DNA-binding transcription factor activity”. TFs of interest were subsequently identified as those that presented an LFC higher than 4 when summed across all time points. The co-expression network was generated in Cytoscape using the CoExpNetViz app ( Tzfadia et al. , 2016 ), using previously identified transcription factors of interest as bait genes. Transcription factors of interest were taken as being differentially expressed in WD plants compared to WW plants (LFC > |1|; p adj -value < 0.05) at any time point, while the networks were generated on all differentially expressed genes irrespective of time point. A stringent filter was applied for significant correlations (p- value ≤ 0.05), with a Pearson’s correlation (r) ≥ |0.95| used to generate the co-expression network. To identify which nodes exerted the most influence on the network - and thus would represent potential transcription factors of interest - we additionally calculated betweenness- centrality between all nodes in Cytoscape. The number of correlated genes for each transcription factor of interest was extracted from the topological overlap matrix, using a weight threshold of 0.4. To visualize the gene expression clusters, z-scores were calculated from the normalized expression values by centering each value on the mean expression and scaling it by the standard deviation of the entire dataset. Results and discussion Ecophysiological responses of tomato leaves to increasing water deficit We monitored ecophysiological changes in plants subjected to a 10d period of water deprivation, while simultaneously taking samples for gene expression analyses at specific time points determined by the monitoring ( Fig. 1 ). Under well-watered (WW) conditions where the leaf water potential (Ψ leaf ) was above -0.4 MPa, the plants showed a leaflet stomatal conductance ( gs ) of 392±18 (average±SE) mmol H 2 O m -2 s -1 and a CO 2 assimilation ( A ) of 8.9±0.4 µmol m -2 s -1 . These leaflets showed 154±5 stomata mm -2 with a length of 18.35±0.2 µm at the abaxial side and an intervessel pit membrane thickness of 0.42±0.01 μm. As the soil water content decreased, Ψ leaf declined as well ( Fig. 1A ) and stomata closed ( Fig. 1B ) at Ψ leaf =-0.7 MPa. The turgor loss point in leaves (Ψ TLP , the key physiological indicator of a plant’s ability to maintain cell turgidity in leaves, defined when the pressure in- and outside the cells is equal) was identified at -1±0.04 MPa (purple line in Fig. 1B ). At a later stage of WD, air embolism appeared and spread in the leaf veins, causing 12% increase of embolized pixels (12% PEP, or P 12 ) at -1.41±0.1 MPa, while the average P 50 and P 88 were detected at -1.75±0.1 and -2.08±0.14 MPa, respectively ( Fig. 1B ). The stomatal safety margin (i.e., the xylem pressure range between stomatal closure and the P 50 ) corresponded to 1.05 MPa. To elucidate the molecular processes underlying the observed physiological changes, we extracted RNA at distinct physiological stages for a comparative transcriptomic analysis through RNA sequencing ( Fig. 1 , 2). To this end, we sampled comparable leaflets from mature leaves for the WW condition and at 4, 5, 8, and 10 days after water deficit (d-WD). These sampling dates corresponded to precise physiological thresholds: at 4d-WD the stomata started to close ( gs 10 , Ψ leaf =-0.41 MPa), at 5d-WD the stomata were completely closed ( gs 90 , Ψ leaf =-0.73 MPa), at 8d-WD the leaves appeared wilted and embolism started to develop in the xylem (4-6% PEP, Ψ leaf =-1.3 MPa), and at 10d-DW, P 12 was reached (Ψ leaf =-1.43 MPa). In addition, at 10d-WD, we also sampled young developing leaves closer to the stem apex that showed no wilting despite being under a lower Ψ leaf compared to mature leaves on the same plant (Ψ youngleaf =-1.61 versus Ψ matureleaf =-1.43 MPa on average). Transcriptomic profiling during increasing water deficit Over 1 billion paired clean reads were successfully mapped across all samples to the SL4.00 genome release, with 93.35% to 95.47% of counted fragments uniquely mapped to genes according to the ITAG4.1 gene annotation release (Suppl. Table S1). A multidimensional scaling plot shows that all the WW samples cluster close together, showing comparable transcriptomic profiles and indicating that our control samples are similar ( Fig. 2A ). In contrast, a clear progressive separation over Dimension 1 (representing 61% of the observed variance) is seen for WD samples, correlated with progressively increasing WD ( Fig. 2A ). Similarly, hierarchical clustering analysis of the top 500 most variable genes across all samples clustered all replicates to their respective time points and conditions ( Fig. 2B ), highlighting the low variability between biological replicates and the overall good quality of the sampling. Download figure Open in new tab Figure 2. Transcriptomic profiling of water deficit (WD) and well-watered (WW) tomato plants. A) Multidimensional scaling plot visualizing (dis)similarity between WD and WW samples. Symbols represent biological replicates, colors represent the different treatments (yellow for WD, green for WW conditions), and different shapes represent the different sampling points (circle for 4d-WD, square for 5d-WD, rhombus for 8d-WD, triangle for 10d- WD, and plus signs for young leaves at 10d-WD). B) Expression heatmap and hierarchical clustering of the top 500 most variable genes. Each column corresponds to one replicate: WW or WD, mature (M) or young (Y) leaves, and different clusters are separated by white lines. All biological replicates cluster to their respective time points and treatments. Color is determined by z-score, ranging from -3 (blue) to 3 (red). C) Volcano plots showing differentially expressed genes between WW and WD conditions in mature leaves for every time point. Differential expression is defined as having a log 2 fold change over WW conditions of at least 1, and a Benjamini-Hochberg FDR-adjusted p -value below 0.05. We compared gene expression in the well-watered and drought-stressed plants for each time point. In total, we found 7179 genes that were differentially expressed (DE) in at least one time point (Suppl. Table S2). At 4 days after water deficit (4d-WD at 10% reduction in gs ), a total of 379 genes were differentially expressed (232 genes were down- and 147 genes were upregulated; Fig. 2C ). At 5d-WD at 90% reduction in gs (i.e., the stomata were fully closed) we found 1637 DE genes (1055 down- and 582 upregulated; Fig. 2C ). At 8d- WD after the leaves had fully wilted and embolism started to develop, we found 5531 DE genes (3115 down- and 2416 upregulated; Fig. 2C ). Lastly, at 10d-WD, when the vulnerability curve started to show an exponential increase in embolism events and the percentage of embolized pixels (PEP) reaches 12% (P 12 ; Fig. 1B ), we found 4637 DE genes (2921 down- and 1716 upregulated; Fig. 2C ). There was a pronounced difference in the number of pre- and post-stomatal closure DE genes (Suppl. Fig. S5). Most DE genes were unique to each time point, highlighting the particular importance of transcription for the specific WD intensity. However, 105 genes were commonly up- or down-regulated during all WD intensities, whereas a much larger portion (1011) was differentially expressed in post-stomatal closure samples (from 5d- to 10d-WD), indicating an important shift in the leaf physiology imposed by (and/or stimulating) stomatal closure. Similarly, a large fraction of DE genes (2514) is shared across the post-wilting samples, when embolism starts to develop in leaf xylem. Next, we examined the clustering and expression patterns of genes via a weighted gene correlation network analysis (WGCNA, Fig. 3 ). Firstly, we restricted our analysis only to the previously identified 7179 genes that showed differential expression in at least one time point (Suppl. Table S2). To effectively characterize gene expression under increasing WD, we used a pseudo-time series, taking 4d-WD well-watered plants as basal condition, following up with WD samples for subsequent time points (4-10-d WD). Nineteen clusters of inter-correlated genes were identified ( Fig. 3A ). We correlated the cluster eigengenes to the physiological traits we measured on the same plants ( Fig. 3B ), showing their change over increasing WD. Three groups of clusters were easily detectable. The first group (clusters i-vi) showed a positive relationship between gene expression and Ψ leaf , gs , and A ( Fig. 3B ), a negative relationship with the increase of embolism (leaf PEP, Fig. 3B ), and a general downregulation during increasing WD ( Fig. 3C ). The second group (clusters vii-xiii) showed the opposite trend: a negative relationship between gene expression and Ψ leaf , gs , and A ( Fig. 3B ), a positive relationship with the increase of embolism (leaf PEP, Fig. 3B ), and a general upregulation during WD ( Fig. 3C ); the third group (clusters xiv-xix) showed an intermediate regulation compared to the two previous groups and a significant relationship with few or no recorded traits. Download figure Open in new tab Figure 3. Weighted gene correlation network analysis (WGCNA). A) Dendrogram denoting WGCNA results for all previously identified DE genes. The y-axis (Height) represents the distance, or dissimilarity between clusters. Nineteen distinct co-expression clusters were identified, with different colors and Roman numerals (i-xix). Gene cluster designations can be found in Suppl. Table S2. B) Cluster-trait relationship heatmap, correlating the module eigengenes (ME) to the measured physiological traits: leaf water potential (Ψ leaf ), stomatal conductance ( gs ), CO 2 assimilation ( A ), and loss of hydraulic conductance in leaves (leaf PEP). Correlations are plotted on a diverging color scale centered on 0 (yellow), with red denoting a positive correlation and green a negative correlation. C) Expression Z-scores for each cluster, using scaled and centered VST normalized counts over increasing water deficit from basal (WW) conditions, consisting of WW plants on the first sampling date, to 10d-WD. Co-expression gene network analyses indicate key regulatory genes involved in the different stages of the water deficit response To further explore the transcriptional regulation of the WD response, we prioritized identifying the most vital transcription factors (TFs) given their pivotal roles in regulating gene expression. We argued that TFs are important for two primary reasons: 1) they could exhibit the highest differential expression at a given time point and 2) they could exert a high influence on the network. To this end, we calculated centrality measures (degree and betweenness centrality) for all previously identified TFs of interest, in addition to examining differential expression per time point ( Fig. 4 ; Suppl. Fig. S6). Two TFs in particular ( Solyc01g095460.3 and Solyc05g050220.3 ; Suppl. Fig. S6) presented the highest number of correlations and betweenness centrality over the four sampled dates. These two TFs belong to the G-box binding factor (GBF) family and their expression increases with increasing WD (cluster viii). They are close homologs of AtGBF3 , an abscisic acid (ABA)-activated Arabidopsis gene that confers increased tolerance to drought when overexpressed (Lu et al. , 1996; Ramegowda et al. , 2017). ABA is a critical hormone in abiotic stress signaling known to accumulate during WD stress ( Waadt et al. , 2022 ). Our data suggest that both genes play a pivotal role in the drought response of tomato leaves, likely by triggering a cascade of events related to the ABA-dependent response to WD. Download figure Open in new tab Figure 4: Transcription factor expression and number of correlated genes over a progressive water deficit. Transcription factors are separated by which sampling time point they first break the differential expression threshold (log 2 fold change > 1; p adj < 0.05), focusing only on upregulated transcription factors. In the columns (from left to right) we indicated: a histogram showing the distribution of correlation strengths (weight) of correlated genes, grouping weights within 0.01 difference together; the gene cluster (see Fig. 3 ); the number of correlated genes from 0 (light yellow) to 2500 (green); the differential expression for each transcription factor across time points is given in log 2 fold change spanning from -8 (blue) to 8 (red) and centered on 0 (white); and the locus ID of the gene. Subsequently, we looked for TFs that drive the WD response at the specific time points. We subdivided the previously identified TFs of interest by the time point at which they first show significant upregulation ( Fig. 4 ). At 4d-WD, at the onset of stomatal closure ( Fig. 1B ), the TFs presenting the highest number of correlations (all belonging to cluster xvi) followed an ABA-independent molecular pathway: Solyc03g124110.2 and Solyc03g026280.3 are two genes presenting a strong difference in expression compared to control plants (4.21 and 2.65 LFC, respectively); they are close homologs of Arabidopsis AtDREB1 (one of the key TFs driving the ABA-independent response to drought; Akhtar et al. , 2012). Alongside, Solyc04g078640.3 , Solyc12g009240.1 , Solyc03g093540.1 , and Solyc03g093550.1 belong to the Ethylene Response Factor (ERF) family, members of which are often, but not always, repressed by ABA ( Müller and Munné-Bosch, 2015 ). It has been shown that hydraulic signals (i.e., changes in water potential) regulate stomatal response to WD before ABA accumulates in the stressed leaves imposing a more stable stomatal closure ( Tombesi et al. , 2015 ; Huber et al. , 2019 ; Abdalla et al. , 2021). Our data corroborates this, as the ABA- independent TFs were driving the co-expression network in tomato leaflets that showed a reduction of 10% in stomatal conductance and a Ψ leaf =-0.41 MPa. At 5d-WD, the time point of full stomata closure, the leaves were not yet wilting, and the onset of embolism did not yet occur in the leaf xylem ( Fig. 1 ). The TFs presenting the highest number of correlations belong to cluster viii (i.e., increasing linearly with WD intensity) and are mostly related to the ABA-dependent pathway ( Fig. 4 ). More specifically, they included Solyc01g095460.3 and Solyc05g050220.3 , close homologs of AtGBF3 described above, and Solyc04g078840.3 ( AREB1 , ABA-responsive element binding protein 1), one of the most important TFs in the ABA-dependent response to drought in many plant species ( Fujita et al. , 2005 ; Barbosa et al. , 2012; Pan et al. , 2023 ). We also retrieved Solyc12g044390.3 , Solyc06g066540.1 , and Solyc03g120840.3 , close homologs of AtDREB3 , known to be involved in abiotic stress tolerance ( Islam and Wang, 2009 ). In addition to the genes involved in the ABA-dependent response to drought, we found two genes involved in the biosynthesis of antioxidant molecules: Solyc10g079150.3 and Solyc10g081840.3 are two SlNF-YA TFs involved in flavonol biosynthesis ( Li et al. , 2016 ; Wang et al. , 2021). Furthermore, Solyc07g040680.3 and Solyc06g053960.3 encoding the two heat shock proteins Sl HsfA2 and Sl HsfA6b, respectively ( Huang et al. , 2016 ), were highly correlated as well. In summary, this shift in expression occurring at 5d-WD likely indicates an accumulation of ABA in the leaves of (rather sensitive) tomato plants that become more drought-stressed. The flavonol biosynthesis TFs point to the beginning of the antioxidant response, probably limited to the most sensitive leaf mesophyll tissues ( Watkins et al. , 2017 ), and the upregulation of heat shock proteins indicates an increase in leaf temperature resulting from fully closed stomata that strongly reduce water loss and leaf cooling. At 8d-WD, the leaves were wilting, and embolism events started to develop in the leaf xylem (4-6% PEP). We observed an increase in the number of correlations and betweenness centrality of additional heat-shock responsive TFs Solyc03g097120.3 , Solyc04g016000.3 , and Solyc06g050520.3 indicating increasing leaf temperature and the pronounced role of heat shock proteins in severe WD ( Fig. 4 ; Mishkra et al. , 2002; Hwang et al. , 2012 ). Three WRKY TFs were also strongly upregulated, probably related to the antioxidant response (see SlWRKY01 , SlWRKY24 , and SlWRKY35 below). In addition, Solyc01g100460.3 , Solyc01g104650.3 , and Solyc02g092090.3 belonging to the bZIP family were highly correlated, probably involved in C and N metabolism reprogramming ( Dietrich et al. , 2011 ). Aside from the analysis of TFs, we observed that genes encoding enzymes for the biosynthesis of trehalose ( Solyc02g071590.3 , Solyc07g006500.3 ) and proline ( Solyc06g019170.3 , Solyc02g068640.3 ) were upregulated from 8d-WD onwards. These are considered two of the most important osmoprotectants, meaning that they regulate the internal cell pressure against the dehydrating apoplast ( Iturriaga et al. , 2009 ; Szabados and Savoure, 2010). This indicates that, in these later stages of stressful WD conditions, retaining residual water as much as possible is vital for leaf mesophyll cells. Interestingly, enzymes involved in proline degradation ( Solyc02g089620.3 ; Solyc02g089630.3 ) were already downregulated at 4d-WD, demonstrating that the osmoprotectant response is one of the earliest molecular responses to WD in tomato plants. At 10d-WD, embolism in mature leaves reached P 12 , corresponding to the point after which embolism started to spread exponentially in the leaf veins. Fewer TFs were upregulated, among which Solyc02g090310.1 and Solyc05g010516.1 presented the highest number of correlations ( Fig. 4 ). Only one of them is well described ( Solyc02g090310.1 ); its Arabidopsis homolog, AtDOF1 (Dof zinc finger protein1), is activated by oxidative stress ( Yanagisawa, 2002 ), indicating the strong presence of reactive oxygen species (ROS) molecules in highly stressed leaves. The role and timing of four transcription factor families during water deficit in mature tomato leaves As an alternative approach to disentangle the transcriptome profiles across the different time points during our tomato drought experiment, we identified genes from four key TF families that presented a significantly higher or lower expression than in control plants ( Table 1 ). View this table: View inline View popup Table 1. Response of main TF families in tomato along with the genes that were differentially expressed in plants under water deficit in at least one time point compared to control plants (i.e., genes presenting log2 fold change > |1| and p adj -value < 0.05). For every gene, its cluster of expression (see Fig. 3 ), its up- (+), down- (-), or no- (=) regulation for different days after drought (4, 5, 8, and 10d-WD) compared to WW plants, and locus ID are indicated, as well as the study describing the function of the different loci. First, the ABF / AREB ( ABRE BINDING FACTOR / ABA-RESPONSIVE ELEMENT BINDING PROTEIN ) genes encode important transcription factors involved in ABA signaling; they are generally upregulated in the response of plants to abiotic stresses ( Pan et al. , 2023 ). We found three genes of this family to be upregulated during WD, confirming their role in the response of tomato to drought ( Table 1 ). Second, the WRKY family of transcription factors has been described to be involved in plants’ growth development and response to a wide range of abiotic and biotic stresses (Chen et al. , 2017). In our experiment, the majority of the differentially regulated WRKY genes were found to be downregulated during WD (i.e. 27 genes belonging to clusters i to iv, Table 1 ), indicating a shift towards WRKY deactivation, possibly reducing cell metabolic activities. For example, SlWRKY81 has been described as a negative regulator of drought tolerance by maintaining open stomata and repressing proline biosynthesis during water shortage ( Ahammed et al. , 2020a ; Ahammed et al. , 2020b ). In our experiment, expression of SlWRKY81 was strongly downregulated at the highest intensities of WD. Conversely, three other WRKY genes ( SlWRKY01 , SlWRKY24 , and SlWRKY35 ) previously described to play a role in biotic stress response and carotenoid biosynthesis in tomato ( Moghaddam et al. , 2019 ; Yuan et al. , 2022 ), were upregulated during increasing WD. Our results also suggest that these three TFs might be involved in triggering and/or intensifying the antioxidant response of leaves during drought. Four other genes in this family ( SlWRKY11 , SlWRKY52 , SlWRKY54 , and SlWRKY55 ) exhibited an intermediate regulation, being upregulated at 5d- WD (at stomatal closure) and downregulated at 8d-WD (leaf wilting and onset of xylem embolism) and 10d-WD (at P 12 ; cluster xvi, Table 1 ). SlWRKY52 was previously described to have a role during drought stress ( Jia et al. , 2023 ), whereby its overexpression generally increased tolerance to drought and osmotic stress in tomato plants. Our results indicate that SlWRKY52 is activated only during stomatal closure and not at higher WD intensities, confirming the hypothesis of Jia et al. (2023) of its potential role in ABA biosynthesis and signal transmission. Third, the NAC transcription factor family included 7 genes that were downregulated with increasing WD intensity, while 15 genes were upregulated during WD ( Table 1 ). Of the latter ones, SlNAC63 / RD26 and SlNAC90 were previously described to be upregulated during WD ( Iovieno et al. , 2016 ), confirming our results showing an expression that increased linearly with increasing stress. Interestingly, three other NAC genes ( SlNAC64 , SlNAC78 , SlNAC88 ) were upregulated only at 5d-WD ( Table 1 ), but downregulated at the other WD intensities, suggesting a specific role during the earlier WD stages at stomatal closure. Fourth, the nuclear transcription factor Y (NF-Y) family encodes genes that operate in a multiprotein complex consisting of three subunits (NF-YA, NF-YB, and NF-YC). While these TFs are mostly described for their role in flavonol biosynthesis during fruit ripening in tomato ( Li et al. , 2016 ; Wang et al. , 2021), they were observed to be involved in growth and flower development in multiple plant species ( Petroni et al. , 2012 ). Interestingly, SlNF-YA9 , SlNF- YB8 , SlNF-YC1 , and SlNF-YC9 showed a progressive increase along the increasing WD stages ( Table 1 ), indicating a possible role of these TFs in the antioxidant response in stressed leaves. Likewise, we also found SlNF-YA10 , which was shown to be involved in increased tolerance to oxidative stress in leaves and fruits of tomato plants ( Chen et al. , 2020 ), to be upregulated during increasing WD. Finally, two more NF-Y TFs ( SlNF-YC04 and SlNF-YC10 ) were upregulated exclusively at 8d-WD, the time point of the onset of xylem embolism in the leaves. Comparative transcriptomic analyses explain stronger water deficit tolerance of young compared to mature leaves Young developing leaves had a more negative water potential than mature and fully expanded leaves, both in WW (Ψ leaf =-0.5 ± 0.05 MPa versus -0.44 ± 0.07 MPa), and in WD (Ψ leaf =-1.6 ± 0.04 MPa versus -1.4 ± 0.07 MPa) conditions. This observation is consistent with the cohesion-tension theory, where an increasingly negative water potential pulls water from the bottom to the top of the plants ( Dixon and Joly, 1895 ). Interestingly, despite this more negative Ψ leaf , the young developing leaves at 10d-WD did not show any sign of wilting and/or yellowing, while the mature leaves were well below their wilting point (Suppl. Fig. S1). From an evolutionary standpoint, natural selection may lead to prioritizing the survival of younger leaves over older ones as older leaves are often less efficient in photosynthesis due to age-related decline ( Procházková & Wilhelmová, 2007 ; Wang et al. , 2014 ). This idea is also corroborated by studies highlighting the allocation of resources to younger tissues during stress- and season-related leaf senescence ( Schippers et al. , 2015 ; Berens et al. , 2019 ). Likewise, it has been shown that leaves produced later in the season have an increased WD tolerance compared to leaves produced in the early spring in several woody species ( Sorek et al. , 2021 ; Sorek et al. , 2022 ). We aimed to elucidate the mechanisms behind this apparent resilience of younger leaves by comparing the transcriptomic profiles between young and mature leaves under WW and 10d-WD conditions in the same plants. We identified 360 differentially regulated genes ( i.e. , genes with an LFC (log 2 fold change) above 2, or below -2, and an FDR-adjusted p -value below 0.05) from the interaction between differentially expressed genes between leaf age and watering regime (Suppl. Table S3). We grouped these genes into clusters of similar patterns of expression via K-means clustering. We identified seven distinct gene expression clusters ( Fig. 5 ), which potentially delineate a functional separation in the response to WD. To avoid confusion with our previous clustering related to the WD response in mature leaves only, we named these clusters 1 to 7. To further dive into the potential biological significance of these expression patterns, we performed a Gene Ontology (GO)-enrichment analysis, a manual examination of the AHRD gene descriptions, as well as an examination of genes commonly found to be involved in the WD response ( Kido et al. , 2019 ; Zhang et al. , 2022 ; Mishra et al. , 2023 ). Below, we limit our discussion to the potential role of the genes in cluster 1 (upregulation of 26 genes in young leaves at 10d-WD) and cluster 5 (downregulation of 97 genes in young leaves at 10d-WD; Fig. 5 ). Download figure Open in new tab Figure 5: Comparative transcriptomic analysis of young and mature leaves under well-watered and water-deficit conditions. A) Heatmap showing z-scores of differentially regulated genes between young and mature leaves under well-watered (WW) and water deficit (WD) conditions (10d-WD). Values are represented by a diverging color scale from low (blue) to high (red), centered on 0 (light yellow). Clusters were determined by K-means clustering and cluster assignments are given on the left. Gene cluster assignments are given in Suppl. Table S3. B) Line plots detailing average (± SE) gene expression in each cluster for young and mature leaves under WW and WD conditions. Values are given in Variance Stabilized Transformed (VST) counts. The average gene expression for mature leaves is shown in yellow, and for young leaves in blue. Cluster 1 contains genes that show an upregulation in WD-young leaves, while they are neither responsive in mature leaves (WW and WD) nor in WW-young leaves ( Fig. 5 ). Thus, this cluster could potentially aid in conferring a WD tolerance specifically to younger leaves. For instance, we found that cluster 1 is enriched for genes involved in cell wall biosynthesis (Suppl. Fig. S7) with multiple xyloglucan endotransglucosylase enzymes (XTH; Solyc03g093120.5, Solyc03g093080.3, Solyc03g093130.3 ). XTHs catalyze the cleavage and (often) re-ligation of xyloglucan molecules to cellulose and can have dual roles in cell wall rigidity ( Kaku et al. , 2004 ; Miedes et al. , 2013 ). In Arabidopsis, constitutive overexpression of hot pepper ( Capsicum annuum ) CaXTH3 enhanced tolerance to salinity and WD, likely through increased cross-linking of xyloglucan molecules, thereby decreasing plasticity and elongation of the cell walls ( Cho et al. , 2006 ; Hayashi & Kaida, 2011 ). It is plausible that increased cell wall deposition and rigidity in younger leaves contribute to the maintenance of turgor pressure in smaller cells, at the expense of cell division and expansion (i.e., ceased leaf growth). This hypothesis is corroborated by the increased expression in young leaves of cellulose synthase ( Solyc02g072240.3 ) and fasciclin-like arabinogalactan protein ( Solyc11g069250.2 ) in cluster 4, and the decreased expression of expansins ( Solyc09g018020.3, Solyc06g049050.3 ) in cluster 5. Cluster 1 also contains a callose synthase ( Solyc11g005980.3 ). Callose deposition is intricately connected to plasmodesmata permeability and is tightly controlled by two antagonistic enzymes: the synthesizing callose synthases/glucan synthase-like protein (CS/GSL), and breakdown β-1,3-glucanases (BGs) ( Wu et al. , 2018 ). Interestingly, one of the closest homologs of Solyc11g005980.3 in Arabidopsis ( AtGSL12 / AtCalS3 ) has been connected to callose deposition in sieve elements in the phloem, where increased activity led to clogging of the plasmodesmata pores and impaired phloem loading ( Vatén et al. , 2011 ; Wu et al. , 2018 ). In our dataset, one BG ( Solyc01g060020.4 ) is strongly downregulated during WD uniquely in young leaves. The increased expression of the callose synthase gene Solyc11g005980.3 , together with the decreased expression of the callose breakdown gene Solyc01g060020.4 in young tomato leaves during WD conditions suggests a combined approach to stimulate callose deposition in the phloem thereby blocking outward sugar transport. Soluble sugar accumulation during WD has been related to an increased osmoprotective power delaying leaf wilting (Ozturk et al. , 2020), as well as a ready source of energy during the antioxidant response ( Sami et al. , 2016 ). Rosa et al. (2009) highlighted a whole chain of gene up- and down-regulations coordinated by soluble sugars during abiotic stresses. Our transcriptomic profile also suggests that young leaves accumulate sugars, which in turn might be related to (if not coordinating) the accelerated response to dehydration in these developing leaves. Cluster 5 contains genes that show a much higher expression in young leaves under WW conditions, while they are downregulated to the same level as in mature leaves during WD ( Fig. 5 ). This cluster may contain genes that are “shut off” during WD in younger leaves. Aside from the previously mentioned expansins, this cluster is enriched with nine genes encoding for chlorophyll a/b binding proteins ( Solyc08g067330.1, Solyc02g070980.1, Solyc06g069730.3, Solyc02g070970.1, Solyc02g070990.1, Solyc02g071000.1, Solyc02g070950.1, and Solyc02g071010.1 ). The increased expression of the chlorophyll a/b binding protein genes under WW conditions is indicative of leaves undergoing active growth, as actively proliferating cells build their photosynthetic machinery ( Cackett et al. , 2022 ). The reduction in expression of these genes during WD coincided with the shifting in cell wall structure (and consequently growth cessation) highlighted in cluster 1, indicating that photosynthesis as well is affected at this WD intensity. Moreover, even under well-watered conditions, the photosynthetic apparatus is one of the major producers of ROS, which is exacerbated by abiotic stress (Asada et al. , 2006; Gill & Tuteja, 2010 ). Therefore, the capacity of young leaves to more effectively reduce photosynthetic activity compared to older leaves may also help to reduce the oxidative pressure under WD. We also found that young leaves maintain an overall higher presence of negative regulators of peptidase activity under WD. Two are upregulated uniquely in young leaves under WD (cluster 1; Solyc09g083435.1, Solyc03g098795.1 ) and three remain highly expressed in young leaves, while they are downregulated in mature leaves (cluster 2), namely one cathepsin D inhibitor protein ( Solyc03g098790.3 ), and two protease inhibitors ( Solyc09g083435.1, Solyc03g098795.1 ). In addition, eleven prote(in)ase inhibitors maintained a significantly higher expression only in WD-young leaves, compared to just three which maintained a higher expression in WD-mature leaves. Proteases are incredibly varied, and their functions may be both detrimental and beneficial for acquiring tolerance to WD ( Jedmowski et al. , 2014 ; Cheng et al. , 2016 ; Malefo et al. , 2020 ; D’Ippólito et al. , 2021 ; Moloi & Ngara, 2023 ). The increased presence of prote(in)ase inhibitors activity in young leaves under WD could denote decreased protease activity. Although increased protease activity may aid in maintaining a functional protein/enzyme pool in the face of oxidative stress by increasing protein turnover, having a high protein turnover is resource-intensive and may not be sustainable in a prolonged drought. Considering the time point of measurement (10d-WD) it would be plausible that there is a focus on maintaining current proteins, specifically in young leaves. Conclusions By disentangling gene expression profiles at precise physiological thresholds during increasing water deficit in mature leaves, we identified key regulatory genes and main TF families that play a role in the drought response of tomato leaves. We highlighted the activation of ABA-independent TFs during the lowest WD intensities (at 4d-WD with partial stomatal closure, i.e. 10% reduction in stomatal conductance ( gs )), followed by a high transcription of ABA-dependent TFs one day later (5d-WD, full stomatal closure, i.e. 90% reduction in gs ). At stronger WD intensities after 8 days, when embolism started to develop in the leaves (4-6% PEP), we observed an activation of HEAT SHOCK PROTEIN genes and a general antioxidant response. Finally, at the strongest intensities after 10 days of WD, when embolism spread is about to exponentially increase (P 12 ), TFs related to oxidative stress were highly upregulated. When transcriptome profiles of young versus mature leaves are compared after 10d-WD in the same tomato plants, we can generalize that the stronger resilience of younger leaves to WD results from a shift in resource allocation from cell proliferation and expansion toward homeostasis maintenance. Compared to mature leaves, young leaves showed increased expression of genes altering the cell wall structure, like XTHs, cellulose synthase, and fasciclin-like arabinogalactan genes, which possibly create smaller cells with more rigid cell walls that can sustain stronger turgor pressures. Moreover, a potential preferential callose deposition at the phloem sieve elements could minimize water, sugar, and nutrient loss from these young leaves, allowing them to maintain osmoprotectant capacity and energy storage. Finally, by downregulating chlorophyll a/b binding protein genes, young leaves might reduce photosynthesis-related ROS production. This combined ecophysiological-molecular dataset in tomato leaves illustrates how gene regulatory pathways and ecophysiological thresholds are tightly intertwined. We anticipate that these outcomes will serve future studies aiming to integrate both disciplines more thoroughly, which will undoubtedly shed new light on how natural and agricultural selection has impacted the underlying mechanisms of drought response across the plant kingdom. Supplementary information Supplementary Figure S1 . Visual change of plant morphology over increasing water deficit. Supplementary Figure S2 . A compound leaf from the 6-8th node, showing leaflets selected for gas exchange, stomatal size, RNA-Seq, leaf PEP, Ψ TLP, and water potential analyses. Supplementary Figure S3 : Clustering accuracy as viewed on the first three principal components. Supplementary Figure S4 . Scale-free topology model fit and mean connectivity for different soft threshold powers. Supplemental Figure S5: Venn diagram of up- and downregulated genes per time point. Supplementary Figure S6: Transcription factors co-expression network. Supplementary Figure S7: Gene ontology enrichment. Supplementary Table S1: Variance Stabilizing Transformation (VST) Counts. Supplementary Table S2: Differentially Expressed (DE) genes in mature leaves. Supplementary Table S3: Differentially Expressed (DE) genes in young leaves at 10 d-WD. Author contributions G.B., T.d.W., B.M.R., S.B., and F.L. designed the experiments. S.B., F.L., and B.M.-R. initiated the research and provided funding. G.B., M.L., and A.T. carried out the physiological measurements. T.d.W. carried out the differential expression and clustering analyses. G.B. and T.d.W. examined genes of interest, identified structural homologs in Arabidopsis, and performed accompanying literature research. G.B. and T.d.W. wrote the manuscript by accepting amendments from all authors. All authors agreed on the final version of the manuscript. Conflict of interest The authors declare that they do not have a conflict of interest. Data availability Raw aligned and mapped counts from all samples, as well as the raw sequencing files are available at https://www.ncbi.nlm.nih.gov/geo/query/acc.cgi?acc=GSE270245 . Differential expression results, as well as the VST-stabilized counts as used for the WGCNA network generation are provided in the Supplementary Data. Figures alt text Figure 1A: Line graph showing leaf water potential variation over time, for both well-watered (green) and water-deficit subjected (yellow) tomato plants during a 10-day drought experiment. Figure 1B: Graph showing stomatal safety margin during leaf water potential decline. The blue curve (right) shows gas exchange, the red curve (left) shows the percentage of embolized pixels. A purple vertical line marks the turgor loss point. Figure 2A: Multidimensional scaling plot showing dissimilarity between water deficit (yellow symbols) and well-watered (green symbols) tomato plant samples. Shapes represent sampling points. Figure 2B: Heatmap and hierarchical clustering of the top 500 most variable genes in water deficit and well-watered tomato plants. Figure 2C: Volcano plots showing differentially expressed genes between well-watered and water deficit conditions in mature leaves at each time point. Differential expression is defined as log2 fold change >1 and FDR-adjusted p-value <0.05. Figure 3A: Dendrogram showing WGCNA results for differentially expressed genes in tomato plants under water deficit. 19 distinct co-expression clusters are identified with different colors and Roman numerals. Figure 3B: Heatmap showing cluster-trait relationships between module eigengenes and physiological traits: leaf water potential (Ψleaf), stomatal conductance (gs), CO2 assimilation (A), and leaf PEP. Figure 3C: Line plot showing expression Z-scores for each WGCNA cluster over increasing water deficit from 0 to 10 days. Expression values are scaled and centered on VST normalized counts. Figure 4: Heatmap showing transcription factor expression and number of correlated genes over progressive water deficit in tomato plants. Transcription factors are separated by the time point they first show differential expression. FIgure 5A: Heatmap showing z-scores of differentially regulated genes between young and mature leaves under well-watered and water deficit conditions at 10 days. Figure 5B: Line plots showing average gene expression (± SE) in each K-means cluster for young and mature leaves under well-watered and water deficit conditions. Supplementary Figure S1: Series of images (A-E) showing visual changes in tomato plant morphology over increasing water deficit from well-watered to 10 days after water deficit. Image F shows a sampled mature compound leaf and three younger leaves from the same plant at 10 days. Supplementary Figure S2: Image of a compound leaf from the 6-8th node, highlighting leaflets selected for various analyses. Supplementary Figure S3: Scatter plot showing clustering accuracy on the first three principal components. Supplemental Figure S4: Two-panel graph showing scale-free topology model fit and mean connectivity for different soft threshold powers. The left panel shows model fit (y-axis) against soft threshold powers (x-axis), and the right panel shows mean connectivity (y-axis) across the same range of powers. Supplemental Figure S5: Venn diagram showing numbers of up- and downregulated genes at each time point and their intersections. Numbers in sections represent gene counts, with percentages indicating the fraction of total differentially expressed genes. Supplemental Figure S6: Network diagram showing transcription factor co-expression network in tomato plants under well-watered and water deficit conditions. Supplemental Figure S7: Dot plot showing significantly enriched gene ontology terms on biological processes within four clusters. Acknowledgments We would like to thank Alex Bos (IBL), Jan Vink (IBL), Bertie-Joan van Hueven (Naturalis), Rob Langelaan (Naturalis), and Régis Burlett (platform Phenobois Bordeaux) for providing technical help. We want to thank Dirk Walther (MPI-MP) for verifying the bioinformatic approaches. B.M.-R. and T.d.W. thank the International Max Planck Research School ‘Primary Metabolism and Plant Growth’ (IMPRS-PMPG; fellowship to T.d.W.) for support. F.L., M.L., and G.B. thank the Dutch Research Council for funding (grant ALWOP.488). Footnotes ↵ * Joint first authors GB: giovanni.bortolami{at}epfl.ch TdW: tobias.dewerk{at}gmail.com ML: maximilien.larter{at}gmail.com AT: ajaree.thonglim{at}gmail.com BMR: bmr{at}uni-potsdam.de SB: s.balazadeh{at}biology.leidenuniv.nl FL: frederic.lens{at}naturalis.nl References 1. Abdalla M , Ahmed MA , Cai G , Wankmüller F , Schwartz N , Litig O , Javaux M , Carminati A . 2022 . Stomatal closure during water deficit is controlled by below-ground hydraulics . Annals of Botany 129 , 161 – 170 . OpenUrl 2. ↵ Ahammed GJ , Li X , Wan H , Zhou G , Cheng Y . 2020a . SlWRKY81 reduces drought tolerance by attenuating proline biosynthesis in tomato . Scientia Horticulturae 270 , 109444 . OpenUrl 3. ↵ Ahammed GJ , Li X , Yang Y , Liu C , Zhou G , Wan H , Cheng Y . 2020b . Tomato WRKY81 acts as a negative regulator for drought tolerance by modulating guard cell H2O2–mediated stomatal closure . Environmental and Experimental Botany 171 , 103960 . OpenUrl 4. Altschul SF , Gish W , Miller W , Myers EW , Lipman DJ . 1990 . Basic Local Alignment Search Tool . Journal of Molecular Biology 215 , 403 – 410 . OpenUrl CrossRef PubMed Web of Science 5. Asada K . 2006 . Production and Scavenging of Reactive Oxygen Species in Chloroplasts and Their Functions . Plant Physiology 141 , 391 – 396 . OpenUrl FREE Full Text 6. ↵ Bailey-Serres J , Parker JE , Ainsworth EA , Oldroyd GED , Schroeder JI . 2019 . Genetic strategies for improving crop yields . Nature 575 , 109 – 118 . OpenUrl CrossRef PubMed 7. Barbosa EGG , Leite JP , Marin SRR , Marinho JP , De Fátima Corrêa Carvalho J , Fuganti-Pagliarini R , Farias JRB , Neumaier N , Marcelino-Guimarães FC , De Oliveira MCN , Yamaguchi-Shinozaki K , Nakashima K , Maruyama K , Kanamori N , Fujita Y , Yoshida T , Nepomuceno AL . 2013 . Overexpression of the ABA-Dependent AREB1 Transcription Factor from Arabidopsis thaliana Improves Soybean Tolerance to Water Deficit . Plant Molecular Biology Report 31 , 719 – 730 . OpenUrl 8. ↵ Berens ML , Wolinska KW , Spaepen S , Ziegler J , Nobori T , Nair A , Krüler V , Winkelmüller TM , Wang Y , Mine A , Becker D , Garrido-Oter R , Schulze-Lefert P , Tsuda K . 2019 . Balancing trade-offs between biotic and abiotic stress responses through leaf age-dependent variation in stress hormone cross-talk . Proceedings of the National Academy of Sciences 116 , 2364 – 2373 . OpenUrl Abstract / FREE Full Text 9. ↵ Bian Z , Wang Y , Zhang X , Grundy S , Hardy K , Yang Q , Lu C . 2021 . A Transcriptome Analysis Revealing the New Insight of Green Light on Tomato Plant Growth and Drought Stress Tolerance . Frontiers of Plant Sciences 12 , 649283 . OpenUrl 10. ↵ Brodribb TJ , Skelton RP , McAdam SAM , Bienaimé D , Lucani CJ , Marmottant P . 2016 . Visual quantification of embolism reveals leaf vulnerability to hydraulic failure . New Phytologist 209 , 1403 – 1409 . OpenUrl CrossRef PubMed 11. ↵ Cackett L , Luginbuehl LH , Schreier TB , Lopez-Juez E , Hibberd JM . 2022 . Chloroplast development in green plant tissues: the interplay between light, hormone, and transcriptional regulation . New Phytologist 233 , 2000 – 2016 . OpenUrl CrossRef 12. ↵ Chase JM , Leibold MA . 2003 . Ecological Niches: Linking Classical and Contemporary Approaches . University of Chicago Press . 13. ↵ Chen W , Hu T , Ye J , Wang B , Liu G , Wang Y , Yuan L , Li J , Li F , Ye Z , Zhang Y. 2020 . A CCAAT-binding factor, SlNFYA10, negatively regulates ascorbate accumulation by modulating the d-mannose/l-galactose pathway in tomato . Horticultural Research 7 , 200 . OpenUrl 14. Cheng C , Krishnakumar V , Chan AP , Thibaud-Nissen F , Schobel S , Town CD . 2017 . Araport11: a complete reannotation of the Arabidopsis thaliana reference genome . The Plant Journal 89 , 789 – 804 . OpenUrl CrossRef PubMed 15. ↵ Cheng L , Wang Y , He Q , Li H , Zhang X , Zhang F . 2016 . Comparative proteomics illustrates the complexity of drought resistance mechanisms in two wheat ( Triticum aestivum L.) cultivars under dehydration and rehydration . BMC Plant Biology , 16 , 188 . OpenUrl 16. ↵ Cho SK , Kim JE , Park J-A , Eom TJ , Kim WT . 2006 . Constitutive expression of abiotic stress- inducible hot pepper CaXTH3, which encodes a xyloglucan endotransglucosylase/hydrolase homolog, improves drought and salt tolerance in transgenic Arabidopsis plants . FEBS Letters , 580 , 3136 – 3144 . OpenUrl CrossRef PubMed Web of Science 17. ↵ Choat B , Brodribb TJ , Brodersen CR , Duursma RA , López R , Medlyn B.E . 2018 . Triggers of tree mortality under drought . Nature 558 , 531 – 539 . OpenUrl CrossRef PubMed 18. ↵ Corso D , Delzon S , Lamarque LJ , Cochard H , Torres-Ruiz JM , King A , Brodribb TJ . 2020 . Neither xylem collapse, cavitation, or changing leaf conductance drive stomatal closure in wheat . Plant Cell and Environment 43 : 854 – 865 . OpenUrl 19. ↵ Da Silva A , Oliveira Silva C , Do Rosario Rosa V , Silva Santos MF , Naomi Kuki K , Dal- Bianco M , Delmond Bueno R , Alves De Oliveira J , Santos Brito D , Costa AC , Ribeiro C . 2022 . Metabolic adjustment and regulation of gene expression are essential for increased resistance to severe water deficit and resilience post-stress in soybean . PeerJ 10 , e13118 . OpenUrl 20. ↵ Dietrich K , Weltmeier F , Ehlert A , Weiste C , Stahl M , Harter K , Dröge-Laser W . 2011 . Heterodimers of the Arabidopsis transcription factors bZIP1 and bZIP53 reprogram amino acid metabolism during low energy stress . Plant Cell 23 , 381 – 95 . OpenUrl Abstract / FREE Full Text 21. ↵ Diouf I , Albert E , Duboscq R , Santoni S , Bitton F , Gricourt J , Causse M . 2020 . Integration of QTL, Transcriptome and Polymorphism Studies Reveals Candidate Genes for Water Stress Response in Tomato . Genes 11 , 900 . OpenUrl 22. ↵ D’Ippólito S , Rey-Burusco MF , Feingold SE , Guevara MG . 2021 . Role of proteases in the response of plants to drought . Plant Physiology and Biochemistry 168 , 1 – 9 . OpenUrl 23. ↵ Dixon HH , Joly J . 1895 . On the Ascent of Sap . Philosophical Transactions of the Royal Society of London . 186 , 563 – 576 . OpenUrl 24. ↵ Dong S , Ling J , Song L , Zhao L , Wang Y , Zhao T . 2023 . Transcriptomic Profiling of Tomato Leaves Identifies Novel Transcription Factors Responding to Dehydration Stress . IJMS 24 , 9725 . OpenUrl 25. ↵ Fujita Y , Fujita M , Satoh R , Maruyama K , Parvez MM , Seki M , Hiratsu K , Ohme-Takagi M , Shinozaki K , Yamaguchi-Shinozaki K . 2005 . AREB1 Is a Transcription Activator of Novel ABRE-Dependent ABA Signaling That Enhances Drought Stress Tolerance in Arabidopsis . The Plant Cell 17 , 3470 – 3488 . OpenUrl Abstract / FREE Full Text 26. Ge SX , Jung D , Yao R . 2020 . ShinyGO: a graphical gene-set enrichment tool for animals and plants . Bioinformatics 36 , 2628 – 2629 . OpenUrl CrossRef PubMed 27. ↵ Gill SS , Tuteja N . 2010 . Reactive oxygen species and antioxidant machinery in abiotic stress tolerance in crop plants . Plant Physiology and Biochemistry , 48 , 909 – 930 . OpenUrl CrossRef PubMed Web of Science 28. ↵ Gong P , Zhang J , Li H , Yang C , Zhang C , Zhang X , Khurram Z , Zhang Y , Wang T , Fei Z , Ye Z . 2010 . Transcriptional profiles of drought-responsive genes in modulating transcription signal transduction, and biochemical pathways in tomato . Journal of Experimental Botany 61 , 3563 – 3575 . OpenUrl CrossRef PubMed Web of Science 29. ↵ Greve P , Gudmundsson L , Seneviratne SI . 2018 . Regional scaling of annual mean precipitation and water availability with global temperature change . Earth System Dynamics 9 , 227 – 240 . OpenUrl 30. ↵ Greve P , Orlowsky B , Mueller B , Sheffield J , Reichstein M , Seneviratne SI . 2014 . Global assessment of trends in wetting and drying over land . Nature Geoscience 7 , 716 – 721 . OpenUrl 31. ↵ Greve P , Roderick ML , Seneviratne SI . 2017 . Simulated changes in aridity from the last glacial maximum to 4xCO . Environmental Research Letters 12 , 114021 . OpenUrl 32. ↵ Guerreiro SB , Dawson RJ , Kilsby C , Lewis E , Ford A . 2018 . Future heat-waves, droughts and floods in 571 European cities . Environmental Research Letters 13 , 034009 . OpenUrl 33. ↵ Harrison Day BL , Carins-Murphy MR , Brodribb TJ . 2022 . Reproductive water supply is prioritized during drought in tomato . Plant, Cell & Environment , 45 , 69 – 79 . OpenUrl CrossRef 34. ↵ Hayashi T , Kaida R . 2011 . Functions of Xyloglucan in Plant Cells . Molecular Plant 4 , 17 – 24 . OpenUrl CrossRef PubMed Web of Science 35. ↵ Haverroth EJ , Rimer IM , Oliveira LA , de Lima LGA , Cesarino I , Martins SCV , McAdam SAM , Cardoso AA . 2024 . Gradients in embolism resistance within stems driven by secondary growth in herbs . Plant, Cell, & Environment , 1 – 13 . 36. Huang S , Gao Y , Liu J , Peng X , Niu X , Fei Z , Cao S , Liu Y . 2012 . Genome-wide analysis of WRKY transcription factors in Solanum lycopersicum . Molecular Genetics and Genomics 287 : 495 – 513 . OpenUrl CrossRef PubMed 37. ↵ Huang Y , Niu C , Yang C , Jinn T . 2016 . The Heat Stress Factor HSFA6b Connects ABA Signaling and ABA-Mediated Heat Responses . Plant Physiology 172 , 1182 – 1199 . OpenUrl Abstract / FREE Full Text 38. ↵ Huber AE , Melcher PJ , Piñeros MA , Setter TL , Bauerle TL . 2019 . Signal coordination before, during and after stomatal closure in response to drought stress . New Phytologist 224 , 675 – 688 . OpenUrl 39. ↵ Hummel I , Pantin F , Sulpice R , Piques M , Rolland G , Dauzat M , Christophe A , Pervent M , Bouteillé M , Stitt M , Gibon Y , Muller B . 2010 . Arabidopsis Plants Acclimate to Water Deficit at Low Cost through Changes of Carbon Usage: An Integrated Perspective Using Growth, Metabolite, Enzyme, and Gene Expression Analysis . Plant Physiology 154 , 357 – 372 . OpenUrl Abstract / FREE Full Text 40. ↵ Hwang I , Sheen J , Müller B . 2012 . Cytokinin signaling networks . Annual Review of Plant Biology 63 , 353 – 80 . OpenUrl CrossRef PubMed Web of Science 41. ↵ Iovieno P , Punzo P , Guida G , Mistretta C , Van Oosten MJ , Nurcato R , Bostan H , Colantuono C , Costa A , Bagnaresi P , Chiusano ML , Albrizio R , Giorio P , Batelli G , Grillo S . 2016 . Transcriptomic Changes Drive Physiological Responses to Progressive Drought Stress and Rehydration in Tomato . Frontiers of Plant Sciences 7 . 42. ↵ Islam MS , Wang M-H . 2009 . Expression of dehydration responsive element-binding protein- 3 (DREB3) under different abiotic stresses in tomato . BMB Reports 42 , 611 – 616 . OpenUrl 43. ↵ Iturriaga G , Suárez R , Nova-Franco B . 2009 . Trehalose Metabolism: From Osmoprotection to Signaling . International Journal of Molucular Sciences 10 , 3793 – 3810 . OpenUrl 44. ↵ Jedmowski C , Ashoub A , Beckhaus T , Berberich T , Karas M , Brüggemann W . 2014 . Comparative Analysis of Sorghum bicolor Proteome in Response to Drought Stress and following Recovery . International Journal of Proteomics 2014 , 395905 . OpenUrl 45. ↵ Jia C , Wang J , Guo B , Li X , Yang T , Yang H , Li N , Wang B , Yu Q . 2023 . A group III WRKY transcription factor, SlWRKY52, positively regulates drought tolerance in tomato . Environmental and Experimental Botany 215 , 105513 . OpenUrl 46. Jin JF , Wang ZQ , He QY , Wang JY , Li PF , Xu JM , Yheng SJ , Fan W , Yang JL . 2020 . Genome-wide identification and expression analysis of the NAC transcription factor family in tomato ( Solanum lycopersicum ) during aluminum stress . BMC Genomics 21 , 288 . OpenUrl 47. ↵ Kaku T , Tabuchi A , Wakabayashi K , Hoson T . 2004 . Xyloglucan Oligosaccharides Cause Cell Wall Loosening by Enhancing Xyloglucan Endotransglucosylase/Hydrolase Activity in Azuki Bean Epicotyls . Plant and Cell Physiology 45 , 77 – 82 . OpenUrl CrossRef PubMed Web of Science 48. ↵ M. A. Hossain , V. Kumar , D. J. Burritt , M. Fujita P. S. A. Mäkelä Kido ÉA , Ferreira-Neto JRC , da Silva MD , Santos VEP , da Silva Filho JLB , Benko-Iseppon A M . 2019 . Osmoprotectant-Related Genes in Plants Under Abiotic Stress: Expression Dynamics, In Silico Genome Mapping, and Biotechnology . In M. A. Hossain , V. Kumar , D. J. Burritt , M. Fujita , & P. S. A. Mäkelä (Eds.), Osmoprotectant-Mediated Abiotic Stress Tolerance in Plants: Recent Advances and Future Perspectives (pp. 1 – 40 ). Springer International Publishing . 49. ↵ Lamarque LJ , Delmas CEL , Charrier G , Burlett R , Dell’Acqua N , Pouzoulet J , Gambetta GA , Delzon S . 2023 . Quantifying the grapevine xylem embolism resistance spectrum to identify varieties and regions at risk in a future dry climate . Scientific Report 13 , 7724 . OpenUrl 50. ↵ Lamarque LJ , Delzon S , Toups H , Gravel A , Corso D , Badel E , Burlett R , Charrier G , Cochard H , Jansen S , King A , Torres-Ruiz JM , Pouzoulet J , Cramer GR , Thompson AJ , Gambetta GA . 2020 . Over-accumulation of abscisic acid in transgenic tomato plants increases the risk of hydraulic failure . Plant Cell Environment 43 , 548 – 562 . OpenUrl 51. ↵ Lamesch P , Berardini TZ , Li D , Swarbreck D , Wilks C , Sasidharan R , Muller R , Dreher K , Alexander DL , Garcia-Hernandez M , Karthikeyan AS , Lee CH , Nelson WD , Ploetz L , Singh S , Wensel A , Huala E . 2012 . The Arabidopsis Information Resource (TAIR): improved gene annotation and new tools . Nucleic Acids Research 40 , D1202 – D1210 . OpenUrl CrossRef PubMed Web of Science 52. ↵ Langfelder P , Horvath S . 2008 . WGCNA: an R package for weighted correlation network analysis . BMC Bioinformatics 9 , 559 OpenUrl CrossRef PubMed 53. Lens F , Gleason SM , Bortolami G , Brodersen C , Delzon S , Jansen S . 2022 . Functional xylem characteristics associated with drought-induced embolism in angiosperms . New Phytologist 236 , 2019 – 2036 . OpenUrl 54. ↵ Lens F , Sperry JS , Christman MA , Choat B , Rabaey D , Jansen S . 2011 . Testing hypotheses that link wood anatomy to cavitation resistance and hydraulic conductivity in the genus Acer . New Phytologist 190 , 709 – 723 . OpenUrl CrossRef PubMed Web of Science 55. ↵ Li S , Li K , Ju Z , Cao D , Fu D , Zhu H , Zhu B , Luo Y . 2016 . Genome-wide analysis of tomato NF-Y factors and their role in fruit ripening . BMC Genomics 17 , 36 . OpenUrl 56. ↵ Lionello P , Scarascia L . 2020 . The relation of climate extremes with global warming in the Mediterranean region and its north versus south contrast . Regional Environmental Change 20 , 31 . OpenUrl 57. ↵ Love MI , Huber W , Anders S . 2014 . Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2 . Genome Biology 15 , 550 . OpenUrl CrossRef PubMed 58. ↵ Mantova M , Herbette S , Cochard H , Torres-Ruiz JM . 2022 . Hydraulic failure and tree mortality: from correlation to causation . Trends Plant in Science 27 , 335 – 345 . OpenUrl 59. Martin-StPaul N , Delzon S , Cochard H . 2017 . Plant resistance to drought depends on timely stomatal closure . Ecology Letters 20 , 1437 – 1447 . OpenUrl CrossRef PubMed 60. ↵ Malefo MB , Mathibela EO , Crampton BG , Makgopa ME . 2020 . Investigating the role of Bowman-Birk serine protease inhibitor in Arabidopsis plants under drought stress . Plant Physiology and Biochemistry , 149 , 286 – 293 . OpenUrl 61. ↵ Miedes E , Suslov D , Vandenbussche F , Kenobi K , Ivakov A , Van Der Straeten D , Lorences EP , Mellerowicz EJ , Verbelen J-P , Vissenberg K . 2013 . Xyloglucan endotransglucosylase/hydrolase (XTH) overexpression affects growth and cell wall mechanics in etiolated Arabidopsis hypocotyls . Journal of Experimental Botany 64 , 2481 – 2497 . OpenUrl CrossRef PubMed Web of Science 62. Mishra SK , Tripp J , Winkelhaus S , Tschiersch B , Theres K , Nover L , Scharf KD . 2002 . In the complex family of heat stress transcription factors, HsfA1 has a unique role as master regulator of thermotolerance in tomato . Genes Development 16 , 1555 – 67 . OpenUrl Abstract / FREE Full Text 63. ↵ Mishra U , Rai A , Kumar R , Singh M , Pandey HP . 2016 . Gene expression analysis of Solanum lycopersicum and Solanum habrochaites under drought conditions . Genomics Data 9 , 40 – 41 . OpenUrl 64. ↵ Mishra N , Jiang C , Chen L , Paul A , Chatterjee A , Shen G . 2023 . Achieving abiotic stress tolerance in plants through antioxidative defense mechanisms . Frontiers in Plant Science , 14 . 65. ↵ Moghaddam GA , Rezayatmand Z , Nasr Esfahani M , Khozaei M . 2019 . Genetic defense analysis of tomatoes in response to early blight disease, Alternaria alternata . Plant Physiology and Biochemistry 142 , 500 – 509 . OpenUrl 66. ↵ Moloi SJ , Ngara R . 2023 . The roles of plant proteases and protease inhibitors in drought response: A review . Frontiers in Plant Science , 14 . 67. ↵ Müller M , Munné-Bosch S . 2015 . Ethylene Response Factors: A Key Regulatory Hub in Hormone and Stress Signaling . Plant Physiology 169 , 32 – 41 . OpenUrl Abstract / FREE Full Text 68. Ozturk M , Turkyilmaz Unal B , García-Caparrós P , Khursheed A , Gul A , Hasanuzzaman M . 2021 . Osmoregulation and its actions during the drought stress in plants . Physiologia Plantarum 172 , 1321 – 1335 . OpenUrl CrossRef 69. ↵ Pan X , Wang C , Liu Z , Gao R , Feng L , Li A , Yao K , Liao W . 2023 . Identification of ABF/AREB gene family in tomato ( Solanum lycopersicum L.) and functional analysis of ABF/AREB in response to ABA and abiotic stresses . PeerJ 11 , e15310 . OpenUrl 70. ↵ Pammenter NW , Van Der Willigen C . 1998 . A mathematical and statistical analysis of the curves illustrating vulnerability of xylem to cavitation . Tree Physiology 18 , 589 – 593 . OpenUrl PubMed Web of Science 71. ↵ Petroni K , Kumimoto RW , Gnesutta N , Calvenzani V , Fornari M , Tonelli C , Holt BF 3rd, Mantovani R . 2012 . The promiscuous life of plant NUCLEAR FACTOR Y transcription factors . Plant Cell 24 , 4777 – 92 . OpenUrl Abstract / FREE Full Text 72. ↵ Pingali PL . 2012 . Green Revolution: Impacts, limits, and the path ahead . Proceedings of the National Academy of Sciences of the United States of America , 109 , 12302 . OpenUrl Abstract / FREE Full Text 73. ↵ Pirona R , Frugis G , Locatelli F , Mattana M , Genga A , Baldoni E . 2023 . Transcriptomic analysis reveals the gene regulatory networks involved in leaf and root response to osmotic stress in tomato . Frontiers of Plant Sciences 14 , 1155797 . OpenUrl 74. ↵ Ray DK , Mueller ND , West , PC , Foley JA . 2013 . Yield Trends Are Insufficient to Double Global Crop Production by 2050 . PLoS ONE , 8 , e66428 . OpenUrl CrossRef PubMed 75. ↵ Rodriguez-Dominguez CM , Forner A , Martorell S , Choat B , Lopez R , Peters JMR. et al. , 2022 . Leaf water potential measurements using the pressure chamber: Synthetic testing of assumptions towards best practices for precision and accuracy . Plant, Cell & Environment 45 , 2037 – 2061 . OpenUrl 76. ↵ Rosa M , Prado C , Podazza G , Interdonato R , González JA , Hilal M , Prado FE . 2009 . Soluble sugars: Metabolism, sensing and abiotic stress: A complex network in the life of plants . Plant Signaling & Behavior 4 , 388 – 393 . OpenUrl 77. ↵ Procházková D , Wilhelmová N . . 2007 . Leaf senescence and activities of the antioxidant enzymes . Biologia plantarum 51 , 401 – 406 . OpenUrl 78. ↵ Raymond F , Ullmann A , Camberlin P , Oueslati B , Drobinski P . 2018 . Atmospheric conditions and weather regimes associated with extreme winter dry spells over the Mediterranean basin . Climate Dynamics 50 , 4437 – 4453 . OpenUrl 79. ↵ Sack L , Pasquet-Kok J , Bartlett M . 2010 . Leaf pressure-volume curve parameters . PROMETHEUS . 80. ↵ Sami F , Yusuf M , Faizan M , Faraz A , Hayat S . 2016 . Role of sugars under abiotic stress . Plant Physiology and Biochemistry 109 , 54 – 61 . OpenUrl CrossRef 81. ↵ Schippers JHM , Schmidt R , Wagstaff C , Jing H-C . 2015 . Living to Die and Dying to Live: The Survival Strategy behind Leaf Senescence . Plant Physiology 169 , 914 – 930 . OpenUrl Abstract / FREE Full Text 82. Schneider CA , Rasband , WS , Eliceiri KW . 2012 . NIH Image to ImageJ: 25 years of image analysis . Nature Methods 9 , 671 – 675 . OpenUrl 83. ↵ Singh PK , Indoliya Y , Agrawal L , Awasthi S , Deeba F , Dwivedi S , Chakrabarty D , Shirke PA , Pandey V , Singh N , Dhankher OP , Barik SK , Tripathi RD . 2022 . Genomic and proteomic responses to drought stress and biotechnological interventions for enhanced drought tolerance in plants . Current Plant Biology 29 , 100239 . OpenUrl 84. Solgenomics Network . 2020 . ITAG4.1_descriptions.txt [Dataset]. Last consulted May 1 st 2023 , from https://solgenomics.net/ftp/tomato_genome/annotation/ITAG4.1_release/ 85. ↵ Sorek Y , Greenstein S , Hochberg U . 2022 . Seasonal adjustment of leaf embolism resistance and its importance for hydraulic safety in deciduous trees . Physiologia Plantarum 174 , e13785 . OpenUrl 86. ↵ Sorek Y , Greenstein S , Netzer Y , Shtein I , Jansen S , Hochberg , U . 2021 . An increase in xylem embolism resistance of grapevine leaves during the growing season is coordinated with stomatal regulation, turgor loss point and intervessel pit membranes . New Phytologist 229 : 1955 – 1969 . OpenUrl CrossRef 87. Szabados L , Savouré A . 2010 . Proline: a multifunctional amino acid . Trends in Plant Science 15 , 89 – 97 . OpenUrl CrossRef PubMed Web of Science 88. ↵ Taiz L , Møller IM , Murphy A , Zeiger E . 2024 . Plant Physiology and Development, 7th edition . Oxford University Press . 89. ↵ Thonglim A , Bortolami G , Delzon S , Larter M , Offringa R , Keurentjes JJB , Smets E , Balazadeh S , Lens F . 2023 . Drought response in Arabidopsis displays synergistic coordination between stems and leaves . Journal of Experimental Botany 74 , 1004 – 1021 . OpenUrl 90. ↵ Thonglim A , Delzon S , Larter M , Karami O , Rahimi A , Offringa R , Keurentjes JJN , Balazadeh S , Smets E , Lens F . 2021 . Intervessel pit membrane thickness best explains variation in embolism resistance amongst stems of Arabidopsis thaliana accessions , Annals of Botany 128 , 171 – 182 . OpenUrl 91. ↵ Tombesi S , Nardini A , Frioni T , Soccolini M , Zadra C , Farinelli D , Poni S , Palliotti A . 2015 . Stomatal closure is induced by hydraulic signals and maintained by ABA in drought-stressed grapevine . Scientific Reports 5 , 12449 . OpenUrl 92. ↵ Treurnicht M , Pagel J , Tonnabel J , Esler KJ , Slingsby JA , Schurr FM . 2020 . Functional traits explain the Hutchinsonian niches of plant species . Global Ecology and Biogeography , 29 , 534 – 545 . OpenUrl 93. ↵ Tzfadia O , Diels T , De Meyer S , Vandepoele K , Aharoni A , Van De Peer Y . 2016 . CoExpNetViz: Comparative Co-Expression Networks Construction and Visualization Tool . Frontiers of Plant Sciences 6 . 94. ↵ Vatén A , Dettmer J , Wu S , Stierhof Y-D , Miyashima S , Yadav SR , Roberts CJ , Campilho A , Bulone V , Lichtenberger R , Lehesranta S , Mähönen AP , Kim J-Y , Jokitalo E , Sauer N , Scheres B , Nakajima K , Carlsbecker A , Gallagher KL , Helariutta Y . 2011 . Callose Biosynthesis Regulates Symplastic Trafficking during Root Development . Developmental Cell 21 , 1144 – 1155 . OpenUrl CrossRef PubMed Web of Science 95. ↵ Waadt R , Seller CA , Hsu P-K , Takahashi Y , Munemasa S , Schroeder JI . 2022 . Plant hormone regulation of abiotic stress responses . Nature Reviews Molecular Cell Biology 23 , 680 – 694 . OpenUrl 96. ↵ Wang R-H , Chang J-C , Li K-T , Lin T-S , Chang L-S . 2014 . Leaf age and light intensity affect gas exchange parameters and photosynthesis within the developing canopy of field net- house-grown papaya trees . Scientia Horticulturae 165 , 365 – 373 . OpenUrl 97. Wang X , Gao L , Jiao C , Stravoravdis S , Hosmani PS , Saha S , Zhang J , Mainiero S , Strickler SR , Catala C , Martin GB , Mueller LA , Vrebalov J , Giovannoni JJ , Wu S , Fei Z . 2020 . Genome of Solanum pimpinellifolium provides insights into structural variants during tomato breeding . Nature Communications 11 , 5817 . OpenUrl 98. ↵ Wartenburger R , Hirschi M , Donat MG , Greve P , Pitman AJ , Seneviratne SI . 2017 . Changes in regional climate extremes as a function of global mean temperature: an interactive plotting framework . Geoscientific Model Development 10 , 3609 – 3634 . OpenUrl 99. ↵ Watkins JM , Chapman JM , Muday GK . 2017 . Abscisic Acid-Induced Reactive Oxygen Species Are Modulated by Flavonols to Control Stomata Aperture . Plant Physiology 175 , 1807 – 1825 . OpenUrl Abstract / FREE Full Text 100. ↵ Wu S-W , Kumar R , Iswanto ABB , Kim J-Y . 2018 . Callose balancing at plasmodesmata . Journal of Experimental Botany 69 , 5325 – 5339 . OpenUrl CrossRef 101. ↵ Yanagisawa S . 2002 . The Dof family of plant transcription factors . Trends in Plant Science 7 , 555 – 560 . OpenUrl CrossRef PubMed Web of Science 102. Yao T , Zhang J , Xie M , Yuan G , Tschaplinski TJ , Muchero W , Chen J-G . 2021 . Transcriptional Regulation of Drought Response in Arabidopsis and Woody Plants . Frontiers of Plant Sciences 11 , 572137 . OpenUrl 103. ↵ Yuan Y , Ren S , Liu X , Su L , Wu Y , Zhang W , Li Y , Jiang Y , Wang H , Fu R , Bouzayen M , Liu M , Zhang Y . 2022 . SlWRKY35 positively regulates carotenoid biosynthesis by activating the MEP pathway in tomato fruit . New Phytologist 234 , 164 – 178 . OpenUrl 104. ↵ Zhang H , Zhu J , Gong Z , Zhu J-K . 2022 . Abiotic stress responses in plants . Nature Reviews Genetics 23 , 104 – 119 . OpenUrl CrossRef 105. ↵ Zhang L , Merlin I , Pascal S , Bert P , Domergue F , Gambetta GA . 2020 . Drought activates MYB41 orthologs and induces suberization of grapevine fine roots . Plant Direct 4 , e00278 . OpenUrl 106. ↵ Zhou R , Yu X , Zhao T , Ottosen C-O , Rosenqvist E , Wu Z . 2019 . Physiological analysis and transcriptome sequencing reveal the effects of combined cold and drought on tomato leaf . BMC Plant Biology 19 , 377 . OpenUrl 107. ↵ Zhu A , Ibrahim JG , Love MI . 2019 . Heavy-tailed prior distributions for sequence count data: removing the noise and preserving large differences . Bioinformatics 35 , 2084 – 2092 . OpenUrl CrossRef PubMed View the discussion thread. Back to top Previous Next Posted July 05, 2024. Download PDF Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Integrating Gene Expression Analysis and Ecophysiological Responses to Water Deficit in Leaves of Tomato Plants Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Integrating Gene Expression Analysis and Ecophysiological Responses to Water Deficit in Leaves of Tomato Plants G Bortolami , T de Werk , M Larter , A Thonglim , B Mueller-Roeber , S. Balazadeh , F. Lens bioRxiv 2024.07.05.602262; doi: https://doi.org/10.1101/2024.07.05.602262 Share This Article: Copy Citation Tools Integrating Gene Expression Analysis and Ecophysiological Responses to Water Deficit in Leaves of Tomato Plants G Bortolami , T de Werk , M Larter , A Thonglim , B Mueller-Roeber , S. Balazadeh , F. Lens bioRxiv 2024.07.05.602262; doi: https://doi.org/10.1101/2024.07.05.602262 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Plant Biology Subject Areas All Articles Animal Behavior and Cognition (7653) Biochemistry (17763) Bioengineering (13944) Bioinformatics (42101) Biophysics (21509) Cancer Biology (18667) Cell Biology (25592) Clinical Trials (138) Developmental Biology (13413) Ecology (19969) Epidemiology (2067) Evolutionary Biology (24393) Genetics (15647) Genomics (22582) Immunology (17791) Microbiology (40524) Molecular Biology (17222) Neuroscience (88860) Paleontology (667) Pathology (2848) Pharmacology and Toxicology (4841) Physiology (7670) Plant Biology (15182) Scientific Communication and Education (2048) Synthetic Biology (4312) Systems Biology (9843) Zoology (2274)