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.