Regulatory Networks of Flowering Genes in Angelica sinensis during Vernalization

Angelica sinensis is a low-temperature and long-day perennial herb that has been widely used for cardio-cerebrovascular diseases in recent years. In commercial cultivation, up to 40% of flowering decreases the officinal yield of roots and accumulation of bioactive compounds. Although the regulatory mechanism of flowering genes during the photoperiod has been revealed, the networks during vernalization have not been mapped. Here, transcriptomics profiles of A. sinensis with uncompleted (T1), completed (T2) and avoided vernalization (T3) were performed using RNA-seq, and genes expression was validated with qRT-PCR. A total of 61,241 isoforms were annotated on KEGG, KOG, Nr and Swiss-Prot databases; 4212 and 5301 differentially expressed genes (DEGs) were observed; and 151 and 155 genes involved in flowering were dug out at T2 vs. T1 and T3 vs. T1, respectively. According to functional annotation, 104 co-expressed genes were classified into six categories: FLC expression (22; e.g., VILs, FCA and FLK), sucrose metabolism (12; e.g., TPSs, SUS3 and SPSs), hormone response (18; e.g., GID1B, RAP2s and IAAs), circadian clock (2; i.e., ELF3 and COR27), downstream floral integrators and meristem identity (15; e.g., SOC1, AGL65 and SPLs) and cold response (35; e.g., PYLs, ERFs and CORs). The expression levels of candidate genes were almost consistent with FPKM values and changes in sugar and hormone contents. Based on their functions, four pathways that regulate flowering during vernalization were mapped, including the vernalization pathway, the autonomic pathway, the age pathway and the GA (hormone) pathway. This transcriptomic analysis provides new insights into the gene-regulatory networks of flowering in A. sinensis.