Transcriptome Analysis of Seedling Development

Total RNA from the root, stem, and leaf of 4-week-old seedlings was extracted with the RNAprep Pure Plant Kit (432; TIANGEN, Sichuan, China) and submitted into ribonucleic acid sequencing using the Illumina HiSeq × Ten system with the paired-end method by Annoroad Gene Technology (Beijing) Co. Ltd. For transcriptome analysis, reference genome and annotation files with a gene model were downloaded from the SOL genomics network (SGN, http://www.sgn.cornell.edu) database directly. Raw reads were trimmed for low quality and length by Fastp v0.20.0 (Chen et al., 2018 (link)). Clean reads were aligned to the reference genome using Hisat v2.0.5 (Pertea et al., 2016 (link)) with the splice-aware method. The gene expression levels from RNA-seq were quantified by feature counts (Liao et al., 2014 (link)) and normalized by fragments per kb per million reads (FPKM) values. Differential expression analysis of two conditions/groups (two biological replicates per condition) was performed using the DESeq2 R package (1.16.1) (Love et al., 2014 (link)). The p-values were adjusted using Benjamini and Hochberg's approach for controlling the false discovery rate. Genes with an adjusted p < 0.05 found by DESeq2 were assigned as differentially expressed. The raw data were deposited in the NCBI SRA database (Accession: PRJNA605912).