Molecular study of drought response in the Mediterranean conifer Pinus pinaster Ait.: Differential transcriptomic profiling reveals constitutive water deficit-independent drought tolerance mechanisms

Maritime Pine Rootstock Genotype Modulates Gene Expression Associated with Stress Tolerance in Grafted Stems

Climate change-induced hazards, such as drought, threaten forest resilience, particularly in vulnerable regions such as the Mediterranean Basin. Maritime pine (Pinus pinaster Aiton), a model species in Western Europe, plays a crucial role in the Mediterranean forest due to its genetic diversity and ecological plasticity. This study characterizes transcriptional profiles of scion and rootstock stems of four P. pinaster graft combinations grown under well-watered conditions. Our grafting scheme combined drought-sensitive and drought-tolerant genotypes for scions (GAL1056: drought-sensitive scion; and Oria6: drought-tolerant scion) and rootstocks (R1S: drought-sensitive rootstock; and R18T: drought-tolerant rootstock). Transcriptomic analysis revealed expression patterns shaped by genotype provenance and graft combination. The accumulation of differentially expressed genes (DEGs) encoding proteins, involved in defense mechanisms and pathogen recognition, was higher in drought-sensitive scion stems and also increased when grafted onto drought-sensitive rootstocks. DEGs involved in drought tolerance mechanisms were identified in drought-tolerant genotypes as well as in drought-sensitive scions grafted onto drought-tolerant rootstocks, suggesting their establishment prior to drought. These mechanisms were associated with ABA metabolism and signaling. They were also involved in the activation of the ROS-scavenging pathways, which included the regulation of flavonoid and terpenoid metabolisms. Our results reveal DEGs potentially associated with the conifer response to drought and point out differences in drought tolerance strategies. These findings suggest genetic trade-offs between pine growth and defense, which could be relevant in selecting more drought-tolerant Pinus pinaster trees.

Under Stress: Searching for Genes Involved in the Response of Abies pinsapo Boiss to Climate Change

Currently, Mediterranean forests are experiencing the deleterious effects of global warming, which mainly include increased temperatures and decreased precipitation in the region. Relict Abies pinsapo fir forests, endemic in the southern Iberian Peninsula, are especially sensitive to these recent environmental disturbances, and identifying the genes involved in the response of this endangered tree species to climate-driven stresses is of paramount importance for mitigating their effects. Genomic resources for A. pinsapo allow for the analysis of candidate genes reacting to warming and aridity in their natural habitats. Several members of the complex gene families encoding late embryogenesis abundant proteins (LEAs) and heat shock proteins (HSPs) have been found to exhibit differential expression patterns between wet and dry seasons when samples from distinct geographical locations and dissimilar exposures to the effects of climate change were analyzed. The observed changes were more perceptible in the roots of trees, particularly in declining forests distributed at lower altitudes in the more vulnerable mountains. These findings align with previous studies and lay the groundwork for further research on the molecular level. Molecular and genomic approaches offer valuable insights for mitigating climate stress and safeguarding this endangered conifer.

Physiological responses to drought stress of three pine species and comparative transcriptome analysis of Pinus yunnanensis var. pygmaea

Drought stress can significantly affect plant growth, development, and yield. Fewer comparative studies have been conducted between different species of pines, particularly involving Pinus yunnanensis var. pygmaea (P. pygmaea). In this study, the physiological indices, photosynthetic pigment and related antioxidant enzyme changes in needles from P. pygmaea, P. elliottii and P. massoniana under drought at 0, 7, 14, 21, 28 and 35 d, as well as 7 days after rehydration, were measured. The PacBio single-molecule real-time (SMRT) and Illumina RNA sequencing were used to uncover the gene expression differences in P. pygmaea under drought and rehydration conditions. The results showed that the total antioxidant capacity (TAOC) of P. pygmaea was significantly higher than P. massoniana and P. elliottii. TAOC showed a continuous increase trend across all species. Soluble sugar (SS), starch content and non-structural carbohydrate (NSC) of all three pines displayed a "W" pattern, declining initially, increasing, and then decreasing again. P. pygmaea exhibits stronger drought tolerance and greater recovery ability under prolonged drought conditions. Through the PacBio SMRT-seq, a total of 50,979 high-quality transcripts were generated, and 6,521 SSR and 5,561 long non-coding RNAs (LncRNAs) were identified. A total of 2310, 1849, 5271, 5947, 7710, and 6854 differentially expressed genes (DEGs) were identified compared to the control (Pp0D) in six pair-wise comparisons of treatment versus control. bHLH, NAC, ERF, MYB_related, C3H transcription factors (TFs) play an important role in drought tolerance of P. pygmaea. KEGG enrichment analysis and Gene set enrichment analysis (GSEA) analysis showed that P. pygmaea may respond to drought by enhancing metabolic processes such as ABA signaling pathway, alpha-linolenic acid. Weighted gene co-expression network analysis (WGCNA) revealed GST, CAT, LEC14B, SEC23 were associated with antioxidant enzyme activity and TAOC. This study provides a basis for further research on drought tolerance differences among coniferous species.

