Stem Growth and Dehydration Responses of Mediterranean Tree Species to Atmospheric and Soil Drought.

Stem growth responses to soil and atmospheric drought are critical to forecasting the tree carbon sink strength. Yet, responses of drought-prone forests remain uncertain despite global aridification trends. Stem diameter variations at an hourly resolution were monitored in five Mediterranean tree species from a mesic and a xeric site for 6 and 12 years. Stem growth and dehydration responses to soil (REW) and atmospheric (VPD) drought were explored at different timescales. Annually, growth was determined by the number of growing days and hours. Seasonally, growth was bimodal (autumn growth ≈ 8%-18% of annual growth), varying among species and sites across the hydrometeorological space, while dehydration consistently responded to REW. Sub-daily, substantial growth occurred during daytime, with nighttime-to-daytime ratios ranging between 1.2 and 3.5 (Arbutus unedo ≈ Quercus faginea < Quercus ilex < Pinus halepensis in the mesic site, and Juniperus thurifera < P. halepensis in the xeric site). Overall, time windows favourable for growth were limited by soil (rather than atmospheric) drought, modulating annual and seasonal growth in Mediterranean species, and stems maintained non-negligible growth during daytime. These patterns contrast with observations from wetter or cooler biomes, demonstrating the growth plasticity of drought-prone species to more arid climate conditions.

Transcription factor PdMYB118 in poplar regulates lignin deposition and xylem differentiation in addition to anthocyanin synthesis through suppressing the expression of PagKNAT2/6b gene.

R2R3-MYB transcription factors function as the master regulators of the phenylpropanoid pathway in which both lignin and anthocyanin are produced. In poplar, R2R3-MYB transcription factor PdMYB118 positively regulates anthocyanin production to change leaf color. However, the molecular mechanism by which it controls different branches of the phenylpropanoid pathway still remains poorly understood. Here, we reported that in addition to anthocyanin synthesis, lignin deposition and xylem differentiation were regulated by PdMYB118 through inhibiting PagKNAT2/6b gene expression. The transgenic poplar plants overexpressing PdMYB118 accumulated more xylem, lignin and anthocyanin. Transcriptome and reverse transcription quantitative PCR analyses revealed that the expression of PagKNAT2/6b gene which inhibited lignin deposition and xylem differentiation was significantly down-regulated in transgenic poplar plants. Subsequent dual-luciferase reporter and yeast-one-hybrid assays demonstrated that PdMYB118 directly inhibited the transcription of PagKNAT2/6b by binding to the AC elements in its promoter region. Further experiments with transgenic poplar plants overexpressing PagKNAT2/6b demonstrated that overexpression of PagKNAT2/6b in the PdMYB118 overexpression background rescued lignin accumulation and xylem width to the same level of wild type plants. The findings in this work suggest that PdMYB118 is involved in the lignin deposition and xylem differentiation via modulating the expression of PagKNAT2/6b, and the PdMYB118- PagKNAT2/6b model can be used for the genetic breeding of new woody tree with high lignin production.

Rice homolog of Arabidopsis Xylem NAC domain 1 (OsXND1), a NAC transcription factor regulates drought stress responsive root system architecture in indica rice.

Rice yield is greatly constrained by drought stress. In Arabidopsis, XYLEM NAC DOMAIN 1 (XND1) gene regulates the xylem formation, efficiency of water transport, and the delicate equilibrium between drought tolerance and resistance to pathogens. However, diversity and the role of rice homologs of OsXND1 is not reported so far. This study hypothesized that the rice homolog of OsXND1 also regulates drought stress tolerance through modulation of root architecture. Initially, phylogenetic analysis identified two OsXND1 homologs (Os02g0555300 and Os04g0437000) in rice. Further, 14 haplotypes were identified in the OsXND1 of which Hap1 and Hap3 were major haplotypes. The association analysis of OsXND1 with 16 different traits, including 10 root traits, showed three SNPs (Chr02:20972728-Promoter variant; Chr02:20972791-5' UTR variant, and Chr02:20973745-3' UTR variant) were significantly associated with root area, root surface area, total root length, and convex hull area only under drought stress in indica rice. Besides, the superior haplotype of OsXND1 increased the root area, root surface area, total root length, and convex hull area by 46%, 40%, 38%, and 42%, respectively, under drought stress conditions. Therefore, the identified superior haplotype of OsXND1 can be utilized in haplotype breeding programs for the improvement of drought tolerance in rice.

