María Dolores Requena-Ramírez
Sergio G. Atienza
Dámaso Hornero-Méndez
Cristina Rodríguez-Suárez- 1Instituto de Agricultura Sostenible ‐ Consejo Superior de Investigaciones Científicas, Córdoba, Spain
- 2Department of Food Phytochemistry, Instituto de la Grasa ‐ Consejo Superior de Investigaciones Científicas, Campus Universidad Pablo de Olavide, Seville, Spain
Carotenoids are essential in human diet, so that the development of programs toward carotenoid enhancement has been promoted in several crops. The cereal tritordeum, the amphiploid derived from the cross between Hordeum chilense Roem. et Schulz. and durum wheat has a remarkable carotenoid content in the endosperm. Besides, a high proportion of these carotenoids are esterified with fatty acids. The identification of the gene(s) responsible for xanthophyll esterification would be useful for breeding as esterified carotenoids show an increased ability to accumulate within plant cells and have a higher stability during post-harvest storage. In this work, we analyzed five genes identified as candidates for coding the xanthophyll acyltransferase (XAT) enzyme responsible for lutein esterification in H. chilense genome. All these genes were expressed during grain development in tritordeum, but only HORCH7HG021460 was highly upregulated. Sequence analysis of HORCH7HG021460 revealed a G-to-T transversion, causing a Glycine to Cysteine substitution in the protein of H290 (the only accession not producing quantifiable amounts of lutein esters, hereinafter referred as zero-ester) of H. chilense compared to the esterifying genotypes. An allele-specific marker was designed for the SNP detection in the H. chilense diversity panel. From the 93 accessions, only H290 showed the T allele and the zero-ester phenotype. Furthermore, HORCH7HG021460 is the orthologue of XAT-7D, which encodes a XAT enzyme responsible for carotenoid esterification in wheat. Thus, HORCH7HG021460 (XAT-7Hch) is a strong candidate for lutein esterification in H. chilense and tritordeum, suggesting a common mechanism of carotenoid esterification in Triticeae species. The transference of XAT-7Hch to wheat may be useful for the enhancement of lutein esters in biofortification programs.
Introduction
Carotenoids pigments are responsible for many of the yellow, orange, and bright red colors of fruits, vegetables, and flowers. They have essential roles in light-harvesting, protection against excess of light energy and oxidative damage (Britton, 2008). Humans are not able to synthesize carotenoids and thus their supply relies on diary intake. The consumption of carotenoid-rich foods has been associated with a lower risk of developing some diseases, including several types of cancer (Nishino et al., 2009). Vitamin A deficiency is a serious concern in developing countries since it leads to increased childhood mortality, growth retardation, and blindness (WHO, 2009). All carotenoids with at least one unsubstituted β-end ring have provitamin A activity, and consequently, β-carotene is the main vitamin A precursor (Giuliano et al., 2008). Lutein has also been suggested to play an important role during neural system development in infant brains (Johnson, 2014). In addition, lutein and zeaxanthin are particularly concentrated in the eye macula, i.e., the central part of the retina (Johnson, 2014), and their consumption alleviates age macular degeneration (AMD; Krinsky et al., 2003).
In wheat and related cereals, lutein is the main carotenoid in grains (Rodríguez-Suárez et al., 2010; Ficco et al., 2014). Tritordeum (×Tritordeum martini A. Pujadas; Pujadas-Salvá, 2016) is the amphiploid derived from the cross between Hordeum chilense Roem. et Schultz. and durum wheat, and it has an unusually high carotenoid content (Martín et al., 1999). On average, tritordeums have over five times more carotenoids in seeds than their corresponding durum wheat parents (Atienza et al., 2007b). Thus, this species is suitable for the development of biofortification programs as in the case of maize (Menkir et al., 2018; Zhu et al., 2018) and rice (Paine et al., 2005). Tritordeum is being commercialized by Agrasys S.L.1 At present, tritordeum products are available in many European countries and Australia.2
The carotenoid biosynthesis pathway in plants is well-known (Rodriguez-Concepcion et al., 2018). The first step is regulated by phytoene synthase (PSY). In grasses, there are three Psy genes, Psy1, Psy2, and Psy3 (Li et al., 2009), but only Psy1 contributes to endosperm carotenoid accumulation (Li et al., 2009). The same happens in tritordeum, with Psy1_Hch being the main responsible for the endosperm carotenoid content (Atienza et al., 2007a; Rodríguez-Suárez et al., 2014; Mattera et al., 2015).
In addition to carotenoid content, the bioavailability of these pigments and its retention through the food chain are of the outmost importance (Fernández-García et al., 2012). Lutein esters have been found to be more stable than free lutein at 60°C (Ahmad et al., 2013). Indeed, collective evidence show that esterification increases carotenoids stability and accumulation (Mellado-Ortega and Hornero-Méndez, 2015b; Meléndez-Martínez et al., 2019); thus, improving carotenoid esterification in cereal endosperm would be useful for biofortification programs.
Tritordeum contains a high proportion of lutein esters in the endosperm, while durum wheat produces low amounts (Atienza et al., 2007b) or no lutein esters at all (Ziegler et al., 2015). Accordingly, this high proportion of lutein esters found in tritordeum was suggested to be caused by genes from the H. chilense genome (Hch). This hypothesis was confirmed by analyzing H. chilense seeds (Mellado-Ortega and Hornero-Méndez, 2015a) and H. chilense-common wheat genetic stocks (Mattera et al., 2015). Furthermore, these results showed the existence of genes for lutein esterification in both chromosomes 7Hch from H. chilense and 7D from common wheat (Mattera et al., 2015). The presence of genes for lutein esterification in the D genome from common wheat is consistent with the diversity observed for carotenoid esterification in this species (Ziegler et al., 2015; Paznocht et al., 2019).
Using a 7D chromosome deletion set in common wheat, along with the reference genome (cv. “Chinese Spring”) and the Wheat Expression Browser (Ramírez-González et al., 2018), we identified the target region harboring the gene(s) for lutein esterification in chromosome 7D (Ávila et al., 2019b). This allowed for the selection of a set of candidate genes with xanthophyll acyltransferase (XAT) function, which could be responsible for the Lute locus associated to lutein esterification in the cross Haruhikari//Sunco/Indis.82 (Ahmad et al., 2015). Due to the high degree of genetic collinearity that H. chilense shows with Triticeae species (Ávila et al., 2019a), the set of candidate genes identified in chromosome 7D is also valid for chromosome 7Hch, since orthologs of these genes are likely present in chromosome 7Hch. Indeed, genetic studies have shown that orthologous genes for traits of interest are in equivalent positions in H. chilense compared to wheat and/or barley, for instance Psy1 (Atienza et al., 2007a), APETALA2-like (Gil-Humanes et al., 2009), carotenoid biosynthesis pathway (Rodríguez-Suárez and Atienza, 2012), Ppo (Rodríguez-Suárez and Atienza, 2014), Wx (Alvarez et al., 2019), and Rfm1 (Rodríguez-Suárez et al., 2020).
The objectives of this work were to identify the orthologue genes of H. chilense corresponding to a set of candidate genes with putative XAT function in chromosome 7D, to investigate their expression profile during grain development in tritordeum, and to identify the gene responsible for lutein esterification using a diversity panel of H. chilense.
Materials and Methods
Candidate Gene Identification and Amplification in Hordeum chilense
Hordeum chilense (H7 and H16), wheat (Triticum aestivum “Chinese Spring”), and barley (Hordeum vulgare “Betzes”) seeds were sown in pots with a mixture 1:4.5 (v:v; sand:premium pflanzerde substrate). Osmocote Exact Mini (1.2 g/L) Everris International B.V. was added. Two-week-old leaves were used for DNA isolation. Genomic DNA was extracted using the CTAB method according to Murray and Thompson (1980) with slight modifications.
A set of five candidate genes located in the homeologous group 7 of the Triticeae species were selected (Table 1). This set of genes was selected from a previous genetic study (Ávila et al., 2019b). The annotation of these genes is related to the esterification process, and all of them are expressed during grain development in common wheat as determined using the Wheat Expression Browser powered by expVIP (Borrill et al., 2016; Ramírez-González et al., 2018). The sequences of these genes in wheat were retrieved from Ensembl Plants website and used to identify the orthologous genes in barley (Howe et al., 2019; Additional File 1, Supplementary Table S1).3
Table 1. Candidate genes for lutein esterification.
To obtain H. chilense genomic sequences, the wheat (cv. Chinese Spring) and barley (cv. Morex) sequences available for each gene at Ensembl Plants were aligned using the multiple sequence alignment tool ClustalW (Madeira et al., 2019). Primers were designed in the conserved regions to ensure amplification in H. chilense templates (Additional File 1, Supplementary Table S2). The NCBI Primer-Blast tool was used for primer design (Ye et al., 2012). PCR amplifications were carried out with MyTaq™ DNA polymerase (Bioline, London, UK) following the manufacturer’s instructions, using H7 and H16 H. chilense accessions, along with barley (cv. Betzes) and wheat (cv. Chinese Spring) samples as controls. Amplification products were resolved in agarose gels and visualized with Safeview™ Nucleic Acid Stain (NBS Biologicals, Ltd., Cambridgeshire, England). The PCR fragments obtained were directly sequenced at STAB VIDA (Portugal).4
Quantitative Real-Time PCR During Grain Development
Tritordeum samples used for quantitative real-time PCR (qPCR) experiments were obtained in a previous work (Rodríguez-Suárez et al., 2014). Briefly, five hexaploid tritordeums [HT240, HT335, HT609, HT621 (Ballesteros et al., 2005), and HT630] were grown at field conditions following a completely randomized design with two replicates. Samples for qPCR analyses were obtained from developing grains of tritordeums. Spikes were marked at anthesis and developing grains were collected with forceps at 8 (Zadoks 71), 14 (Zadoks 73), 18 (Zakoks 77), and 25 (Zakoks 83) days post-anthesis (dpa; Zadoks et al., 1974) and immediately frozen in liquid nitrogen at −80°C. From each block, two independent biological replicates were collected (total of four biological replicates, two samples per block and two blocks in the field trial). From each biological replicate, four grains were randomly selected for RNA extraction. Thus, eighty samples were analyzed (five genotypes × four developmental stages × four biological samples).
Total RNA was isolated from whole grains using the TRIzol® Reagent (Invitrogen, Carlsbad, CA) according to manufacturer’s instructions with slight modifications, and cDNA was obtained as described previously (Rodríguez-Suárez et al., 2011).
Once H. chilense partial sequences were obtained for the five genes, RT-qPCR primers were designed to specifically amplify H. chilense gene copy in tritordeum. For this purpose, sequence sets for each gene were aligned and primers were designed as mentioned above considering polymorphic regions. Primer pair amplification specificity was checked by the standard PCR.
Real-time qPCR reactions for target genes were carried out using four biological replicates with two technical duplicates each from cDNA samples, which is consisted in a 1:4 dilution of each sample. Amplicon sizes ranged from 99 to 155 bp. RT-qPCR reactions were carried out using the SYBR® Green on an Applied Biosystems™ 7500 Real-Time PCR System (Applied Biosystems, Foster City, CA, USA). The PCR conditions were 95°C for 10 min followed by 40 cycles at 95°C for 15 s and 60°C for 1 min. qPCR reactions were performed using 5 μl of diluted cDNA (10 ng/μl), 5 μl of iTaqTM Universal SYBR® Green Supermix (Bio-Rad, Hercules, California), and a primer pair concentration of 0.22 μM each. No template controls were included in the reactions.
Homeolog-specific amplification was confirmed by the presence of a single peak in the qPCR dissociation curve analysis. PCR efficiency of each primer pair was determined by the LinRegPCR quantitative PCR data analysis program (version 11.0; Ruijter et al., 2009) using raw normalized fluorescence as input data. Expression of target genes for each sample (N0) was determined using the equation N0 = 0.2/ECq, where E is the PCR efficiency for each primer and Cq is the number of cycles needed to reach 0.2 arbitrary units of fluorescence.
Geometric mean of reference genes ADP-RF(m), RLI(a), and CDC(a) was used for normalization (Giménez et al., 2011). Expression stability of these genes and normalization factors for each sample were assessed using GeNorm (Vandesompele et al., 2002).
SNP Detection in HORCH7HG021460 Gene
Three H. chilense genotypes were used to obtain the complete sequence of this gene, H7, H16, and H290, the only H. chilense accession not producing quantifiable amounts of lutein esters (hereinafter referred as zero-ester genotype; Ávila et al., 2019a). For the complete DNA genomic sequencing of HORCH7HG021460, three overlapping fragments were amplified, sequenced directly from PCR products and assembled (Additional File 2, Supplementary Figure S1).
SignalP 4.0 software was used for signal peptide and cellular location prediction of the expected proteins (Petersen et al., 2011). Search for conserved domain analysis was carried out by using the Conserved Domain Database (Lu et al., 2019).
A SNP was identified between H290 (zero-esters) and H7, H16. Primer 1 web service was used for primer design, following a Tetra-Primer ARMS (amplification refractory mutation system) strategy (Collins and Ke, 2012) for SNP detection in HORCH7HG021460: Fw_out GCTTACTGCTCCTGGCCATCTTCCTC; Rv_out TTGAACCCAATCTTTTCAGCTGCAACAA; Fw_inn_H290 CATGTCGCCACAGGGAGGTTCTGCAGCT; Rv_inn_wild AACCTATGAAATCGATAAGCAGCTTACC.
A diversity panel of 93 accessions of H. chilense, previously characterized for carotenoid content and profile (Ávila et al., 2019a), was used for the validation of the SNP detected in this work. The carotenoid content and profile of these accessions are shown in Additional File 3. The origin of these accessions has been previously reported (Vaz Patto et al., 2001).
Results
Candidate Genes Identification in Hordeum chilense
To identify H. chilense orthologs of the candidate genes for lutein esterification in the D genome (Table 1), wheat homeologues, along with H. vulgare orthologs, were identified using the wheat reference genome (Additional File 1, Supplementary Table S1). Gene sequences were aligned and used as a template for primer design considering the conserved regions. Amplification products obtained in H. chilense accessions H7 and H16 were directly sequenced. The identity of the H. chilense sequences obtained was confirmed by BLASTn, showing high homologies with the corresponding orthologue genes (Additional File 1, Supplementary Table S3). The five genes were then named based on H. vulgare gene’s ID: HORCH7HG001330, HORCH7HG007590, HORCH7HG021310, HORCH7HG021460, and HORCH7HG024220. Partial sequences of these genes were used along with those from barley and wheat to design gene-specific primers for RT-qPCR (Additional File 1, Supplementary Table S2). These primers only amplified H. chilense copies as confirmed by conventional PCR.
Transcript Profiling of Hordeum chilense Candidate Genes for Lutein Esterification in Tritordeum
Tritordeum samples from four grain developing stages (St) were used for transcript profiling: St1 [Zadoks 71, 8 days post-anthesis(dpa)]; St2 (Zadoks 73, 14 dpa); St3 (Zadoks 77, 18 dpa) and St4 (Zadoks 83, 25 dpa; Zadoks et al., 1974). These developmental stages correspond to the caryopsis water ripe (Z71), early milk (Z73), late milk (Z77) and early dough (Z83; Zadoks et al., 1974). Several primer pairs were designed for each candidate gene, selecting those with specific amplification (Additional File 1, Supplementary Table S2). For each gene, the change in expression was calculated relative to the expression at St1 (Figure 1) considering the mean values and standard errors calculated with all tritordeum samples (five genotypes × two biological replicates × two field plots). The five H. chilense genes were expressed during grain development but they showed great differences in their expression patterns (Figure 1). Indeed, HORCH7HG021460 increased remarkably its relative abundance during grain development with a 38.6-fold up-regulation in St4 relative to St1 (Figure 1). All the other genes were in the range 0.7–2.0 relative to St1 with the exception of HORCH7HG021310 which showed a 2.5 fold increase in St3.
Figure 1. Fold change expression of candidate genes during grain development in tritordeum. For each gene (represented by a different color), the fold change of expression at each developmental stage was calculated relative to stage 1 (Zadoks 71). Ratios above 1 indicate upregulation of the gene compared to the first stage (Zadoks 71). The gene HORCH7HG001330 has ratios below 1, which indicates a down regulation of this gene compared to the first stage (Zadoks 71). Bars were calculated as the mean of five different tritordeum genotypes ± SE.
Sequence Analysis in a Hordeum chilense Accession Unable to Synthesize Lutein Esters
This high increase in HORCH7HG021460 expression at the end of the grain development may indicate a putative role in lutein esterification. Recently we identified the genotype H290, a H. chilense accession without quantifiable amounts of lutein esters (zero-ester genotype; Ávila et al., 2019a). We obtained the complete gene sequence of HORCH7HG021460 in H290 and compared it to two H. chilense esterifying genotypes H7 and H16. For that purpose, a set of three primers pairs were designed (Additional File 1, Supplementary Table S2) to amplify overlapping fragments and obtain the gene sequence in the selected genotypes. The intron/exon structure was determined based on that of H. vulgare gene HORVU7Hr1G021460. The H. chilense gene was predicted to have a five exon structure, with ORFs of 1,462 bp in H290 and 1,460 bp in H7 and H16. A 22 amino acid-length signal peptide was found in the sequences of H7, H16, and H290, predicting an extracellular protein location. HORCH7HG021460 encodes for a 352 amino acid-length GDSL protein (characterized by the conserved GDSL motif), where the active sites at positions 38 (S), 115 (T), 169 (C), 327 (D), and 330 (H) have been identified (Figure 2). After comparing the coding sequences, a single SNP was found at the end of exon 1: a G in H7 and H16 and a T in H290 (Supplementary Figure S1; Additional File 2). In the predicted proteins, this transversion causes a Glycine (Gly) to Cysteine (Cys) amino acid change at position 75 in zero-ester genotype (Figure 2).
Figure 2. Alignment of XAT proteins. Protein alignment of HORCH7HG021460 of Hordeum chilense accessions H16 (identical to H7) and H290, HORVU7Hr1G021460 of Hordeum vulgare and TraesCS7D094000 of bread wheat. The predicted active sites are shown with black triangles. Black arrows show the mutations described in Gladius TILLING population (Watkins et al., 2019). The mutation detected in the zero-ester genotype H290 is pointed with a red arrow, within the conserved consensus sequence of block IIB (Akoh et al., 2004). Signal peptide in positions 1–22 is also highlighted in HORCH7HG021460 proteins.
Validation of the SNP for Zero-Ester in the Hordeum chilense Diversity Panel
A tetra-Primer ARMS strategy was optimized to detect the G/T SNP in gene HORCH7HG021460. The outer primer pair amplified a common fragment of 344 bp in all genotypes. A 163 bp fragment amplified in the allele present in H290 whereas a 237 bp fragment was obtained in H7 and H16 genotypes. The SNP was scored in the H. chilense diversity panel composed by 93 accessions, previously characterized for their carotenoid and esterification profile (Ávila et al., 2019a). All lines producing esters (92) showed the 237 bp allele, and only H290 showed the 163 bp allele and the zero-ester phenotype.
Discussion
The genes studied in this work were initially identified as candidates to explain the underlying variation for lutein esterification in chromosome 7D of common wheat (Ávila et al., 2019b). They were selected considering their annotated function, their position in the “Chinese Spring” genome, and their expression profile retrieved from the Wheat Expression Browser. Hordeum chilense shows a good degree of collinearity with H. vulgare (Ávila et al., 2019a) and, thus, with all the members of the Triticeae tribe (Mayer et al., 2011). Indeed, the position of H. chilense genes identified using candidate gene approaches in previous works is in agreement with their position in other Triticeae species (Atienza et al., 2007a; Gil-Humanes et al., 2009; Rodríguez-Suárez and Atienza, 2012; Alvarez et al., 2019; Rodríguez-Suárez et al., 2020). Thus, these candidate genes in chromosome 7D were a good starting point for identifying the gene involved in lutein esterification in Hch genome.
One of the hypotheses to be tested was if the orthologue genes from H. chilense were expressed during grain development, as it happens with the candidate genes from D genome in wheat. In fact, all the candidate genes tested were expressed during grain development with different expression patterns (Figure 1), among which HORCH7HG021460 showed a remarkable upregulation. Interestingly, its wheat orthologue TRAESCS7D02g094000 (XAT-7D) has been recently identified as the gene coding for the XAT enzyme responsible for lutein esterification in common wheat (Watkins et al., 2019). As observed for XAT-7D, HORCH7HG021460 encodes for a GDSL protein, predicted to have an extracellular (apoplastic) location based on the conserved domain and signal peptide analyses.
To clear up whether HORCH7HG021460 may be associated with variation for lutein esterification in H. chilense, we studied the accession H290, the only zero-ester member within the species diversity panel (Ávila et al., 2019a). Genomic sequence of HORCH7HG021460 in H290 was compared with those of the two selected esterifying accessions, H7 and H16. This analysis revealed a mutation leading to an amino acid change from Gly to Cys at position 75 in the H290 protein (Figure 2). This position takes part of the conserved domain for block IIB (TGRF-NG----D--) defined for GDSL proteins by Akoh et al. (2004). Although this amino acid position is not identified as an active site, it is highly conserved among distant species (Additional File 4, Supplementary Figure S2), which may suggest its importance in the carotenoid esterification protein activity. Watkins et al. (2019) analyzed a TILLING population derived from EMS-mutagenized wheat cultivar Glaudius and found six independent mutations leading to the absence of lutein esters. None of these mutations correspond to the active sites defined for GDSL proteins (Figure 2), and they are different from that found in H290. Thus, not only mutations in the active sites may lead to non-functional GDSL proteins but the mutation found in accession H290 may also explain the zero-ester phenotype. Additionally, all accessions from the diversity panel harbor the “G” allele (237 bp fragment) of HORCH7HG021460 and are able to produce esters, and only H290 shows both the zero-ester phenotype and the “T” allele (163 bp fragment). Thus, we propose HORCH7HG021460 (XAT-7Hch) as the gene responsible for lutein esterification in H. chilense and tritordeum. Although a functional analysis would be necessary for further confirmation, the findings in the diversity panel and the functional validation of the wheat ortholog XAT-7D are strong evidences supporting the role of XAT-7Hch in the esterification of lutein in tritordeum.
Here, we provide evidence that the same XAT mechanism may be acting in both D and Hch genomes. This is relevant for Triticeae species since it points out to a general mechanism of carotenoid esterification useful for cereals biofortification. Indeed, the simultaneous presence of chromosomes 7D and 7Hch (carrying XAT-7D and XAT-7Hch, respectively) in H. chilense-wheat chromosome substitution lines results in higher amounts of lutein esters (Mattera et al., 2015). This may indicate an additive effect of XAT-7D and XAT-7Hch, which is in agreement with previous results, suggesting that the amount of lutein esters may be limited by the existence of specific XAT enzymes and not by substrate availability (Mellado-Ortega and Hornero-Méndez, 2018).
Our collective data suggest some similarities and some differences between the XAT enzymes in bread wheat and H. chilense. Watkins et al. (2019) predicted XAT to be localized to the apoplast, while carotenoids are produced in plastids (Shumskaya et al., 2012). As a consequence, these authors hypothesized that this subcellular partitioning of substrates and enzymes would constitute a mechanism of regulation of lutein esterification in wheat grain throughout time. The signal peptide of XAT-7Hch also predicts an apoplastic location, which is in agreement with our observations. In a previous work, we investigated the carotenoid profile from 20 to 40 dpa in developing seeds of tritordeum and in three common wheat-H. chilense chromosome substitution lines carrying chromosome 7Hch (Mattera et al., 2020). Despite the high amounts of carotenoids, no lutein esters were detected until 36 dpa (Mattera et al., 2020), when desiccation had been initiated and endosperm cells would be under their programmed cell death, resulting in the loss of membrane integrity (Brown and Lemmon, 2007). From this point onward, carotenoids may be accessed by XAT-7Hch, resulting in the production of carotenoid esters.
As it has been found for most fruits, the formation of xanthophyll esters during fruit ripening is characterized for a lack of correlation between the fatty acid profile of the total lipid pool and the major fatty acids involved in the esterification, suggesting that the responsible enzyme is rather specific toward the acyl group (Hornero-Méndez, 2019). Regarding XAT-7D specificity, Watkins et al. (2019) reported a good correlation between the fatty acid abundance in the endosperm and the fatty acids esterified to lutein in endosperm in wheat, which would indicate that XAT from bread wheat has low discrimination for the fatty acids. In the case of tritordeum, palmitic acid (16:0) and linoleic acid (18:2) are the fatty acids involved in the esterification of lutein, with a general preference for palmitic acid (16:0) being observed (Mellado-Ortega and Hornero-Mendez, 2012; Mellado-Ortega et al., 2015; Mattera et al., 2017) despite linoleic acid (18:2) being the most abundant in tritordeum endosperm (Mellado-Ortega and Hornero-Méndez, 2018) as it happens in bread wheat (Ziegler et al., 2015) and other Triticeae species. This suggests a higher specificity of XAT-7Hch compared to XAT-7D from bread wheat. Indeed, in a wide survey for carotenoid esterification patterns in H. chilense, the majority of accessions showed preference toward palmitic acid or similar preference toward palmitic and linoleic acids, while only three accessions used linoleic as the preferred fatty acid for lutein esterification (Ávila et al., 2019a).
In conclusion, a set of H. chilense genes putatively involved in xanthophylls esterification were successfully identified, and its expression during grain development in tritordeum was confirmed. Among them, HORCH7HG021460, referred as XAT-7Hch, seems to be responsible for the esterification of lutein in H. chilense. Thus, lutein esterification in H. chilense and tritordeum would be regulated by the same enzymatic system than in bread wheat suggesting the existence of a common mechanism for lutein esterification in the Triticeae. The transference of XAT-7Hch to durum or bread wheat may allow the improvement of their nutritional value by increasing the carotenoid accumulation and improving their stability and retention during post-harvest storage of grains. The SNP marker developed in this work is useful for a marker-assisted selection (MAS) strategy with this purpose. Further studies will be needed for determining the specificity of XAT-7Hch and its potential application in tritordeum breeding programs.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: GenBank accessions MT880902, MT880903, and MT880904.
Author Contributions
SA and DH-M conceived the idea. CR-S and SA developed the field experiment and collected the field samples. MR-R and CR-S developed the molecular studies. All authors contributed to the analysis and interpretation of the results, the redaction of the manuscript draft and approved the final version.
Funding
This research was funded by grant AGL2017-85368-P from Ministerio de Ciencia, Innovación y Universidades, and co-funded by FEDER. MR-R was recipient of a fellowship from Ministerio de Ciencia, Innovación y Universidades (PRE2018-084037). DH-M is member of CaRed Network, funded by Ministerio de Economía y Competitividad (BIO2017-90877-REDT) and participated in EUROCAROTEN COSTAction (CA15136). SA is member of FiRCMe Network, funded by Ministerio de Economía, Industria y Competitividad (AGL2016-81855-REDT).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
We acknowledge the support of the publication fee by the CSIC Open Access Publication Support Initiative through its Unit of Information Resources for Research (URICI).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2020.592515/full#supplementary-material
Footnotes
2. https://www.tritordeum.com/ww/?lang=en, last accessed 11/11/2020.
References
Ahmad, F. T., Asenstorfer, R. E., Soriano, I. R., and Mares, D. J. (2013). Effect of temperature on lutein esterification and lutein stability in wheat grain. J. Cereal Sci. 58, 408–413. doi: 10.1016/j.jcs.2013.08.004
Ahmad, F. T., Mather, D. E., Law, H. -Y., Li, M., Yousif, S. -J., Chalmers, K. J., et al. (2015). Genetic control of lutein esterification in wheat (Triticum aestivum L.) grain. J. Cereal Sci. 64, 109–115. doi: 10.1016/j.jcs.2015.05.007
Akoh, C. C., Lee, G. -C., Liaw, Y. -C., Huang, T. -H., and Shaw, J. -F. (2004). GDSL family of serine esterases/lipases. Prog. Lipid Res. 43, 534–552. doi: 10.1016/j.plipres.2004.09.002
Alvarez, J. B., Castellano, L., Recio, R., and Cabrera, A. (2019). Wx gene in Hordeum chilense: chromosomal location and characterisation of the allelic variation in the two main ecotypes of the species. Agronomy 9:261. doi: 10.3390/agronomy9050261
Atienza, S. G., Avila, C. M., and Martín, A. (2007a). The development of a PCR-based marker for PSY1 from Hordeum chilense, a candidate gene for carotenoid content accumulation in tritordeum seeds. Aust. J. Agric. Res. 58, 767–773. doi: 10.1071/AR06338
Atienza, S. G., Ballesteros, J., Martín, A., and Hornero-Méndez, D. (2007b). Genetic variability of carotenoid concentration and degree of esterification among tritordeum (xTritordeum Ascherson et Graebner) and durum wheat accessions. J. Agric. Food Chem. 55, 4244–4251. doi: 10.1021/jf070342p
Ávila, C. M., Mattera, M. G., Rodríguez-Suárez, C., Palomino, C., Ramírez, M. C., Martin, A., et al. (2019a). Diversification of seed carotenoid content and profile in wild barley (Hordeum chilense Roem. et Schultz.) and Hordeum vulgare L.–H. chilense synteny as revealed by DArTSeq markers. Euphytica 215:45. doi: 10.1007/s10681-019-2369-6
Ávila, C. M., Palomino, M. C., Hornero-Méndez, D., and Atienza, S. G. (2019b). Identification of candidate genes for lutein esterification in common wheat (Triticum aestivum) using physical mapping and genomics tools. Crop Pasture Sci. 70, 567–574. doi: 10.1071/CP18531
Ballesteros, J., Ramirez, M. C., Martinez, C., Atienza, S. G., and Martin, A. (2005). Registration of HT621, a high carotenoid content tritordeum germplasm line. Crop Sci. 45, 2662–2663. doi: 10.2135/cropsci2004.0641
Borrill, P., Ramirez-Gonzalez, R., and Uauy, C. (2016). expVIP: a customizable RNA-seq data analysis and visualization platform. Plant Physiol. 170:2172. doi: 10.1104/pp.15.01667
Britton, G. (2008). “Functions of intact carotenoids” in Carotenoids. eds. G. Britton, S. Liaaen-Jensen, and H. Pfander (Basel, Boston, Berlin: Birkhäuser Verlag), 189–212.
Brown, R. C., and Lemmon, B. E. (2007). “The developmental biology of cereal endosperm” in Endosperm: Developmental and molecular biology. ed. O. -A. Olsen (Berlin, Heidelberg: Springer Berlin Heidelberg), 1–20.
Collins, A., and Ke, X. (2012). Primer1: primer design web service for tetra-Primer ARMS-PCR. Open Bioinforma. J. 6, 55–58. doi: 10.2174/1875036201206010055
Fernández-García, E., Carvajal-Lérida, I., Jarén-Galán, M., Garrido-Fernández, J., Pérez-Gálvez, A., and Hornero-Méndez, D. (2012). Carotenoids bioavailability from foods: from plant pigments to efficient biological activities. Food Res. Int. 46, 438–450. doi: 10.1016/j.foodres.2011.06.007
Ficco, D. B. M., Mastrangelo, A. M., Trono, D., Borrelli, G. M., De Vita, P., Fares, C., et al. (2014). The colours of durum wheat: a review. Crop Pasture Sci. 65, 1–15. doi: 10.1071/CP13293
Gil-Humanes, J., Piston, F., Martin, A., and Barro, F. (2009). Comparative genomic analysis and expression of the APETALA2-like genes from barley, wheat, and barley-wheat amphiploids. BMC Plant Biol. 9:66. doi: 10.1186/1471-2229-9-66
Giménez, M. J., Pistón, F., and Atienza, S. G. (2011). Identification of suitable reference genes for normalization of qPCR data in comparative transcriptomics analyses in the Triticeae. Planta 233, 163–173. doi: 10.1007/s00425-010-1290-y
Giuliano, G., Tavazza, R., Diretto, G., Beyer, P., and Taylor, M. A. (2008). Metabolic engineering of carotenoid biosynthesis in plants. Trends Biotechnol. 26, 139–145. doi: 10.1016/j.tibtech.2007.12.003
Hornero-Méndez, D. (2019). “Occurrence of carotenoids esters in foods” in Carotenoid Esters in Foods: Physical, Chemical and Biological Properties ed. A. Z. Mercadante (Royal Society of Chemistry). 182–284.
Howe, K. L., Contreras-Moreira, B., De Silva, N., Maslen, G., Akanni, W., Allen, J., et al. (2019). Ensembl Genomes 2020—enabling non-vertebrate genomic research. Nucleic Acids Res. 48, D689–D695. doi: 10.1093/nar/gkz890
Johnson, E. J. (2014). Role of lutein and zeaxanthin in visual and cognitive function throughout the lifespan. Nutr. Rev. 72, 605–612. doi: 10.1111/nure.12133
Krinsky, N. I., Landrum, J., and Bone, R. A. (2003). Biologic mechanisms of the protective role of lutein and zeaxanthin in the eye. Ann. Rev. Nutr. 23, 171–201. doi: 10.1146/annurev.nutr.23.011702.073307
Li, F., Tsfadia, O., and Wurtzel, E. T. (2009). The phytoene synthase gene family in the grasses: subfunctionalization provides tissue-specific control of carotenogenesis. Plant Signal. Behav. 4, 208–211. doi: 10.4161/psb.4.3.7798
Lu, S., Wang, J., Chitsaz, F., Derbyshire, M. K., Geer, R. C., Gonzales, N. R., et al. (2019). CDD/SPARCLE: the conserved domain database in 2020. Nucleic Acids Res. 48, D265–D268. doi: 10.1093/nar/gkz991
Madeira, F., Park, Y. M., Lee, J., Buso, N., Gur, T., Madhusoodanan, N., et al. (2019). The EMBL-EBI search and sequence analysis tools APIs in 2019. Nucleic Acids Res. 47, W636–W641. doi: 10.1093/nar/gkz268
Martín, A., Alvarez, J. B., Martín, L. M., Barro, F., and Ballesteros, J. (1999). The development of tritordeum: a novel cereal for food processing. J. Cereal Sci. 30, 85–95. doi: 10.1006/jcrs.1998.0235
Mattera, M. G., Cabrera, A., Hornero-Méndez, D., and Atienza, S. G. (2015). Lutein esterification in wheat endosperm is controlled by the homoeologous group 7, and is increased by the simultaneous presence of chromosomes 7D and 7Hch from Hordeum chilense. Crop Pasture Sci. 66, 912–921. doi: 10.1071/CP15091
Mattera, M. G., Hornero-Méndez, D., and Atienza, S. G. (2017). Lutein ester profile in wheat and tritordeum can be modulated by temperature: evidences for regioselectivity and fatty acid preferential of enzymes encoded by genes on chromosomes 7D and 7Hch. Food Chem. 219, 199–206. doi: 10.1016/j.foodchem.2016.09.133
Mattera, M. G., Hornero-Méndez, D., and Atienza, S. G. (2020). Carotenoid content in tritordeum is not primarily associated with esterification during grain development. Food Chem. 310:125847. doi: 10.1016/j.foodchem.2019.125847
Mayer, K. F. X., Martis, M., Hedley, P. E., Šimková, H., Liu, H., Morris, J. A., et al. (2011). Unlocking the barley genome by chromosomal and comparative genomics. Plant Cell 23, 1249–1263. doi: 10.1105/tpc.110.082537
Meléndez-Martínez, A. J., Mapelli-Brahm, P., Hornero-Méndez, D., and Vicario, I. M. (2019). “CHAPTER 1: structures, nomenclature and general chemistry of carotenoids and their esters” in Food chemistry, function and analysis. ed. A. Z. Mercadante (London: Royal Society of Chemistry), 3–50.
Mellado-Ortega, E., Atienza, S. G., and Hornero-Méndez, D. (2015). Carotenoid evolution during postharvest storage of durum wheat (Triticum turgidum conv. durum) and tritordeum (×Tritordeum Ascherson et Graebner) grains. J. Cereal Sci. 62, 134–142. doi: 10.1016/j.jcs.2015.01.006
Mellado-Ortega, E., and Hornero-Mendez, D. (2012). Isolation and identification of lutein esters, including their regioisomers, in tritordeum (xTritordeum Ascherson et Graebner) grains: evidence for a preferential xanthophyll acyltransferase activity. Food Chem. 135, 1344–1352. doi: 10.1016/j.foodchem.2012.05.046
Mellado-Ortega, E., and Hornero-Méndez, D. (2015a). Carotenoid profiling of Hordeum chilense grains: the parental proof for the origin of the high carotenoid content and esterification pattern of tritordeum. J. Cereal Sci. 62, 15–21. doi: 10.1016/j.jcs.2014.12.005
Mellado-Ortega, E., and Hornero-Méndez, D. (2015b). Carotenoids in cereals: an ancient resource with present and future applications. Phytochem. Rev. 15, 873–890. doi: 10.1007/s11101-015-9423-3
Mellado-Ortega, E., and Hornero-Méndez, D. (2018). Effect of lutein esterification on the differential distribution of carotenoids in germ and endosperm fractions from tritordeum grains. J. Cereal Sci. 79, 462–468. doi: 10.1016/j.jcs.2017.12.006
Menkir, A., Palacios-Rojas, N., Alamu, O., Dias Paes, M. C., Dhliwayo, T., Maziya-dixon, B., et al. (2018). Vitamin A-biofortified maize: exploiting native genetic variation for nutrient enrichment. Sci. Br. Biofortification No. 2 (February 2018). CIMMYT, IITA, EMBRAPA, Harvest. Crop Turst. Bonn, Ger.
Murray, Y. H. G., and Thompson, W. F. (1980). Rapid isolation of high molecular weight plant DNA. Nucleic Acids Res. 8, 4321–4326. doi: 10.1093/nar/8.19.4321
Nishino, H., Murakoshi, M., Tokuda, H., and Satomi, Y. (2009). Cancer prevention by carotenoids. Arch. Biochem. Biophys. 483, 165–168. doi: 10.1016/j.abb.2008.09.011
Paine, J. A., Shipton, C. A., Chaggar, S., Howells, R. M., Kennedy, M. J., Vernon, G., et al. (2005). Improving the nutritional value of Golden Rice through increased pro-vitamin A content. Nat. Biotechnol. 23, 482–487. doi: 10.1038/nbt1082
Paznocht, L., Kotíková, Z., Orsák, M., Lachman, J., and Martinek, P. (2019). Carotenoid changes of colored-grain wheat flours during bun-making. Food Chem. 277, 725–734. doi: 10.1016/j.foodchem.2018.11.019
Petersen, T. N., Brunak, S., von Heijne, G., and Nielsen, H. (2011). SignalP 4.0: discriminating signal peptides from transmembrane regions. Nat. Methods 8, 785–786. doi: 10.1038/nmeth.1701
Pujadas-Salvá, A. J. (2016). × Tritordeum martinii A. Pujadas (Poaceae) nothosp. nov. × Tritordeum martinii A. Pujadas (Poaceae) nothosp. nov. Acta Bot. Malacit. 41, 325–328. doi: 10.24310/abm.v41i0.2486
Ramírez-González, R. H., Borrill, P., Lang, D., Harrington, S. A., Brinton, J., Venturini, L., et al. (2018). The transcriptional landscape of polyploid wheat. Science 361:eaar6089. doi: 10.1126/science.aar6089
Rodriguez-Concepcion, M., Avalos, J., Bonet, M. L., Boronat, A., Gomez-Gomez, L., Hornero-Mendez, D., et al. (2018). A global perspective on carotenoids: metabolism, biotechnology, and benefits for nutrition and health. Prog. Lipid Res. 70, 62–93. doi: 10.1016/j.plipres.2018.04.004
Rodríguez-Suárez, C., and Atienza, S. G. G. (2012). Hordeum chilense genome, a useful tool to investigate the endosperm yellow pigment content in the Triticeae. BMC Plant Biol. 12:200. doi: 10.1186/1471-2229-12-200
Rodríguez-Suárez, C., and Atienza, S. G. (2014). Polyphenol oxidase genes in Hordeum chilense and implications in tritordeum breeding. Mol. Breed. 34, 1867–1877. doi: 10.1007/s11032-014-0145-9
Rodríguez-Suárez, C., Atienza, S. G., and Pistón, F. (2011). Allelic variation, alternative splicing and expression analysis of Psy1 gene in Hordeum chilense Roem. et Schult. PLoS One 6:e19885. doi: 10.1371/journal.pone.0019885
Rodríguez-Suárez, C., Bagnaresi, P., Cattivelli, L., Pistón, F., Castillo, A., Martín, A. C., et al. (2020). Transcriptomics, chromosome engineering and mapping identify a restorer-of-fertility region in the CMS wheat system msH1. Theor. Appl. Genet. 133, 283–295. doi: 10.1007/s00122-019-03457-3
Rodríguez-Suárez, C., Giménez, M. J., and Atienza, S. G. (2010). Progress and perspectives for carotenoid accumulation in selected Triticeae species. Crop Pasture Sci. 61, 743–751. doi: 10.1071/CP10025
Rodríguez-Suárez, C., Mellado-Ortega, E., Hornero-Méndez, D., and Atienza, S. G. (2014). Increase in transcript accumulation of Psy1 and e-Lcy genes in grain development is associated with differences in seed carotenoid content between durum wheat and tritordeum. Plant Mol. Biol. 84, 659–673. doi: 10.1007/s11103-013-0160-y
Ruijter, J. M., Ramakers, C., Hoogaars, W. M. H., Karlen, Y., Bakker, O., van den Hoff, M. J. B., et al. (2009). Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. Nucleic Acids Res. 37:e45. doi: 10.1093/nar/gkp045
Shumskaya, M., Bradbury, L. M. T., Monaco, R. R., and Wurtzel, E. T. (2012). Plastid localization of the key carotenoid enzyme phytoene synthase is altered by isozyme, allelic variation, and activity. Plant Cell 24, 3725–3741. doi: 10.1105/tpc.112.104174
Vandesompele, J., De Preter, K., Pattyn, F., Poppe, B., Van Roy, N., De Paepe, A., et al. (2002). Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 3:research0034. doi: 10.1186/gb-2002-3-7-research0034
Vaz Patto, M. C., Aardse, A., Buntjer, J., Rubiales, D., Martín, A., and Niks, R. E. (2001). Morphology and AFLP markers suggest three Hordeum chilense ecotypes that differ in avoidance to rust fungi. Can. J. Bot. 79, 204–213. doi: 10.1139/cjb-79-2-204
Watkins, J. L., Li, M., McQuinn, R. P., Chan, K. X., McFarlane, H. E., Ermakova, M., et al. (2019). A GDSL esterase/lipase catalyzes the esterification of lutein in bread wheat. Plant Cell 31, 3092–3112. doi: 10.1105/tpc.19.00272
WHO (2009). Global prevalence of vitamin A deficiency in populations at risk 1995–2005. WHO Glob. database Vitam. A Defic.
Ye, J., Coulouris, G., Zaretskaya, I., Cutcutache, I., Rozen, S., and Madden, T. L. (2012). Primer-BLAST: a tool to design target-specific primers for polymerase chain reaction. BMC Bioinformatics 13:134. doi: 10.1186/1471-2105-13-134
Zadoks, J. C., Chang, T. T., and Konzak, C. F. (1974). A decimal code for the growth stages of cereals. Weed Res. 14, 415–421. doi: 10.1111/j.1365-3180.1974.tb01084.x
Zhu, C., Farré, G., Zanga, D., Lloveras, J., Michelena, A., Ferrio, J. P., et al. (2018). High-carotenoid maize: development of plant biotechnology prototypes for human and animal health and nutrition. Phytochem. Rev. 17, 195–209. doi: 10.1007/s11101-017-9506-4
Keywords: carotenoids, H. chilense, lutein esters, tritordeum, XAT, xanthophyll esterification
Citation: Requena-Ramírez MD, Atienza SG, Hornero-Méndez D and Rodríguez-Suárez C (2020) Mediation of a GDSL Esterase/Lipase in Carotenoid Esterification in Tritordeum Suggests a Common Mechanism of Carotenoid Esterification in Triticeae Species. Front. Plant Sci. 11:592515. doi: 10.3389/fpls.2020.592515
Edited by:
Petr Smýkal, Palacký University, Olomouc, CzechiaReviewed by:
Andrés Ricardo Schwember, Pontificia Universidad Católica de Chile, ChileR. Varma Penmetsa, University of California, Davis, United States
Copyright © 2020 Requena-Ramírez, Atienza, Hornero-Méndez and Rodríguez-Suárez. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Sergio G. Atienza, c2dhdGllbnphQGlhcy5jc2ljLmVz, orcid.org/0000-0001-5175-8303