Warming appears as the main risk of non-adaptedness for western Mediterranean relict fir forests under expected climate change scenarios

Circum-Mediterranean firs are considered among the most drought-sensitive species to climate change. Understanding the genetic basis of trees' adaptive capacity and intra-specific variability to drought avoidance is mandatory to define conservation measures, thus potentially preventing their extinction. We focus here on Abies pinsapo and Abies marocana, both relict tree species, endemic from south Spain and north Morocco, respectively. A total of 607 samples were collected from eight nuclei: six from Spanish fir and two from Moroccan fir. A genotyping by sequencing technique called double digestion restriction site-associated DNA sequencing (ddRAD-seq) was performed to obtain a genetic matrix based on single-nucleotide polymorphisms (SNPs). This matrix was utilized to study the genetic structure of A. pinsapo populations and to carry out selection signature studies. In order to understand how Spanish fir and Moroccan fir cope with climate change, genotype-environment associations (GEAs) were identified. Further, the vulnerability of these species to climate variations was estimated by the risk of non-adaptedness (RONA). The filtering of the de novo assembly of A. pinsapo provided 3,982 SNPs from 504 out of 509 trees sequenced. Principal component analysis (PCA) genetically separated Grazalema from the rest of the Spanish populations. However, FST values showed significant differences among the sampling points. We found 51 loci potentially under selection. Homolog sequences were found for some proteins related to abiotic stress response, such as dehydration-responsive element binding transcription factor, regulation of abscisic acid signaling, and methylation pathway. A total of 15 associations with 11 different loci were observed in the GEA studies, with the maximum temperature of the warmest month being the variable with the highest number of associated loci. This temperature sensitivity was also supported by the risk of non-adaptedness, which yielded a higher risk for both A. pinsapo and A. marocana under the high emission scenario (Representative Concentration Pathway (RCP) 8.5). This study sheds light on the response to climate change of these two endemic species.

Heat-Priming during Somatic Embryogenesis Increased Resilience to Drought Stress in the Generated Maritime Pine (Pinus pinaster) Plants

Drought stress is becoming the most important factor of global warming in forests, hampering the production of reproductive material with improved resilience. Previously, we reported that heat-priming maritime pine (Pinus pinaster) megagametophytes during SE produced epigenetic changes that generated plants better adapted to subsequent heat stress. In this work, we tested, in an experiment performed under greenhouse conditions, whether heat-priming will produce cross-tolerance to mild drought stress (30 days) in 3-year-old priming-derived plants. We found that they maintain constitutive physiological differences as compared to controls, such as higher proline, abscisic acid, starch, and reduced glutathione and total protein contents, as well as higher ΦPSII yield. Primed plants also displayed a constitutive upregulation of the WRKY transcription factor and the Responsive to Dehydration 22 (RD22) genes, as well as of those coding for antioxidant enzymes (APX, SOD, and GST) and for proteins that avoid cell damage (HSP70 and DHNs). Furthermore, osmoprotectants as total soluble sugars and proteins were early accumulated in primed plants during the stress. Prolongated water withdrawal increased ABA accumulation and negatively affected photosynthesis in all plants but primed-derived plants recovered faster than controls. We concluded that high temperature pulses during somatic embryogenesis resulted in transcriptomic and physiological changes in maritime pine plants that can increase their resilience to drought stress, since heat-primed plants exhibit permanent activation of mechanisms for cell protection and overexpression of stress pathways that pre-adapt them to respond more efficiently to soil water deficit.

Effects of drought stress memory on the accumulation of stress-protective compounds in naturally grown pine and spruce