Embolism propagation does not rely on pressure only: time-based shifts in xylem vulnerability curves of angiosperms determine the accuracy of the flow-centrifuge method.

Centrifuges provide a fast approach to quantify embolism resistance of xylem in vulnerability curves (VCs). Since embolism formation is assumingly driven by pressure only, spin time is not standardised for flow centrifuge experiments. Here, we explore to what extent embolism resistance could be spin-time dependent, and hypothesise that changes in hydraulic conductivity (Kh) would shift VCs towards higher water potential (Ψ) values over time. We quantified time-based shifts in flow-centrifuge VCs and their parameter estimations for six angiosperm species by measuring Kh over 15 minutes of spinning at a particular speed, before a higher speed was applied to the same sample. We compared various VCs per sample based on cumulative spin time, and modelled the relationship between Kh, Ψ, and spin-time. Time-based changes of Kh showed considerable increases and decreases at low and high centrifuge speeds, respectively, which generally shifted VCs towards more positive Ψ values. Values corresponding to 50% loss of hydraulic conductivity (P50) became less negative by up to 0.72 MPa in Acer pseudoplatanus, and on average by 8.5% for all six species compared to VCs that did not consider spin-time. By employing an asymptotic exponential model, we estimated time-stable Kh, which improved the statistical significance of VCs in 5 of the 6 species studied. This model also revealed the instability of VCs at short spin times with embolism formation in flow-centrifuges following a saturating exponential growth curve. Although pressure remains the major determinant of embolism formation, spin-time should be considered in flow-centrifuge VCs because not considering the time-dependent stability of Kh overestimates embolism resistance. This spin-time artefact is species-specific, and likely based on relatively slow gas diffusion that is associated with embolism propagation. The accuracy of VCs is improved by determining time-stable Kh values for each centrifuge speed, without considerably extending the experimental time to construct VCs.

An early cladoxylopsid with complex vascular architecture: Paracladoxylon kespekianum gen. et sp. nov. from the Lower Devonian (Emsian) of Quebec, Canada.

PREMISE: Cladoxylopsids, one of the first lineages with complex organization to rise from the plexus of structurally simple plants that comprised the earliest euphyllophyte floras, are moniliformopsid euphyllophytes. They formed Earth's earliest forests by the Middle Devonian and are thought to have given rise to the equisetopsids and probably some fern lineages. The Lower Devonian (Emsian) Battery Point Formation (Quebec, Canada) contains previously unrecognized cladoxylopsids preserved anatomically. One of these provides new data on structural evolution among euphyllophytes and is described here. METHODS: The anatomy and morphology of permineralized axes of the new plant were studied with light and electron microscopy on sections produced using the cellulose acetate peel technique. Morphological comparisons and phylogenetic analysis were used for taxonomic placement of the plant. RESULTS: The plant represents a new species, Paracladoxylon kespekianum Chu et Tomescu, gen. et sp. nov., that has tracheids with modern-looking bordered pits and the complex cauline vascular architecture characteristic of the genus Cladoxylon. Its dissected ultimate appendages have complex regular taxis and a pattern of vascularization that suggests bilateral symmetry. CONCLUSIONS: Paracladoxylon kespekianum is one of the largest Early Devonian euphyllophytes, among the oldest representatives of the cladoxylopsid group, and older than any species of the closely related Cladoxylon by at least 35 million years. It is also one of the oldest anatomically preserved representatives of the cladoxylopsid group. Its anatomical organization pushes the rise of complex vascular architecture among moniliformopsid euphyllophytes deeper in time than previously recognized.

Productive Poplar Genotypes Exhibited Temporally Stable Low Stem Embolism Resistance and Hydraulic Resistance Segmentation at the Stem-Leaf Transition.

Breeding tree genotypes that are both productive and drought-resistant is a primary goal in forestry. However, the relationships between plant hydraulics and yield at the genotype level, and their temporal stabilities, remain unclear. We selected six poplar genotypes from I-101 (Populus alba) × 84 K (P. alba × Popolus tremula var. glandulosa) for experiments in the first and fourth years after planting in a common garden. Measurements included stem embolism resistance, shoot hydraulic resistance and its partitioning between stems and leaves, vessel- and pit-level anatomy, leaf carbon acquisition capacity, carbon allocation to leaves, and aboveground biomass (yield proxy). Significant genetic variations in hydraulic properties and yield were found among genotypes in both years. Productive genotypes had wide vessels, large thin pit membranes, small pit apertures, and shallow pit chambers. Hydraulic resistance was negatively correlated with yield, enabling high stomatal conductance and assimilation rates. Productive genotypes allocated less aboveground carbon and hydraulic resistance to leaves. Temporally stable trade-offs between stem embolism resistance and yield, and between hydraulic segmentation and yield, were identified. These findings highlight the tight link between hydraulic function and yield and suggest that stable trade-offs may challenge breeding poplar genotypes that are both productive and drought-resistant.

