SRX6463491 GSM3953310 SRP215231 SRS5114672 GSM3953310 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from the tissues using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. About 60 mg of tissues were ground into powder by liquid nitrogen in a 2 mL tube, followed by being homogenized for 2 minutes and rested horizontally for 5 minutes. The mix was centrifuged for 5 minutes at 12,000×g at 4°C, then the supernatant was transferred into a new EP tube with 0.3 mL chloroform/isoamyl alcohol (24:1). The mix was shacked vigorously for 15s, and then centrifuged at 12,000×g for 10 minutes at 4°C. After centrifugation, the upper aqueous phase where RNA remained was transferred into a new tube with equal volume of supernatant of isopropyl alcohol, then centrifuged at 13,600 rpm for 20 minutes at 4°C. After deserting the supernatant, the RNA pellet was washed twice with 1 mL 75% ethanol, then the mix was centrifuged at 13,600 rpm for 3 minutes at 4°C to collect residual ethanol, followed by the pellet air dry for 5-10 minutes in the biosafety cabinet. Finally, 25µL~100µL of DEPC-treated water was added to dissolve the RNA. Subsequently, total RNA was qualified and quantified using a Nano Drop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, MA, USA). DNase I was used to digest double-stranded and single-stranded DNA in total RNA, then, magnetic beads were purified to recover the reaction products, RNase H or Ribo-Zero method (human, mouse, plants) (Illumina, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads, afterwards, A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplyfied by PCR, and the products were purified by Ampure XP Beads. Library was validating on the Agilent Technologies 2100 bioanalyzer for quality control. The double stranded PCR products above were heated denatured and circularized by the splint oligo sequence. The single strand circle DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29(Thermo Fisher Scientific, MA, USA) to make DNA nanoball (DNB) which had more than 300 copies of one molecular, DNBs were loaded into the patterned nanoarray and single end 50 bases reads were generated on BGISEQ500 platform (BGI Shenzhen, China). BGISEQ-500 gds 303953310 GSM3953310 GEO Accession GSM3953310 SRX6463492 GSM3953311 SRP215231 SRS5114671 GSM3953311 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from the tissues using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. About 60 mg of tissues were ground into powder by liquid nitrogen in a 2 mL tube, followed by being homogenized for 2 minutes and rested horizontally for 5 minutes. The mix was centrifuged for 5 minutes at 12,000×g at 4°C, then the supernatant was transferred into a new EP tube with 0.3 mL chloroform/isoamyl alcohol (24:1). The mix was shacked vigorously for 15s, and then centrifuged at 12,000×g for 10 minutes at 4°C. After centrifugation, the upper aqueous phase where RNA remained was transferred into a new tube with equal volume of supernatant of isopropyl alcohol, then centrifuged at 13,600 rpm for 20 minutes at 4°C. After deserting the supernatant, the RNA pellet was washed twice with 1 mL 75% ethanol, then the mix was centrifuged at 13,600 rpm for 3 minutes at 4°C to collect residual ethanol, followed by the pellet air dry for 5-10 minutes in the biosafety cabinet. Finally, 25µL~100µL of DEPC-treated water was added to dissolve the RNA. Subsequently, total RNA was qualified and quantified using a Nano Drop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, MA, USA). DNase I was used to digest double-stranded and single-stranded DNA in total RNA, then, magnetic beads were purified to recover the reaction products, RNase H or Ribo-Zero method (human, mouse, plants) (Illumina, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads, afterwards, A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplyfied by PCR, and the products were purified by Ampure XP Beads. Library was validating on the Agilent Technologies 2100 bioanalyzer for quality control. The double stranded PCR products above were heated denatured and circularized by the splint oligo sequence. The single strand circle DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29(Thermo Fisher Scientific, MA, USA) to make DNA nanoball (DNB) which had more than 300 copies of one molecular, DNBs were loaded into the patterned nanoarray and single end 50 bases reads were generated on BGISEQ500 platform (BGI Shenzhen, China). BGISEQ-500 gds 303953311 GSM3953311 GEO Accession GSM3953311 SRX6463493 GSM3953312 SRP215231 SRS5114673 GSM3953312 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from the tissues using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. About 60 mg of tissues were ground into powder by liquid nitrogen in a 2 mL tube, followed by being homogenized for 2 minutes and rested horizontally for 5 minutes. The mix was centrifuged for 5 minutes at 12,000×g at 4°C, then the supernatant was transferred into a new EP tube with 0.3 mL chloroform/isoamyl alcohol (24:1). The mix was shacked vigorously for 15s, and then centrifuged at 12,000×g for 10 minutes at 4°C. After centrifugation, the upper aqueous phase where RNA remained was transferred into a new tube with equal volume of supernatant of isopropyl alcohol, then centrifuged at 13,600 rpm for 20 minutes at 4°C. After deserting the supernatant, the RNA pellet was washed twice with 1 mL 75% ethanol, then the mix was centrifuged at 13,600 rpm for 3 minutes at 4°C to collect residual ethanol, followed by the pellet air dry for 5-10 minutes in the biosafety cabinet. Finally, 25µL~100µL of DEPC-treated water was added to dissolve the RNA. Subsequently, total RNA was qualified and quantified using a Nano Drop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, MA, USA). DNase I was used to digest double-stranded and single-stranded DNA in total RNA, then, magnetic beads were purified to recover the reaction products, RNase H or Ribo-Zero method (human, mouse, plants) (Illumina, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads, afterwards, A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplyfied by PCR, and the products were purified by Ampure XP Beads. Library was validating on the Agilent Technologies 2100 bioanalyzer for quality control. The double stranded PCR products above were heated denatured and circularized by the splint oligo sequence. The single strand circle DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29(Thermo Fisher Scientific, MA, USA) to make DNA nanoball (DNB) which had more than 300 copies of one molecular, DNBs were loaded into the patterned nanoarray and single end 50 bases reads were generated on BGISEQ500 platform (BGI Shenzhen, China). BGISEQ-500 gds 303953312 GSM3953312 GEO Accession GSM3953312 SRX6463494 GSM3953313 SRP215231 SRS5114674 GSM3953313 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from the tissues using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. About 60 mg of tissues were ground into powder by liquid nitrogen in a 2 mL tube, followed by being homogenized for 2 minutes and rested horizontally for 5 minutes. The mix was centrifuged for 5 minutes at 12,000×g at 4°C, then the supernatant was transferred into a new EP tube with 0.3 mL chloroform/isoamyl alcohol (24:1). The mix was shacked vigorously for 15s, and then centrifuged at 12,000×g for 10 minutes at 4°C. After centrifugation, the upper aqueous phase where RNA remained was transferred into a new tube with equal volume of supernatant of isopropyl alcohol, then centrifuged at 13,600 rpm for 20 minutes at 4°C. After deserting the supernatant, the RNA pellet was washed twice with 1 mL 75% ethanol, then the mix was centrifuged at 13,600 rpm for 3 minutes at 4°C to collect residual ethanol, followed by the pellet air dry for 5-10 minutes in the biosafety cabinet. Finally, 25µL~100µL of DEPC-treated water was added to dissolve the RNA. Subsequently, total RNA was qualified and quantified using a Nano Drop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, MA, USA). DNase I was used to digest double-stranded and single-stranded DNA in total RNA, then, magnetic beads were purified to recover the reaction products, RNase H or Ribo-Zero method (human, mouse, plants) (Illumina, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads, afterwards, A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplyfied by PCR, and the products were purified by Ampure XP Beads. Library was validating on the Agilent Technologies 2100 bioanalyzer for quality control. The double stranded PCR products above were heated denatured and circularized by the splint oligo sequence. The single strand circle DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29(Thermo Fisher Scientific, MA, USA) to make DNA nanoball (DNB) which had more than 300 copies of one molecular, DNBs were loaded into the patterned nanoarray and single end 50 bases reads were generated on BGISEQ500 platform (BGI Shenzhen, China). BGISEQ-500 gds 303953313 GSM3953313 GEO Accession GSM3953313 SRX6463495 GSM3953314 SRP215231 SRS5114676 GSM3953314 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from the tissues using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. About 60 mg of tissues were ground into powder by liquid nitrogen in a 2 mL tube, followed by being homogenized for 2 minutes and rested horizontally for 5 minutes. The mix was centrifuged for 5 minutes at 12,000×g at 4°C, then the supernatant was transferred into a new EP tube with 0.3 mL chloroform/isoamyl alcohol (24:1). The mix was shacked vigorously for 15s, and then centrifuged at 12,000×g for 10 minutes at 4°C. After centrifugation, the upper aqueous phase where RNA remained was transferred into a new tube with equal volume of supernatant of isopropyl alcohol, then centrifuged at 13,600 rpm for 20 minutes at 4°C. After deserting the supernatant, the RNA pellet was washed twice with 1 mL 75% ethanol, then the mix was centrifuged at 13,600 rpm for 3 minutes at 4°C to collect residual ethanol, followed by the pellet air dry for 5-10 minutes in the biosafety cabinet. Finally, 25µL~100µL of DEPC-treated water was added to dissolve the RNA. Subsequently, total RNA was qualified and quantified using a Nano Drop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, MA, USA). DNase I was used to digest double-stranded and single-stranded DNA in total RNA, then, magnetic beads were purified to recover the reaction products, RNase H or Ribo-Zero method (human, mouse, plants) (Illumina, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads, afterwards, A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplyfied by PCR, and the products were purified by Ampure XP Beads. Library was validating on the Agilent Technologies 2100 bioanalyzer for quality control. The double stranded PCR products above were heated denatured and circularized by the splint oligo sequence. The single strand circle DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29(Thermo Fisher Scientific, MA, USA) to make DNA nanoball (DNB) which had more than 300 copies of one molecular, DNBs were loaded into the patterned nanoarray and single end 50 bases reads were generated on BGISEQ500 platform (BGI Shenzhen, China). BGISEQ-500 gds 303953314 GSM3953314 GEO Accession GSM3953314 SRX6463496 GSM3953315 SRP215231 SRS5114675 GSM3953315 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from the tissues using Trizol (Invitrogen, Carlsbad, CA, USA) according to manual instruction. About 60 mg of tissues were ground into powder by liquid nitrogen in a 2 mL tube, followed by being homogenized for 2 minutes and rested horizontally for 5 minutes. The mix was centrifuged for 5 minutes at 12,000×g at 4°C, then the supernatant was transferred into a new EP tube with 0.3 mL chloroform/isoamyl alcohol (24:1). The mix was shacked vigorously for 15s, and then centrifuged at 12,000×g for 10 minutes at 4°C. After centrifugation, the upper aqueous phase where RNA remained was transferred into a new tube with equal volume of supernatant of isopropyl alcohol, then centrifuged at 13,600 rpm for 20 minutes at 4°C. After deserting the supernatant, the RNA pellet was washed twice with 1 mL 75% ethanol, then the mix was centrifuged at 13,600 rpm for 3 minutes at 4°C to collect residual ethanol, followed by the pellet air dry for 5-10 minutes in the biosafety cabinet. Finally, 25µL~100µL of DEPC-treated water was added to dissolve the RNA. Subsequently, total RNA was qualified and quantified using a Nano Drop and Agilent 2100 bioanalyzer (Thermo Fisher Scientific, MA, USA). DNase I was used to digest double-stranded and single-stranded DNA in total RNA, then, magnetic beads were purified to recover the reaction products, RNase H or Ribo-Zero method (human, mouse, plants) (Illumina, USA) was used to remove the rRNA. Purified mRNA from previous steps was fragmented into small pieces with fragment buffer at appropriate temperature. Then, First-strand cDNA was generated in First Strand Reaction System by PCR, and the second-strand cDNA was generated as well. The reaction product was purified by magnetic beads, afterwards, A-Tailing Mix and RNA Index Adapters were added by incubating to carry out end repair. The cDNA fragments with adapters were amplyfied by PCR, and the products were purified by Ampure XP Beads. Library was validating on the Agilent Technologies 2100 bioanalyzer for quality control. The double stranded PCR products above were heated denatured and circularized by the splint oligo sequence. The single strand circle DNA (ssCir DNA) was formatted as the final library. The final library was amplified with phi29(Thermo Fisher Scientific, MA, USA) to make DNA nanoball (DNB) which had more than 300 copies of one molecular, DNBs were loaded into the patterned nanoarray and single end 50 bases reads were generated on BGISEQ500 platform (BGI Shenzhen, China). BGISEQ-500 gds 303953315 GSM3953315 GEO Accession GSM3953315