Tracing 100 million years of grass genome evolutionary plasticity.

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

Recent progress in molecular genetics and omics-driven research in seed biology.

Elucidating the mechanisms that control seed development, metabolism, and physiology is a fundamental issue in biology. Michel Caboche had long been a catalyst for seed biology research in France up until his untimely passing away last year. To honour his memory, we have updated a review written under his coordination in 2010 entitled "Arabidopsis seed secrets unravelled after a decade of genetic and omics-driven research". This review encompassed different molecular aspects of seed development, reserve accumulation, dormancy and germination, that are studied in the lab created by M. Caboche. We have extended the scope of this review to highlight original experimental approaches implemented in the field over the past decade such as omics approaches aimed at investigating the control of gene expression, protein modifications, primary and specialized metabolites at the tissue or even cellular level, as well as seed biodiversity and the impact of the environment on seed quality.

L'élucidation des mécanismes qui contrôlent le développement, le métabolisme et la physiologie des graines est une question fondamentale en biologie. Michel Caboche a longtemps été un catalyseur de la recherche en biologie des graines en France jusqu'à son décès prématuré l'année dernière. Pour honorer sa mémoire, nous avons mis à jour une revue écrite sous sa coordination en 2010 intitulée « Arabidopsis seed secrets unravelled after a decade of genetic and omics-driven research ¼. Cette revue englobait différents aspects moléculaires du développement des graines, de l'accumulation des réserves, de la dormance et de la germination, qui sont étudiés dans le laboratoire créé par M. Caboche. Nous avons étendu la portée de cette revue pour mettre en évidence des approches expérimentales originales mises en œuvre dans le domaine au cours de la dernière décennie, telles que les approches omiques visant à étudier le contrôle de l'expression des gènes, les modifications des protéines, les métabolites primaires et spécialisés au niveau des tissus ou même des cellules, tout en tenant compte de la biodiversité des graines et de l'impact de l'environnement sur leur qualité.

Custom methods to identify conserved genetic modules applied to novel transcriptomic data from Amborella trichopoda.

We have devised a procedure for the inter-species comparison of transcriptomic data and used this procedure to reconstruct the expression dynamics of major genetic modules that were present at least 149 million years ago in the most recent common ancestor of living angiosperms. We began by using laser-assisted microdissection to generate novel transcriptomic data from female flower tissues of Amborella trichopoda, the likely sister to all other living angiosperms. We then employed a gene-expression clustering method, followed by a custom procedure to compare genetic modules on the basis of gene orthology between Amborella and the molecular-genetic model angiosperm Arabidopsis thaliana. Using this protocol, we succeeded in identifying nine major genetic modules that appear to have conserved their expression dynamics from an early stage in angiosperm evolution. The genes of these modules, representing over 5000 orthogroups, include around one third of those known to control female reproductive development in Arabidopsis. Our study constitutes a proof of concept for the comparison of transcriptomic data between widely diverged plant species and represents a first step in the large-scale analysis of gene expression dynamics in a macro-evolutionary context.

Genetic and Molecular Control of Somatic Embryogenesis.

Somatic embryogenesis is a method of asexual reproduction that can occur naturally in various plant species and is widely used for clonal propagation, transformation and regeneration of different crops. Somatic embryogenesis shares some developmental and physiological similarities with zygotic embryogenesis as it involves common actors of hormonal, transcriptional, developmental and epigenetic controls. Here, we provide an overview of the main signaling pathways involved in the induction and regulation of somatic embryogenesis with a focus on the master regulators of seed development, LEAFY COTYLEDON 1 and 2, ABSCISIC ACID INSENSITIVE 3 and FUSCA 3 transcription factors whose precise role during both zygotic and somatic embryogenesis remains to be fully elucidated.

The Seed Development Factors TT2 and MYB5 Regulate Heat Stress Response in Arabidopsis.

