Were identified in Rt vs. St, such as 1594 (66.0 ) up-regulated genes and 820 (34.0 ) down-regulated genes, as well as the log2 fold-change of most DEGs was roughly + 1 to + five. In St vs. Sck and Rt vs. Rck, 2286 and 1068 DEGs were detected, respectively. From the 2286 DEGs within the S line, 245 (10.7 ) had been up-regulated and 2041 (89.3 ) were down-regulated, as well as the log2 fold-change of most DEGs ranged from – five to – 1. The 1068 DEGs from the R line integrated 458 (42.9 ) up-regulated genes and 610 (57.1 ) down-regulated genes. The log2 fold-change was amongst – two and three.Fig. 2 FPKM density distribution of genes in the four simplesWang et al. BMC Genomics(2021) 22:Page 4 ofFig. three Venn diagram of the number of DEGs detected in 4 simples. a. Venn diagram indicated the number of up-regulated DEGs. b. Venn diagram indicated the number of down-regulated DEGsEnrichment analysis of DEGs in Rt vs. St, St vs. Sck and Rt vs. RckThe DEGs in Rt vs. St, St vs. Sck and Rt vs. Rck were annotated into 19, 17 and 14 significant GO terms, respectively (Fig. 5). Below biological processes, oxidationreduction reactions were MAP3K8 Synonyms overrepresented in Rt vs. St, St vs. Sck and Rt vs. Rck. DEGs inside the S and R lines had been annotated for responses to oxidative tension. Under cellular components, ubiquitin ligase complex, extracellular region, and apoplast had been HDAC5 Gene ID essentially the most abundant terms in Rt vs. St; and DEGs in the S and R lines were mainlyannotated to the extracellular region and membranes, respectively. As for molecular functions, the DEGs within the 3 groups have been primarily associated with oxidoreductase activity. Furthermore, DEGs in Rt vs. St have been also involved in transcriptional regulation and DNA binding, and DEGs in the S and R lines participated in catalytic activity. KEGG enrichment was completed to identify in which metabolic pathways the DEGs were involved. As shown in Table 1, the DEGs in Rt vs. St were considerably enriched in phenylpropanoid biosynthesis, cysteine andFig. 4 log2fold transform within the DEGs detected in Rck VS Sck, Rt VS St, St VS Sck and Rt VS Rck. a. Variety of genes using a log2fold adjust -5. b. Variety of genes with -5 log2fold alter -3; c. Quantity of genes with -3 log2fold adjust -2. d. Variety of genes with -2 log2fold transform -1. e. Quantity of genes with 1 log2fold change three; f. Variety of genes with three log2fold adjust 5; g. Number of genes with log2fold changeWang et al. BMC Genomics(2021) 22:Web page five ofFig. five GO classification of DEGs. a. GO classification of DEGs in Rt VS St. b. GO classification of DEGs in St VS Sck. c. GO classification of DEGs in Rt VS Rck. BP: biological method; MF: molecular function; CC: cellular component. The x-axis represents essentially the most abundant categories of every group, and the y-axis represents the number of the total genes in every single categorymethionine metabolism, plant-pathogen interaction, MAPK signaling, alpha-linolenic acid metabolism, and linoleic acid metabolism. The DEGs within the S and R lines were considerably enriched in 18 and 9 metabolic pathways, respectively and five pathways had been shared by both S and R lines, like phenylpropanoid biosynthesis, alpha-linolenic acid metabolism, tyrosine metabolism, plant hormone signal transduction, cysteine, and methionine metabolism. There had been 13 unique pathways inside the S line, which includes plant-pathogen interactions, glucosinolate biosynthesis, and MAPK signaling, although 4 exclusive pathways such as valine, leucine and isoleucine degradation have been located within the R line.Functional class.