Next Article in Journal
Cryopreservation of 13 Commercial Cannabis sativa Genotypes Using In Vitro Nodal Explants
Next Article in Special Issue
Genomic-Wide Identification and Characterization of the Uridine Diphosphate Glycosyltransferase Family in Eucommia ulmoides Oliver
Previous Article in Journal
Comparative Analysis of Detectors and Feature Descriptors for Multispectral Image Matching in Rice Crops
Previous Article in Special Issue
Genome-Wide Association Study Reveals the Genetic Basis of Chilling Tolerance in Rice at the Reproductive Stage
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Plants
Volume 10
Issue 9
10.3390/plants10091792
plants-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Biosynthetic Pathway of Proanthocyanidins in Major Cash Crops

by
Insu Lim
and
Jungmin Ha
*
Department of Plant Science, Gangneung-Wonju National University, Gangneung 25457, Korea
*
Author to whom correspondence should be addressed.
Plants 2021, 10(9), 1792; https://doi.org/10.3390/plants10091792
Submission received: 6 July 2021 / Revised: 17 August 2021 / Accepted: 18 August 2021 / Published: 28 August 2021
(This article belongs to the Collection Exploration and Application of Useful Agricultural Genes)
Browse Figures
Versions Notes

Abstract

:
Proanthocyanidins (PAs) are a group of oligomers or polymers composed of monomeric flavanols. They offer many benefits for human fitness, such as antioxidant, anticancer, and anti-inflammatory activities. To date, three types of PA have been observed in nature: procyanidins, propelargonidins, and prodelphinidins. These are synthesized as some of the end-products of the flavonoid pathway by different consecutive enzymatic activities, from the same precursor—naringenin. Although the general biosynthetic pathways of PAs have been reported in a few model plant species, little is known about the species-specific pathways in major crops containing different types of PA. In the present study, we identified the species-specific pathways in 10 major crops, based on the presence/absence of flavanol-based intermediates in the metabolic pathway, and found 202 orthologous genes in the reference genomic database of each species, which may encode for key enzymes involved in the biosynthetic pathways of PAs. Parallel enzymatic reactions in the pathway are responsible for the ratio between PAs and anthocyanins, as well as among the three types of PAs. Our study suggests a promising strategy for molecular breeding, to regulate the content of PAs and anthocyanins and improve the nutritional quality of food sources globally.

1. Introduction

Proanthocyanidins (PAs), or condensed tannins, are oligomeric or polymeric end-products of flavonoid metabolism, starting from the central phenylpropanoid pathway [1]. PAs are brown-pigmented and present in the seed coats or seeds, fruits, bark, and leaves of a wide range of plant species, including important cash crops, such as apples, grapes, soybeans, common beans, cereals, and most berries [2]. These phytochemicals are brown-pigmented, increase plant resistance to herbivory, and protect plants from biotic and abiotic stresses, such as pathogens, insect attacks, and ultraviolet (UV)-B radiation [3,4].
In plants, PAs play important roles in resistance to several abiotic and biotic stresses [3,4]; several studies have reported that PAs increase tolerance to severe environmental conditions such as low temperature [5], drought [6], and UV-B radiation [7]. Furthermore, PAs impart astringency and bitterness to young leaves and fruits, deterring herbivory [8]. Accumulation of PAs also enhances tolerance to infection by biotrophic fungi [9,10] and other plant pathogens [7]. More PAs were accumulated after mechanical wounding and attack by herbivores [8].
PAs have attracted much attention because of their biological and therapeutic potential in humans. PAs provide unique flavors and acerbity to many foods and drinks, such as chocolate, fruit juice, tea, and wine [11]. Tea, which is the most commonly consumed beverage worldwide, is rich in PAs [12]. Cacao is a major food source with high PA content in the confectionery industry [13]. Pharmacological studies have shown beneficial effects of PAs in humans, such as antimicrobial, antidiabetic, antiaging, antioxidant, anticancer, and anti-inflammatory effects [14,15]. In addition to bioavailability, PAs have been reported to improve eyesight and neuroprotective functions, and promote flexibility in joints and blood circulation [16]. With these various benefits for human health, they have been considered to be important food-derived bioactive compounds in the pharmacological and cosmetic industries. Barks of common cinnamons are abundant with PAs, used as a fold medicine or supplement [17]. In soybeans, a cultivar containing high PA levels has been used as an ingredient for cosmetic products [18]. In grapes, because PAs are intensively accumulated in their seeds, seed oil products have been used as supplements for health promotion [19].
In the flavonoid biosynthetic pathway, PAs are the end-products of a branch of the anthocyanidin biosynthetic pathway. PAs and anthocyanins are derived from the same precursor—anthocyanidin—and share a common biosynthetic process for the conversion of phenylalanine to anthocyanidin [3]. PAs are a group of oligomers and polymers composed of flavan-3-ols—the most common subclass of flavonoids. Depending on the composition of the monomer precursor, the type of PA varies, including catechin, epicatechin, gallocatechin, epigallocatechin, afzelechin, and epiafzelechin, which are commonly found in the plant kingdom [2,20,21,22].
However, because of their beneficial properties, the accumulation of PAs has become a target of breeding and genetic engineering in a few model species; yet, the key genes involved in the pathway remain unclear in most cash crops [18]. In this study, species-specific biosynthetic pathways of PAs are reported based on metabolic intermediates, and orthologous genes of the key enzymes involved in each pathway are identified in 10 major cash crops.

2. Biosynthetic Pathways of PAs and Anthocyanins

2.1. Chemical Structure of PAs

PAs are oligomers or polymers composed of flavan-3-ol units. The structure of PAs varies depending on the features of flavan-3-ols (as a starter or extension unit), the degree of polymerization [3,21], the position of the interflavan linkage between the monomeric units [4,21,23], and esterification of the 3-hydoxyl group [3,20,24,25]. Among these factors, the nature of flavan-3-ols (stereochemistry at the chiral center and hydroxyl patterns) is a key concept for understanding the structural variability of PAs. Flavan-3-ol units contain C6–C3–C6 basic flavonoid skeletons with two benzene rings, referred to as A and B rings, and one heterocyclic benzopyran ring, referred to as the C ring (Figure 1A) [26]. Flavan-3-ols have two chiral centers—the R and S configurations—on 2,3 of the C ring. Therefore, there are four diastereoisomers of flavan-3-ol with the same hydroxyl pattern: RR, RS, SR, and SS. The most common types of flavan-3-ol that serve as the building blocks of PAs include 2,3-trans-(+)-catechin and 2,3-cis-(−)-epicatechin (epimer) [27].
Another key factor that determines the type of flavan-3-ol is the hydroxyl pattern. At the early stage of the PA biosynthetic pathway, the formation of the 3′-hydroxyl and 3′5′-hydroxyl groups on the B-ring of flavan-3-ols is catalyzed by flavonoid 3′-hydroxylase (F3′H) and flavonoid 3′5′-hydroxylase (F3′5′H), respectively (Figure 1B). Depending on the hydroxyl pattern of flavan-3-ols, PAs can be categorized as procyanidins, propelargonidins, or prodelphinidins, which are homo-oligomeric PAs, with catechin and/or epicatechin as building units, and a 3′4′-dihydroxyl pattern. Propelargonidins and prodelphinidins are oligomers and polymers composed of extension units of afzelechin and epiafzelechin, respectively, with a 4′-hydroxyl pattern, and gallocatechin and epigallocatechin, respectively, with 4′5′6′-trihydroxyl groups [28,29]. Although, theoretically, there can be more types of PA depending on the differences in the hydroxyl pattern on the B-ring, except for the three types of PAs (procyanidins, propelargonidins, and prodelphinidins), no other types have been found in nature; therefore, these will not be discussed further here [30].

2.2. Biosynthetic Pathway of PAs

