Association Analysis of Transcriptome and Targeted Metabolites Identifies Key Genes Involved in Iris germanica Anthocyanin Biosynthesis.

The anthocyanin biosynthetic pathway is the main pathway regulating floral coloration in Iris germanica , a well-known ornamental plant. We investigated the transcriptome profiles and targeted metabolites to elucidate the relationship between genes and metabolites in anthocyanin biosynthesis in the bitone flower cultivar 'Clarence', which has a deep blue outer perianth and nearly white inner perianth. In this study, delphinidin-, pelargonidin-, and cyanidin-based anthocyanins were detected in the flowers. The content of delphinidin-based anthocyanins increased with the development of the flower. At full bloom (stage 3), delphinidin-based anthocyanins accounted for most of the total anthocyanin metabolites, whereas the content of pelargonidin- and cyanidin-based anthocyanins was relatively low. Based on functional annotations, a number of novel genes in the anthocyanin pathway were identified, which included early biosynthetic genes IgCHS , IgCHI, and IgF3H and late biosynthetic genes Ig F3'5'H , IgANS, and IgDFR . The expression of key structural genes encoding enzymes, such as IgF3H , Ig F3'5'H , IgANS, and IgDFR , was significantly upregulated in the outer perianth compared to the inner perianth. In addition, most structural genes exhibited their highest expression at the half-color stage rather than at the full-bloom stage, which indicates that these genes function ahead of anthocyanins synthesis. Moreover, transcription factors (TFs) of plant R2R3-myeloblastosis (R2R3-MYB) related to the regulation of anthocyanin biosynthesis were identified. Among 56 R2R3-MYB genes, 2 members belonged to subgroup 4, with them regulating the expression of late biosynthetic genes in the anthocyanin biosynthetic pathway, and 4 members belonged to subgroup 7, with them regulating the expression of early biosynthetic genes in the anthocyanin biosynthetic pathway. Quantitative real-time PCR (qRT-PCR) analysis was used to validate the data of RNA sequencing (RNA-Seq). The relative expression profiles of most candidate genes were consistent with the FPKM of RNA-seq. This study identified the key structural genes encoding enzymes and TFs that affect anthocyanin biosynthesis, which provides a basis and reference for the regulation of plant anthocyanin biosynthesis in I. germanica.