HEAT SHOCK FACTOR A2 (HSFA2) is a regulator of multiple environmental stress responses required for stress acclimation. We analyzed HSFA2 co-regulated genes and identified 43 genes strongly co-regulated with HSFA2 during multiple stresses. Motif enrichment analysis revealed an over-representation of the site II element (SIIE) in the promoters of these genes. In a yeast 1-hybrid screen with the SIIE, we identified the closely related R2R3-MYB transcription factors TT2 and MYB5. We found overexpression of MYB5 or TT2 rendered plants heat stress tolerant. In contrast, tt2, myb5, and tt2/myb5 loss of function mutants showed heat stress hypersensitivity. Transient expression assays confirmed that MYB5 and TT2 can regulate the HSFA2 promoter together with the other members of the MBW complex, TT8 and TRANSPARENT TESTA GLABRA 1 (TTG1) and that the SIIE was involved in this regulation. Transcriptomic analysis revealed that TT2/MYB5 target promoters were enriched in SIIE. Overall, we report a new function of TT2 and MYB5 in stress resistance and a role in SIIE-mediated HSFA2 regulation.

Variation in Expression of the HECT E3 Ligase UPL3 Modulates LEC2 Levels, Seed Size, and Crop Yields in Brassica napus.

Identifying genetic variation that increases crop yields is a primary objective in plant breeding. We used association analyses of oilseed rape/canola (Brassica napus) accessions to identify genetic variation that influences seed size, lipid content, and final crop yield. Variation in the promoter region of the HECT E3 ligase gene BnaUPL3 C03 made a major contribution to variation in seed weight per pod, with accessions exhibiting high seed weight per pod having lower levels of BnaUPL3 C03 expression. We defined a mechanism in which UPL3 mediated the proteasomal degradation of LEC2, a master transcriptional regulator of seed maturation. Accessions with reduced UPL3 expression had increased LEC2 protein levels, larger seeds, and prolonged expression of lipid biosynthetic genes during seed maturation. Natural variation in BnaUPL3 C03 expression appears not to have been exploited in current B napus breeding lines and could therefore be used as a new approach to maximize future yields in this important oil crop.

Regulation of FUSCA3 Expression During Seed Development in Arabidopsis.

FUSCA3 (FUS3) is a master regulator of seed development important in establishing and maintaining embryonic identity whose expression is tightly regulated at genetic and epigenetic levels. Despite this prominent role, the control of FUS3 expression remains poorly understood. Promoter and functional complementation analyses provided insight into the regulation of FUS3. W-boxes present in the promoter proximal to the start of transcription are recognized by WRKY type-1 factors which are necessary for the activation of FUS3 expression. The RY motif, the binding site of B3 factors, is important for the activation of FUS3 in the embryo proper but not in the suspensor. The loss of a negative regulatory sequence (NRS) leads to preferential expression of FUS3 in the vasculature of vegetative tissues. Since the NRS includes the RY motif, mechanisms of activation and repression target adjacent or overlapping regions. These findings discriminate the regulation of FUS3 from that of LEAFY COTYLEDON2 by the control exerted by WRKY factors and by the presence of the RY motif, yet also confirm conservation of certain regulatory elements, thereby implicating potential regulation by BASIC PENTACYSTEINE (BPC) factors and POLYCOMB REPRESSIVE COMPLEX2 (PRC2).

Combining laser-assisted microdissection (LAM) and RNA-seq allows to perform a comprehensive transcriptomic analysis of epidermal cells of Arabidopsis embryo.

BACKGROUND: Genome-wide characterization of tissue- or cell-specific gene expression is a recurrent bottleneck in biology. We have developed a sensitive approach based on ultra-low RNA sequencing coupled to laser assisted microdissection for analyzing different tissues of the small Arabidopsis embryo. METHODS AND RESULTS: We first characterized the number of genes detected according to the quantity of tissue yield and total RNA extracted. Our results revealed that as low as 0.02 mm2 of tissue and 50 pg of total RNA can be used without compromising the number of genes detected. The optimised protocol was used to compare the epidermal versus mesophyll cell transcriptomes of cotyledons at the torpedo-shaped stage of embryo development. The approach was validated by the recovery of well-known epidermal genes such AtML1 or AtPDF2 and genes involved in flavonoid and cuticular waxes pathways. Moreover, the interest and sensitivity of this approach were highlighted by the characterization of several transcription factors preferentially expressed in epidermal cells. CONCLUSION: This technical advance unlocks some current limitations of transcriptomic analyses and allows to investigate further and efficiently new biological questions for which only a very small amounts of cells need to be isolated. For instance, it paves the way to increasing the spatial accuracy of regulatory networks in developing small embryo of Arabidopsis or other plant tissues.