PAs are synthesized as one of the end-products of the flavonoid pathway, which shares the biosynthetic pathway with anthocyanins, from phenylalanine to leucoanthocyanidins (Figure 1B). According to the reference pathway of phenylpropanoid and flavonoid biosynthesis in the KEGG pathway database (www.genome.jp, accessed on 26 February 2021), the biosynthetic pathway begins with the conversion of phenylalanine to naringenin chalcone through cis-trans formation and cyclization by consecutive reactions mediated by phenylalanine ammonia-lyase, cinnamic acid 4-hydroxylase, and chalcone synthase (CHS). The cyclization of naringenin chalcone catalyzes the formation of naringenin by chalcone isomerase (CHI) [42,43,44,45].
In the second step, the type of flavan-3-ol (based on the B-ring on the hydroxyl pattern) is determined by two enzymes: F3′H and F3′5′H [46]. F3′H and F3′5′H can catalyze the conversion of naringenin to eriodictyol and pentahydroxy flavanone, respectively. Subsequently, naringenin and its two derivatives (eriodictyol and pentahydroxy flavanone) are catalyzed to produce dihydroflavanols by flavanone 3′-hydroxylase (F3H) [47]. DFR is a key enzyme that determines the type of flavan-3-ol by catalysis of three dihydroflavanols—2,3-dihydroquercetin, dihydrokaempferol, and 2,3-dihydromyricetin—into leucoanthocyanidins [48].
In the last steps of the biosynthetic pathway, the end-products are determined by the expressional competition and/or the presence/absence of ANS, ANR, and leucoanthocyanidin reductase (LAR) in each species [49,50,51]. Leucoanthocyanidins are synthesized into anthocyanidins by ANS, or into (2R, 3S)-flavan-3-ols—such as catechin, afzelechin, and gallocatechin—by LAR. Anthocyanidins (cyanidin, pelargonidin, and delphinidin) can be converted into epicatechin, epiafzelechin, and epigallocatechin, respectively, by ANR.
In the investigation of the biosynthetic pathways to produce different types of PAs, the constituents of flavan-3-ols were identified through a literature review to identify species-specific pathways in 10 major crops and the model plant species Arabidopsis thaliana. Catechin and epicatechin are commonly detected in apples, almonds, blueberries, cacao, common beans, grapes, strawberries, soybeans, tea plants, and peanuts [21,22,52,53,54,55,56]. In contrast, A. thaliana has been reported to have only epicatechin [50,57,58,59]. Afzelechin and epiafzelechin are mainly detected in almonds and strawberries [52,60]. Epigallocatechin and gallocatechin are found in grapevines, blueberries, and tea plants [53,54,61]. Other crops—including apples, cacao, common beans, soybeans, and peanuts—have been reported to contain only (epi)catechin or trace amounts of (epi)afzelechin and (epi)gallocatechin [21,22].

2.3. Regulatory Mechanisms of Flavonoids in Model Plants

Two stereoisomers—the flavan-3-ols (+)-catechin (2,3-trans) and (–)-epicatechin (2,3-cis)—are the building blocks of PAs. (–)-Epicatechin serves as a more common extender unit in the plant kingdom [1,3]. These play a critical role in seed dormancy, germination, and longevity [62].
The regulatory mechanisms underlying anthocyanin and PA biosynthesis have been intensively studied in A. thaliana The contents of anthocyanin and PAs have been reported to be responsive to plant hormones [63], temperature [64], UV irradiation [65], and drought [66]. Using A. thaliana, transparent testa (symbolized tt) mutants—TT2, TT8, and TTG1—have been reported to play key roles in regulating PA accumulation in seeds [67]. Flowering locus C and flowering locus T regulate PA biosynthesis to modulate seed dormancy [68]. In vegetative tissues, MBW (MYB-bHLH-WD repeat) complexes have been reported to regulate the anthocyanin biosynthetic pathway: MYB proteins including production of anthocyanin pigment 1 (PAP1)/MYB75, PAP2/MYB90, MYB113, and MYB114; and bHLH proteins including GLABRA3 (GL3), enhancer of GLABRA3 (EGL3), TT8, and WD-repeat protein TTG1 [67,69]. Leucoanthocyanidins—the common precursors of flavonoids in A. thaliana—are synthesized by CHS, CHI, F3H, F3′H, and DFR. Leucoanthocyanidin dioxygenase/anthocyanidin synthase (LDOX/ANS) and BANYULS/anthocyanidin reductase (BAN/ANR), TT19, TT12, and AHA10 are involved in the biosynthesis of PAs and anthocyanins from leucoanthocyanidins [70]. DFR, LDOX, BAN, TT19, TT12, and AHA10 are directly targeted by the MBW complexes, and MYB11, MYB12, and MYB111 transcription factors regulate the transcription of CHS, CHI, and F3H [67,71].
M. truncatula and tobacco have also been well studied for bioengineering of PAs and anthocyanins [72]. When the MtANR gene was overexpressed, the levels of PAs were significantly increased, while the levels of anthocyanins were reduced in M. truncatula [72]. In V. vinifera, the VvLAR1 gene was functionally characterized, and overexpression of VvLAR1 in tobacco resulted in an increase in PA levels [73]. Similarly, the expression of the CHS, CHI, F3H, and DFR genes was reported to be upregulated by the VvMYB5a gene in tobacco [74].
Because the end-products of flavonoids are accumulated in the testa as pigments, and the levels of each flavonoid are significantly associated with seed coat color, studies on flavonoids have been intensively conducted in model legume species, where seed coat color varies and genomic information has been relatively well constructed. In soybeans, the color of seed coats ranges from yellow, green, buff, and brown to black, resulting from a combination of anthocyanin and PA compounds [75]. At least six loci—I, T, Wp, W1, R, and O—have been reported to be involved in this effect, and all of the genes corresponding to each locus have been characterized, including the CHS, F3′H, F3H, F3′5′H, ANR, and R2R3 MYB transcription factor genes [76,77,78,79,80,81]. Measurements of flavonoid compounds in various soybean germplasms have indicated that there are large variations in anthocyanin and PA levels, and that wild soybeans (Glycine soja), in general, have higher PA levels than cultivated soybeans (G. max) [82]. Black and red soybeans have been reported to possess catechin, epicatechin, and cyanidin 3-glucoside, whereas no PAs have been detected in yellow and buff soybeans. Brown soybeans have only PAs, but no anthocyanins [38,75,83,84]. IT109098—a G. max genotype with a brown seed coat—has been reported to have PA levels as high as those of G. soja [18]. It has been reported that the differences in the expression levels of genes encoding delphinidin 3-O-glucoside 2″-O-glucosyltransferases and ANR result in variations in the levels of anthocyanin and PA between the three genotypes with different seed coat colors: cultivated soybean (yellow), landrace (brown), and wild soybean (black).

2.4. Major Cash Crops with High PA Contents

In this study, a total of 10 high-PA-content cash crops with currently available genomic databases were selected: almond (Prunus dulcis), apple (Malus domestica), blueberry (Vaccinium corymbosum), cacao bean (Theobroma cacao), common bean (Phaseolus vulgaris), grape (Vitis vinifera), peanut (Arachis hypogaea), soybean (Glycine max), strawberry (Fragaria × ananassa), and tea tree (Camellia sinensis) (Table 1). These crops play a critical role in the global food market and cosmetic industry; according to statistical data from the Food and Agriculture Organization of the United Nations (www.fao.org, accessed on 5 February 2021), apples, grapes, and soybeans were ranked in the top 20 for commodity production in the USA in 2019. Almonds and blueberries were also produced at capacities of about 193.7 kt and 30.9 kt, respectively, in 2019 in the USA. Peanuts and strawberries are important commercial crops in several countries, including China, India, Indonesia, and the USA. In China, which has the largest market, the trade volumes of peanuts and strawberries were USD 471 million and 75 million, respectively, in 2019. Cacao beans are consumed worldwide as an ingredient in a variety of processed products—such as chocolate, cocoa powder, and cocoa butter—at a capacity of 3991 kt (in 2016, with the latest available data), as reported by Statista (www.statista.com/, accessed on 6 February 2021).
The PA content ranged from 145.0 to 3532.2 mg/100 g among the 10 species (Table 1). Grape seeds have been reported to have the highest level of PAs (3532.2 g/100 g), followed by cacao and common beans (1460.0 g/100 g and 1000.1 g/100 g, respectively) [2,14,34]. The levels of PAs detected in the 10 species can vary depending on the cultivars or measuring conditions used in each study, and the amount of actual absorption by humans can change depending on the manner of intake as food; for example, fresh grapes have much lower levels of PAs (0.05 g/100 g) than grape seeds (3532.2 g/100 g) [2]; in soybeans, seed coats contained significantly higher contents of PAs than the embryos [38]. In common beans, the levels of anthocyanin and PAs vary according to the color features of the bean coats, depending on the genotypes [34,85].

2.5. Identification of Orthologous Genes Involved in PA Biosynthesis in Major Cash Crops

Despite the nutritional importance of flavonoids, key enzymes involved in these pathways have not been identified in many major food crops with high flavonoid contents (Table 1). To characterize the PA biosynthetic pathways in the selected species, the orthologous genes encoding key enzymes involved in all pathways were searched in the latest reference genomic databases of each species (Table 1). N50 and genome coverage are two of the most important factors for evaluating the quality of genome assembly, where N50 is defined as the length of the contig, scaffold, super-scaffold, and pseudomolecule together being shorter than or equal to 50% of the total genome assembly length, and coverage is calculated from the percentage of the total assembly size over the reference genome. In this study, contig N50 varied widely, ranging from 79 kb (strawberry) to 1900 kb (common bean), and the highest, lowest, and mean genome coverage were 102% (blueberry and grape), 80% (common bean and cacao), and 94%, respectively (Table 1).
To identify orthologous genes involved in the PA pathway in the major crops, we collected a list of key enzymes (F3′5′H, F3′H, F3H, DFR, LAR, ANS, and ANR) based on the EC number from the KEGG database. Of the eight species, 92 genes were identified, including two wild relative species of strawberry and peanut (Arachis duranensis: 6, A. thaliana: 6, Fragaria vesca: 4, G. max: 20, M. domestica: 12, P. vulgaris: 12, T. cacao: 9, V. vinifera: 23) (Table 2). The protein sequences were downloaded from the NCBI GenBank database (www.ncbi.nlm.nih.gov, accessed on 4 March 2021) and cross-checked using the reference genomic database of each species (Table 2). We found that the sequences of three genes—XP_014632702.1 (G. max), XP_007138248 (P. vulgaris), and XP_017985307.1 (T. cacao)—were not available in the latest reference annotation data (Supplementary Table S1). In total, 83 genes were identified in the reference genome, and 14 paralogs, which had not been reported to be involved in the pathways in the KEGG database, were additionally detected in two species (A. hypogea: 5 and F. ananassa: 9; Table 2). Using these 97 genes, orthologous genes were identified in the eight crop species, and 105 genes were newly identified as candidate genes involved in the flavonoid biosynthetic pathway from five species (A. hypogea: 1, C. sinensis: 19, F. ananassa: 7, P. dulcis: 4, and V. corymbosum: 74) (Table 2, Supplementary Table S2). The annotation data currently available for the newly identified orthologous and paralogous genes agreed well with the EC numbers from the KEGG database. In C. sinensis, P. dulcis, and V. corymbosum, a total of 97 orthologs were first found in our study, and they have not been reported to be involved in the phenylpropanoid pathway according to the KEGG pathway.

