SRX7726254 e1_BHI_1 SRP249507 PRJNA600825 For RNA-seq analysis, raw reads were assessed for quality by FASTQC v0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and filtered to remove low-quality reads with Trimmomatic v0.36. Filtered reads were then aligned against to the corresponding reference genome of E. spnifera using STAR v2.5.3a. Read counts for each gene were calculated using FeatureCounts v1.5.2. Differentially expressed genes were detected by R Bioconductor package DESeq2 with the Benjamini-Hochberg adjusted p-value < 0.05 and fold change > 3. SRS6147714 e1BHI1 e1_BHI_1 RNA-Seq TRANSCRIPTOMIC other HiSeq X Ten SRX7726255 e1_BHI_2 SRP249507 PRJNA600825 For RNA-seq analysis, raw reads were assessed for quality by FASTQC v0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and filtered to remove low-quality reads with Trimmomatic v0.36. Filtered reads were then aligned against to the corresponding reference genome of E. spnifera using STAR v2.5.3a. Read counts for each gene were calculated using FeatureCounts v1.5.2. Differentially expressed genes were detected by R Bioconductor package DESeq2 with the Benjamini-Hochberg adjusted p-value < 0.05 and fold change > 3. SRS6147715 e1BHI2 e1_BHI_2 RNA-Seq TRANSCRIPTOMIC other HiSeq X Ten SRX7726256 e1_SDB_1 SRP249507 PRJNA600825 For RNA-seq analysis, raw reads were assessed for quality by FASTQC v0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and filtered to remove low-quality reads with Trimmomatic v0.36. Filtered reads were then aligned against to the corresponding reference genome of E. spnifera using STAR v2.5.3a. Read counts for each gene were calculated using FeatureCounts v1.5.2. Differentially expressed genes were detected by R Bioconductor package DESeq2 with the Benjamini-Hochberg adjusted p-value < 0.05 and fold change > 3. SRS6147716 e1SDB1 e1_SDB_1 RNA-Seq TRANSCRIPTOMIC other HiSeq X Ten SRX7726257 e1_SDB_2 SRP249507 PRJNA600825 For RNA-seq analysis, raw reads were assessed for quality by FASTQC v0.11.5 (http://www.bioinformatics.babraham.ac.uk/projects/fastqc/) and filtered to remove low-quality reads with Trimmomatic v0.36. Filtered reads were then aligned against to the corresponding reference genome of E. spnifera using STAR v2.5.3a. Read counts for each gene were calculated using FeatureCounts v1.5.2. Differentially expressed genes were detected by R Bioconductor package DESeq2 with the Benjamini-Hochberg adjusted p-value < 0.05 and fold change > 3. SRS6147717 e1SDB2 e1_SDB_2 RNA-Seq TRANSCRIPTOMIC other HiSeq X Ten