Main conclusion Abscisic acid is involved in the drought response of Ilex paraguariensis. Acclimation includes root growth stimulation, stomatal closure, osmotic adjustment, photoprotection, and regulation of nonstructural carbohydrates and amino acid metabolisms.

Reverse transcription of RNA followed by real-time quantitative PCR (qPCR) is to date the most reliable method for gene expression studies. However, to control the errors introduced along the numerous experimental procedures, it requires a normalization using internal reference genes with stable expression. To address this issue, nine candidate reference genes were investigated in Ilex paraguariensis leaves subjected to water stress. To facilitate the selection, we analysed the real-time qPCR data with three different software programs. The obtained results support the conclusion that RNA polymerase associated protein rtf1 homolog (RTF) combined with any of the following pairs is the most suitable triad of genes to compute a normalization factor: elongation factor 1-alpha + tubulin alpha chain (EF1a + α-Tub), actin + cyclophilin 38 (ACT + CYP38), or cyclophilin 38 + vacuolar protein sorting-associated protein 18 homologs (CYP38 + VPS). Our analysis constitutes the first in-depth study to identify the appropriate reference genes for the quantification of transcription in Ilex paraguariensis leaves during drought and provides essential information for further gene expression studies in this tree species.

The appearance of water stress episodes triggers leaf abscission and decreases Ilex paraguariensis yield. To explore the mechanisms that allow it to overcome dehydration, we investigated how the root gene expression varied between water-stressed and non-stressed plants and how the modulation of gene expression was linked to metabolite composition and physiological status. After water deprivation, 5160 differentially expressed transcripts were obtained through RNA-seq. The functional enrichment of induced transcripts revealed significant transcriptional remodelling of stress-related perception, signalling, transcription, and metabolism. Simultaneously, the induction of the enzyme 9-cis-expoxycarotenoid dioxygenase (NCED) transcripts reflected the central role of the hormone abscisic acid in this response. Consequently, the total content of amino acids and soluble sugars increased, and that of starch decreased. Likewise, osmotic adjustment and radical growth were significantly promoted to preserve cell membranes and water uptake. This study provides a valuable resource for future research to understand the molecular adaptation of I. paraguariensis plants under drought conditions and facilitates the exploration of drought-tolerant candidate genes.