Identification of BpEXP family genes and functional characterization of the BpEXPA1 gene in the stems development of Betula platyphylla.

Expansins (EXPs) are unique plant cell wall proteins with the ability to induce cell wall expansion and play potential roles in xylem development. In the present study, a total of 25 BpEXP genes were identified in Betula platyphylla. Results of bioinformatics analysis described that BpEXP gene family was highly conserved in the process of evolution. All these genes were clustered into four groups, EXPA (Expansin A), EXPB (Expansin B), EXLA (Expansin-like A) and EXLB (Expansin-like B), according to phylogenetic analysis and BpEXPA1 was highly homologous to PttEXP1 and PttEXP2. The results of RT-qPCR showed that BpEXPA1 was expressed higher in stems and preferentially expressed in the first internodes, followed by apical buds and the third internodes, promoter expression analysis with GUS assay demonstrated that it was expressed in developing xylem, suggesting that BpEXPA1 might be involved in the development of the primary stems of birch. Overexpression of BpEXPA1 can promote cortex cell expansion and then enlarge the cortex cell area and layer, however inhibit the secondary cell wall deposition and result in the thinner cell wall and larger lumens of xylem fiber in transgenic plants. This study will provide information for investigating the regulation mechanism of BpEXP family genes and gene resources for birch genetics improvement.

Interorgan, intraorgan and interplant communication mediated by nitric oxide and related species.

Plant survival to a potential plethora of diverse environmental insults is underpinned by coordinated communication amongst organs to help shape effective responses to these environmental challenges at the whole plant level. This interorgan communication is supported by a complex signal network that regulates growth, development and environmental responses. Nitric oxide (NO) has emerged as a key signalling molecule in plants. However, its potential role in interorgan communication has only recently started to come into view. Direct and indirect evidence has emerged supporting that NO and related species (S-nitrosoglutathione, nitro-linolenic acid) are mobile interorgan signals transmitting responses to stresses such as hypoxia and heat. Beyond their role as mobile signals, NO and related species are involved in mediating xylem development, thus contributing to efficient root-shoot communication. Moreover, NO and related species are regulators in intraorgan systemic defence responses aiming an effective, coordinated defence against pathogens. Beyond its in planta signalling role, NO and related species may act as ex planta signals coordinating external leaf-to-leaf, root-to-leaf but also plant-to-plant communication. Here, we discuss these exciting developments and emphasise how their manipulation may provide novel strategies for crop improvement.

Hydrogen Peroxide Is Involved in Methane-Alleviated Cadmium Toxicity in Alfalfa (Medicago sativa L.) Seedlings by Enhancing Cadmium Chelation onto Root Cell Walls.

Although previous studies have demonstrated that methane (CH4) can mitigate the toxicity of cadmium (Cd) in alfalfa seedlings, the CH4-rich water used in these studies may create hypoxic conditions, potentially influencing the experimental outcomes. Therefore, this study aimed to investigate whether CH4 can reduce Cd toxicity in alfalfa seedlings without the interference of hypoxia and to analyze its underlying mechanisms. Here, it was observed that supplementing oxygen with saturated CH4-rich water can significantly alleviate the inhibition of 75 µM CdCl2 on the growth of alfalfa (Medicago sativa L.) seedlings. Less Cd accumulation was also observed in both root and shoot parts, which could be explained by the CH4-altered cell wall components in alfalfa seedling roots, including covalent and ionic soluble pectin, and the degree of demethylation in pectin, thus enabling a higher proportion of Cd binding to the cell walls and reducing the entry of Cd into the cells. The above actions of CH4 were accompanied by an increase in hydrogen peroxide (H2O2) content and NADPH oxidase activity, which could be blocked by the addition of the NADPH oxidase inhibitor diphenylene iodonium (DPI). Taken together, these results implied that exogenously applied CH4 could alleviate Cd toxicity in alfalfa seedlings by enhancing Cd chelation onto the root cell walls, which might be closely associated with NADPH oxidase-dependent H2O2 signals. These findings could provide insight into the mechanism through which CH4 alleviates Cd toxicity in alfalfa plants.