2.6. Phylogenetic Analysis of Orthologs

Using the three orthologous genes (ANR, ANS, and DFR), we conducted a phylogenetic analysis of the 12 species, including 10 cash crops, A. thaliana, and Oryza sativa as an outgroup. As expected, O. sativa was positioned outermost in the gene tree, and A. thaliana was also clearly distinguished from the other 10 major crops. In total, these 12 species could be categorized into 7 orders, and the phylogenetic tree agreed well with the taxonomy (Figure 2).

3. Regulatory Mechanisms for the Biosynthesis of PAs and Anthocyanins

PAs and anthocyanins are the two main classes of flavonoids, both of which have numerous beneficial effects on human health [15,16] and contribute to plant coloration, development, and responses to biotic and abiotic environments [86,87]. Over the last few decades, with advances in technologies to detect secondary metabolites, various types of flavanol-based metabolite have been identified in a large number of plant species [21,22]. Because of their diverse benefits, many studies have been performed to identify key enzymes involved in PAs in model plants [47,49,51,88,89]; however, the biosynthesis leading to the diversity of PAs has not been fully understood in most of the major crops cultivated worldwide.

3.1. Competitions between Parallel Pathways in the Flavonoid Pathway

In the flavonoid biosynthetic pathway, PAs and anthocyanins share a common biosynthetic pathway for the conversion of phenylalanine to anthocyanidin [3]. PAs and anthocyanins are derived from the same precursor—anthocyanidin—by the activity of ANR and anthocyanin 3′-O-beta-glucosyltransferase (3GT), respectively, [51,90]. PA and anthocyanin contents have been reported to be determined by the competition between these parallel pathways (Figure 1). Expression of VvANR increases the levels of PAs in grape fruit [73]. More PAs have been found to accumulate in the leaves of MtANR-overexpressing plants, as compared to those of the wild types of tobacco and M. truncatula [91]. In soybeans, the expression levels of ANR and 3GT have been shown to be associated with the levels of PAs and anthocyanins, respectively [18]. In blackberries, the expression levels of ANS and ANR have been found to correlate with the changing patterns of PA and anthocyanin levels [92]. We identified 71 orthologs of LAR (31), ANS (21), and ANR (19) (Table 3), thus curating a list of candidate genes that may play a critical role in regulating the levels of PAs or anthocyanins in major crop species.
The types of PA are also determined by another competition between parallel pathways. Different flavan-3-ol monomers are required as building blocks to polymerize the three types of PAs. Procyanidins, propelargonidins, and prodelphinidins are composed of (epi)catechin, (epi)afzelechin, and (epi)gallocatechin, respectively (Figure 1B) [93,94,95,96]. These monomeric building blocks are converted from the same precursor—naringenin—by consecutive enzymatic activities. The types of building blocks are determined by the location and number of hydroxyl groups on the B ring of flavan-3-ols. (Leuco)cyanidins and (leuco)delphinidins are catalyzed by eriodictyol and pentahydroxy flavone, respectively, which are hydroxylated by F3′H and F3′5′H from naringenin, respectively (Figure 1B). (Leuco)pelargonidins are catalyzed directly from naringenin, without hydroxylation. We presumed the enzymatic reactions occurring in each target crop based on the presence/absence of their metabolic intermediates, and detected the orthologous genes encoding key enzymes involved in the biosynthetic pathways of PAs.

3.2. Species-Specific Pathways for the Biosynthesis of Catechin and Epicatechin

Among PAs, procyanidin is the most abundant in the plant kingdom, and its monomeric building blocks—catechin and epicatechin—have also been reported as among the most common PAs [93,94]. These are also found in all of our target crops (Table 1), except for catechin, which has not been identified in A. thaliana to date [58,59,97]. This is consistent with our results, which showed that orthologous genes of LAR were not detected in A. thaliana, which catalyzed catechin from leucoanthocyanidins. In the present study, although both catechin and epicatechin were detected in V. vinifera and P. dulcis [21,22,52,53,73], orthologous genes of ANR were not detected in V. vinifera, while F3H and F3′H were not detected in P. dulcis. The complexity of their genome structure and/or the quality of the reference genome assemblies of V. vinifera and P. dulcis may impede the search for orthologous genes based on sequence similarity (Table 1). Alternatively, there might be species-specific pathways for the biosynthesis of catechin and epicatechin, which have not been reported even in model species.
Both procyanidins and prodelphinidins have been reported to accumulate in C. sinensis, V. corymbosum, and V. vinifera [2,53,61,98]. In contrast to Arabidopsis, the presence of LAR and the high activity of F3′5′H in C. sinensis, Medicago truncatula, V. corymbosum, and V. vinifera can promote the formation of dihydromyricetin, resulting in high PA content in leaf tissues [99]. In grapevines, the accumulation of procyanidins and prodelphinidins has been reported to depend on the presence of the two genes VvF3′H and VvF3′5′H [53,98]. Our analysis was also consistent with previous metabolic studies, which identified that most of the 49 orthologous genes of F3′5′H are present in the three species C. sinensis (5), V. corymbosum (37), and V. vinifera (7). In tobacco, the overexpression of F3′5′H increased delphinidin derivatives in transgenic plants, indicating that the existence of F3′5′H might be responsible for the high PA content in plants such as tea plants, blueberries, and grapes (Table 3).
Propelargonidins are the rarest group among PAs, because these are synthesized only when naringenin is not hydroxylated by either F3′H or F3′5′H (Figure 1). Among the 10 species, the orthologous genes of F3′5′H have not been detected in four species—A. hypogea, F. ananassa, M. domestica, and P. dulcis—where prodelphinidins have not been reported [21,22]. In P. dulcis, only propelargonidin and its monomer, afzelechin, have been detected, where both the orthologous genes of F3′H and F3′5′H have not been identified [21,22,52]. In F. ananassa, although orthologs of F3′H have been identified [21,22], both propelargonidins and procyanidins have been identified where the activity of the F3′H enzyme may be relatively low in the competition of parallel pathways.
In phylogenetic analysis, the three species with prodelphinidins—C. sinensis, V. corymbosum, and V. vinifera—were grouped closely, while the two species containing propelargonidins—F. ananassa, and P. dulcis—were also classified into the same clade as the Rosaceae family (Figure 2).

4. Conclusions

Although different types of PAs are present in many cash crops, PAs have been relatively underinvestigated compared to other flavonoids, such as anthocyanins. This study identified species-specific biosynthesis pathways for PAs and a list of responsive orthologs in the reference genomic data of each species. The competition between parallel pathways may represent a significant regulatory mechanism not only for PA content, but also for the types of PAs, depending on the species. Our results will play a role in molecular breeding, to improve the nutritional quality of dietary food sources. In addition to the food market, PAs have been receiving attention in the pharmacological and cosmetic industries, with therapeutic potential for humans. Furthermore, molecular engineering of the parallel pathways could be a prominent approach to regulate the PA and anthocyanin levels, depending on the desired purpose.

5. Materials and Methods

5.1. Sequence Analysis

The protein sequences of each enzyme corresponding to the enzyme commission number (EC number) in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were downloaded from the National Center for Biotechnology Information (NCBI) GenBank database (www.ncbi.nlm.nih.gov, accessed on 4 March 2021). Because these sequences had been uploaded by different research projects and, therefore, had different names even in the same species, these were blasted against the latest reference genomic databases of each species using BLASTP version 2.9.0, downloaded from the NCBI FTP server (ncbi.nlm.nih.gov/blast/executables/, accessed on 4 March 2021). A total of 92 genes from 8 species (A. duranensis, A. thaliana, F. vesca G. max, M. domestica, P. vulgaris, T. cacao, and V. vinifera) were used as queries to search for orthologous genes. Blast hits with percentage identity below 65% were removed (Table 2).

5.2. Phylogenetic Analysis

