4.7. Expression Analysis of BnaFBA Genes

The publically available RNA-seq dataset of 12 different tissues (sepal, pistil, stamen, ovule, pericarp, blossomy pistil, wilting pistil, root, flower, leaf, silique wall, and stem) collected from different stages of the growth (BioProject ID: PRJNA394926) [36], RNA-seq dataset of leaf and root tissues under drought stress (BioProject ID: PRJNA256233) [74], RNA-seq dataset of seeds across four phases of the development (BioProject ID: PRJNA311067) [56], RNA-seq dataset of stems and leaves after Sclerotinia sclerotiorum infection (BioProject ID: PRJNA321917) [75], and time-series RNA-seq dataset of roots under a synthetic analog of strigolactones (rac-GR24) treatments (BioProject ID: PRJNA484313) [76] were downloaded from the NCBI SRA database, and further used as main sources to perform gene expression profiling of BnaFBA genes in B. napus. The transcriptome reads were mapped to B. napus var. Darmor-bzh reference genome using HISAT2 (version 2.1.0, Baltimore, MD, USA) with the default settings [77]. The read counts per gene were generated by featureCounts [78]. Fragments per kilobase of exon per million fragments mapped (FPKM) was used for the quantification of gene expression. The clustered heatmaps were visualized with expression levels (log2) of BnaFBA genes by R software using the pheatmap function package (https://cran.r-project.org/web/packages/pheatmap/).