Tree Rings Mercury Controlled by Atmospheric Gaseous Elemental Mercury and Tree Physiology.

Tree rings are an emerging atmospheric mercury (Hg) archive. Questions have arisen, though, regarding their mechanistic controls and reliability. Here, we report contrasting tree-ring Hg records in three collocated conifer species: Norway spruce (Picea abies), Scots pine (Pinus sylvestris), and European larch (Larix decidua), which are from a remote boreal forest. Centennial atmospheric Hg trends at the site, derived from varved lake sediments, peats, and atmospheric monitoring, indicated a steady rise from the 1800s, peaking in the 1970s, and then declining. Prior to ca. 2005, larch and spruce tree rings reproduced the peak in the atmospheric Hg trend, while pine tree rings peaked in the 1930s, likely due to the prolonged sapwood period and ambiguity in the heartwood-sapwood boundary of pine. Since ca. 2005, tree rings from all species showed increasing Hg concentrations in the physiologically active outer rings despite declining atmospheric Hg concentrations. The good agreement between Hg and nitrogen concentrations in active tree-ring cells indicates a similar transport mechanism and cautions against their applicability as atmospheric Hg archives. Our results suggest that tree-ring Hg records are controlled by atmospheric Hg and tree physiology. We provide recommendations for using tree-ring Hg archives that take tree physiology into account.

A one-way ticket: wheat roots do not functionally refill xylem emboli following rehydration.

Understanding xylem embolism spread in roots is essential for predicting the loss of function across root systems during drought. However, the lasting relevance of root embolism to plant recovery depends on whether roots can refill xylem emboli and resume function after rehydration. Using MicroCT and optical and dye staining methods, we investigated embolism repair in rehydrated intact roots of wheat (Triticum aestivum L.'Krichauff') exposed to a severe water deficit of -3.5 MPa, known to cause approximately 30 % total root network embolism in this species. Air emboli in the xylem vessels of intact roots remained clearly observable using MicroCT after overnight rehydration. This result was verified by xylem staining of the root system and optical quantification of emboli, both of which indicated a lack of functional root xylem recovery 60 hours following soil re-saturation. The absence of root xylem refilling in wheat has substantial implications for how we understand plant recovery after drought. Our findings suggest that xylem embolism causes irreversible damage to the soil-root hydraulic connection in affected parts of the root network.

De novo Whole-Genome Assembly of the 10-Gigabase Fokienia Hodginsii Genome to Reveal Differential Epigenetic Events Between Callus and Xylem.

Fokienia hodginsii (F. hodginsii), belonging to the genus Fokienia of the Cupressaceae. F. hodginsii has significant application value due to its wood properties and great research value in evolutionary studies as a gymnosperm. However, the genome of F. hodginsii remains unknown due to the large size of gymnosperms genome. Pacific Bioscience sequencing, Hi-C mapping, whole-genome Bisulfite Sequencing (BS-Seq), long-read isoform sequencing (Iso-Seq), direct RNA sequencing (DRS), quantitative proteomics, and metabonomics analysis are employed to facilitate genome assembly, gene annotation, and investigation into epigenetic mechanisms. In this study, the 10G F. hodginsii genome is assembled into 11 chromosomes. Furthermore, 50 521 protein-coding genes are annotated and determined that 65% of F. hodginsii genome comprises repetitive sequences. It is discovered that transposable element (TE)-including introns is associated with higher expression. The DNA methylome of F. hodginsii reveals that xylem has a higher DNA methylation level compared to callus. Moreover, DRS reveals the significant alterations in RNA full-length ratio, which potentially associated with poly(A) length (PAL) and alternative polyadenylation (APA). Finally, the morphology measurement and metabonomics analysis revealed the difference of 14 cultivars. In summary, the genomes and epigenetics datasets provide a molecular basis for callus formation in the gymnosperm family.

PagHAM4a-PagSCL21 and PagHAM4b-PagTCP20 modules positively regulate cambial activity and its differentiation into secondary xylem in poplar.

Stem secondary xylem produced by cambial division and differentiation is the main source of tree biomass. The secondary xylem formation involves a complex transcriptional regulatory network, however, the underlying mechanism is still being explored. Here, we report that PagHAM4a and PagHAM4b are positive regulators of cambial differentiation into secondary xylem in hybrid poplar (Populus alba × Populus glandulosa clone 84K). Overexpression of PagHAM4a and PagHAM4b enhanced cambial activity and increased the number of secondary xylem cells in the stems of poplar. By contrast, single or double mutations of PagHAM4a and PagHAM4b by CRISPR-Cas9 decreased cambial activity, leading to a significant reduction of secondary xylem. Neither overexpression nor mutation of the two genes affected the size of vessels and fibers in xylem. Both PagHAM4a- and PagHAM4b-regulated gene networks were mainly centered in the stage when cambium had just initiated secondary growth, but the molecular networks regulated by the two genes were distinct. Further analysis revealed that PagSCL21 and PagTCP20 are direct targets of PagHAM4a and PagHAM4b, respectively, and their overexpression also promoted cambial differentiation into secondary xylem. Taken together, we identified two novel key regulatory modules in poplar, PagHAM4a-PagSCL21 and PagHAM4b-PagTCP20, which provide new insights into the mechanism of secondary xylem formation in trees.

Optimizing 'red Fuji' apple quality: Auxin-mediated calcium distribution via fruit-stalk in bagging practices.

In apples, a bottleneck effect in calcium (Ca) transport within fruit stalk has been observed. To elucidate that how auxin affects Ca forms and distribution in the apple fruit stalk, we investigated the effects of different concentrations of auxin treatment (0, 10, 20, and 30 mg·L-1) on Ca content, forms, distribution, and fruit quality during later stages of fruit expansion. The results showed that auxin treatment led to a dramatic reduction in total Ca content in stalk, while an approximately 30 % increase in fruit. Furthermore, auxin treatment effectively enhanced the functionality of xylem vessels in vascular bundles of the stalk in bagged apples. Finally, TOPSIS method was used to assess fruit quality, with treatments ranked as follows: IAA20 > NAA20 > IAA30 > IAA10 > CK > NPA. The findings lay a foundation for further studies on the bottleneck in Ca transport within stalk, uneven distribution of Ca in fruit, and provide insights into Ca utilization efficiency in bagged apples.

An appreciation of apex-to-base variation in xylem traits will lead to more precise understanding of xylem phenotypic plasticity.

Xylem air embolism is the primary cause of drought-related tree mortality. Phenotypic plasticity of xylem traits is key for species acclimation to environmental variability and evolution. It is widely believed that plants increase xylem embolism resistance in response to drought. However, I argue that this hypothesis, based on extensive literature, relies on sampling methods that overlook predictable anatomical patterns, potentially biasing our understanding of acclimation and adaptation strategies.

Weak link or strong foundation? Vulnerability of fine root networks and stems to xylem embolism.

Resolving the position of roots in the whole-plant hierarchy of drought-induced xylem embolism resistance is fundamental for predicting when species become isolated from soil water resources. Published research generally suggests that roots are the most vulnerable organ of the plant vascular system, although estimates vary significantly. However, our knowledge of root embolism excludes the fine roots (< 2 mm diameter) that form the bulk of total absorptive surface area of the root network for water and nutrient uptake. We measured fine root and stem xylem vulnerability in 10 vascular plant species from the major land plant clades (five angiosperms, three conifers, a fern and lycophyte), using standardised in situ methods (Optical Methods and MicroCT). Mean fine root embolism resistance across the network matched or exceeded stems in all study species. In six of these species (one fern, one lycophyte, three conifers and one angiosperm), fine roots were significantly more embolism resistant than stems. No clear relationship was found between root xylem conduit diameter and vulnerability. These results provide insight into the resistance of the plant hydraulic pathway at the site of water and nutrient uptake, and challenge the long-standing assumption that fine roots are more vulnerable than stems.

Growing at the arid edge: Anatomical variations in leaves are more extensive than in stems of five Mediterranean species across contrasting moisture regimes.