Forest trees are subjected to multiple stressors during their long lifetime and therefore require effective and finely regulated stress-protective systems. Stressors can induce protective systems either directly or with the involvement of stress memory mechanisms. Stress memory has only begun to be uncovered in model plants and is unexplored in coniferous species. Therefore, we studied the possible role of stress memory in the regulation of the accumulation of stress-protective compounds (heat shock proteins, dehydrins, proline) in the needles of naturally grown Scots pine and Norway spruce trees subjected to the subsequent action of long-term (multiyear) and short-term (seasonal) water shortages. Although the water deficit was relatively mild, it significantly influenced the pattern of expression of stress memory-related heat shock factor (HSF) and SWI/SNF genes, indicating the formation of stress memory in both species. In spruce, dehydrin accumulation was increased by water shortage in a manner compatible with Type II stress memory. The accumulation of HSP40 in spruce needles was positively influenced by long-term water shortage, but this increase was unlikely to be of biological importance due to the concomitant decrease in HSP70, HSP90 and HSP101 accumulation. Finally, proline accumulation was negatively influenced by short-term water deficit in spruce. In pine, no one protective compound accumulated in response to water stress. Taken together, the results indicate that the accumulation of stress-protective compounds was generally independent of stress memory effects both in pine and in spruce.

Differences in heat tolerance, water use efficiency and growth among Douglas-fir families and varieties evidenced by GWAS and common garden studies

Severe and frequent heat and drought events challenge the survival and development of long-generation trees. In this study, we investigated the genomic basis of heat tolerance, water use efficiency and growth by performing genome-wide association studies in coastal Douglas-fir (Pseudotsuga menziesii) and intervarietal (menziesii × glauca) hybrid seedlings. GWAS results identified 32 candidate genes involved in primary and secondary metabolism, abiotic stress and signaling, among other functions. Water use efficiency (inferred from carbon isotope discrimination), photosynthetic capacity (inferred from %N), height and heat tolerance (inferred from electrolyte leakage in a heat stress experiment) were significantly different among Douglas-fir families and varieties. High-elevation seed sources had increased water use efficiency, which could be a result of higher photosynthetic capacity. Similarly, families with greater heat tolerance also had higher water use efficiency and slower growth, suggesting a conservative growth strategy. Intervarietal hybrids showed increased heat tolerance (lower electrolyte leakage at 50 and 55 °C) and higher water use efficiency compared with coastal families, suggesting that hybridization might be a source of pre-adapted alleles to warming climates and should be considered for large-scale reforestation projects under increasingly arid conditions.

De novo transcriptome sequencing and gene co-expression reveal a genomic basis for drought sensitivity and evidence of a rapid local adaptation on Atlas cedar (Cedrus atlantica)

Introduction

Understanding the adaptive capacity to current climate change of drought-sensitive tree species is mandatory, given their limited prospect of migration and adaptation as long-lived, sessile organisms. Knowledge about the molecular and eco-physiological mechanisms that control drought resilience is thus key, since water shortage appears as one of the main abiotic factors threatening forests ecosystems. However, our current background is scarce, especially in conifers, due to their huge and complex genomes.

Methods

Here we investigated the eco-physiological and transcriptomic basis of drought response of the climate change-threatened conifer Cedrus atlantica. We studied C. atlantica seedlings from two locations with contrasting drought conditions to investigate a local adaptation. Seedlings were subjected to experimental drought conditions, and were monitored at immediate (24 hours) and extended (20 days) times. In addition, post-drought recovery was investigated, depicting two contrasting responses in both locations (drought resilient and non-resilient). Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of drought resilience and investigate a rapid local adaptation of C. atlantica.

Results

De novo transcriptome assembly was performed for the first time in this species, providing differences in gene expression between the immediate and extended treatments, as well as among the post-drought recovery phenotypes. Weighted gene co-expression network analysis showed a regulation of stomatal closing and photosynthetic activity during the immediate drought, consistent with an isohydric dynamic. During the extended drought, growth and flavonoid biosynthesis inhibition mechanisms prevailed, probably to increase root-to-shoot ratio and to limit the energy-intensive biosynthesis of secondary metabolites. Drought sensitive individuals failed in metabolism and photosynthesis regulation under drought stress, and in limiting secondary metabolite production. Moreover, genomic differences (SNPs) were found between drought resilient and sensitive seedlings, and between the two studied locations, which were mostly related to transposable elements.

Discussion

This work provides novel insights into the transcriptomic basis of drought response of C. atlantica, a set of candidate genes mechanistically involved in its drought sensitivity and evidence of a rapid local adaptation. Our results may help guide conservation programs for this threatened conifer, contribute to advance drought-resilience research and shed light on trees' adaptive potential to current climate change.

Correlating the above- and belowground genotype of Pinus pinaster trees and rhizosphere bacterial communities under drought conditions

