Deciphering the genomic structure, function and evolution of carotenogenesis related phytoene synthases in grasses
© Dibari et al; licensee BioMed Central Ltd. 2012
Received: 9 February 2012
Accepted: 6 June 2012
Published: 6 June 2012
Carotenoids are isoprenoid pigments, essential for photosynthesis and photoprotection in plants. The enzyme phytoene synthase (PSY) plays an essential role in mediating condensation of two geranylgeranyl diphosphate molecules, the first committed step in carotenogenesis. PSY are nuclear enzymes encoded by a small gene family consisting of three paralogous genes (PSY 1-3) that have been widely characterized in rice, maize and sorghum.
In wheat, for which yellow pigment content is extremely important for flour colour, only PSY 1 has been extensively studied because of its association with QTLs reported for yellow pigment whereas PSY 2 has been partially characterized. Here, we report the isolation of bread wheat PSY 3 genes from a Renan BAC library using Brachypodium as a model genome for the Triticeae to develop Conserved Orthologous Set markers prior to gene cloning and sequencing. Wheat PSY 3 homoeologous genes were sequenced and annotated, unravelling their novel structure associated with intron-loss events and consequent exonic fusions. A wheat PSY 3 promoter region was also investigated for the presence of cis-acting elements involved in the response to abscisic acid (ABA), since carotenoids also play an important role as precursors of signalling molecules devoted to plant development and biotic/abiotic stress responses. Expression of wheat PSYs in leaves and roots was investigated during ABA treatment to confirm the up-regulation of PSY 3 during abiotic stress.
We investigated the structural and functional determinisms of PSY genes in wheat. More generally, among eudicots and monocots, the PSY gene family was found to be associated with differences in gene copy numbers, allowing us to propose an evolutionary model for the entire PSY gene family in Grasses.
KeywordsCarotenoids Phytoene synthase Wheat Intron loss Abiotic stress Evolution
Carotenoids represent a diverse group of pigments found in bacteria, fungi and plants . In plants, carotenoids have several important functions, essential for plant development as: (i) they serve as accessory pigments to harvest light for photosynthesis and constitute the basic structural units of photosynthetic apparatus and (ii) they also act as photoprotectors for plants to adapt to high light stress . In addition, oxidative cleavage of carotenoids produces apocarotenoids considered as (i) signals in plant development, (ii) antifungal agents, (iii) contributors to flower and fruit flavours and aromas . The phyto-hormone abscisic acid, the most well known apocarotenoid, plays an important role in plant response to stress conditions . In addition, carotenoids have long been considered as essential nutrients as well as important health beneficial compounds. In animals including humans, unable to synthesize carotenoids de novo, dietary carotenoids are essential precursors of vitamin A and retinoid compounds needed in development [4, 5].
PSY genes are of primary importance in wheat due to their impact on flour colour, so that numerous genetics studies have been conducted on yellow pigment (YP) content in wheat grain . A high YP concentration is a desirable trait in durum wheat and is the target of breeding programs worldwide . However, low YP is a breeding target in hexaploid wheat, for which bright white flour is preferred. The inheritance of YP is complex, controlled by additive gene effects and highly heritable . Several QTLs for YP have been identified in both durum and bread wheat on chromosomes 1A , 1B [17, 18], 3A , 3B [16, 20], 4A and 5A , 2A, 4B and 6B , 5B , and 6A . However, the majority of studies agreed that chromosome group 7 has the greatest impact in controlling YP content in bread wheat [13, 16–19, 22–24]. Elouafi et al. identified two minor QTLs on chromosome 7A and a third major QTL on chromosome 7BL. This last QTL has been corroborated in two other durum wheat populations [13, 24] and is ortholog to the one identified in hexaploid wheat on the same chromosome . A QTL for YP in durum wheat was identified by Patil et al. on chromosome 7A that accounted for up to 55% of the trait variation. Mares and Campbell  identified QTLs on chromosome 7A in two hexaploid wheat populations that had a large impact on flour yellowness, and correspond to loci also associated with xanthophyll concentration. Cloning of genes controlling flour pigmentation would be of great help in breeding for genotypes with appropriate levels of yellow pigments in durum and bread wheat. Using a PCR-based approach Cenci et al. identified and mapped clones from a durum wheat BAC library thought to contain genes coding for important enzymes involved in carotenogenesis. In their study, putative phytoene synthase, phytoene desaturase and ζ-carotene desaturase clones were assigned to the chromosome groups 5, 4 and 2, respectively. To date, four PSY genes have been identified in durum wheat, forming two paralogous series on chromosome groups 5 (Psy 2) and 7 (Psy 1) . One of these, Psy 1-B1, co-segregates with a QTL for YP on 7BL [13, 24]. Psy 1-A1 was localized on the distal end of the chromosome 7AL in one hexaploid wheat mapping population, associated with variation for flour colour . Wheat PSY 1 genomic sequences show similar exon-intron characteristic to PSY 1 known in other Grass species, with six exons separated by five introns . A partial sequence of PSY 2 was also obtained and it was mapped on the short arm of the chromosome group 5 , with no clear data regarding the impact of PSY 2 on carotenoid accumulation.
In this article, we report evidence that the wheat PSY gene family also consists of three paralogous genes and that PSY 3 is characterized by a novel gene structure, consequence of intron loss events. Since the wheat PSY 3 promoter region is characterized by cis- acting elements involved in plant response to abiotic stress, we studied changes in wheat PSY expression levels in leaves and roots during exogenous treatment with ABA. Finally, after comparative genomic investigation of PSY loci, we reconstructed the paleo-history for this gene family in Grasses identifying evolutionary events that have shaped PSY gene copies in modern species.
PSY3 gene sequencing and mapping in bread wheat
Identification and characterization of PSY3 loci
BAC sequencing and annotation (cf Materials and Methods) made it possible to identify PSY 3 orthologous regions. We selected ungapped contigs of 25 000 bp, 53 411 bp and 46 803 bp harbouring, respectively, PSY 3A, PSY 3B and PSY 3D genes ( Additional file 1: S1S2S3). Figure 2b illustrates BAC annotation results as well as micro-colinearity with rice, Brachypodium, maize and sorghum genome sequences, confirming that PSY 3 and Bradi4g37530.1 genes are at orthologous positions. In the PSY 3 loci we found the insertion of a non-syntenic gene (orthologous to Bradi3g60840). According to paleogenomic data in monocots , this gene is located on the syntenic chromosomes 2, 4 and 3 of, rice, sorghum and Brachypodium genomes respectively. In wheat, this gene would be expected to map on the orthologous chromosome group 6 and not group 5, where it was found based on the PSY sequenced region assignation. We can thus suppose that this gene has been specifically transposed in wheat from an ancestral donor region (group 6) to the modern current acceptor region (group 5).
Despite intron/exon shuffling events reported previously, it was possible to predict protein structure. PSY 3A encodes a protein of 413 amino acids, longer than those coded by PSY 3B and PSY 3D homoeologs (401 and 403 amino acids, respectively) ( Additional file 1: S1S2S3). The NCBI Conserved Domain Search software allowed us to identify the trans-Isoprenyl Diphosphate Synthase (Trans-IPPS) domain which includes squalene and phytoene synthases. According to the previous domain characterization, wheat PSY 3 was classed in the IPSS superfamily and Class I terpene cyclise. The conserved trans-IPPS domain was identified at positions: aa123-386, aa112-375 and aa113-376, respectively, for PSY 3A, B and D copies. Protein motifs identified based on ScanProsite are available in Additional file 1: Figure S1. ChloroP1.1  identified transit peptides with 56 and 57 residues for PSY 3B and PSY 3D, but only 6-residues for PSY 3A and a signal peptide of two amino acids for Brachypodium PSY 3. However, according to ProtComp 9.0  prediction, wheat PSY 3 may have a putative chloroplast localization. We also analyzed sequence homology between wheat PSY 3 proteins: PSY 3A and PSY 3B showed 91%, PSY 3A and PSY 3D 90% and PSY 3B and PSY 3D 96% amino acid sequence identity.
Evolutionary history of PSY gene family
PSY genes expression profile and function investigation in bread wheat
Since there is no evidence of PSY 3 expression in the literature, we used the public wheat Affymetrix chip data  based on array hybridization with RNA extracted from wheat grains at different development stages (100 Degree Days (DD), 200DD, 250DD, 300DD, 500DD), from cellular division to starch filling, to unravel the expression of PSY genes in wheat. We blast-searched the PSY 3 sequences obtained against Affymetrix chip sequences and found two probes, Ta.20776.1.S1_at and Ta.18880.1.S1_at, homoeologs of PSY 1 and PSY 2, respectively. The expression level of PSY 1 in grains was twice that of PSY 2 during all developmental stages except at 500DD when PSY 1 levels decreased ( Additional file 1: Figure S2). PSY 2 expression was constant during grain development, whereas PSY 1 expression exhibited significant differences between 350DD and 500DD (p-value = 2.22E-10, using the Bonferroni test). Our data showed that PSY 1 increased at 200DD while its level remained constant at 300DD and then decreased at 500DD. No evidence was found for PSY 3 expression in grain at these developmental stages. To gain insight into PSY gene expression during grain development, we used Single-Strand Conformation Polymorphism (SSCP) analysis  to screen wheat RNAs (cv Recital) extracted from leaves, roots and grains at different development stages. SSCP offered in this experiment the opportunity to (i) detect minor sequence changes in polymerase chain reaction-amplified samples and (ii) evaluate the expression level of three homoeologous copies of each PSY gene in the same running time . In the same PCR reaction, we amplified specifically one paralogous gene and its homoeologous copies that were successfully separated, by capillary electrophoresis according to their sequence differences and assigned to each genome, using amplification profiling approaches on nullitetrasomic lines as references . PSY 1 expression was evaluated based on the primer pair PSY 1_F1R1 (amplification product of 218 bp) designed on the fourth exon. Its expression in developing grain (100DD, 200DD, 250DD, 300DD and 500DD) was constant and PSY 1-B showed the highest transcript level. PSY 1 expression profile did not change in leaves and roots where PSY 1-B was always the copy most expressed. We also amplified PSY 2 in each tissue using the primer pair PSY 2F1R1, designed on the fourth exon (amplification product of 141 bp). PSY 2 showed a constant expression pattern in developing grains, the B and D copies showing the highest transcript levels. In leaves, all three homoeolog copies presented the same expression profiles; whereas, in roots, the D transcript decreased compared to B and A, whose levels remain unchanged. Using this SSCP approach we were able to detect PSY 3 expression using the primer pair PSY 3F2R1, designed on the third exon (amplification product of 216 bp). During grain development at 100DD the D copy showed the highest transcript level, while, at 200DD, amplification intensity of the B copy increased in contrast with the two other homoeologs, which decreased. For the remaining stages, expression levels of the three homoeologous copies were constant. In leaves, A and D copies had the highest expression levels, while in roots, the D copy showed a greater expression than the A and B copies ( Additional file 1: Figure S3).
To investigate changes in PSY expression during stress, we determined transcript levels in leaf and root samples after stress treatment. PSY responses to the different ABA concentrations tested were constant. PSY 1 mRNA levels, in stressed leaves, underwent to a minimum reduction during the first 2 h to then returned to the initial level at 5 h. PSY 2 transcripts slightly decreased in stressed leaves after 2 h and then remained constant. PSY 3 expression levels increased about 4-fold after 5 h of stress treatment and then returned to initial levels. In stressed roots, differences in expression levels became more obvious: PSY 1 mRNA levels doubled after 5 h of treatment, while PSY 2 transcript levels doubled after only 2 h of stress and then decreased gradually after 5 and 8 h of treatments to initial expression levels (Figure 5c-d-e). PSY 3 expression levels started to increase just 1 h after the ABA application to 18-fold after 2 h and then returned to lower levels (see Additional file 1: Table S2 for ANOVA analysis). Overall, Figure 5f clearly shows the effectiveness of ABA treatment on PSY 3 up-regulation in roots. Normalized expression data of PSY genes obtained in this experiment are available as Additional file 1: Figure S4.
Carotenoids are widely recognized as essential pigments in plants because of their role in development, photosynthesis, photoprotection and in response to stress since they are essential precursors of apocarotenoids, such as the abscisic acid (ABA). The PSY gene family, playing a central role in carotenoid biosynthesis, is characterized by duplicated gene copies. Their characterization in Grasses and the recent identification of a third copy in rice, maize and sorghum raises questions about PSY gene family organization in wheat. The role and function of PSY 1 have been largely investigated because of its association with major QTLs for carotenoid accumulation in wheat grain. Partial sequences of PSY 2 from several wheat varieties are available but little is known about its function and, up to now, there was no evidence of a PSY 3 homolog and associated activity in wheat. Our current analysis has demonstrated that wheat PSY gene family is characterized by the existence of PSY 1 to 3 gene copies. We report here the precise nature of gene family structure, evolution and function.
Because of the polyploid nature of wheat and the lack of a whole genome sequence, we identified orthologous regions in reference genomes such as Brachypodium and developed COS markers to isolate wheat PSY 3 from a BAC library. Wheat PSY 3 is characterized by a novel gene structure due to intron loss events leading to four exons and three introns. Rice and maize PSY 3, whose gene structure is highly conserved, are characterized by a third intron starting with a non-canonical 5′splicing site GC instead of GT. This type of intron is deleted by the same spliceosome mechanism that removes GT-AG introns, mainly involved in alternative splicing rather than in introns loss events . Brachypodium and sorghum show differences in PSY gene structure, characterized by five exons and four introns. The third intron has been deleted but the causal mechanism is not yet clear. The classical model to explain intron loss is based on recombination of a genomic copy of a gene with a reverse-transcribed copy of a spliced mRNA transcript, deleting one or more adjacent introns. A key aspect of this theory is that only multicellular eukaryotes that have genes transcribed in germ cells would be susceptible to mRNA-mediated intron loss [41–45]. We identified inverted and repeated short motifs at deletion breakpoints in the wheat PSY 3 sequence, corresponding to the 5′ and the 3′ splicing sites, which may have driven illegitimate recombination leading to the observed intron loss events. Brachypodium and sorghum PSY 3 sequences are not characterized by such repeated motifs at deletion breakpoint, suggesting that intron loss occurred by simple genomic deletion. Maize PSY3 intron 2 and rice intron 3 and maize intron 3 and rice intron 2 are similar in size, leading to the hypothesis that these introns may have been inverted but this could not be confirmed due to limited sequence similarities between introns. Howitt et al. described the presence within the intron #3 in wheat PSY-D1 of an insertion of 1200 bp of an inverted repeat motif. The characterization of either repeated motifs and/or TEs in this intron in Grasses may have driven intron loss or shuffling events leading to the distinct PSY 3 gene structures observed.
PSY duplicated genes exist throughout monocots and eudicots. Gallagher et al. proposed an evolutionary model for the Grass PSY gene family in which monocots and eudicots share the same ancestral gene and the same duplication event, giving rise to genes PSY 1 and PSY 2. Welsch and collaborators  proposed that rice PSY 3 originated from PSY 2 duplication or vice versa, because of their close similarity in the 5′UTR region. We analyzed homology between PSY genes in reference genomes and showed that PSY 3 had greater homology with PSY 1 sequence than PSY 2, suggesting that PSY 3 originated from a duplication of PSY 1. Furthermore, public paleogenomics data in Grasses indicated that PSY genes are not located on known duplicated chromosomes originating from the characterized ancestral WGD event reported by Murat et al.. We suggest, according to Jaillon et al. and Tang et al.[36, 37], that a more ancient WGD event, specific to monocots has shaped the PSY 3 copy as a duplication of PSY 1.
Based on a semi-quantitative approach, such as SSCP strategy, we observed PSY expression in wheat grains, leaves and roots. PSY 1 and PSY 2 homoeologous copies showed stable expression profiles in the different tissues. PSY-B1 was the copy most expressed, while, PSY-B2 and PSY-D2 showed the largest transcript levels in developing grains and, in leaves and roots, PSY-A2 transcripts increased. PSY 3 homoeologous copies showed a variable pattern of expression during the first stage of cellular division during early grain development and we found PSY-A3 and PSY-D3 most expressed in leaves and PSY-D3 in roots. However, with this approach, a statistically basis comparison was not possible between the expression levels of the three paralogous genes because it is semi-quantitative and only the broad transcript levels of homoeologous copies of each gene can be compared when amplified with the same primer pairs and in the same conditions. Access to recently published Affymetrix data  provided the opportunity to investigate PSY gene expression pattern and showed that PSY 1 expression level was higher than PSY 2 in developing grains. PSY 1 transcripts increased at 200DD to remain constant until 500DD, when they were reduced. In developing maize endosperm, Li et al. showed that during endosperm carotenogenesis, only PSY 1 transcript increased significantly, starting after 12DAP until 22DAP. Basis on this evidence, we can speculate that wheat PSY 1 in the endosperm follows the same kinetics as the maize PSY 1 ortholog and that apparently there is a correspondence in timing of expression between wheat and maize PSY 1 orthologs. Our results are consistent with the evidence that PSY 1 plays an important role in carotenoid accumulation in the endosperm and with its association with QTLs for YP on chromosome group 7 in wheat [13, 18, 24].
Gene expression is closely associated with cis- regulatory DNA sequences (known as cis-elements) located in its upstream regions . In rice and maize, PSY 3 expression is suggested to be strongly inducible in roots by drought and salt stress and application of ABA, since their promoter regions harbour cis-elements involved in plant response to abiotic stress [11, 39]. We identified several regulatory elements in the wheat PSY 3 promoter region, (i) the ABA-responsive element (ABRE) whose core is represented by the G-box, and (ii) a Coupling Element (CE1 or CE3) that is necessary to obtain ABRE functional [48, 49]. Through a Real-Time PCR approach, we quantified wheat PSY transcript levels in leaves and roots before and after the exogenous application of ABA, observing that PSY 2 represents the most abundant transcript in unstressed leaves, while PSY 1 transcript levels were 3-fold lower than in carotenoid accumulation in grains [13, 18, 24]. Wheat PSY 2 may thus play an important and complementary role in YP accumulation in photosynthetic tissues. PSY 3 expression levels were low in non-stressed leaves compared to PSY 1 and PSY 2 while its levels increased in non-stressed roots. PSY 1 and PSY 2 transcript levels did not show significant changes in leaves after ABA treatment. In stressed roots, PSY 1 and PSY 2 levels doubled after 5 h and 2 h of treatment respectively. PSY 3 showed the most intense response to ABA as in leaves, after 5 h of ABA application, PSY 3 transcript was 4-fold higher, while in stressed roots after 2 h its transcript levels increase about 18-fold. PSY 3 was the only gene that showed significant changes of transcript levels in the studied conditions. Carotenoids provide precursors for ABA which is the main apocarotenoid that plants can accumulate under drought conditions, ABA applications having been shown to affect plant growth and development, mimicking the effect of water stress and improving capacity to escape or tolerate such stress . Since PSY genes are considered as rate-limiting in carotenogenesis, we can explain the up-regulation of PSY 3 by ABA as the necessity of plants to control carotenoid flux to roots during stress conditions to provide ABA precursors and thus activate stress responses. In a previous study, Singh et al., studying the association of PSY 1-A1 allelic variation with YP colour in durum wheat, speculated about the possible existence of a third PSY on the chromosome 7AL to explain the presence of a second QTL near the marker Xwmc116. We identified the chromosomal bin localization for each homoeolog and the expected syntenic relationships were confirmed assigning PSY 3 to chromosome group 5 which has the highest concentration of QTLs and major loci controlling plant adaptation to environment, particularly, those controlling heading date, frost and salt tolerance, whereas a region with a crucial role in drought tolerance is located on chromosome group 7. Multiple-stress QTLs and linked markers have also been detected, suggesting the existence of common mechanisms for different stresses or of gene clusters controlling different stress tolerance processes . For example, in wheat a significant QTL for ABA accumulation in droughted leaves was found on chromosome 5A co-located with the major vernalization responsive gene, Vrn1 , and in the same chromosomal bin we mapped PSY-3A. Aprile et al. studied different physiological reactions to water stress and the substantially different molecular responses of bread and durum wheat genotypes demonstrating that there is a difference between durum and bread wheat response to water stress and that deletion of bin 5AL-10 from Chinese Spring determined a lower response to such stresses. They also showed that the different response by durum and bread wheat can be associated with the absence of D genome, where genes or factors that modulate expression of cluster of genes may be localized. PSY 3 localization and function in wheat are consistent with these data, especially the PSY 3D genomic structure, function and expression reported in the current study. Finally, we also showed here that PSY 3s are characterized by intron losses, inverted repeats and up-regulation by abiotic stresses. It is important to note that biotic/abiotic stresses could also activate transposons and retrotransposons  whose expression may have played a role in PSY 3 gene structure modification.
PSY gene family plays an important role in plant development (photoprotection and resistance to stress) as well as is a key contributor to flower and fruit flavours and aromas. The current analysis provides relevant structural (sequence, primers), functional (expression patterns) and evolutionary (gene conservation and specificity) informations that can be used as applied tools in wheat breeding programs.
PSY gene orthologs identification
The methodology used to reassess synteny between wheat/rice/Brachypodium/sorghum/maize genomes as well as the identification of intra-chromosomal duplications in wheat was previously described [28, 29, 54, 55]. Briefly, wheat [30, 31], rice (41046 genes), Brachypodium (25504 genes), sorghum (34008 genes) and maize (32540 genes) genomes were aligned with the publicly available PSY gene sequence  to identify orthologs. Three parameters were used to increase stringency and significance of BLAST sequence alignment by parsing BLASTN results and rebuilding HSPs (high scoring pairs) or pairwise sequence alignments. The first parameter, AL (aligned length), corresponds to the sum of all HSP lengths. The second, CIP (cumulative identity percentage) corresponds to the cumulative percent of sequence identity obtained for all the HSPs (CIP = ∑ nb ID by HSP/AL) x 100). The third parameter, CALP is the cumulative alignment length percentage. It represents the sum of the HSP lengths (AL) for all the HSPs divided by the length of the query sequence (CALP = AL/Query length). These parameters allow identification of the best alignment, with the highest cumulative percentage of identity in the longest cumulative length, taking into account all HSPs obtained for any pairwise alignment, regarding PSY sequences. These parameters were applied to all the BLAST alignments that were performed in the present study. Based on the genome-wide synteny analysis, gene relationships between species were then referenced as COS (for conserved gene pairs), CNV (for tandem duplicated genes) and PAV (for non –conserved genes).
BAC clone screening, sequencing and annotation
A BAC library from T. aestivum ‘Renan’  was used in this study, the library has a 6.9x genome coverage. It consists of 812,544 clones organized in 2592 384-well plates pooled in two steps. The first step aims at mixing the 384 clones of each plate; the second step organizes plate pools in two dimensional matrices. For each 384-well plate, 16-line pools and 24-column pools and a total of 40 pools were available. Overall, 274 line and column pools were available to identify BAC clones of interest through PCR screening using PSY primers. Each amplicon (or BAC of interest) in a 384-wells plate was identified by coordinates which identified positive line and column pools in a matrix. By crossing the coordinates of positive pools we determined the plate harbouring the clone of interest. The individual BAC clone of interest was then identified on specific plates by PCR screening. BAC sequencing was performed on the GS Junior (Roche) using the 454 sequencing technology according the procedure described by the manufacturer at the genotyping platform GENTYANE (UMR-1095 GDEC, Clermont-Ferrand).
Genes and repeated elements (TEs) were identified by computing and integrating results based on BLAST algorithms , predictor programs and software described below. Gene structure and putative functions were identified by combining results of BLASTN and BLASTX alignments against dbEST  and SwissProt databases  with the results of predictor program, FgeneSH  with default parameters. Known genes were named based on BLASTX results against protein with known functions (SwissProt). Promoter analysis was performed using both TSSP database and Toucan platform  to search for cis-elements. PolyA Miner Signal  was used in terminator analysis. Characterization of protein functional domain was performed using InterPro  and SMART database [64, 65], while protein motifs were analysed with ScanProsite . ChloroP1.1 and ProtComp 9.0 were used, respectively, to identify signal peptide and to predict protein subcellular localization.
Transposable elements (TEs) were detected by comparison with two databases of repetitive elements: TREP  and Repbase . TEs boundaries were identified with the REPET package . Insertion profile of TEs is identified using a modified version of svg_ltr.pl script .
Stress treatment experimental design and quantitative real-time PCR
Recital seeds were surface sterilized and embryos at 16 days after pollination (DAP) were collected and grown on Murashige and Skoog agar medium with 2% saccharose and indole-3-acetic acid (IAA). In each pot, three embryos were sown in order to obtain three biological replications for each condition. Plants were grown during the first 3 days in the dark at 25 °C and were transferred to continuous light at 22 °C at the two-leaf stage.
To carry out stress experiments, plants were treated with solution of ABA at different concentration: 50 μM, 100 μM and 150 μM. Leaf and root samples were collected after 0, 1, 2, 5 and 8 hours (h) of ABA treatment and were stored at −80 °C. RNA extraction was performed according to TRIzol® Reagent protocol. DNase set (Qiagen) and RNeasy Minelute Cleanup (Qiagen) kits were used for RNA purification. Integrity and quality of mRNA (contamination assessment of genomic DNA) were checked with the A260/A280 ratio with Nanodrop (Thermo Scientific) and on agarose with a classical PCR using a reference gene (see primers below). cDNA synthesis was carried out with a Transcriptor first strand cDNA synthesis kit (Roche) following the procedure described by the manufacturer with 1 μg of RNA (cycling conditions: 10′ at 25 °C/30′ at 55 °C/5′ at 85 °C°). Four primer pairs used for Real-Time PCR analysis ( Additional file 1: Table S1) were designed with primer3 in order to obtain exonic specific amplification of 100pb. Real-time PCR was performed using the LightCycler® 480 (Roche Diagnostics) and all reactions were performed using the LightCycler® 480 DNA SYBR Green I Master (Roche Diagnostics) in 10 μl with 4 ng of cDNA and 55 PCR cycles at 60 °C (hybridation temperature). All samples were prepared in three technical replicates and a negative control using water as a template was included. Quantification cycles (Cq) were analysed using the Light Cycler software version 1.5.0 and normalized with two reference genes (wheat spastin (primers: TGCCACTGCCTGATCCAAAT and AGCAAGCCTCTCCAGATCATG) and RNase L inhibitor-like (primers: TTGAGCAACTCATGGACCAG and GCTTTCCAAGGCACAAACAT) ) using ‘Advanced Relative Quantification’ module to obtain normalized ratio Et-CqT/Er-CqR (with CqT/CqR : Cycle number at target/reference detection threshold (crossing point) and ET/ER : Efficiency of target/reference amplification (10-1/slope)). Specificity of amplification was confirmed via melting curve analysis of final PCR products by increasing the temperature from 65 °C to 95 °C. PCR efficiency was calculated for each gene using a standard curve of serial dilutions and used in relative expression analysis. To calculate transcript abundance under stress conditions, the transcript levels were obtained for the controls, at each considered time point and ABA concentration and tissues. All observations were expressed as mean ± SD; the significance (P ≤ 0.05) of variables studied was assessed by one-way analysis of variance (ANOVA) using R software version 2.13.1.
This work was supported by grants from the Agence Nationale de la Recherche, Program ANRjc-PaleoCereal (ref: ANR-09-JCJC-0058-01) and Program ANR-Blanc-PAGE (ref: ANR-2011-BSV6-00801). The authors would like also to thank Felicity Vear (INRA Clermont-Ferrand, UMR 1095, France) for her fruitful contribution during the manuscript preparation.
- Britton G, Liaaen-Jensen S, Pfander H: Carotenoids Handbook. 2004, Birkhäuser Verlag, BaselView ArticleGoogle Scholar
- Lu S, Li L: Carotenoid Metabolism: Biosynthesis, Regulation and Beyond. J of Integrative Plant Biology. 2008, 50: 778-785. 10.1111/j.1744-7909.2008.00708.x.View ArticleGoogle Scholar
- Auldrige ME, McCarty DR, Klee HJ: Plant carotenoid cleavage oxygenases and their apocarotenoid products. Curr Opin Plant Biol. 2006, 9: 315-321. 10.1016/j.pbi.2006.03.005.View ArticleGoogle Scholar
- Fraser PD, Bramley PM: The biosynthesis and nutrional uses of carotenoids. Prog Lipid Res. 2004, 43: 228-265. 10.1016/j.plipres.2003.10.002.View ArticlePubMedGoogle Scholar
- Gallagher CE, Matthews PD, Li F, Wurtzel ET: Gene duplication in the carotenoid biosynthetic pathway preceded evolution of the Grasses. Plant Physiol. 2004, 135: 1776-1783. 10.1104/pp.104.039818.PubMed CentralView ArticlePubMedGoogle Scholar
- Cunningham FX, Gantt E: Gene and enzymes of carotenoid biosynthesis in plants. Annu RevPlant Physiol Plant Mol Biol. 1998, 49: 557-583. 10.1146/annurev.arplant.49.1.557.View ArticleGoogle Scholar
- Hirschberg J: Carotenoid biosynthesis in flowering plants. Curr Opin Plant Biol. 2001, 4: 210-218. 10.1016/S1369-5266(00)00163-1.View ArticlePubMedGoogle Scholar
- Ye X, Al-Babili S, Klöti A, Zhang J, Lucca P, Beyer P, Potrykus I: Engineering the pro-vitamin A (β-carotene) biosynthetic pathway into (carotenoid-free) rice endosperm. Science. 2000, 287: 303-305. 10.1126/science.287.5451.303.View ArticlePubMedGoogle Scholar
- Paine JA, Shipton CA, Chaggar S, Howells RM, Kennedy MJ, Vernon G, Wright SY, Hinchliffe E, Adams JL, Silverstone AL, Drake R: Improving the nutritional value of Golden Rice through increased pro-vitamin A content. Nat Biotechnol. 2005, 23: 429-430. 10.1038/nbt0405-429.View ArticleGoogle Scholar
- Al-Babili S, Hoa TT, Schaub P: Exploring the potential of the bacterial carotene desaturase CrtI to increase the beta-carotene content in Golden Rice. J Exp Bot. 2006, 57: 1007-1014. 10.1093/jxb/erj086.View ArticlePubMedGoogle Scholar
- Welsh R, Wüst F, Bär C, Al-babili S, Beyer P: A Third Phytoene Synthase Is Devoted to Abiotic Stress-Induced Abscisic Acid Formation in Rice and Defines Functional Diversification of Phytoene Synthase Genes. Plant Physiol. 2008, 147: 367-380. 10.1104/pp.108.117028.View ArticleGoogle Scholar
- Li F, Vallabhaneni R, Wurtzel ET: PSY3, a new member of the Phytoene Synthase gene famlily conserved in the Poaceae and regulator of abiotic stress – induced root carotenogenesis. Plant Physiol. 2008, 146: 1333-1345. 10.1104/pp.107.111120.PubMed CentralView ArticlePubMedGoogle Scholar
- Pozniak CJ, Knox RE, Clarke FR, Clarke JM: Identification of QTL and association of a phytoene synthase gene with endosperm colour in durum wheat. Theor Appl Genet. 2007, 114: 525-537. 10.1007/s00122-006-0453-5.View ArticlePubMedGoogle Scholar
- Troccoli A, Borrelli GM, De Vita P, Fares C, Di Fonzo N: Durum wheat quality: a multidisciplinary concept. J Cereal Sci. 2000, 32: 99-113. 10.1006/jcrs.2000.0322.View ArticleGoogle Scholar
- Singh A, Reimer S, Pozniak CJ, Clarke FR, Clarke JM, Knox RE, Singh AK: Allelic variation at Psy1-A1 and association with yellow pigment in durum wheat grain. Theor Appl Genet. 2009, 118: 1539-1548. 10.1007/s00122-009-1001-x.View ArticlePubMedGoogle Scholar
- Patil RM, Oak MD, Tamhankar SA, Sourdille P, Rao VS: Mapping and validation of a major QTL for yellow pigment content on 7AL in durum wheat. Mol Breed. 2008, 21: 485-496. 10.1007/s11032-007-9147-1.View ArticleGoogle Scholar
- Cervigni G, Zhang W, Picca A, Carrera A, Helguera M, Manthey F, Miranda R, Dubcosky J, Echenique V: QTL mapping for LOX activity and quality traits in durum wheat.In Proceedings of the 7thinternational wheat conference SAGPyA/INTA: 27 November-2 December 2005. 2005, Mar del Plata, Argentina, -Google Scholar
- He XY, Zhang YL, He ZH, Wu YP, Xiao YG, Ma CX, Xia XC: Characterization of Phytoene synthase 1 gene (Psy1) located on common wheat chromosome 7A and development of a functional marker. Theor Appl Genet. 2008, 116: 213-221. 10.1007/s00122-007-0660-8.View ArticlePubMedGoogle Scholar
- Parker GD, Chalmers KJ, Rathjien AJ, Langridge P: Mapping loci associated with flour colour in wheat. Theor Appl Genet. 1998, 97: 238-245. 10.1007/s001220050891.View ArticleGoogle Scholar
- Mares DJ, Campbell AW: Mapping components of flour and noodle colour in Australian wheat. Aust J Agric Res. 2001, 52: 1297-1309. 10.1071/AR01048.View ArticleGoogle Scholar
- Hessler TG, Thomson MJ, Benscher D, Nachit MM, Sorrells ME: Association of a lipoxygenase locus, Lpx-B1, with variation in lipoxygenase activity in durum wheat seeds. Crop Sci. 2002, 42: 1695-1700. 10.2135/cropsci2002.1695.View ArticleGoogle Scholar
- Elouafi I, Nachit MM, Martin LM: Identification of a microsatellite on chromosome 7B showing a strong linkage with yellow pigment in durum wheat. Hereditas. 2001, 135: 255-261.View ArticlePubMedGoogle Scholar
- Atienza SG, Ballesteros J, Martin A, Hornero-Mendes D: Genetic variability of carotenoid concentration and degree of esterification among Tritordeum (xTritordeum Ascherson et Graebner) and durum wheat accessions. Agric Food Chem. 2007, 55: 4244-4251. 10.1021/jf070342p.View ArticleGoogle Scholar
- Zhang W, Dubcovsky J: Association between allelic variation at the Phytoene synthase 1 gene and yellow pigment content in the wheat grain. Theor Appl Genet. 2008, 116: 635-645. 10.1007/s00122-007-0697-8.View ArticlePubMedGoogle Scholar
- Kuchel H, Langridge P, Mosionek L, Williams K, Jefferies SP: The genetic control of milling yield, dough rheology and baking quality of bread. Theor Appl Genet. 2006, 112: 1487-1495. 10.1007/s00122-006-0252-z.View ArticlePubMedGoogle Scholar
- Cenci A, Somma S, Chantret N, Dubrovsky J, Blanco A: PCR identification of durum wheat BAC clones containing genes coding for carotenoid biosynthesis enzymes and their chromosome localization. Genome. 2004, 47: 911-917. 10.1139/g04-033.View ArticlePubMedGoogle Scholar
- Quraishi UM, Abrouck M, Bolot S, Pont C, Throude M, Guilhot N, Confolent C, Bortolini F, Praud S, Murigneux A, Chrmet G, Salse J: Genomics in cereals: from genome-wide conserved orthologous set (COS) sequences to candidate genes for trait dissection. Funct Integr Genomics. 2009, 9: 473-484. 10.1007/s10142-009-0129-8.View ArticlePubMedGoogle Scholar
- Abrouk et al: Paleogenomics of plants: synteny – based modelling of extinct ancestors. Trends Plant Sci. 2010, 15: 479-487. 10.1016/j.tplants.2010.06.001.View ArticleGoogle Scholar
- Murat F, Xu JH, Tannier E, Abrouk M, Guilhot N, Pont C, Messing J, Salse J: Ancestral grass karyotype reconstruction unravels new mechanisms of genome shuffling as a source of plant evolution. Genome Res. 2010, 20: 1545-1557. 10.1101/gr.109744.110.PubMed CentralView ArticlePubMedGoogle Scholar
- Affymetrix: [http://www.affymetrix.com/estore/]
- Wheat Genomics. [http://www.cerealsdb.uk.net/]
- INRA – CNRGV. [http://cnrgv.toulouse.inra.fr/]
- Sears ER: Nullisomic – tetrasomic combination in hexaploid wheat. Chromosomes manipulations and plant genetics. Edited by: Riley R, Lewis KR. 1966, Oliver and Boyd, Edimburgh, 29-45.View ArticleGoogle Scholar
- ChloroP1.1 Server. [http://www.cbs.dtu.dk/services/ChloroP/]
- Softberry. [http://linux1.softberry.com/all.htm]
- Jaillon O: The grapevine genome sequence suggests ancestral hexaploidization in major angiosperm phyla. Nature. 2007, 449: 463-467. 10.1038/nature06148.View ArticlePubMedGoogle Scholar
- Tang H, Bowers JE, Wang X, Paterson AH: Angiosperm genome comparisons reveal early polyploidy in the monocot lineage. PNAS. 2010, 107: 472-477. 10.1073/pnas.0908007107.PubMed CentralView ArticlePubMedGoogle Scholar
- Pont C, Murat F, Confolent C, Balzergue S, Salse J: RNA-seq in grain unveils fate of neo-and paleopolyploidization events in bread wheat (Triticum aestivum L.). Genome Biol. 2011, 12 (12): R119-10.1186/gb-2011-12-12-r119.PubMed CentralView ArticlePubMedGoogle Scholar
- Li F, Tzfadia O, Wurtzel ET: The phytoene synthase gene family in the Grasses. Plant Signaling & Behavior. 2009, 3: 208-211.View ArticleGoogle Scholar
- Farrer T, Brock Roller A, Kent WJ, Zahler AM: Analysis of the role of Caenorhabditis elegans GC-AG introns in regulated splicing. Nucleic Acids Res. 2002, 30: 3360-3367. 10.1093/nar/gkf465.PubMed CentralView ArticlePubMedGoogle Scholar
- Roy SW, Gilbert W: The evolution of spliceosomal introns: patterns, puzzles and progress. Nature Review Genetics. 2006, 7: 211-221.Google Scholar
- Deng-Ke N, Wen-Ru H, Shu-Wei L: mRNA-Mediated Intron Losses: Evidence from Extraordinarily Large Exons. Mol Biol Evol. 2005, 22: 1475-1481. 10.1093/molbev/msi138.View ArticleGoogle Scholar
- Lecharny A, Boudet N, Gy I, Aubourg S, Kreis M: Introns in, intron out in plant gene families: a genomic approach of the dynamics of gene structure. J Structural and Functional Genomics. 2003, 3: 111-116. 10.1023/A:1022614001371.View ArticleGoogle Scholar
- Drouin G, Moniz de Sá: Loss of introns in the pollen-specific actin gene subfamily members of potato and tomato. J Mol Evol. 1997, 45: 509-513. 10.1007/PL00006255.View ArticleGoogle Scholar
- Frugoli JA, McPeek MA, Thomas TL, McClung CR: Intron loss and gain during evolution of the catalase gene family in Angiosperms. Genetics. 1998, 149: 355-365.PubMed CentralPubMedGoogle Scholar
- Howitt CA, Cavanagh CR, Bowerman AF, Cazzonelli C, Rampling L, Mimica JL, Pogson BJ: Alternative splicing, activation of cryptic exons and amino acid substitutions in carotenoid biosynthetic genes are associated with lutein accumulation in wheat endosperm. Funct Integr Genomics. 2009, 9: 363-376. 10.1007/s10142-009-0121-3.View ArticlePubMedGoogle Scholar
- Brivanlou A, Darnell J: Signal transduction and control of gene expression. Science. 2002, 295: 813-818. 10.1126/science.1066355.View ArticlePubMedGoogle Scholar
- Shen Q, Ho TH: Functional dissection of an abscissic acid (ABA) – inducible gene reveals two independent ABA – responsive complexes each containing a G – box and a novel cis - acting element. Plant Cell. 1995, 7: 295-307.PubMed CentralView ArticlePubMedGoogle Scholar
- Shen Q: Modular nature of abscissic acid (ABA) response complexes: composite promoter units that are necessary and sufficient for ABA induction of gene expression in barley. Plant Cell. 1996, 8: 1107-1119.PubMed CentralView ArticlePubMedGoogle Scholar
- Quarrie SA, Gulli M, Calestani C, Steed A, Marmiroli N: Location of a gene regulation drought-induced abscisic acid production on the long arm of chromosome 5A of wheat. Theor Appl Genetic. 1994, 89: 794-800.View ArticleGoogle Scholar
- Cattivelli L, Baldi P, Crosatti C, Di Fonzo N, Faccioli P, Grossi M, Mastrangelo AM, Pecchioni N, Stanca AM: Chromosome regions and stress-related sequences involved in resistance to abiotic stress in Triticeae. Plant Mol Biol. 2002, 48: 649-665. 10.1023/A:1014824404623.View ArticleGoogle Scholar
- Aprile A, Mastrangelo AM, De Leonardis AM, Galiba G, Roncaglia E, Ferrari F, De Bellis L, Turchi L, Giuliano G, Cattivelli L: Transciptional profiling in response to terminal frought stress reveals differential responses along the wheat genome. BMC Genomics. 2009, 10: 279-10.1186/1471-2164-10-279.PubMed CentralView ArticlePubMedGoogle Scholar
- Rapp RA, Wendel JF: Epigenetics and plant evolution. New Phytol. 2005, 168: 81-91. 10.1111/j.1469-8137.2005.01491.x.View ArticlePubMedGoogle Scholar
- Salse J, Abrouk M, Murat F, Quraishi UM, Feuillet C: Improved criteria and comparative genomics tool to provide new insights into grass paleogenomics. Briefings in Bioinformatics. 2009, 6: 619-630.View ArticleGoogle Scholar
- Bolot S, Abrouk M, Quraishi UM, Stein N, Messing J, Feuillet C, Salse J: The ‘inner circle’ of the cereal genomes. Current Opinion in Plant Biology. 2009, 12: 119-125. 10.1016/j.pbi.2008.10.011.View ArticlePubMedGoogle Scholar
- Phytozome. [http://www.phytozome.net/]
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic local alignment search tool. J Mol Biol. 1990, 215: 403-410.View ArticlePubMedGoogle Scholar
- NCBI. [http://www.ncbi.nlm.nih.gov/]
- Expasy. [http://expasy.org/sprot/]
- Salamov A, Solovyev V: Ab initio gene finding in Drosophila genomic DNA. Genome Res. 2000, 10: 516-522. 10.1101/gr.10.4.516.PubMed CentralView ArticlePubMedGoogle Scholar
- Toucan 2. [http://homes.esat.kuleuven.be/~saerts/software/toucan.php#screenshots]
- PolyA Miner Signal. [http://dnafsminer.bic.nus.edu.sg/PolyA.html]
- Hunter S, Apweiler R, Attwood TK, Bairoch A, Bateman A, Binns D, Bork P, Das U, Daugherty L, Duquenne L, Finn RD, Gough J, Haft D, Hulo N, Kahn D, Kelly E, Laugraud A, Letunic I, Lonsdale D, Lopez R, Madera M, Maslen J, McAnulla C, McDowall J, Mistry J, Mitchell A, Mulder N, Natale D, Orengo C, Quinn AF, Selengut JD, Sigrist CJ, Thimma M, Thomas PD, Valentin F, Wilson D, Wu CH, Yeast C: InterPro: the integrative protein signature database. Nucleic Acid Res. 2009, 37: 224-228. 10.1093/nar/gkn833.View ArticleGoogle Scholar
- Schultz J, Milpets F, Bork P, Ponting CP: SMART, a simple modular architecture research tool: identification of signaling domain. Proc Natl Acad Sci USA. 1998, 95: 5857-5864. 10.1073/pnas.95.11.5857.PubMed CentralView ArticlePubMedGoogle Scholar
- Letunic I, Doerks T, Bork P: SMART 6: recent uploads and developments. Nucleic Acid Res. 2009, 37: 229-232.View ArticleGoogle Scholar
- De Castro E, Sigrist CJA, Gattiker A, Bulliard V, Langendijk-Genevaux PS, Gasteiger E, Bairoch A, Hulo N: ScanProsite: detection of PROSITE signature matches and ProRule-associated functional residues in proteins. Nucleic Acids Res. 2006, 34: W362-W365. 10.1093/nar/gkl124.PubMed CentralView ArticlePubMedGoogle Scholar
- Graingenes: A database for Triticeae and Avena. [http://wheat.pw.usda.gov./ITMI/Repeats/]
- Jurka J: Repbase update: a database and an electronic journal of repetitive elements. Trends Genet. 2000, 9: 418-420.View ArticleGoogle Scholar
- Flutre T, Duprat E, Feuillet C, Quesneville H: Considering transposable element diversification in de novo annotation approaches. PLoS One. 2011, 6 (1): e16526-10.1371/journal.pone.0016526.PubMed CentralView ArticlePubMedGoogle Scholar
- Kronmiller BA, Wise RP: TEnest: Automated Chronological Annotation and Visualization of Nested Plant Transposable Elements. Plant Physiol. 2008, 146: 45-59.PubMed CentralView ArticlePubMedGoogle Scholar
- Giménez MJ, Pistón F, Atienza SG: Identification of suitable reference genes for normalization of qPCR data in comparative transcriptomics analyses in the Triticeae. Planta. 2011, 233: 163-173. 10.1007/s00425-010-1290-y.View ArticlePubMedGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.