The Physcomitrella patens System for Transient Gene Expression Assays.

Transient expression assays are valuable techniques to study in vivo the transcriptional regulation of gene expression. These methods allow to assess the transcriptional properties of a given transcription factor (TF) or a complex of regulatory proteins against specific DNA motifs, called cis-regulatory elements. Here, we describe a fast, efficient, and reliable method based on the use of Physcomitrella patens protoplasts that allows the study of gene expression in a qualitative and quantitative manner by combining the advantage of GFP (green fluorescent protein) as a marker of promoter activity with flow cytometry for accurate measurement of fluorescence in individual cells.

Deciphering and modifying LAFL transcriptional regulatory network in seed for improving yield and quality of storage compounds.

Increasing yield and quality of seed storage compounds in a sustainable way is a key challenge for our societies. Genome-wide analyses conducted in both monocot and dicot angiosperms emphasized drastic transcriptional switches that occur during seed development. In Arabidopsis thaliana, a reference species, genetic and molecular analyses have demonstrated the key role of LAFL (LEC1, ABI3, FUS3, and LEC2) transcription factors (TFs), in controlling gene expression programs essential to accomplish seed maturation and the accumulation of storage compounds. Here, we summarize recent progress obtained in the characterization of these LAFL proteins, their regulation, partners and target genes. Moreover, we illustrate how these evolutionary conserved TFs can be used to engineer new crops with altered seed compositions and point out the current limitations. Last, we discuss about the interest of investigating further the environmental and epigenetic regulation of this network for the coming years.

Deciphering the Molecular Mechanisms Underpinning the Transcriptional Control of Gene Expression by Master Transcriptional Regulators in Arabidopsis Seed.

In Arabidopsis (Arabidopsis thaliana), transcriptional control of seed maturation involves three related regulators with a B3 domain, namely LEAFY COTYLEDON2 (LEC2), ABSCISIC ACID INSENSITIVE3 (ABI3), and FUSCA3 (ABI3/FUS3/LEC2 [AFLs]). Although genetic analyses have demonstrated partially overlapping functions of these regulators, the underlying molecular mechanisms remained elusive. The results presented here confirmed that the three proteins bind RY DNA elements (with a 5'-CATG-3' core sequence) but with different specificities for flanking nucleotides. In planta as in the moss Physcomitrella patens protoplasts, the presence of RY-like (RYL) elements is necessary but not sufficient for the regulation of the OLEOSIN1 (OLE1) promoter by the B3 AFLs. G box-like domains, located in the vicinity of the RYL elements, also are required for proper activation of the promoter, suggesting that several proteins are involved. Consistent with this idea, LEC2 and ABI3 showed synergistic effects on the activation of the OLE1 promoter. What is more, LEC1 (a homolog of the NF-YB subunit of the CCAAT-binding complex) further enhanced the activation of this target promoter in the presence of LEC2 and ABI3. Finally, recombinant LEC1 and LEC2 proteins produced in Arabidopsis protoplasts could form a ternary complex with NF-YC2 in vitro, providing a molecular explanation for their functional interactions. Taken together, these results allow us to propose a molecular model for the transcriptional regulation of seed genes by the L-AFL proteins, based on the formation of regulatory multiprotein complexes between NF-YBs, which carry a specific aspartate-55 residue, and B3 transcription factors.

Analysis of the DNA-Binding Activities of the Arabidopsis R2R3-MYB Transcription Factor Family by One-Hybrid Experiments in Yeast.