The protein sequences of three genes—anthocyanidin synthase (ANS), anthocyanidin reductase (ANR), and dihydroflavonol 4-reductase (DFR)—were used for phylogenetic analysis between 12 species, including O. sativa as an outgroup (A. hypogea, A. thaliana, F. ananassa, G. max, C. sinensis, M. domestica, P. vulgaris, T. cacao, V. vinifera, P. dulcis, and V. corymbosum). Orthologous genes were selected based on the ANR and ANS of cacao and DFR of peanuts, which are single genes. The sequences were concatenated and aligned using the MUSCLE algorithm, and conserved regions in the alignments were selected using Gblocks software (version 0.91) [100,101]. The analysis was then performed using alignments of a set of sequences, with the maximum likelihood method and MEGA X software (https://www.megasoftware.net/download_form, accessed on 9 April 2021). The phylogenetic tree was visualized using FigTree software version 1.4.4 (http://tree.bio.ed.ac.uk/software/figtree/, accessed on 14 April 2021) [102].

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/plants10091792/s1. Table S1: List of genes involved in the biosynthetic pathways of PAs; Table S2: Identification of orthologous genes involved in PAs biosynthetic pathways.

Author Contributions

Conceptualization, J.H. and I.L.; methodology, J.H.; software, I.L.; validation, J.H. and I.L.; formal analysis, J.H. and I.L.; investigation, J.H. and I.L.; resources, J.H. and I.L.; data curation, J.H. and I.L.; writing—original draft preparation, I.L.; writing—review and editing, J.H.; visualization, J.H. and I.L.; supervision, J.H.; project administration, J.H.; funding acquisition, J.H. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by a grant from the National Research Foundation of Korea (NRF), funded by the Korean government (MSIT; no. 2021R1C1C1004233).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing not applicable.

Acknowledgments

This work was supported by a grant from the National Research Foundation of Korea (NRF), funded by the Korean government (MSIT; no. 2021R1C1C1004233).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Lepiniec, L.; Debeaujon, I.; Routaboul, J.-M.; Baudry, A.; Pourcel, L.; Nesi, N.; Caboche, M. Genetics and Biochemistry of Seed Flavonoids. Annu. Rev. Plant. Biol. 2006, 57, 405–430. [Google Scholar] [CrossRef]
  2. Hellström, J.K.; Törrönen, A.R.; Mattila, P.H. Proanthocyanidins in Common Food Products of Plant Origin. J. Agric. Food Chem. 2009, 57, 7899–7906. [Google Scholar] [CrossRef]
  3. Dixon, R.A.; Xie, D.; Sharma, S.B. Proanthocyanidins—A Final Frontier in Flavonoid Research? New Phytol. 2005, 165, 9–28. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Xie, D.-Y.; Dixon, R.A. Proanthocyanidin Biosynthesis–Still More Questions than Answers? Phytochemistry 2005, 66, 2127–2144. [Google Scholar] [CrossRef] [Green Version]
  5. An, J.-P.; Li, R.; Qu, F.-J.; You, C.-X.; Wang, X.-F.; Hao, Y.-J. R2R3-MYB Transcription Factor MdMYB23 Is Involved in the Cold Tolerance and Proanthocyanidin Accumulation in Apple. Plant J. 2018, 96, 562–577. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Malisch, C.S.; Salminen, J.-P.; Kölliker, R.; Engström, M.T.; Suter, D.; Studer, B.; Lüscher, A. Drought Effects on Proanthocyanidins in Sainfoin (Onobrychis Viciifolia Scop.) Are Dependent on the Plant’s Ontogenetic Stage. J. Agric. Food Chem. 2016, 64, 9307–9316. [Google Scholar] [CrossRef]
  7. Mellway, R.D.; Tran, L.T.; Prouse, M.B.; Campbell, M.M.; Constabel, C.P. The Wound-, Pathogen-, and Ultraviolet B-Responsive MYB134 Gene Encodes an R2R3 MYB Transcription Factor That Regulates Proanthocyanidin Synthesis in Poplar. Plant Physiol. 2009, 150, 924–941. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Yu, D.; Huang, T.; Tian, B.; Zhan, J. Advances in Biosynthesis and Biological Functions of Proanthocyanidins in Horticultural Plants. Foods 2020, 9, 1774. [Google Scholar] [CrossRef]
  9. Ullah, C.; Unsicker, S.B.; Fellenberg, C.; Constabel, C.P.; Schmidt, A.; Gershenzon, J.; Hammerbacher, A. Flavan-3-Ols Are an Effective Chemical Defense against Rust Infection. Plant Physiol. 2017, 175, 1560–1578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Yuan, L.; Wang, L.; Han, Z.; Jiang, Y.; Zhao, L.; Liu, H.; Yang, L.; Luo, K. Molecular Cloning and Characterization of PtrLAR3, a Gene Encoding Leucoanthocyanidin Reductase from Populus Trichocarpa, and Its Constitutive Expression Enhances Fungal Resistance in Transgenic Plants. J. Exp. Bot. 2012, 63, 2513–2524. [Google Scholar] [CrossRef] [PubMed]
  11. Santos-Buelga, C.; Scalbert, A. Proanthocyanidins and Tannin-like Compounds–Nature, Occurrence, Dietary Intake and Effects on Nutrition and Health. J. Sci. Food Agric. 2000, 80, 1094–1117. [Google Scholar] [CrossRef]
  12. Engelhardt, U.H.; Lakenbrink, C.; Pokorny, O. Proanthocyanidins, Bisflavanols, and Hydrolyzable Tannins in Green and Black Teas. In Nutraceutical Beverages; (ACS Symposium Series); American Chemical Society: Washington, DC, USA, 2003; Volume 871, pp. 254–264. ISBN 978-0-8412-3823-7. [Google Scholar]
  13. Schroeter, H.; Heiss, C.; Balzer, J.; Kleinbongard, P.; Keen, C.L.; Hollenberg, N.K.; Sies, H.; Kwik-Uribe, C.; Schmitz, H.H.; Kelm, M. (–)-Epicatechin Mediates Beneficial Effects of Flavanol-Rich Cocoa on Vascular Function in Humans. Proc. Natl. Acad. Sci. USA 2006, 103, 1024–1029. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Prior, R.L.; Gu, L. Occurrence and Biological Significance of Proanthocyanidins in the American Diet. Phytochemistry 2005, 66, 2264–2280. [Google Scholar] [CrossRef] [PubMed]
  15. Zhao, M.; Yang, B.; Wang, J.; Liu, Y.; Yu, L.; Jiang, Y. Immunomodulatory and Anticancer Activities of Flavonoids Extracted from Litchi (Litchi Chinensis Sonn.) Pericarp. Int. Immunopharmacol. 2007, 7, 162–166. [Google Scholar] [CrossRef] [PubMed]
  16. Rauf, A.; Imran, M.; Abu-Izneid, T.; Patel, S.; Pan, X.; Naz, S.; Sanches Silva, A.; Saeed, F.; Rasul Suleria, H.A. Proanthocyanidins: A Comprehensive Review. Biomed. Pharmacother. 2019, 116, 108999. [Google Scholar] [CrossRef]
  17. Lin, G.-M.; Lin, H.-Y.; Hsu, C.-Y.; Chang, S.-T. Structural Characterization and Bioactivity of Proanthocyanidins from Indigenous Cinnamon (Cinnamomum Osmophloeum). J. Sci. Food Agric. 2016, 96, 4749–4759. [Google Scholar] [CrossRef]
  18. Ha, J.; Kim, M.; Kim, M.Y.; Lee, T.; Yoon, M.Y.; Lee, J.; Lee, Y.-H.; Kang, Y.-G.; Park, J.S.; Lee, J.H.; et al. Transcriptomic Variation in Proanthocyanidin Biosynthesis Pathway Genes in Soybean (Glycine Spp.). J. Sci. Food Agric. 2018, 98, 2138–2146. [Google Scholar] [CrossRef]
  19. Rodríguez-Pérez, C.; García-Villanova, B.; Guerra-Hernández, E.; Verardo, V. Grape Seeds Proanthocyanidins: An Overview of In Vivo Bioactivity in Animal Models. Nutrients 2019, 11, 2435. [Google Scholar] [CrossRef] [Green Version]
  20. Buzzini, P.; Turchetti, B.; Ieri, F.; Goretti, M.; Branda, E.; Mulinacci, N.; Romani, A. Catechins and Proanthocyanidins: Naturally Occurring O-Heterocycleswith Antimicrobial Activity. In Bioactive Heterocycles IV; Khan, M.T.H., Ed.; Topics in Heterocyclic Chemistry; Springer: Berlin/Heidelber, Germany, 2007; pp. 239–263. ISBN 978-3-540-73404-8. [Google Scholar]
  21. Gu, L.; Kelm, M.A.; Hammerstone, J.F.; Beecher, G.; Holden, J.; Haytowitz, D.; Prior, R.L. Screening of Foods Containing Proanthocyanidins and Their Structural Characterization Using LC-MS/MS and Thiolytic Degradation. J. Agric. Food Chem. 2003, 51, 7513–7521. [Google Scholar] [CrossRef]
  22. Gu, L.; Kelm, M.A.; Hammerstone, J.F.; Beecher, G.; Holden, J.; Haytowitz, D.; Gebhardt, S.; Prior, R.L. Concentrations of Proanthocyanidins in Common Foods and Estimations of Normal Consumption. J. Nutr. 2004, 134, 613–617. [Google Scholar] [CrossRef]
  23. Kimura, H.; Ogawa, S.; Akihiro, T.; Yokota, K. Structural Analysis of A-Type or B-Type Highly Polymeric Proanthocyanidins by Thiolytic Degradation and the Implication in Their Inhibitory Effects on Pancreatic Lipase. J. Chromatogr. A 2011, 1218, 7704–7712. [Google Scholar] [CrossRef]
  24. Forrest, G.I.; Bendall, D.S. The Distribution of Polyphenols in the Tea Plant (Camellia Sinensis L.). Biochem. J. 1969, 113, 741–755. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Huo, Q.; Kong, X.; Yang, X.; Wang, Y.; Ma, L.; Luo, M.; Xu, D. Molecular Modification of Proanthocyanidins. Bioengineered 2016, 7, 274–281. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Alianza SIDALC. Available online: http://www.sidalc.net/cgi-bin/wxis.exe/?IsisScript=COLPOS.xis&method=post&formato=2&cantidad=1&expresion=mfn=039519 (accessed on 23 March 2021).
  27. Foo, L.Y.; Porter, L.J. The Phytochemistry of Proanthocyanidin Polymers. Phytochemistry 1980, 19, 1747–1754. [Google Scholar] [CrossRef]
  28. Ferreira, D.; Brandt, E.V.; Coetzee, J.; Malan, E. Condensed Tannins. In Fortschritte der Chemie Organischer Naturstoffe/Progress in the Chemistry of Organic Natural Products; Ayer, W.A., Brandt, E.V., Coetzee, J., Daniewski, W.M., Ferreira, D., Malan, E., Trifonov, L.S., Vidari, G., Eds.; Springer: Vienna, Austria, 1999; pp. 21–67. ISBN 978-3-7091-6366-5. [Google Scholar]
  29. Porter, L.J. Flavans and proanthocyanidins. In The Flavonoids: Advances in Research since 1980; Harborne, J.B., Ed.; Springer: Boston, MA, USA, 1988; pp. 21–62. ISBN 978-1-4899-2913-6. [Google Scholar]
  30. Ferreira, D.; Li, X.-C. Oligomeric Proanthocyanidins: Naturally Occurring O-Heterocycles. Nat. Prod. Rep. 2000, 17, 193–212. [Google Scholar] [CrossRef]
  31. Alioto, T.; Alexiou, K.G.; Bardil, A.; Barteri, F.; Castanera, R.; Cruz, F.; Dhingra, A.; Duval, H.; Fernández I Martí, Á.; Frias, L.; et al. Transposons Played a Major Role in the Diversification between the Closely Related Almond and Peach Genomes: Results from the Almond Genome Sequence. Plant J. Cell Mol. Biol. 2020, 101, 455–472. [Google Scholar] [CrossRef] [Green Version]
  32. Daccord, N.; Celton, J.-M.; Linsmith, G.; Becker, C.; Choisne, N.; Schijlen, E.; van de Geest, H.; Bianco, L.; Micheletti, D.; Velasco, R.; et al. High-Quality de Novo Assembly of the Apple Genome and Methylome Dynamics of Early Fruit Development. Nat. Genet. 2017, 49, 1099–1106. [Google Scholar] [CrossRef]
  33. Colle, M.; Leisner, C.P.; Wai, C.M.; Ou, S.; Bird, K.A.; Wang, J.; Wisecaver, J.H.; Yocca, A.E.; Alger, E.I.; Tang, H.; et al. Haplotype-Phased Genome and Evolution of Phytonutrient Pathways of Tetraploid Blueberry. GigaScience 2019, 8, giz012. [Google Scholar] [CrossRef] [Green Version]
  34. Kan, L.; Nie, S.; Hu, J.; Liu, Z.; Xie, M. Antioxidant Activities and Anthocyanins Composition of Seed Coats from Twenty-Six Kidney Bean Cultivars. J. Funct. Foods 2016, 26, 622–631. [Google Scholar] [CrossRef]
  35. Schmutz, J.; McClean, P.E.; Mamidi, S.; Wu, G.A.; Cannon, S.B.; Grimwood, J.; Jenkins, J.; Shu, S.; Song, Q.; Chavarro, C.; et al. A Reference Genome for Common Bean and Genome-Wide Analysis of Dual Domestications. Nat. Genet. 2014, 46, 707–713. [Google Scholar] [CrossRef] [Green Version]
  36. Jaillon, O.; Aury, J.-M.; Noel, B.; Policriti, A.; Clepet, C.; Casagrande, A.; Choisne, N.; Aubourg, S.; Vitulo, N.; Jubin, C.; et al. The Grapevine Genome Sequence Suggests Ancestral Hexaploidization in Major Angiosperm Phyla. Nature 2007, 449, 463–467. [Google Scholar] [CrossRef]
  37. Zhuang, W.; Chen, H.; Yang, M.; Wang, J.; Pandey, M.K.; Zhang, C.; Chang, W.-C.; Zhang, L.; Zhang, X.; Tang, R.; et al. The Genome of Cultivated Peanut Provides Insight into Legume Karyotypes, Polyploid Evolution and Crop Domestication. Nat. Genet. 2019, 51, 865–876. [Google Scholar] [CrossRef]
  38. Lee, C.C.; Dudonné, S.; Dubé, P.; Desjardins, Y.; Kim, J.H.; Kim, J.S.; Kim, J.-E.; Park, J.H.Y.; Lee, K.W.; Lee, C.Y. Comprehensive Phenolic Composition Analysis and Evaluation of Yak-Kong Soybean (Glycine Max) for the Prevention of Atherosclerosis. Food Chem. 2017, 234, 486–493. [Google Scholar] [CrossRef] [PubMed]
  39. Schmutz, J.; Cannon, S.B.; Schlueter, J.; Ma, J.; Mitros, T.; Nelson, W.; Hyten, D.L.; Song, Q.; Thelen, J.J.; Cheng, J.; et al. Genome Sequence of the Palaeopolyploid Soybean. Nature 2010, 463, 178–183. [Google Scholar] [CrossRef] [Green Version]
  40. Edger, P.P.; Poorten, T.J.; VanBuren, R.; Hardigan, M.A.; Colle, M.; McKain, M.R.; Smith, R.D.; Teresi, S.J.; Nelson, A.D.L.; Wai, C.M.; et al. Origin and Evolution of the Octoploid Strawberry Genome. Nat. Genet. 2019, 51, 541–547. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Xia, E.; Li, F.; Tong, W.; Yang, H.; Wang, S.; Zhao, J.; Liu, C.; Gao, L.; Tai, Y.; She, G.; et al. The Tea Plant Reference Genome and Improved Gene Annotation Using Long-Read and Paired-End Sequencing Data. Sci. Data 2019, 6, 122. [Google Scholar] [CrossRef] [Green Version]
  42. Bell-Lelong, D.A.; Cusumano, J.C.; Meyer, K.; Chapple, C. Cinnamate-4-Hydroxylase Expression in Arabidopsis (Regulation in Response to Development and the Environment). Plant Physiol. 1997, 113, 729–738. [Google Scholar] [CrossRef] [Green Version]
  43. Cochrane, F.C.; Davin, L.B.; Lewis, N.G. The Arabidopsis Phenylalanine Ammonia Lyase Gene Family: Kinetic Characterization of the Four PAL Isoforms. Phytochemistry 2004, 65, 1557–1564. [Google Scholar] [CrossRef]
  44. Feinbaum, R.L.; Ausubel, F.M. Transcriptional Regulation of the Arabidopsis Thaliana Chalcone Synthase Gene. Mol. Cell. Biol. 1988, 8, 1985–1992. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Shirley, B.W.; Hanley, S.; Goodman, H.M. Effects of Ionizing Radiation on a Plant Genome: Analysis of Two Arabidopsis Transparent Testa Mutations. Plant Cell 1992, 4, 333–347. [Google Scholar] [CrossRef] [PubMed]
  46. Seitz, C.; Ameres, S.; Forkmann, G. Identification of the Molecular Basis for the Functional Difference between Flavonoid 3′-Hydroxylase and Flavonoid 3′,5′-Hydroxylase. FEBS Lett. 2007, 581, 3429–3434. [Google Scholar] [CrossRef] [PubMed]
  47. Charrier, B.; Coronado, C.; Kondorosi, A.; Ratet, P. Molecular Characterization and Expression of Alfalfa (Medicago Sativa L.) Flavanone-3-Hydroxylase and Dihydroflavonol-4-Reductase Encoding Genes. Plant Mol. Biol. 1995, 29, 773–786. [Google Scholar] [CrossRef]
  48. Fischer, T.C.; Halbwirth, H.; Meisel, B.; Stich, K.; Forkmann, G. Molecular Cloning, Substrate Specificity of the Functionally Expressed Dihydroflavonol 4-Reductases from Malus Domestica and Pyrus Communis Cultivars and the Consequences for Flavonoid Metabolism. Arch. Biochem. Biophys. 2003, 412, 223–230. [Google Scholar] [CrossRef]
  49. Saito, K.; Kobayashi, M.; Gong, Z.; Tanaka, Y.; Yamazaki, M. Direct Evidence for Anthocyanidin Synthase as a 2-Oxoglutarate-Dependent Oxygenase: Molecular Cloning and Functional Expression of CDNA from a Red Forma of Perilla Frutescens. Plant J. Cell Mol. Biol. 1999, 17, 181–189. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Tanner, G.J.; Francki, K.T.; Abrahams, S.; Watson, J.M.; Larkin, P.J.; Ashton, A.R. Proanthocyanidin Biosynthesis in Plants. Purification of Legume Leucoanthocyanidin Reductase and Molecular Cloning of Its CDNA. J. Biol. Chem. 2003, 278, 31647–31656. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Xie, D.-Y.; Sharma, S.B.; Paiva, N.L.; Ferreira, D.; Dixon, R.A. Role of Anthocyanidin Reductase, Encoded by BANYULS in Plant Flavonoid Biosynthesis. Science 2003, 299, 396–399. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Bolling, B.W. Almond Polyphenols: Methods of Analysis, Contribution to Food Quality, and Health Promotion. Compr. Rev. Food Sci. Food Saf. 2017, 16, 346–368. [Google Scholar] [CrossRef] [Green Version]
  53. Jeong, S.T.; Goto-Yamamoto, N.; Hashizume, K.; Esaka, M. Expression of the Flavonoid 3′-Hydroxylase and Flavonoid 3′,5′-Hydroxylase Genes and Flavonoid Composition in Grape (Vitis Vinifera). Plant Sci. 2006, 170, 61–69. [Google Scholar] [CrossRef]
  54. Liu, M.; Tian, H.; Wu, J.-H.; Cang, R.-R.; Wang, R.-X.; Qi, X.-H.; Xu, Q.; Chen, X.-H. Relationship between Gene Expression and the Accumulation of Catechin during Spring and Autumn in Tea Plants (Camellia Sinensis L.). Hortic. Res. 2015, 2, 1–8. [Google Scholar] [CrossRef] [Green Version]
  55. Yang, Q.-Q.; Gan, R.-Y.; Ge, Y.-Y.; Zhang, D.; Corke, H. Polyphenols in Common Beans (Phaseolus Vulgaris L.): Chemistry, Analysis, and Factors Affecting Composition. Compr. Rev. Food Sci. Food Saf. 2018, 17, 1518–1539. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  56. Yoshioka, Y.; Li, X.; Zhang, T.; Mitani, T.; Yasuda, M.; Nanba, F.; Toda, T.; Yamashita, Y.; Ashida, H. Black Soybean Seed Coat Polyphenols Prevent AAPH-Induced Oxidative DNA-Damage in HepG2 Cells. J. Clin. Biochem. Nutr. 2017, 60, 108–114. [Google Scholar] [CrossRef] [Green Version]
  57. Abrahams, S.; Lee, E.; Walker, A.R.; Tanner, G.J.; Larkin, P.J.; Ashton, A.R. The Arabidopsis TDS4 Gene Encodes Leucoanthocyanidin Dioxygenase (LDOX) and Is Essential for Proanthocyanidin Synthesis and Vacuole Development. Plant J. 2003, 35, 624–636. [Google Scholar] [CrossRef] [PubMed]
  58. Kitamura, S.; Matsuda, F.; Tohge, T.; Yonekura-Sakakibara, K.; Yamazaki, M.; Saito, K.; Narumi, I. Metabolic Profiling and Cytological Analysis of Proanthocyanidins in Immature Seeds of Arabidopsis Thaliana Flavonoid Accumulation Mutants. Plant J. 2010, 62, 549–559. [Google Scholar] [CrossRef]
  59. Routaboul, J.-M.; Kerhoas, L.; Debeaujon, I.; Pourcel, L.; Caboche, M.; Einhorn, J.; Lepiniec, L. Flavonoid Diversity and Biosynthesis in Seed of Arabidopsis Thaliana. Planta 2006, 224, 96–107. [Google Scholar] [CrossRef] [PubMed]
  60. Aaby, K.; Ekeberg, D.; Skrede, G. Characterization of Phenolic Compounds in Strawberry (Fragaria × Ananassa) Fruits by Different HPLC Detectors and Contribution of Individual Compounds to Total Antioxidant Capacity. J. Agric. Food Chem. 2007, 55, 4395–4406. [Google Scholar] [CrossRef]
  61. Riihinen, K.; Jaakola, L.; Kärenlampi, S.; Hohtola, A. Organ-Specific Distribution of Phenolic Compounds in Bilberry (Vaccinium Myrtillus) and “northblue” Blueberry (Vaccinium Corymbosum × V. Angustifolium). Food Chem. 2008, 110, 156–160. [Google Scholar] [CrossRef] [PubMed]
  62. Debeaujon, I.; Léon-Kloosterziel, K.M.; Koornneef, M. Influence of the Testa on Seed Dormancy, Germination, and Longevity in Arabidopsis. Plant Physiol. 2000, 122, 403–414. [Google Scholar] [CrossRef] [Green Version]
  63. Qi, T.; Song, S.; Ren, Q.; Wu, D.; Huang, H.; Chen, Y.; Fan, M.; Peng, W.; Ren, C.; Xie, D. The Jasmonate-ZIM-Domain Proteins Interact with the WD-Repeat/BHLH/MYB Complexes to Regulate Jasmonate-Mediated Anthocyanin Accumulation and Trichome Initiation in Arabidopsis Thaliana. Plant Cell 2011, 23, 1795–1814. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Catalá, R.; Medina, J.; Salinas, J. Integration of Low Temperature and Light Signaling during Cold Acclimation Response in Arabidopsis. Proc. Natl. Acad. Sci. USA 2011, 108, 16475–16480. [Google Scholar] [CrossRef] [Green Version]
  65. Cominelli, E.; Gusmaroli, G.; Allegra, D.; Galbiati, M.; Wade, H.K.; Jenkins, G.I.; Tonelli, C. Expression Analysis of Anthocyanin Regulatory Genes in Response to Different Light Qualities in Arabidopsis Thaliana. J. Plant Physiol. 2008, 165, 886–894. [Google Scholar] [CrossRef]
  66. Li, P.; Li, Y.-J.; Zhang, F.-J.; Zhang, G.-Z.; Jiang, X.-Y.; Yu, H.-M.; Hou, B.-K. The Arabidopsis UDP-Glycosyltransferases UGT79B2 and UGT79B3, Contribute to Cold, Salt and Drought Stress Tolerance via Modulating Anthocyanin Accumulation. Plant J. 2017, 89, 85–103. [Google Scholar] [CrossRef] [Green Version]
  67. Xu, W.; Grain, D.; Bobet, S.; Gourrierec, J.L.; Thévenin, J.; Kelemen, Z.; Lepiniec, L.; Dubos, C. Complexity and Robustness of the Flavonoid Transcriptional Regulatory Network Revealed by Comprehensive Analyses of MYB–BHLH–WDR Complexes and Their Targets in Arabidopsis Seed. New Phytol. 2014, 202, 132–144. [Google Scholar] [CrossRef]
  68. Chen, M.; Penfield, S. Feedback Regulation of COOLAIR Expression Controls Seed Dormancy and Flowering Time. Science 2018, 360, 1014–1017. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  69. Gonzalez, A.; Zhao, M.; Leavitt, J.M.; Lloyd, A.M. Regulation of the Anthocyanin Biosynthetic Pathway by the TTG1/BHLH/Myb Transcriptional Complex in Arabidopsis Seedlings. Plant J. 2008, 53, 814–827. [Google Scholar] [CrossRef] [PubMed]
  70. Routaboul, J.-M.; Dubos, C.; Beck, G.; Marquis, C.; Bidzinski, P.; Loudet, O.; Lepiniec, L. Metabolite Profiling and Quantitative Genetics of Natural Variation for Flavonoids in Arabidopsis. J. Exp. Bot. 2012, 63, 3749–3764. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Stracke, R.; Ishihara, H.; Huep, G.; Barsch, A.; Mehrtens, F.; Niehaus, K.; Weisshaar, B. Differential Regulation of Closely Related R2R3-MYB Transcription Factors Controls Flavonol Accumulation in Different Parts of the Arabidopsis Thaliana Seedling. Plant J. 2007, 50, 660–677. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  72. Tian, L.; Pang, Y.; Dixon, R.A. Biosynthesis and Genetic Engineering of Proanthocyanidins and (Iso)Flavonoids. Phytochem. Rev. 2008, 7, 445–465. [Google Scholar] [CrossRef]
  73. Bogs, J.; Downey, M.O.; Harvey, J.S.; Ashton, A.R.; Tanner, G.J.; Robinson, S.P. Proanthocyanidin Synthesis and Expression of Genes Encoding Leucoanthocyanidin Reductase and Anthocyanidin Reductase in Developing Grape Berries and Grapevine Leaves. Plant Physiol. 2005, 139, 652–663. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  74. Bogs, J.; Jaffé, F.W.; Takos, A.M.; Walker, A.R.; Robinson, S.P. The Grapevine Transcription Factor VvMYBPA1 Regulates Proanthocyanidin Synthesis during Fruit Development. Plant Physiol. 2007, 143, 1347–1361. [Google Scholar] [CrossRef] [Green Version]
  75. Todd, J.J.; Vodkin, L.O. Pigmented Soybean (Glycine Max) Seed Coats Accumulate Proanthocyanidins during Development. Plant Physiol. 1993, 102, 663–670. [Google Scholar] [CrossRef] [Green Version]
  76. Tuteja, J.H.; Zabala, G.; Varala, K.; Hudson, M.; Vodkin, L.O. Endogenous, Tissue-Specific Short Interfering RNAs Silence the Chalcone Synthase Gene Family in Glycine Max Seed Coats. Plant Cell 2009, 21, 3063–3077. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Zabala, G.; Vodkin, L. Cloning of the Pleiotropic T Locus in Soybean and Two Recessive Alleles That Differentially Affect Structure and Expression of the Encoded Flavonoid 3′ Hydroxylase. Genetics 2003, 163, 295–309. [Google Scholar] [CrossRef] [PubMed]
  78. Zabala, G.; Vodkin, L.O. The Wp Mutation of Glycine Max Carries a Gene-Fragment-Rich Transposon of the CACTA Superfamily. Plant Cell 2005, 17, 2619–2632. [Google Scholar] [CrossRef] [Green Version]
  79. Zabala, G.; Vodkin, L.O. A Rearrangement Resulting in Small Tandem Repeats in the F3′ 5′ H Gene of White Flower Genotypes Is Associated with the Soybean W1 Locus. Crop Sci. 2007, 47, S-113. [Google Scholar] [CrossRef]
  80. Gillman, J.D.; Tetlow, A.; Lee, J.-D.; Shannon, J.G.; Bilyeu, K. Loss-of-Function Mutations Affecting a Specific Glycine Max R2R3 MYB Transcription Factor Result in Brown Hilum and Brown Seed Coats. BMC Plant Biol. 2011, 11, 1–12. [Google Scholar] [CrossRef] [Green Version]
  81. Yang, K.; Jeong, N.; Moon, J.-K.; Lee, Y.-H.; Lee, S.-H.; Kim, H.M.; Hwang, C.H.; Back, K.; Palmer, R.G.; Jeong, S.-C. Genetic Analysis of Genes Controlling Natural Variation of Seed Coat and Flower Colors in Soybean. J. Hered. 2010, 101, 757–768. [Google Scholar] [CrossRef] [Green Version]
  82. Zhou, S.; Sekizaki, H.; Yang, Z.; Sawa, S.; Pan, J. Phenolics in the Seed Coat of Wild Soybean (Glycine Soja) and Their Significance for Seed Hardness and Seed Germination. J. Agric. Food Chem. 2010, 58, 10972–10978. [Google Scholar] [CrossRef] [PubMed]
  83. Ito, C.; Oki, T.; Yoshida, T.; Nanba, F.; Yamada, K.; Toda, T. Characterisation of Proanthocyanidins from Black Soybeans: Isolation and Characterisation of Proanthocyanidin Oligomers from Black Soybean Seed Coats. Food Chem. 2013, 141, 2507–2512. [Google Scholar] [CrossRef]
  84. Kovinich, N.; Saleem, A.; Arnason, J.T.; Miki, B. Combined Analysis of Transcriptome and Metabolite Data Reveals Extensive Differences between Black and Brown Nearly-Isogenic Soybean (Glycine Max) Seed Coats Enabling the Identification of Pigment Isogenes. BMC Genomics 2011, 12, 381. [Google Scholar] [CrossRef] [Green Version]
  85. Bursać, M.; Krstonošić, M.A.; Miladinović, J.; Malenčić, Đ.; Gvozdenović, L.; Cvejić, J.H. Isoflavone Composition, Total Phenolic Content and Antioxidant Capacity of Soybeans with Colored Seed Coat. Nat. Prod. Commun. 2017, 12, 1934578X1701200417. [Google Scholar] [CrossRef] [Green Version]
  86. Maloney, G.S.; DiNapoli, K.T.; Muday, G.K. The Anthocyanin Reduced Tomato Mutant Demonstrates the Role of Flavonols in Tomato Lateral Root and Root Hair Development. Plant Physiol. 2014, 166, 614–631. [Google Scholar] [CrossRef] [Green Version]
  87. Buer, C.S.; Imin, N.; Djordjevic, M.A. Flavonoids: New Roles for Old Molecules. J. Integr. Plant Biol. 2010, 52, 98–111. [Google Scholar] [CrossRef] [PubMed]
  88. Baudry, A.; Heim, M.A.; Dubreucq, B.; Caboche, M.; Weisshaar, B.; Lepiniec, L. TT2, TT8, and TTG1 Synergistically Specify the Expression of BANYULS and Proanthocyanidin Biosynthesis in Arabidopsis Thaliana. Plant J. 2004, 39, 366–380. [Google Scholar] [CrossRef] [PubMed]
  89. Feuillet, C.; Schachermayr, G.; Keller, B. Molecular Cloning of a New Receptor-like Kinase Gene Encoded at the Lr10 Disease Resistance Locus of Wheat. Plant J. 1997, 11, 45–52. [Google Scholar] [CrossRef]
  90. Tohge, T.; Nishiyama, Y.; Hirai, M.Y.; Yano, M.; Nakajima, J.; Awazuhara, M.; Inoue, E.; Takahashi, H.; Goodenowe, D.B.; Kitayama, M. Functional Genomics by Integrated Analysis of Metabolome and Transcriptome of Arabidopsis Plants Over-expressing an MYB Transcription Factor. Plant J. 2005, 42, 218–235. [Google Scholar] [CrossRef] [PubMed]
  91. Xie, D.-Y.; Sharma, S.B.; Wright, E.; Wang, Z.-Y.; Dixon, R.A. Metabolic Engineering of Proanthocyanidins through Co-Expression of Anthocyanidin Reductase and the PAP1 MYB Transcription Factor. Plant J. 2006, 45, 895–907. [Google Scholar] [CrossRef] [Green Version]
  92. Chen, Q.; Yu, H.; Tang, H.; Wang, X. Identification and Expression Analysis of Genes Involved in Anthocyanin and Proanthocyanidin Biosynthesis in the Fruit of Blackberry. Sci. Hortic. 2012, 141, 61–68. [Google Scholar] [CrossRef]
  93. De Pascual-Teresa, S.; Moreno, D.A.; García-Viguera, C. Flavanols and Anthocyanins in Cardiovascular Health: A Review of Current Evidence. Int. J. Mol. Sci. 2010, 11, 1679–1703. [Google Scholar] [CrossRef] [Green Version]
  94. Hammerstone, J.F.; Lazarus, S.A.; Schmitz, H.H. Procyanidin Content and Variation in Some Commonly Consumed Foods. J. Nutr. 2000, 130, 2086S–2092S. [Google Scholar] [CrossRef]
  95. Ferreira, D.; Nel, R.J.J.; Bekker, R. 3.19—Condensed Tannins. In Comprehensive Natural Products Chemistry; Barton, S.D., Nakanishi, K., Meth-Cohn, O., Eds.; Pergamon: Oxford, UK, 1999; ISBN 978-0-08-091283-7. [Google Scholar]
  96. Ferreira, D.; Marais, J.P.J.; Coleman, C.M.; Slade, D. 6.18—Proanthocyanidins: Chemistry and Biology. In Comprehensive Natural Products II; Liu, H.-W., Mander, L., Eds.; Elsevier: Oxford, UK, 2010; pp. 605–661. ISBN 978-0-08-045382-8. [Google Scholar]
  97. Dixon, R.A.; Sarnala, S. Proanthocyanidin Biosynthesis—A Matter of Protection. Plant Physiol. 2020, 184, 579–591. [Google Scholar] [CrossRef]
  98. Kobayashi, H.; Suzuki, S.; Tanzawa, F.; Takayanagi, T. Low Expression of Flavonoid 3′,5′-Hydroxylase (F3′,5′H) Associated with Cyanidin-Based Anthocyanins in Grape Leaf. Am. J. Enol. Vitic. 2009, 60, 362–367. [Google Scholar]
  99. Dixon, R.A.; Liu, C.; Jun, J.H. Metabolic Engineering of Anthocyanins and Condensed Tannins in Plants. Curr. Opin. Biotechnol. 2013, 24, 329–335. [Google Scholar] [CrossRef] [PubMed]
  100. Edgar, R.C. MUSCLE: Multiple Sequence Alignment with High Accuracy and High Throughput. Nucleic Acids Res. 2004, 32, 1792–1797. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  101. Castresana, J. Selection of Conserved Blocks from Multiple Alignments for Their Use in Phylogenetic Analysis. Mol. Biol. Evol. 2000, 17, 540–552. [Google Scholar] [CrossRef] [Green Version]
  102. Tamura, K.; Stecher, G.; Kumar, S. MEGA11: Molecular Evolutionary Genetics Analysis Version 11. Mol. Biol. Evol. 2021. [Google Scholar] [CrossRef] [PubMed]
Plants 10 01792 g001 550
Figure 1. Biosynthetic pathway of PAs: (A) Basic chemical structure of flavan-3-ols. Flavan-3-ols have two benzene rings and one heterocyclic benzopyran ring. (B) A schematic diagram with chemical structures of flavanol compounds. Triangles indicate the species-specific enzymes involved in each biosynthetic pathway. The green triangle shows the enzymes that exist in all of our target species (Table 1); exceptionally, the LAR* gene, which catalyzes (2R, 3S)-flavanols such as catechin, was not identified in A. thaliana. A blue or red triangle indicates the enzyme that is present in grape, blueberry, and tea plant, or strawberry and almond, respectively. The names of the metabolites have been colored based on their hydroxyl pattern on the B-ring, and the black, orange, and pink designate 3′4′-dihydroxyl, 4′-hydroxyl, and 4′5′6′-trihydroxyl patterns, respectively. PAs: proanthocyanidins; ANR: anthocyanidin reductase; ANS: anthocyanidin synthase; LAR: leucoanthocyanidin reductase; F3′H: flavonoid 3′-hydroxylase; F3′5′H: flavonoid 3′,5′-hydroxylase; DFR: dihydroflavonol.
Figure 1. Biosynthetic pathway of PAs: (A) Basic chemical structure of flavan-3-ols. Flavan-3-ols have two benzene rings and one heterocyclic benzopyran ring. (B) A schematic diagram with chemical structures of flavanol compounds. Triangles indicate the species-specific enzymes involved in each biosynthetic pathway. The green triangle shows the enzymes that exist in all of our target species (Table 1); exceptionally, the LAR* gene, which catalyzes (2R, 3S)-flavanols such as catechin, was not identified in A. thaliana. A blue or red triangle indicates the enzyme that is present in grape, blueberry, and tea plant, or strawberry and almond, respectively. The names of the metabolites have been colored based on their hydroxyl pattern on the B-ring, and the black, orange, and pink designate 3′4′-dihydroxyl, 4′-hydroxyl, and 4′5′6′-trihydroxyl patterns, respectively. PAs: proanthocyanidins; ANR: anthocyanidin reductase; ANS: anthocyanidin synthase; LAR: leucoanthocyanidin reductase; F3′H: flavonoid 3′-hydroxylase; F3′5′H: flavonoid 3′,5′-hydroxylase; DFR: dihydroflavonol.
Plants 10 01792 g001
Plants 10 01792 g002 550
Figure 2. The phylogenetic tree was constructed using three genes: ANR, ANS, and DFR. The protein sequences of the three genes were obtained from the reference genomic data of the 12 species (including Oryza sativa as an outgroup). The protein sequences were concatenated and aligned using the MUSCLE algorithm, and the analysis was performed using the maximum likelihood method. The numbers between nodes are bootstrap values indicating the reliability of the genetic tree, and branches have been colored with their family indicated on the right-hand side of the figure. ANR: anthocyanidin reductase; ANS: anthocyanidin synthase; DFR: dihydroflavonol 4-reductase.
Figure 2. The phylogenetic tree was constructed using three genes: ANR, ANS, and DFR. The protein sequences of the three genes were obtained from the reference genomic data of the 12 species (including Oryza sativa as an outgroup). The protein sequences were concatenated and aligned using the MUSCLE algorithm, and the analysis was performed using the maximum likelihood method. The numbers between nodes are bootstrap values indicating the reliability of the genetic tree, and branches have been colored with their family indicated on the right-hand side of the figure. ANR: anthocyanidin reductase; ANS: anthocyanidin synthase; DFR: dihydroflavonol 4-reductase.
Plants 10 01792 g002
Table
Table 1. List of 10 major crops containing high levels of proanthocyanidins and their reference genomic information. Among major crops with high contents of PAs, the species with currently available genomic databases were selected for this study.
Table 1. List of 10 major crops containing high levels of proanthocyanidins and their reference genomic information. Among major crops with high contents of PAs, the species with currently available genomic databases were selected for this study.
SpeciesProanthocyanidins Content
(mg/100 g)		ReferenceGenome DatabaseReferenceAssembly Size
(Mb)		Coverage
(%)		Contig N50
(Kb)		Number of Genes Predicted
Almond
(Prunus dulcis)		184Prior & Gu, 2005 [14]Prunus dulcis Lauranne Genome v1.0
(http://rosaceae.org/, accessed on 12 February 2021)		Alioto et al., 2020 [31]2279510327,969
Apple tree
(Malus domestica)		162Hellström et al., 2009 [2](iris.angers.inra.fr/gddh13/, accessed on 12 February 2021)Daccord et al., 2017 [32]64310062042,140
Blueberry
(Vaccinium corymbosum)		255Prior & Gu, 2005 [14]V_corymbosum v1.0
(http://gigadb.org/, accessed on 12 February 2021)		Colle et al., 2019 [33]16801021532,140
Cacao bean
(Theobroma cacao)		1460Hellström et al., 2009 [2]Cacao Matina1-6 Genome v2.1
(http://cacaogenomedb.org, accessed on 12 February 2021)		Publication in progress (http://cacaogenomedb.org)34680108027,379
Common bean
(Phaseolus vulgaris)		1000Kan et al., 2016 [34]Phaseolus vulgaris v2.1
(http://phytozome.jgi.doe.gov/, accessed on 12 February 2021)		Schmutz et al., 2014 [35]60080190027,433
Grape
(Vitis vinifera)		3532Prior & Gu, 2005 [14]Vitis vinifera v2.1
(http://phytozome.jgi.doe.gov/, accessed on 12 February 2021)		Jaillon et al., 2007 [36]48710256626,346
Peanut
(Arachis hypogaea)		186Hellström et al., 2009 [2](http://peanutgr.fafu.edu.cn/, accessed on 12 February 2021)Zhuang et al., 2019 [37]253894150983,709
Soybean
(Glycine max)		300Lee et al., 2017 [38]Glycine max Wm82.a4
(http://www.soybase.org, accessed on 12 February 2021)		Schmutz et al., 2010 [39]115095149246,430
Strawberry
(Fragaria × ananassa)		145Prior & Gu, 2005 [14](https://datadryad.org/, accessed on 12 February 2021)Edger et al., 2019 [40]813 *9979108,087
Tea tree
(Camellia sinensis)		189Engelhardt et al., 2003 [12](http://tpia.teaplant.org/, accessed on 12 February 2021)Xia et al., 2019 [41]289095 6753,512
* This assembly genome size (Fragaria × ananassa) is the haploid genome size.
Table
Table 2. Numbers of genes involved in proanthocyanidins production in the 11 species. Based on sequence similarity, orthologous genes were searched using BLASTP.
Table 2. Numbers of genes involved in proanthocyanidins production in the 11 species. Based on sequence similarity, orthologous genes were searched using BLASTP.
Speciesa Number of Genes from KEGG Pathwayb Number of Genes Confirmed from Reference Databasec Number of Genes Newly Identified Using aNumber of Orthologous Genes Newly Identified Using b+c
Arachis hypogea6 *651
Arabidopsis thaliana66--
Fragaria × ananassa4 *497
Glycine max2019--
Camelia sinensis---19
Malus domestica1212--
Phaseolus vulgaris1211--
Theobroma cacao98--
Vitis vinifera2317--
Prunus dulcis---4
Vaccinium corymbosum---74
Total928314105
* The genes involved in the biosynthetic pathways have been reported in wild relative species of strawberry (Fragaria vesca) and peanut (Arachis duranensis). a the number of genes identified from KEGG database, b the number ofgenes confirmed from reference database, c the number ofgenes newly identified from reference genome using a.
Table
Table 3. Numbers of the seven orthologs encoding key enzymes that synthesize PAs identified in the reference genome of each crop.
Table 3. Numbers of the seven orthologs encoding key enzymes that synthesize PAs identified in the reference genome of each crop.
Arachis hypogeaArabidopsis thalianaFragaria × ananassaGlycine maxMalus domesticaPhaseolus vulgarisTheobroma cacaoVitis viniferaCamelia sinensisPrunus dulcisVaccinium corymbosum
F3′H21352112104
F3H21442112208
F3′5′H000102175037
DFR21432313316
LAR204221223112
ANS22322111313
ANR21222210214
Total1262019121181719474
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Lim, I.; Ha, J. Biosynthetic Pathway of Proanthocyanidins in Major Cash Crops. Plants 2021, 10, 1792. https://doi.org/10.3390/plants10091792

AMA Style

Lim I, Ha J. Biosynthetic Pathway of Proanthocyanidins in Major Cash Crops. Plants. 2021; 10(9):1792. https://doi.org/10.3390/plants10091792

Chicago/Turabian Style

Lim, Insu, and Jungmin Ha. 2021. "Biosynthetic Pathway of Proanthocyanidins in Major Cash Crops" Plants 10, no. 9: 1792. https://doi.org/10.3390/plants10091792

APA Style

Lim, I., & Ha, J. (2021). Biosynthetic Pathway of Proanthocyanidins in Major Cash Crops. Plants, 10(9), 1792. https://doi.org/10.3390/plants10091792

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop