Cold adaptation drives variability in needle structure and anatomy in P inus sylvestris L. along a 1,900 km temperate–boreal transect

Physiological trait coordination and variability across and within three Pinus species

Studies have explored how traits separate plants ecologically and the trade-offs that underpin this separation. However, uncertainty remains as to the taxonomic scale at which traits can predictably separate species. We studied how physiological traits separated three Pinus (Pinus banksiana, Pinus resinosa, and Pinus strobus) species across three sites. We collected traits from four common leaf and branch measurements (light-response curves, CO2-response curves, pressure-volume curves, and hydraulic vulnerability curves) across each species and site. While common, these measurements are not typically measured together due to logistical constraints. Few traits varied across species and sites as expected given the ecological preferences of the species and environmental site characteristics. Some trait trade-offs present at broad taxonomic scales were observed across the three species, but most were absent within species. Certain trade-offs contrasted expectations observed at broader scales but followed expectations given the species' ecological preferences. We emphasize the need to both clarify why certain traits are being studied, as variation in unexpected but ecologically meaningful ways often occurs and certain traits might not vary substantially within a given lineage (e.g. hydraulic vulnerability in Pinus), highlighting the role a trait selection in trait ecology.

Phylotranscriptomic and ecological analyses reveal the evolution and morphological adaptation of Abies

Coniferous forests are under severe threat of the rapid anthropogenic climate warming. Abies (firs), the fourth-largest conifer genus, is a keystone component of the boreal and temperate dark-coniferous forests and harbors a remarkably large number of relict taxa. However, the uncertainty of the phylogenetic and biogeographic history of Abies significantly impedes our prediction of future dynamics and efficient conservation of firs. In this study, using 1,533 nuclear genes generated from transcriptome sequencing and a complete sampling of all widely recognized species, we have successfully reconstructed a robust phylogeny of global firs, in which four clades are strongly supported and all intersectional relationships are resolved, although phylogenetic discordance caused mainly by incomplete lineage sorting and hybridization was detected. Molecular dating and ancestral area reconstruction suggest a Northern Hemisphere high-latitude origin of Abies during the Late Cretaceous, but all extant firs diversified during the Miocene to the Pleistocene, and multiple continental and intercontinental dispersals took place in response to the late Neogene climate cooling and orogenic movements. Notably, four critically endangered firs endemic to subtropical mountains of China, including A. beshanzuensis, A. ziyuanensis, A. fanjingshanensis and A. yuanbaoshanensis from east to west, have different origins and evolutionary histories. Moreover, three hotspots of species richness, including western North America, central Japan, and the Hengduan Mountains, were identified in Abies. Elevation and precipitation, particularly precipitation of the coldest quarter, are the most significant environmental factors driving the global distribution pattern of fir species diversity. Some morphological traits are evolutionarily constrained, and those linked to elevational variation (e.g., purple cone) and cold resistance (e.g., pubescent branch and resinous bud) may have contributed to the diversification of global firs. Our study sheds new light on the spatiotemporal evolution of global firs, which will be of great help to forest management and species conservation in a warming world.

Exploring Intraspecific Trait Variation in a Xerophytic Moss Species Indusiella thianschanica (Ptychomitriaceae) across Environmental Gradients on the Tibetan Plateau

Investigating intraspecific trait variability is crucial for understanding plant adaptation to various environments, yet research on lithophytic mosses in extreme environments remains scarce. This study focuses on Indusiella thianschanica Broth. Hal., a unique lithophytic moss species in the extreme environments of the Tibetan Plateau, aiming to uncover its adaptation and response mechanisms to environmental changes. Specimens were collected from 26 sites across elevations ranging from 3642 m to 5528 m, and the relationships between 23 morphological traits and 15 environmental factors were analyzed. Results indicated that coefficients of variation (CV) ranged from 5.91% to 36.11%, with gametophyte height (GH) and basal cell transverse wall thickness (STW) showing the highest and lowest variations, respectively. Temperature, elevation, and potential evapo-transpiration (PET) emerged as primary environmental drivers. Leaf traits, especially those of the leaf sheath, exhibited a more pronounced response to the environment. The traits exhibited apparent covariation in response to environmental challenges and indicated flexible adaptive strategies. This study revealed the adaptation and response patterns of different morphological traits of I. thianschanica to environmental changes on the Tibetan Plateau, emphasizing the significant effect of temperature on trait variation. Our findings deepen the understanding of the ecology and adaptive strategies of lithophytic mosses.

Anatomical and morphological changes in Pinus sylvestris and Larix sibirica needles under impact of emissions from a large aluminum enterprise

Species-specific anatomical and morphological characteristics of Pinus sylvestris and Larix sibirica needles were studied at different levels of tree stand pollution by aluminum smelter emissions. The anatomical characteristics of the needle were studied using light microscopy. The level of tree stand pollution was determined using the cluster analysis outcomes of the pollutant elements content (fluorine, sulfur, and heavy metals) in the needles. Four levels of tree stand pollution were separated: low, moderate, high, and critical, as well as background tree stand in unpolluted areas. It was found that the state of tree phytomass deteriorated with increasing levels of pollution (from low to critical): pine crown defoliation increased to 85%, and larch defoliation increased to 65%. The life span of pine needles was reduced to 2-3 years, with a background value of 6-7 years. The change of morphological parameters was more pronounced in P. sylvestris: the weight and length of the 2-year-old shoot decreased by 2.7-3.1 times compared to the background values; the weight of needles on the shoot and the number of needle pairs on the shoot-by 1.9-2.1 times. The length of the needle and shoot and the number of L. sibirica brachyblasts decreased by 1.8-1.9 times. The anatomical parameters of the needle also changed to a greater extent in P. sylvestris. Up to the high level of tree pollution, we observed a decrease in the cross-sectional area of the needle, central cylinder, vascular bundle, area and thickness of mesophyll, number and diameter of resin ducts by 18-66% compared to background values. At the critical pollution level, when the content of pollutant elements in pine needles reached maximum values, the anatomical parameters of the remaining few green needles were close to background values. In our opinion, this may be due to the activation of mechanisms aimed at maintaining the viability of trees. A reduction in thickness and area of assimilation tissue in the L. sibirica needle was detected only at the critical pollution level. An upward trend in these parameters was found at low, medium, and high pollution levels of tree stand, which may indicate an adaptive nature. The results suggested that at a similar pollution level of trees, the greatest amount of negative anatomical and morphological changes were recorded in pine needles, which indicates a greater sensitivity of this species to technogenic emissions.

Long term effects of ionising radiation in the Chernobyl Exclusion zone on DNA integrity and chemical defence systems of Scots pine (Pinus sylvestris)

The Chernobyl Nuclear Power Plant (ChNPP) accident in 1986 resulted in extremely high levels of acute ionising radiation, that killed or damaged Scots pine (Pinus sylvestris) trees in the surrounding areas. Dead trees were cleared and buried, and new plantations established a few years later. Today, more than three decades later, gamma and beta-radiation near the ChNPP is still elevated compared with ambient levels but have decreased by a factor of 300 and 100, respectively. In the present work, Scots pine-trees growing at High (220 μGy h-1), Medium (11 μGy h-1), and Low (0.2 μGy h-1) total (internal + external) dose rates of chronically elevated ionising radiation in the Chernobyl Exclusion zone were investigated with respect to possible damage to DNA, cells and organelles, as well as potentially increased levels of phenolic and terpenoid antioxidants. Scots pine from the High and Medium radiation sites had elevated levels of DNA damage in shoot tips and needles as shown by the COMET assay, as well as increased numbers of resin ducts and subcellular abnormalities in needles. Needles from the High radiation site showed elevated levels of monoterpenes and condensed tannins compared with those from the other sites. In conclusion, more than three decades after the ChNPP accident substantial DNA damage and (sub)cellular effects, but also mobilisation of stress-protective substances possessing antioxidant activity were observed in Scots pine trees growing at elevated levels of ionising radiation. This demonstrates that the radiation levels in the Red Forest still significantly impact the plant community.

The energy-water limitation threshold explains divergent drought responses in tree growth, needle length, and stable isotope ratios

Predicted increases in extreme droughts will likely cause major shifts in carbon sequestration and forest composition. Although growth declines during drought are widely documented, an increasing number of studies have reported both positive and negative responses to the same drought. These divergent growth patterns may reflect thresholds (i.e., nonlinear responses) promoted by changes in the dominant climatic constraints on tree growth. Here we tested whether stemwood growth exhibited linear or nonlinear responses to temperature and precipitation and whether stemwood growth thresholds co-occurred with multiple thresholds in source and sink processes that limit tree growth. We extracted 772 tree cores, 1398 needle length records, and 1075 stable isotope samples from 27 sites across whitebark pine's (Pinus albicaulis Engelm.) climatic niche in the Sierra Nevada. Our results indicated that a temperature threshold in stemwood growth occurred at 8.4°C (7.12-9.51°C; estimated using fall-spring maximum temperature). This threshold was significantly correlated with thresholds in foliar growth, as well as carbon (δ13 C) and nitrogen (δ15 N) stable isotope ratios, that emerged during drought. These co-occurring thresholds reflected the transition between energy- and water-limited tree growth (i.e., the E-W limitation threshold). This transition likely mediated carbon and nutrient cycling, as well as important differences in growth-defense trade-offs and drought adaptations. Furthermore, whitebark pine growing in energy-limited regions may continue to experience elevated growth in response to climate change. The positive effect of warming, however, may be offset by growth declines in water-limited regions, threatening the long-term sustainability of the recently listed whitebark pine species in the Sierra Nevada.

Miniature Noninvasive Sensor Based on Impedance-Change Detection in Branches for Measuring Branch Ice Content in Overwintering Woody Plants

Advancements in detection instruments have enabled the real-time acquisition of water information during plant growth; however, the real-time monitoring of freeze-thaw information during plant overwintering remains a challenge. Based on the relationship between the change in the water-ice ratio and branch impedance during freezing, a miniature noninvasive branch volume ice content (BVIC) sensor was developed for monitoring real-time changes in volumetric ice content and the ice freeze-thaw rate of woody plant branches during the overwintering period. The results of the performance analysis of the impedance measurement circuit show that the circuit has a lateral sensitivity range, measurement range, resolution, measurement accuracy, and power consumption of 0-35 mm, 0-100%, 0.05%, ±1.76%, and 0.25 W, respectively. The dynamic response time was 0.296 s. The maximum allowable error by the output voltage fluctuation, owing to the ambient temperature and humidity, was only ±0.635%, which meets the actual use requirements. The calibration curve fit coefficients were >0.98, indicating a significant correlation. The ice content of plant branches under cold stress was measured for indoor and field environments, and the sensors could effectively monitor changes in the branch ice content in plants exposed to cold stress. Additionally, they can differentiate between plants with different cold resistances, indicating the reliability of the BVIC sensor.

At least it is a dry cold: the global distribution of freeze-thaw and drought stress and the traits that may impart poly-tolerance in conifers

Conifers inhabit some of the most challenging landscapes where multiple abiotic stressors (e.g., aridity, freezing temperatures) often co-occur. Physiological tolerance to multiple stressors ('poly-tolerance') is thought to be rare because exposure to one stress generally limits responses to another through functional trade-offs. However, the capacity to exhibit poly-tolerance may be greater when combined abiotic stressors have similar physiological impacts, such as the disruption of hydraulic function imposed by drought or freezing. Here, we reviewed empirical data in light of theoretical expectations for conifer adaptations to drought and freeze-thaw cycles with particular attention to hydraulic traits of the stem and leaf. Additionally, we examined the commonality and spatial distribution of poly-stress along indices of these combined stressors. We found that locations with the highest values of our poly-stress index (PSi) are characterized by moderate drought and moderate freeze-thaw, and most of the global conifer distribution occupies areas of moderate poly-stress. Among traits examined, we found diverse responses to the stressors. Turgor loss point did not correlate with freeze-thaw or drought stress individually, but did with the PSi, albeit inverse to what was hypothesized. Leaf mass per area was more strongly linked with drought stress than the poly-stress and not at all with freeze-thaw stress. In stems, the water potential causing 50% loss of hydraulic conductivity became more negative with increasing drought stress and poly-stress but did not correlate with freeze-thaw stress. For these traits, we identified a striking lack of coverage for substantial portions of species ranges, particularly at the upper boundaries of their respective PSis, demonstrating a critical gap in our understanding of trait prevalence and plasticity along these stress gradients. Future research should investigate traits that confer tolerance to both freeze-thaw and drought stress in a wide range of species across broad geographic scales.

How much biology is in the product? Role and relevance of biological evolution and function for bio-inspired design

Bio-inspired design (BID) means the concept of transferring functional principles from biology to technology. The core idea driving BID-related work is that evolution has shaped functional attributes, which are termed “adaptations” in biology, to a high functional performance by relentless selective pressure. For current methods and tools, such as data bases, it is implicitly supposed that the considered biological models are adaptations and their functions already clarified. Often, however, the identification of adaptations and their functional features is a difficult task which is not yet accomplished for numerous biological structures, as happens to be the case also for various organismic features from which successful BID developments were derived. This appears to question the relevance of the much stressed importance of evolution for BID. While it is obviously possible to derive an attractive technical principle from an observed biological effect without knowing its original functionality, this kind of BID (“analog BID”) has no further ties to biology. In contrast, a BID based on an adaptation and its function (“homolog BID”) is deeply embedded in biology. It is suggested that a serious and honest clarification of the functional background of a biological structure is an essential first step in devising a BID project, to recognize possible problems and pitfalls as well as to evaluate the need for further biological analysis.

Stomatal density in Pinus sylvestris as an indicator of temperature rather than CO2 : Evidence from a pan-European transect

The commonly observed negative relationship between stomatal density (SD) and atmospheric CO₂ has led to SD being proposed as an indicator of atmospheric CO₂ concentration. The use of SD as a proxy for CO₂ , however, has been hampered by an insufficient understanding of the intraspecific variation of this trait. We hypothesized that SD in Pinus sylvestris, a widely distributed conifer, varies geographically and that this variation is determined by major climatic variables. By sampling needles from naturally growing trees along a latitudinal range of 32.25°, equivalent to 13.7°C gradient of mean annual temperature (MAT) across Europe, we found that SD decreased from the warmest southern sites to the coldest sites in the north at a rate of 4 stomata per mm² for each 1°C, with MAT explaining 44% of the variation. Additionally, samples from a provenance trial exhibited a positive relationship between SD and the MAT of the original localities, suggesting that high SD is an adaptation to warm temperature. Our study revealed one of the strongest intraspecific relationships between SD and climate in any woody species, supporting the utility of SD as a temperature, rather than direct CO₂ , proxy. In addition, our results predict the response of SD to climate warming.

Higher biomass partitioning to absorptive roots improves needle nutrition but does not alleviate stomatal limitation of northern Scots pine

Harsh environmental conditions affect both leaf structure and root traits. However, shoot growth in high-latitude systems is predominately under photoperiod control while root growth may occur for as long as thermal conditions are favorable. The different sensitivities of these organs may alter functional relationships above- and belowground along environmental gradients. We examined the relationship between absorptive root and foliar traits of Scots pine trees growing in situ along a temperate-boreal transect and in trees grown in a long-term common garden at a temperate latitude. We related changes in foliar nitrogen, phosphorus, specific leaf area, needle mass and ¹³ C signatures to geographic trends in absorptive root biomass to better understand patterns of altered tree nutrition and water balance. Increased allocation to absorptive fine roots was associated with greater uptake of soil nutrients and subsequently higher needle nutrient contents in the northern provenances compared with more southern provenances when grown together in a common garden setting. In contrast, the leaf δ¹³ C in northern and southern provenances were similar within the common garden suggesting that higher absorptive root biomass fractions could not adequately increase water supply in warmer climates. These results highlight the importance of allocation within the fine-root system and its impacts on needle nutrition while also suggesting increasing stomatal limitation of photosynthesis in the context of anticipated climatic changes.

Non-adaptive evolutionary processes governed the diversification of a temperate conifer lineage after its migration into the tropics

Constructing phylogenetic relationships among closely related species is a recurrent challenge in evolutionary biology, particularly for long-lived taxa with large effective population sizes and uncomplete reproductive isolation, like conifers. Conifers further have slow evolutionary rates, which raises the question of whether adaptive or non/adaptive processes were predominantly involved when they rapidly diversified after migrating from temperate regions into the tropical mountains. Indeed, fine-scale phylogenetic relationships within several conifer genus remain under debate. Here, we studied the phylogenetic relationships of endemic firs (Abies, Pinaceae) discontinuously distributed in the montane forests from the Southwestern United States to Guatemala, and addressed several hypotheses related to adaptive and non-adaptive radiations. We derived over 80 K SNPs from genotyping by sequencing (GBS) for 45 individuals of nine Mesoamerican species to perform phylogenetic analyses. Both Maximum Likelihood and quartets-inference phylogenies resulted in a well-resolved topology, showing a single fir lineage divided in four subgroups that coincided with the main mountain ranges of Mesoamerica; thus having important taxonomic implications. Such subdivision fitted a North-South isolation by distance framework, in which non-adaptive allopatric processes seemed the rule. Interestingly, several reticulations were observed within subgroups, especially in the central-south region, which may explain past difficulties for generating infrageneric phylogenies. Further evidence for non-adaptive processes was obtained from analyses of 21 candidate-gene regions, which exhibited diminishing values of π_a/π_s and K_a/K_s with latitude, thus indicating reduced efficiency of purifying selection towards the Equator. Our study indicates that non-adaptive allopatric processes may be key generators of species diversity and endemism in the tropics.

Divergent patterns between phenotypic and genetic variation in Scots pine

In boreal forests, autumn frost tolerance in seedlings is a critical fitness component because it determines survival rates during regeneration. To understand the forces that drive local adaptation in this trait, we conducted freezing tests in a common garden setting for 54 Pinus sylvestris (Scots pine) populations (>5000 seedlings) collected across Scandinavia into western Russia, and genotyped 24 of these populations (>900 seedlings) at >10 000 SNPs. Variation in cold hardiness among populations, as measured by Q_ST , was above 80% and followed a distinct cline along latitude and longitude, demonstrating significant adaptation to climate at origin. In contrast, the genetic differentiation was very weak (mean F_ST 0.37%). Despite even allele frequency distribution in the vast majority of SNPs among all populations, a few rare alleles appeared at very high or at fixation in marginal populations restricted to northwestern Fennoscandia. Genotype-environment associations showed that climate variables explained 2.9% of the genetic differentiation, while genotype-phenotype associations revealed a high marker-estimated heritability of frost hardiness of 0.56, but identified no major loci. Very extensive gene flow, strong local adaptation, and signals of complex demographic history across markers are interesting topics of forthcoming studies on this species to better clarify signatures of selection and demography.

The flattened and needlelike leaves of the pine family (Pinaceae) share a conserved genetic network for adaxial-abaxial polarity but have diverged for photosynthetic adaptation

Background

Leaves have highly diverse morphologies. However, with an evolutionary history of approximately 200 million years, leaves of the pine family are relatively monotonous and often collectively called “needles”, although they vary in length, width and cross-section shapes. It would be of great interest to determine whether Pinaceae leaves share similar morpho-physiological features and even consistent developmental and adaptive mechanisms.

Results

Based on a detailed morpho-anatomical study of leaves from all 11 Pinaceae genera, we particularly investigated the expression patterns of adaxial-abaxial polarity genes in two types of leaves (needlelike and flattened) and compared their photosynthetic capacities. We found that the two types of leaves share conserved spatial patterning of vasculatures and genetic networks for adaxial-abaxial polarity, although they display different anatomical structures in the mesophyll tissue differentiation and distribution direction. In addition, the species with needlelike leaves exhibited better photosynthetic capacity than the species with flattened leaves.

Conclusions

Our study provides the first evidence for the existence of a conserved genetic module controlling adaxial-abaxial polarity in the development of different Pinaceae leaves.

A fingerprint of climate change across pine forests of Sweden

Climate change has likely altered high-latitude forests globally, but direct evidence remains rare. Here we show that throughout a ≈1000-km transect in Scots pine (Pinus sylvestris L.) forests in Sweden, mature trees in ≈2015 had longer needles with shorter lifetimes than did trees in ≈1915. These century-scale shifts in needle traits were detected by sampling needles at 74 sites from 2012 to 2017 along the same transect where needle traits had been assessed at 57 sites in 1914-1915. Climate warming of ≈1 °C all along the transect in the past century has driven this temporal shift in foliage traits known to be physiologically critical to growth and carbon cycling processes. These century-scale changes in Scandinavian Scots pine forests represent a fingerprint of climate change on a fundamental biological element, the leaf, with repercussions for productivity and sensitivity to future climate, which are likely to be mirrored by similar changes for evergreen conifers across the boreal biome.

Does climate-related in situ variability of Scots pine (Pinus sylvestris L.) needles have a genetic basis? Evidence from common garden experiments

The correlations of phenotypic traits with environmental drivers suggest that variability of these traits is a result of natural selection, especially if such trait correlations are based on genetic variability. We hypothesized that in situ correlations of structural needle traits of Scots pine (Pinus sylvestris L) with minimal winter temperature (Tmin) reported previously from a temperate/boreal transect would be conserved when plants are cultivated under common conditions. We tested this hypothesis by analyzing needles from two common gardens located in the temperate zone, one including adult trees and the other juvenile seedlings. The majority of adult needle traits for which correlations with Tmin were found in the field turned out to be under environmental influence. In contrast, the majority of traits studied in juvenile needles were correlated with the original Tmin suggesting the role of past natural selection in shaping their variability. Juvenile needles thus appeared to be inherently less plastic than adult needles, perhaps reflecting the stronger selective pressure acting during juvenile, as compared with adult, ontogenetic stage. Genetically based cold-climate adaptation in either juvenile or adult needles, or both, involved an increase in leaf mass per area and leaf density, decrease in needle length, reduction in the amount of xylem and phloem, increase in thickness of epidermis, decrease in tracheid diameter and increase in tracheid density, and increase in diameter and volume fraction of resin ducts. We also show that at least some traits, such as transverse xylem and phloem areas and number of fibers, scale with needle length, suggesting that climate-related trait variation may also be mediated by changes in needle length. Moreover, slopes of these allometric relationships may themselves be plastically modified. The phenotypic syndrome typical of needles from cold environments may thus be under environmental, genetic and allometric control.

Is Resistance to Mountain Pine Beetle Associated with Genetic Resistance to White Pine Blister Rust in Limber Pine?

Limber pine (Pinus flexilis James) co-evolved with the mountain pine beetle (Dendroctonus ponderosae Hopkins; MPB) and is now also challenged by the non-native pathogen Cronartium ribicola (J.C. Fisch.) that causes the lethal disease white pine blister rust (WPBR). Previous research suggests that trees infected with WPBR can be preferred hosts for MPB. Using resin duct traits associated with MPB resistance, we tested for a relationship between resistance to MPB and WPBR in limber pine, in the absence of either biological agent. These analyses will help evaluate if MPB historically may have contributed to natural selection for WPBR resistance in advance of WPBR invasion, and could help explain the unusually high frequency of the dominant Cr4 allele for complete resistance to WPBR in limber pine populations of the Southern Rocky Mountains. Resin duct production, density and relative duct area did not differ between healthy trees previously inferred to carry the dominant Cr4 allele and trees that lack it at 22 sites, though some duct traits varied with elevation. MPB resistance does not appear to have played an evolutionary role in contributing to the high frequency of Cr4 in naïve populations, however, MPB may affect the future evolution of resistance to WPBR in the pines where the two pests coincide and WPBR will affect forest recovery after MPB epidemics. MPB-WPBR interactions in a changing climate will affect the future trajectory of limber pine.

Comparisons of three plant species in accumulating polycyclic aromatic hydrocarbons (PAHs) from the atmosphere: a review

Plant leaves play a key role in the accumulation of PAHs, as they are able to capture PAHs from the air. In this paper, the mechanism, including absorption and adsorption, for plants to scavenge PAHs from the air was reviewed. Moreover, the differences of PAHs accumulating capability are mainly compared among three representative plant species, including pine needles, Holm oak leaves, and moss. On the whole, it is shown that oak leaves present the strongest PAHs accumulating capability for total PAHs among three plants species. Oak leaves and pine needles show higher accumulating tendency for light and medium molecular weight PAHs, whereas moss presents stronger accumulating tendency for heavy molecular weight PAHs. Environmental factors (i.e., temperature, seasonality, and photolysis) also account for the process of PAHs transferred from air to plants. With the temperature climbing, the concentration of PAHs in the air will increase. Due to the meteorological conditions and the human activities changed with seasons, it was shown that the PAHs were greatly accumulated in leaf surface in winter than in summer. Photolysis was also able to influence the PAHs on leaf surface, which are significant to this process. In conclusion, oak, pine, and moss can be used to filter PAHs when considering urban landscaping. Besides combining the traditional analytical methods with in situ determination, there might be able to provide a novel method to further study the specific absorption mechanisms. The accumulation of PAHs in crop leaf surface related to the application of surfactants is also worth studying.