The control of growth and development of all living organisms is a complex and dynamic process that requires the harmonious expression of numerous genes. Gene expression is mainly controlled by the activity of sequence-specific DNA binding proteins called transcription factors (TFs). Amongst the various classes of eukaryotic TFs, the MYB superfamily is one of the largest and most diverse, and it has considerably expanded in the plant kingdom. R2R3-MYBs have been extensively studied over the last 15 years. However, DNA-binding specificity has been characterized for only a small subset of these proteins. Therefore, one of the remaining challenges is the exhaustive characterization of the DNA-binding specificity of all R2R3-MYB proteins. In this study, we have developed a library of Arabidopsis thaliana R2R3-MYB open reading frames, whose DNA-binding activities were assayed in vivo (yeast one-hybrid experiments) with a pool of selected cis-regulatory elements. Altogether 1904 interactions were assayed leading to the discovery of specific patterns of interactions between the various R2R3-MYB subgroups and their DNA target sequences and to the identification of key features that govern these interactions. The present work provides a comprehensive in vivo analysis of R2R3-MYB binding activities that should help in predicting new DNA motifs and identifying new putative target genes for each member of this very large family of TFs. In a broader perspective, the generated data will help to better understand how TF interact with their target DNA sequences.

ZmZHOUPI, an endosperm-specific basic helix-loop-helix transcription factor involved in maize seed development.

In angiosperm seeds the embryo is embedded within the endosperm, which is in turn enveloped by the seed coat, making inter-compartmental communication essential for coordinated seed growth. In this context the basic helix-loop-helix domain transcription factor AtZHOUPI (AtZOU) fulfils a key role in both the lysis of the transient endosperm and in embryo cuticle formation in Arabidopsis thaliana. In maize (Zea mays), a cereal with a persistent endosperm, a single gene, ZmZOU, falls into the same phylogenetic clade as AtZOU. Its expression is limited to the endosperm where it peaks during the filling stage. In ZmZOU-RNA interference knock-down lines embryo size is slightly reduced and the embryonic suspensor and the adjacent embryo surrounding region show retarded breakdown. Ectopic expression of ZmZOU reduces stomatal number, possibly due to inappropriate protein interactions. ZmZOU forms functional heterodimers with AtICE/AtSCREAM and the closely related maize proteins ZmICEb and ZmICEc, but its interaction is more efficient with the ZmICEa protein, which shows sequence divergence and only has close homologues in other monocotyledonous species. Consistent with the observation that these complexes can trans-activate target gene promoters from Arabidopsis, ZmZOU partially complements the Atzou-4 mutant. However, structural, trans-activation and gene expression data support the hypothesis that ZmZOU and ZmICEa may have coevolved to form a functional complex unique to monocot seeds. This divergence may explain the reduced functionality of ZmZOU in Arabidopsis, and reflect functional specificities which are unique to the monocotyledon lineage.

Role of B3 domain transcription factors of the AFL family in maize kernel filling.

In the dicot Arabidopsis thaliana, the B3 transcription factors, ABA-INSENSITIVE 3 (ABI3), FUSCA 3 (FUS3) and LEAFY COTYLEDON 2 (LEC2) are key regulators of seed maturation. This raises the question of the role of ABI3/FUS3/LEC2 (AFL) proteins in cereals, where not only the embryo but also the persistent endosperm accumulates reserve substances. Among the five ZmAFL genes identified in the maize genome, ZmAFL2 and ZmAFL3/ZmVp1 closely resemble FUS3 and ABI3, respectively, in terms of their sequences, domain structure and gene activity profiles. Of the three genes that fall into the LEC2 phylogenetic sub-clade, ZmAFL5 and ZmAFL6 have constitutive gene activity, whereas ZmAFL4, like LEC2, has preferential gene activity in pollen and seed, although its seed gene activity is restricted to the endosperm during reserve accumulation. Knock down of ZmAFL4 gene activity perturbs carbon metabolism and reduces starch content in the developing endosperm at 20 DAP. ZmAFL4 and ZmAFL3/ZmVp1 trans-activate a maize oleosin promoter in a heterologous moss system. In conclusion our results suggest, based on gene activity profiles, that the functions of FUS3 and ABI3 could be conserved between dicot and monocot species. In contrast, LEC2 function may have partially diverged in cereals where our findings provide first evidence of the specialization of ZmAFL4 for roles in the endosperm.

The structural organization of seed oil bodies could explain the contrasted oil extractability observed in two rapeseed genotypes.

