Approximately 100mg of tissue was used for RNA extraction using Spectrum Plant Total RNA Kit (Sigma-Aldrich, USA) according to the manufacturer instructions, followed by DNase treatment (On-Column DNase I Digestion Set, Sigma-Aldrich, USA). 100ng of RNA for samples were aliquoted for sequencing MIT BioMicro Center performed RNA-seq library construction and sequencing following NEB Ultra II Directional RNA with Poly(A) Selection (New England Biolabs, Ipswich, MA) 150 + 150 bases pair-end run with 8 + 8 nucleotide indexes