SRX5796515 GSM3753196 SRP195576 PRJNA541294 SRS4727715 SAMN11581951 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753196 GSM3753196 GEO Accession GSM3753196 SRX5796516 GSM3753197 SRP195576 PRJNA541294 SRS4727716 SAMN11581962 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753197 GSM3753197 GEO Accession GSM3753197 SRX5796517 GSM3753198 SRP195576 PRJNA541294 SRS4727717 SAMN11581961 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753198 GSM3753198 GEO Accession GSM3753198 SRX5796518 GSM3753199 SRP195576 PRJNA541294 SRS4727718 SAMN11581960 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753199 GSM3753199 GEO Accession GSM3753199 SRX5796519 GSM3753200 SRP195576 PRJNA541294 SRS4727719 SAMN11581959 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753200 GSM3753200 GEO Accession GSM3753200 SRX5796520 GSM3753201 SRP195576 PRJNA541294 SRS4727720 SAMN11581958 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753201 GSM3753201 GEO Accession GSM3753201 SRX5796521 GSM3753202 SRP195576 PRJNA541294 SRS4727721 SAMN11581957 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753202 GSM3753202 GEO Accession GSM3753202 SRX5796522 GSM3753203 SRP195576 PRJNA541294 SRS4727722 SAMN11581956 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753203 GSM3753203 GEO Accession GSM3753203 SRX5796523 GSM3753204 SRP195576 PRJNA541294 SRS4727723 SAMN11581955 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753204 GSM3753204 GEO Accession GSM3753204 SRX5796524 GSM3753205 SRP195576 PRJNA541294 SRS4727724 SAMN11581954 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753205 GSM3753205 GEO Accession GSM3753205 SRX5796525 GSM3753206 SRP195576 PRJNA541294 SRS4727725 SAMN11581953 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753206 GSM3753206 GEO Accession GSM3753206 SRX5796526 GSM3753195 SRP195576 PRJNA541294 SRS4727726 SAMN11581952 RNA-Seq TRANSCRIPTOMIC cDNA Total RNA was extracted from frozen samples using the Qiagen RNA plant mini-kit (Qiagen, Hilden, Germany) following the manufacturer's instructions. The quality and quantity of RNA were determined by Nanodrop and Qubit, respectively, and Agilent 2100 was used to evaluate the RNA Integrity. High quality RNA samples with Integrity Number >8 were used for constructing cDNA libraries and sequencing using Illumina HiSeq 2000 platform (Beijing Novogene Bioinformatics Technology Co. Ltd., Beijing, China). A total amount of 3μg RNA per sample was used as input material for the RNA sample preparations. Sequencing libraries were generated using Illumina TruseqTM RNA sample prep Kit for Illumina® (NEB, USA) following manufacturer's recommendations and index codes were added to attribute sequences to each sample. Briefly, mRNA was purified from total RNA using poly-T oligo-attached magnetic beads. Fragmentation was carried out using divalent cations under elevated temperature in NEBNext First Strand Synthesis Reaction Buffer (5X). First strand cDNA was synthesized using random hexamer primer and M-MuLV Reverse Transcriptase. Second strand cDNA synthesis was subsequently performed using DNA Polymerase I and RNase H. Remaining overhangs were converted into blunt ends via exonuclease/polymerase activities. After adenylation of 3' ends of DNA fragments, NEBNext Adaptor with hairpin loop structure were ligated to prepare for hybridization. In order to select cDNA fragments of preferentially 150~200 bp in length, the library fragments were purified with AMPure XP system (Beckman Coulter, Beverly, USA). Then 3 μl USER Enzyme (NEB, USA) was used with size-selected, adaptor-ligated cDNA at 37oC for 15 min followed by 5 min at 95 oC before PCR. Then PCR was performed with Phusion High-Fidelity DNA polymerase, Universal PCR primers and Index (X) Primer. At last, PCR products were purified (AMPure XP system) and library quality was assessed on the Agilent Bioanalyzer 2100 system. Illumina HiSeq 2000 gds 303753195 GSM3753195 GEO Accession GSM3753195