MAIN CONCLUSION: The protein, phospholipid and sterol composition of the oil body surface from the seeds of two rapeseed genotypes was compared in order to explain their contrasted oil extractability. In the mature seeds of oleaginous plants, storage lipids accumulate in specialized structures called oil bodies (OBs). These organelles consist of a core of neutral lipids surrounded by a phospholipid monolayer in which structural proteins are embedded. The physical stability of OBs is a consequence of the interactions between proteins and phospholipids. A detailed study of OB characteristics in mature seeds as well as throughout seed development was carried out on two contrasting rapeseed genotypes Amber and Warzanwski. These two accessions were chosen because they differ dramatically in (1) crushing ability, (2) oil extraction yield and, (3) the stability of purified OBs. Warzanwski has higher crushing ability, better oil extraction yield and less stable purified OBs than Amber. OB morphology was investigated in situ using fluorescence microscopy, transmission electron microscopy and pulsed field gradient NMR. During seed development, OB diameter first increased and then decreased 30 days after pollination in both Amber and Warzanwski embryos. In mature seeds, Amber OBs were significantly smaller. The protein, phospholipid and sterol composition of the hemi-membrane was compared between the two accessions. Amber OBs were enriched with H-oleosins and steroleosins, suggesting increased coverage of the OB surface consistent with their higher stability. The nature and composition of phospholipids and sterols in Amber OBs suggest that the hemi-membrane would have a more rigid structure than that of Warzanwski OBs.

MYB118 represses endosperm maturation in seeds of Arabidopsis.

In the exalbuminous species Arabidopsis thaliana, seed maturation is accompanied by the deposition of oil and storage proteins and the reduction of the endosperm to one cell layer. Here, we consider reserve partitioning between embryo and endosperm compartments. The pattern of deposition, final amount, and composition of these reserves differ between the two compartments, with the embryo representing the principal storage tissue in mature seeds. Complex regulatory mechanisms are known to prevent activation of maturation-related programs during embryo morphogenesis and, later, during vegetative growth. Here, we describe a regulator that represses the expression of maturation-related genes during maturation within the endosperm. MYB118 is transcriptionally induced in the maturing endosperm, and seeds of myb118 mutants exhibit an endosperm-specific derepression of maturation-related genes associated with a partial relocation of storage compounds from the embryo to the endosperm. Moreover, MYB118 activates endosperm-induced genes through the recognition of TAACGG elements. These results demonstrate that the differential partitioning of reserves between the embryo and endosperm in exalbuminous Arabidopsis seeds does not only result from developmental programs that establish the embryo as the preponderant tissue within seeds. This differential partitioning is also regulated by MYB118, which regulates the biosynthesis of reserves at the spatial level during maturation.

Function and localization of the Arabidopsis thaliana diacylglycerol acyltransferase DGAT2 expressed in yeast.

Diacylglycerol acyltransferases (DGATs) catalyze the final and only committed step of triacylglycerol synthesis. DGAT activity is rate limiting for triacylglycerol accumulation in mammals, plants and microbes. DGATs belong to three different evolutionary classes. In Arabidopsis thaliana, DGAT1, encoded by At2g19450, is the major DGAT enzyme involved in triacylglycerol accumulation in seeds. Until recently, the function of DGAT2 (At3g51520) has remained elusive. Previous attempts to characterize its enzymatic function by heterologous expression in yeast were unsuccessful. In the present report we demonstrate that expression of a codon-optimized version of the DGAT2 gene is able to restore neutral lipid accumulation in the Saccharomyces cerevisiae mutant strain (H1246), which is defective in triacylglycerol biosynthesis. Heterologous expression of codon-optimized DGAT2 and DGAT1 induced the biogenesis of subcellular lipid droplets containing triacylglycerols and squalene. Both DGAT proteins were found to be associated with these lipid droplets. The fatty acid composition was affected by the nature of the acyltransferase expressed. DGAT2 preferentially incorporated C16:1 fatty acids whereas DGAT1 displayed preference for C16:0, strongly suggesting that these enzymes have contrasting substrate specificities.

Endosperm breakdown in Arabidopsis requires heterodimers of the basic helix-loop-helix proteins ZHOUPI and INDUCER OF CBP EXPRESSION 1.