Increasing temperatures along with severe droughts are factors that may jeopardize the survival of the forests in the Mediterranean basin. In this region, Pinus pinaster is a common conifer species, that has been used as a model species in evolutionary studies due to its adaptive response to changing environments. Although its drought tolerance mechanisms are already known, knowledge about the dynamics of its root microbiota is still scarce. We aimed to decipher the structural (bacterial abundance), compositional, functional and associative changes of the P. pinaster rhizosphere bacterial communities in spring and summer, at DNA and RNA level (environmental DNA, live and dead cells, and those synthesizing proteins). A fundamental aspect of root microbiome-based approaches is to guarantee the correct origin of the samples. Thus, we assessed the genotype of host needles and roots from which rhizosphere samples were obtained. For more than 50% of the selected trees, genotype discrepancies were found and in three cases the plant species could not be determined. Rhizosphere bacterial communities were homogeneous with respect to diversity and structural levels regardless of the host genotype in both seasons. Nonetheless, significant changes were seen in the taxonomic profiles depending on the season. Seasonal changes were also evident in the bacterial co-occurrence patterns, both in DNA and RNA libraries. While spring communities switched to more complex networks, summer populations resulted in more compartmentalized networks, suggesting that these communities were facing a disturbance. These results may mirror the future status of bacterial communities in a context of climate change. A keystone hub was ascribed to the genus Phenylobacterium in the functional network calculated for summer. Overall, it is important to validate the origin and identity of plant samples in any plant-microbiota study so that more reliable ecological analyses are performed.

miR160 Interacts in vivo With Pinus pinaster AUXIN RESPONSE FACTOR 18 Target Site and Negatively Regulates Its Expression During Conifer Somatic Embryo Development

MicroRNAs (miRNAs) are key regulators of several plant developmental processes including embryogenesis. Most miRNA families are conserved across major groups of plant species, but their regulatory roles have been studied mainly in model species like Arabidopsis and other angiosperms. In gymnosperms, miRNA-dependent regulation has been less studied since functional approaches in these species are often difficult to establish. Given the fundamental roles of auxin signaling in somatic embryogenesis (SE) induction and embryo development, we investigated a previously predicted interaction between miR160 and a putative target encoding AUXIN RESPONSE FACTOR 18 in Pinus pinaster (PpARF18) embryonic tissues. Phylogenetic analysis of AUXIN RESPONSE FACTOR 18 (ARF18) from Pinus pinaster and Picea abies, used here as a model system of conifer embryogenesis, showed their close relatedness to AUXIN RESPONSE FACTOR (ARF) genes known to be targeted by miR160 in other species, including Arabidopsis ARF10 and ARF16. By using a luciferase (LUC) reporter system for miRNA activity in Arabidopsis protoplasts, we have confirmed that P. pinaster miR160 (ppi-miR160) interacts in vivo with PpARF18 target site. When the primary miR160 from P. pinaster was overexpressed in protoplasts under non-limiting levels of ARGONAUTE1, a significant increase of miR160 target cleavage activity was observed. In contrast, co-expression of the primary miRNA and the target mimic MIM160 led to a decrease of miR160 activity. Our results further support that this interaction is functional during consecutive stages of SE in the conifer model P. abies. Expression analyses conducted in five stages of development, from proembryogenic masses (PEMs) to the mature embryo, show that conifer ARF18 is negatively regulated by miR160 toward the fully developed mature embryo when miR160 reached its highest expression level. This study reports the first in vivo validation of a predicted target site of a conifer miRNA supporting the conservation of miR160 interaction with ARF targets in gymnosperms. The approach used here should be useful for future characterization of miRNA functions in conifer embryogenesis.

Candidate regulators and target genes of drought stress in needles and roots of Norway spruce

Drought stress impacts seedling establishment, survival and whole-plant productivity. Molecular responses to drought stress have been most extensively studied in herbaceous species, mostly considering only aboveground tissues. Coniferous tree species dominate boreal forests, which are predicted to be exposed to more frequent and acute drought as a result of ongoing climate change. The associated impact at all stages of the forest tree life cycle is expected to have large-scale ecological and economic impacts. However, the molecular response to drought has not been comprehensively profiled for coniferous species. We assayed the physiological and transcriptional response of Picea abies (L.) H. Karst seedling needles and roots after exposure to mild and severe drought. Shoots and needles showed an extensive reversible plasticity for physiological measures indicative of drought-response mechanisms, including changes in stomatal conductance (gs), shoot water potential and abscisic acid (ABA). In both tissues, the most commonly observed expression profiles in response to drought were highly correlated with the ABA levels. Still, root and needle transcriptional responses contrasted, with extensive root-specific down-regulation of growth. Comparison between previously characterized Arabidopsis thaliana L. drought-response genes and P. abies revealed both conservation and divergence of transcriptional response to drought. In P. abies, transcription factors belonging to the ABA responsive element(ABRE) binding/ABRE binding factors ABA-dependent pathway had a more limited role. These results highlight the importance of profiling both above- and belowground tissues, and provide a comprehensive framework to advance the understanding of the drought response of P. abies. The results demonstrate that a short-term, severe drought induces severe physiological responses coupled to extensive transcriptome modulation and highlight the susceptibility of Norway spruce seedlings to such drought events.

