SRX23388622 ACX_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249778 ACX_mRNA ACX_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388623 BD_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249779 BD_mRNA BD_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388624 LDS_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249780 LDS_mRNA LDS_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388625 LSD_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249781 LSD_mRNA LSD_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388626 LTH_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249782 LTH_mRNA LTH_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388627 LWS_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249783 LWS_mRNA LWS_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388628 MG_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249784 MG_mRNA MG_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388629 MLS_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249785 MLS_mRNA MLS_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388630 MYX_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249786 MYX_mRNA MYX_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388631 MZ_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249787 MZ_mRNA MZ_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388632 NYGK_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249788 NYGK_mRNA NYGK_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388633 QNDX_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249789 QNDX_mRNA QNDX_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388634 BGGC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249790 BGGC_mRNA BGGC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388635 SWC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249791 SWC_mRNA SWC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388636 SZX_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249792 SZX_mRNA SZX_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388637 WQC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249793 WQC_mRNA WQC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388638 XMX_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249794 XMX_mRNA XMX_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388639 YBJ_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249795 YBJ_mRNA YBJ_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388640 YH_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249796 YH_mRNA YH_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388641 ACX_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249797 ACX_ncRNA ACX_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388642 CLL_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249798 CLL_ncRNA CLL_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388643 DDC_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249799 DDC_ncRNA DDC_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388644 DQC_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249800 DQC_ncRNA DQC_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388645 CLL_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249801 CLL_mRNA CLL_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388646 GJ_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249802 GJ_ncRNA GJ_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388647 WQC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249803 WQC_WGS WQC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388648 XMX_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249804 XMX_WGS XMX_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388649 YBJ_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249805 YBJ_WGS YBJ_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388650 DXX_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249806 DXX_mRNA DXX_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388651 YH_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249807 YH_WGS YH_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388652 GJ_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249808 GJ_mRNA GJ_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388653 GRC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249809 GRC_mRNA GRC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388654 LBC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249810 LBC_mRNA LBC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388655 GRC_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249811 GRC_ncRNA GRC_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388656 LWS_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249812 LWS_ncRNA LWS_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388657 MLS_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249813 MLS_ncRNA MLS_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388658 YBJ_ncRNA SRP485848 bp0 Total RNA was extracted from the tissue using TRIzol Reagent (Plant RNA Purification Reagent for plant tissue) according the manufacturers instructions (Invitrogen) and genomic DNA was removed using DNase I (TaKara). RNA degradation and contamination was monitored on 1%agarose gels. Then RNA quality was determined by 2100 Bioanalyser (Agilent Technologies) and quantified using the ND-2000 (NanoDrop Technologies). Only high-quality RNA sample (OD260/280=1.8~2.2, OD260/2302.0, RIN8.0, 28S:18S1.0, >1g) was used to construct sequencing library. RNA purification, reverse transcription, library construction and sequencing were performed at Shanghai Majorbio Bio-pharm Biotechnology Co., Ltd. (Shanghai, China) according to the manufacturers instructions (Illumina, San Diego, CA). RNA-seq transcriptome strand library was prepared following TruSeqTM Stranded Total RNA Library Prep Kit from Illumina (San Diego, CA) using 1g of total RNA. Shortly, ribosomal RNA (rRNA) depletion instead of poly(A) purification was performed by Ribo-off rRNA Depletion Kit and then fragmented by fragmentation buffer firstly. Secondly first-stranded cDNA was synthesized with random hexamer primers. Then we removed the RNA template and synthesized a replacement strand, incorporating dUTP in place of dTTP to generate ds cDNA. The incorporation of dUTP quenched the second strand during amplification, because the polymerase did not incorporate past this nucleotide. AMPure XP beads were used to separate the ds cDNA from the second strand reaction mix. A single A nucleotide was added to the 3 ends of these blunt fragments to prevent them from ligating to one another during the adapter ligation reaction. Lastly, multiple indexing adapters were ligated to the ends of the ds cDNA. Libraries were size selected for cDNA target fragments of 200300 bp on 2% Low Range Ultra Agarose followed by PCR amplified using Phusion DNA polymerase (NEB) for 15PCR cycles. After quantified by TBS380, paired-end RNA-seq sequencing library was sequenced with the Illumina Novaseq 6000 sequencer (2150 bp read length). SRS20249814 YBJ_ncRNA YBJ_ncRNA ncRNA-Seq TRANSCRIPTOMIC Inverse rRNA Illumina NovaSeq 6000 SRX23388659 ACX_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249815 ACX_WGS ACX_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388660 BD_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249816 BD_WGS BD_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388661 BGGC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249817 BGGC_WGS BGGC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388662 CLL_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249818 CLL_WGS CLL_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388663 DDC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249819 DDC_WGS DDC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388664 DDC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249820 DDC_mRNA DDC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388665 DQC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249821 DQC_WGS DQC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388666 DXX_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249822 DXX_WGS DXX_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388667 GJ_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249823 GJ_WGS GJ_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388668 GRC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249824 GRC_WGS GRC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388669 LBC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249825 LBC_WGS LBC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388670 LDS_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249826 LDS_WGS LDS_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388671 LSD_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249827 LSD_WGS LSD_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388672 LTH_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249828 LTH_WGS LTH_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388673 LWS_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249829 LWS_WGS LWS_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388674 MG_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249830 MG_WGS MG_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388675 DQC_mRNA SRP485848 bp0 Total RNA was separately extracted from young leaves of each individual using the QIAGEN RNeasy Plant Mini Kit (QIAGEN, Germany). RNA samples were then treated with DNaseI and purified using Dynabeads Oligo (dT)25 (Life Technologies, Carlsbad, CA, USA) and 150-bp PE libraries were constructed following the manufacturers instructions for the NEBNext Ultra TM RNA Library Prep Kit for Illumina (NEB, Ipswich, MA, USA). The generated cDNA libraries were amplified using the TruSeq PE Cluster kit (Illumina, San Diego, CA, USA) and sequenced on the HiSeq6000 platform at Majorbio Bio-Tech (Shanghai, China). SRS20249831 DQC_mRNA DQC_mRNA RNA-Seq TRANSCRIPTOMIC Oligo-dT Illumina NovaSeq 6000 SRX23388676 MLS_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249832 MLS_WGS MLS_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388677 MYX_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249833 MYX_WGS MYX_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388678 MZ_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249834 MZ_WGS MZ_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388679 NYGK_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249835 NYGK_WGS NYGK_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388680 QNDX_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249836 QNDX_WGS QNDX_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388681 SWC_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249837 SWC_WGS SWC_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000 SRX23388682 SZX_WGS SRP485848 bp0 DNA was extracted from the leaf tissue using CTAB method. Then DNA quality and quantity were determined by the ND-2000 (NanoDrop Technologies). Only high-quality DNA sample (OD260/280=1.8~2.0, OD260/2302.0) was used to construct sequencing library. DNA-seq libraries were prepared following TruSeqTM Nano DNA Library Prep Kit from Illumina (San Diego, CA) using 1g of DNA. Shortly, and then fragmented by fragmentation buffer firstly. Secondly, fragmentated DNA was subjected to end-repairing, phosphorylation and A base addition according to Illuminas library construction protocol. After quantified by TBS380, paired-end DNA-seq sequencing library was sequenced with the Illumina NovaSeq 6000 sequencer (PE150). SRS20249838 SZX_WGS SZX_WGS WGS GENOMIC cDNA Illumina NovaSeq 6000