In Arabidopsis seeds, embryo growth is coordinated with endosperm breakdown. Mutants in the endosperm-specific gene ZHOUPI (ZOU), which encodes a unique basic helix-loop-helix (bHLH) transcription factor, have an abnormal endosperm that persists throughout seed development, significantly impeding embryo growth. Here we show that loss of function of the bHLH-encoding gene INDUCER OF CBP EXPRESSION 1 (ICE1) causes an identical endosperm persistence phenotype. We show that ZOU and ICE1 are co-expressed in the endosperm and interact in yeast via their bHLH domains. We show both genetically and in a heterologous plant system that, despite the fact that both ZOU and ICE1 can form homodimers in yeast, their role in endosperm breakdown requires their heterodimerization. Consistent with this conclusion, we confirm that ZOU and ICE1 regulate the expression of common target genes in the developing endosperm. Finally, we show that heterodimerization of ZOU and ICE1 is likely to be necessary for their binding to specific targets, rather than for their nuclear localization in the endosperm. By comparing our results with paradigms of bHLH function and evolution in animal systems we propose that the ZOU/ICE1 complex might have ancient origins, acquiring novel megagametophyte-specific functions in heterosporous land plants that were conserved in the angiosperm endosperm.

Specialization of oleosins in oil body dynamics during seed development in Arabidopsis seeds.

Oil bodies (OBs) are seed-specific lipid storage organelles that allow the accumulation of neutral lipids that sustain plantlet development after the onset of germination. OBs are covered with specific proteins embedded in a single layer of phospholipids. Using fluorescent dyes and confocal microscopy, we monitored the dynamics of OBs in living Arabidopsis (Arabidopsis thaliana) embryos at different stages of development. Analyses were carried out with different genotypes: the wild type and three mutants affected in the accumulation of various oleosins (OLE1, OLE2, and OLE4), three major OB proteins. Image acquisition was followed by a detailed statistical analysis of OB size and distribution during seed development in the four dimensions (x, y, z, and t). Our results indicate that OB size increases sharply during seed maturation, in part by OB fusion, and then decreases until the end of the maturation process. In single, double, and triple mutant backgrounds, the size and spatial distribution of OBs are modified, affecting in turn the total lipid content, which suggests that the oleosins studied have specific functions in the dynamics of lipid accumulation.

Expression variation in connected recombinant populations of Arabidopsis thaliana highlights distinct transcriptome architectures.

BACKGROUND: Expression traits can vary quantitatively between individuals and have a complex inheritance. Identification of the genetics underlying transcript variation can help in the understanding of phenotypic variation due to genetic factors regulating transcript abundance and shed light into divergence patterns. So far, only a limited number of studies have addressed this subject in Arabidopsis, with contrasting results due to dissimilar statistical power. Here, we present the transcriptome architecture in leaf tissue of two RIL sets obtained from a connected-cross design involving 3 commonly used accessions. We also present the transcriptome architecture observed in developing seeds of a third independent cross. RESULTS: The utilisation of the novel R/eqtl package (which goal is to automatize and extend functions from the R/qtl package) allowed us to map 4,290 and 6,534 eQTLs in the Cvi-0 × Col-0 and Bur-0 × Col-0 recombinant populations respectively. In agreement with previous studies, we observed a larger phenotypic variance explained by eQTLs in linkage with the controlled gene (potentially cis-acting), compared to distant loci (acting necessarily indirectly or in trans). Distant eQTLs hotspots were essentially not conserved between crosses, but instead, cross-specific. Accounting for confounding factors using a probabilistic approach (VBQTL) increased the mapping resolution and the number of significant associations. Moreover, using local eQTLs obtained from this approach, we detected evidence for a directional allelic effect in genes with related function, where significantly more eQTLs than expected by chance were up-regulated from one of the accessions. Primary experimental data, analysis parameters, eQTL results and visualisation of LOD score curves presented here are stored and accessible through the QTLstore service database http://qtlstore.versailles.inra.fr/. CONCLUSIONS: Our results demonstrate the extensive diversity and moderately conserved eQTL landscape between crosses and validate the utilisation of expression traits to explore for candidates behind phenotypic variation among accessions. Furthermore, this stresses the need for a wider spectrum of diversity to fully understand expression trait variation within a species.