Rootstock effects on scion gene expression in maritime pine

Pines are the dominant conifers in Mediterranean forests. As long-lived sessile organisms that seasonally have to cope with drought periods, they have developed a variety of adaptive responses. However, during last decades, highly intense and long-lasting drought events could have contributed to decay and mortality of the most susceptible trees. Among conifer species, Pinus pinaster Ait. shows remarkable ability to adapt to different environments. Previous molecular analysis of a full-sib family designed to study drought response led us to find active transcriptional activity of stress-responding genes even without water deprivation in tolerant genotypes. To improve our knowledge about communication between above- and below-ground organs of maritime pine, we have analyzed four graft-type constructions using two siblings as rootstocks and their progenitors, Gal 1056 and Oria 6, as scions. Transcriptomic profiles of needles from both scions were modified by the rootstock they were grafted on. However, the most significant differential gene expression was observed in drought-sensitive Gal 1056, while in drought-tolerant Oria 6, differential gene expression was very much lower. Furthermore, both scions grafted onto drought-tolerant rootstocks showed activation of genes involved in tolerance to abiotic stress, and is most remarkable in Oria 6 grafts where higher accumulation of transcripts involved in phytohormone action, transcriptional regulation, photosynthesis and signaling has been found. Additionally, processes, such as those related to secondary metabolism, were mainly associated with the scion genotype. This study provides pioneering information about rootstock effects on scion gene expression in conifers.

Comprehensive analysis of the isomiRome in the vegetative organs of the conifer Pinus pinaster under contrasting water availability

An increasing number of microRNAs (miRNAs) and miRNA-related sequences produced during miRNA biogenesis, comprising the isomiRome, have been recently highlighted in different species as critical mediators of environmental stress responses. Conifers have some of the largest known genomes but an extensive characterization of the isomiRome from any conifer species has been lacking. We provide here a comprehensive overview of the Pinus pinaster isomiRome expressed in roots, stem and needles under well-watered and drought conditions. From the 13,441 unique small RNA sequences identified, 2,980 were annotated as canonical miRNAs or miRNA* and the remaining were classified as isomiRNA or miRNA-like sequences. A survey of their expression patterns highlighted roots as the most responsive organ under drought, where specific sequences of which a 24-nt novel miRNA stood out, were strongly down-regulated. Given the putative roles of the miRNA-targeted transcripts validated specifically in root tissues, some of the miRNAs, conserved and novel, are shortlisted as potential regulators of drought response. These results provide a valuable resource for comparative studies between gymnosperms and angiosperms. Furthermore, it evidences high transferability of the isomiRome between pine species being a useful basis for further molecular regulation and physiological studies, and especially those focused on adaptation to drought conditions.

Short De-Etiolation Increases the Rooting of VC801 Avocado Rootstock

Dark-grown (etiolated) branches of many recalcitrant plant species root better than their green counterparts. Here it was hypothesized that changes in cell-wall properties and hormones occurring during etiolation contribute to rooting efficiency. Measurements of chlorophyll, carbohydrate and auxin contents, as well as tissue compression, histological analysis and gene-expression profiles were determined in etiolated and de-etiolated branches of the avocado rootstock VC801. Differences in chlorophyll content and tissue rigidity, and changes in xyloglucan and pectin in cambium and parenchyma cells were found. Interestingly, lignin and sugar contents were similar, suggesting that de-etiolated branches resemble the etiolated ones in this respect. Surprisingly, the branches that underwent short de-etiolation rooted better than the etiolated ones, and only a slight difference in IAA content between the two was observed. Gene-expression profiles revealed an increase in ethylene-responsive transcripts in the etiolated branches, which correlated with enrichment in xyloglucan hydrolases. In contrast, transcripts encoding pectin methylesterase and pectolyases were enriched in the de-etiolated branches. Taken together, it seems that the short de-etiolation period led to fine tuning of the conditions favoring adventitious root formation in terms of auxin-ethylene balance and cell-wall properties.