Transcriptome datasets of β-Aminobutyric acid (BABA)-primed mono- and dicotyledonous plants, Hordeum vulgare and Arabidopsis thaliana

The non-protein amino acid β-Aminobutyric acid (BABA) may trigger the immune responses of plants to various biotic and abiotic stresses leading to a long-term resistance (primed state). We present RNA-seq datasets of BABA - primed mono- and dicotyledonous plants, such as Arabidopsis thaliana and Hordeum vulgare. Illumina NextSeq550 sequencing were carried out after 72 h of BABA exposure. 87 bp long sequence reads were preprocessed of treated and control samples and deposited in the NCBI SRA database. Transcriptome datasets were de novo assembled of each species and deposited in the NCBI TSA database. These SRA and TSA depositions are under the Bioproject accession: PRJNA791573. Pairwise differential expression with enrichment analyses were performed and the most specific DEGs were determined and annotated in both plants.


β-Aminobutyric acid
Arabidopsis thaliana Hordeum vulgare Transcriptome Illumina sequencing a b s t r a c t The non-protein amino acid β-Aminobutyric acid (BABA) may trigger the immune responses of plants to various biotic and abiotic stresses leading to a long-term resistance (primed state). We present RNA-seq datasets of BABA -primed monoand dicotyledonous plants, such as Arabidopsis thaliana and Hordeum vulgare . Illumina NextSeq550 sequencing were carried out after 72 h of BABA exposure. 87 bp long sequence reads were preprocessed of treated and control samples and deposited in the NCBI SRA database. Transcriptome datasets were de novo assembled of each species and deposited in the NCBI TSA database. These SRA and TSA depositions are under the Bioproject accession: PRJNA791573. Pairwise differential expression with enrichment analyses were performed and the most specific DEGs were determined and annotated in both plants.

Value of the Data
• These data contribute to the knowledge of genetic background of pathogen-free immune priming of plants. Immune response stimulation may differ in mono-and dicotyledonous plants, the presented data represent information for both classes. • There is an emerging role of sustainable plant protection which may benefit from these data.
The effect of priming-active elicitor, BABA on protection and signalling pathways has already been demonstrated, so its potential use in the agriculture is well-founded and highlighted. • Illumina GEx sequencing represents the whole transcriptomic gene expression profiling.
There is a rare data on the direct exposure and comparison of BABA for taxonomically distant monocotyledonous and dicotyledonous plants in an experiment. Our datasets may base the understanding the underlying differences of plant physiological processes helping fundamental and applied research as well.

Data Description
Plant defence mechanisms against pathogens can be triggered by various inducers, of which BABA has been shown to be highly effective against viruses, bacteria, fungi, oomycetes, nematodes arthropods and abiotic stresses in several studies [1][2][3][4][5]

Plant materials
H. vulgare cv. Nure and A. thaliana cv. Columbia plants were cultured in phytotrons. BABA pretreatment was performed on 14-day-old plants. Fresh leaves of 17-day old plant samples were collected. Samples were stored in DNA/RNA Shield (Zymo research) at -25 °C until sequencing.

Pre-processing and assembly
Reads were pre-processed using Trimmomatic software [6] . During this step adapters and contamination sequences were removed, low quality bases, short and low-quality reads were filtered out. Full-length transcriptome assembly of cleaned and combined read sets ( A. thaliana combined, H. vulgare combined) from shallow RNA-Seq data were performed by using Trinity and Bowtie2 [7 , 8] . Output statistics of transcriptomes are summarized in Table 1 .

Gene level quantification
To estimate gene expression from RNA-sequencing CountTable was created. To count how many reads map to each feature of interest (genes) each sample reads were aligned to the combined transcriptomes ( Table 1 ). Count Table creation was performed with OmixBox.BioBam ( https://www.biobam.com/omicsbox/ ) using the HTseq package [9] and Bowtie2 [8] . Based on the data of CountTable DEGs were determined as response to BABA treatment with both species (Supplemental Tables 1 and 2).

Pairwise differential expression analysis
Numerical analysis of DEGs in a pairwise comparison of two different experimental conditions was carried out using OmixBox.BioBam. The used application is based on the RSEM and edgeR program implementing quantitative statistical methods to evaluate the significance of individual genes between two experimental conditions [10 , 11] . TMM (Weighted trimmed mean of M-values) normalization method was performed. Distribution of DEGs is shown in the Fig. 1 .

Functional annotation
Annotation of pairwise differential expression with most specific DEGs were determined by Fisher's Exact Test. Overexpressed GO categories (biological process, molecular function and cellular component) as the response to BABA treatment are indicated in bar chart of both plants ( Figs. 2 and 3 ). Gene set enrichment analysis with GO names, GO categories and statistics are detailed in Supplemental Table 3 and 4.

Declaration of Competing Interest
The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.