PREMISE: Increasing aridity in the Mediterranean region affects ecosystems and plant life. Various anatomical changes in plants help them cope with dry conditions. This study focused on anatomical differences in leaves and xylem of five co-occurring Mediterranean plant species namely Quercus calliprinos, Pistacia palaestina, Pistacia lentiscus, Rhamnus lycioides, and Phillyrea latifolia in wet and dry sites. METHODS: Stomatal density, stomatal length, leaf mass area, lamina composition, percentage of intercellular air spaces, and mesophyll cell area in leaves of plants in wet and dry sites were analyzed. Xylem anatomy was assessed through vessel length and area in branches. RESULTS: In the dry site, three species had increased stomatal density and decreased stomatal length. Four species had increased palisade mesophyll and reduced air space volume. In contrast, phenotypic changes in the xylem were less pronounced; vessel length was unaffected by site conditions, but vessel diameter decreased in two species. Intercellular air spaces proved to be the most dynamic anatomical feature. Quercus calliprinos had the most extensive anatomical changes; Rhamnus lycioides had only minor changes. All these changes were observed in comparison to the species in the wet site. CONCLUSIONS: This study elucidated variations in anatomical responses in leaves among co-occurring Mediterranean plant species and identified the most dynamic traits. Understanding these adaptations provides valuable insights into the ability of plants to thrive under changing climate conditions.

Tree drought-mortality risk depends more on intrinsic species resistance than on stand species diversity.

Increasing tree diversity is considered a key management option to adapt forests to climate change. However, the effect of species diversity on a forest's ability to cope with extreme drought remains elusive. In this study, we assessed drought tolerance (xylem vulnerability to cavitation) and water stress (water potential), and combined them into a metric of drought-mortality risk (hydraulic safety margin) during extreme 2021 or 2022 summer droughts in five European tree diversity experiments encompassing different biomes. Overall, we found that drought-mortality risk was primarily driven by species identity (56.7% of the total variability), while tree diversity had a much lower effect (8% of the total variability). This result remained valid at the local scale (i.e within experiment) and across the studied European biomes. Tree diversity effect on drought-mortality risk was mediated by changes in water stress intensity, not by changes in xylem vulnerability to cavitation. Significant diversity effects were observed in all experiments, but those effects often varied from positive to negative across mixtures for a given species. Indeed, we found that the composition of the mixtures (i.e., the identities of the species mixed), but not the species richness of the mixture per se, is a driver of tree drought-mortality risk. This calls for a better understanding of the underlying mechanisms before tree diversity can be considered an operational adaption tool to extreme drought. Forest diversification should be considered jointly with management strategies focussed on favouring drought-tolerant species.

Nearly instantaneous stem diameter response to fluctuations in the atmospheric water demand.

Changes in vapour pressure deficit (VPD) can lead to the depletion and replenishment of stem water pools to buffer water potential variations in the xylem. Yet, the precise velocity at which stem water pools track environmental cues remains poorly explored. Nine eucalyptus seedlings grown in a glasshouse experienced high-frequency environmental oscillations, and their stem radial variations (ΔR) were monitored at a 30-second temporal resolution in upper and lower stem locations and on the bark and xylem. The stem ΔR response to VPD changes was nearly instantaneous (< 1 minute), while temperature lagged behind stem ΔR. No temporal differences in the stem ΔR response were observed between locations. Punctual gravimetric measurements confirmed the synchrony between transpiration and stem ΔR dynamics. These results indicate (i) that stem-stored water can respond to the atmospheric evaporative demand much faster than commonly assumed and (ii) that the origin of the water released to the transpiration stream seems critical in determining time lags in stem water pool dynamics. Near-zero time lags may be explained by the high elasticity of eucalyptus woody tissues and the predominant water use from the xylem, circumventing the hydraulic radial barriers to water flow from/to the outer tissues.

Spatial and Temporal Freezing Dynamics of Leaves Revealed by Time-Lapse Imaging.

Freezing air temperatures kill most leaves, yet the leaves of some species can survive these events. Tracking the temporal and spatial dynamics of freezing remains an impediment to characterizing frost tolerance. Here we deploye time-lapse imaging and image subtraction analysis, coupled with fine wire thermocouples, to discern the in situ spatial dynamics of freezing and thawing. Our method of analysis of pixel brightness reveals that ice formation in leaves exposed to natural frosts initiates in mesophyll before spreading to veins, and that while ex situ xylem sap freezes near 0°C, in situ xylem sap has a freezing point of -2°C in our model freezing-resistant species of Lonicera. Photosynthetic rates in leaves that have been exposed to a rapid freeze or thaw do not recover, but leaves exposed to a slow, natural freezing and thawing to -10°C do recover. Using this method, we are able to quantify the spatial formation and timing of freezing events in leaves, and suggest that in situ and ex situ freezing points for xylem sap can differ by more than 4°C depending on the rate of temperature decline.