Grain carbon isotope composition is a marker for allocation and harvest index in wheat.

The natural 13 C abundance (δ13 C) in plant leaves has been used for decades with great success in agronomy to monitor water-use efficiency and select modern cultivars adapted to dry conditions. However, in wheat, it is also important to find genotypes with high carbon allocation to spikes and grains, and thus with a high harvest index (HI) and/or low carbon losses via respiration. Finding isotope-based markers of carbon partitioning to grains would be extremely useful since isotope analyses are inexpensive and can be performed routinely at high throughput. Here, we took the advantage of a set of field trials made of more than 600 plots with several wheat cultivars and measured agronomic parameters as well as δ13 C values in leaves and grains. We find a linear relationship between the apparent isotope discrimination between leaves and grain (denoted as Δδcorr ), and the respiration use efficiency-to-HI ratio. It means that overall, efficient carbon allocation to grains is associated with a small isotopic difference between leaves and grains. This effect is explained by postphotosynthetic isotope fractionations, and we show that this can be modelled by equations describing the carbon isotope composition in grains along the wheat growth cycle. Our results show that 13 C natural abundance in grains could be useful to find genotypes with better carbon allocation properties and assist current wheat breeding technologies.

Breeding for Economically and Environmentally Sustainable Wheat Varieties: An Integrated Approach from Genomics to Selection.

There is currently a strong societal demand for sustainability, quality, and safety in bread wheat production. To address these challenges, new and innovative knowledge, resources, tools, and methods to facilitate breeding are needed. This starts with the development of high throughput genomic tools including single nucleotide polymorphism (SNP) arrays, high density molecular marker maps, and full genome sequences. Such powerful tools are essential to perform genome-wide association studies (GWAS), to implement genomic and phenomic selection, and to characterize the worldwide diversity. This is also useful to breeders to broaden the genetic basis of elite varieties through the introduction of novel sources of genetic diversity. Improvement in varieties particularly relies on the detection of genomic regions involved in agronomical traits including tolerance to biotic (diseases and pests) and abiotic (drought, nutrient deficiency, high temperature) stresses. When enough resolution is achieved, this can result in the identification of candidate genes that could further be characterized to identify relevant alleles. Breeding must also now be approached through in silico modeling to simulate plant development, investigate genotype × environment interactions, and introduce marker-trait linkage information in the models to better implement genomic selection. Breeders must be aware of new developments and the information must be made available to the world wheat community to develop new high-yielding varieties that can meet the challenge of higher wheat production in a sustainable and fluctuating agricultural context. In this review, we compiled all knowledge and tools produced during the BREEDWHEAT project to show how they may contribute to face this challenge in the coming years.

Imprints of natural selection along environmental gradients in phenology-related genes of Quercus petraea.

We explored single nucleotide polymorphism (SNP) variation in candidate genes for bud burst from Quercus petraea populations sampled along gradients of latitude and altitude in Western Europe. SNP diversity was monitored for 106 candidate genes, in 758 individuals from 32 natural populations. We investigated whether SNP variation reflected the clinal pattern of bud burst observed in common garden experiments. We used different methods to detect imprints of natural selection (FST outlier, clinal variation at allelic frequencies, association tests) and compared the results obtained for the two gradients. FST outlier SNPs were found in 15 genes, 5 of which were common to both gradients. The type of selection differed between the two gradients (directional or balancing) for 3 of these 5. Clinal variations were observed for six SNPs, and one cline was conserved across both gradients. Association tests between the phenotypic or breeding values of trees and SNP genotypes identified 14 significant associations, involving 12 genes. The results of outlier detection on the basis of population differentiation or clinal variation were not very consistent with the results of association tests. The discrepancies between these approaches may reflect the different hierarchical levels of selection considered (inter- and intrapopulation selection). Finally, we obtained evidence for convergent selection (similar for gradients) and clinal variation for a few genes, suggesting that comparisons between parallel gradients could be used to screen for major candidate genes responding to natural selection in trees.

Transcriptional profiling of bud dormancy induction and release in oak by next-generation sequencing.

BACKGROUND: In temperate regions, the time lag between vegetative bud burst and bud set determines the duration of the growing season of trees (i.e. the duration of wood biomass production). Dormancy, the period during which the plant is not growing, allows trees to avoid cold injury resulting from exposure to low temperatures. An understanding of the molecular machinery controlling the shift between these two phenological states is of key importance in the context of climatic change. The objective of this study was to identify genes upregulated during endo- and ecodormancy, the two main stages of bud dormancy. Sessile oak is a widely distributed European white oak species. A forcing test on young trees was first carried out to identify the period most likely to correspond to these two stages. Total RNA was then extracted from apical buds displaying endo- and ecodormancy. This RNA was used for the generation of cDNA libraries, and in-depth transcriptome characterization was performed with 454 FLX pyrosequencing technology. RESULTS: Pyrosequencing produced a total of 495,915 reads. The data were cleaned, duplicated reads removed, and sequences were mapped onto the oak UniGene data. Digital gene expression analysis was performed, with both R statistics and the R-Bioconductor packages (edgeR and DESeq), on 6,471 contigs with read numbers ≥ 5 within any contigs. The number of sequences displaying significant differences in expression level (read abundance) between endo- and ecodormancy conditions ranged from 75 to 161, depending on the algorithm used. 13 genes displaying significant differences between conditions were selected for further analysis, and 11 of these genes, including those for glutathione-S-transferase (GST) and dehydrin xero2 (XERO2) were validated by quantitative PCR. CONCLUSIONS: The identification and functional annotation of differentially expressed genes involved in the "response to abscisic acid", "response to cold stress" and "response to oxidative stress" categories constitutes a major step towards characterization of the molecular network underlying vegetative bud dormancy, an important life history trait of long-lived organisms.

Bioinformatic analysis of ESTs collected by Sanger and pyrosequencing methods for a keystone forest tree species: oak.

BACKGROUND: The Fagaceae family comprises about 1,000 woody species worldwide. About half belong to the Quercus family. These oaks are often a source of raw material for biomass wood and fiber. Pedunculate and sessile oaks, are among the most important deciduous forest tree species in Europe. Despite their ecological and economical importance, very few genomic resources have yet been generated for these species. Here, we describe the development of an EST catalogue that will support ecosystem genomics studies, where geneticists, ecophysiologists, molecular biologists and ecologists join their efforts for understanding, monitoring and predicting functional genetic diversity. RESULTS: We generated 145,827 sequence reads from 20 cDNA libraries using the Sanger method. Unexploitable chromatograms and quality checking lead us to eliminate 19,941 sequences. Finally a total of 125,925 ESTs were retained from 111,361 cDNA clones. Pyrosequencing was also conducted for 14 libraries, generating 1,948,579 reads, from which 370,566 sequences (19.0%) were eliminated, resulting in 1,578,192 sequences. Following clustering and assembly using TGICL pipeline, 1,704,117 EST sequences collapsed into 69,154 tentative contigs and 153,517 singletons, providing 222,671 non-redundant sequences (including alternative transcripts). We also assembled the sequences using MIRA and PartiGene software and compared the three unigene sets. Gene ontology annotation was then assigned to 29,303 unigene elements. Blast search against the SWISS-PROT database revealed putative homologs for 32,810 (14.7%) unigene elements, but more extensive search with Pfam, Refseq_protein, Refseq_RNA and eight gene indices revealed homology for 67.4% of them. The EST catalogue was examined for putative homologs of candidate genes involved in bud phenology, cuticle formation, phenylpropanoids biosynthesis and cell wall formation. Our results suggest a good coverage of genes involved in these traits. Comparative orthologous sequences (COS) with other plant gene models were identified and allow to unravel the oak paleo-history. Simple sequence repeats (SSRs) and single nucleotide polymorphisms (SNPs) were searched, resulting in 52,834 SSRs and 36,411 SNPs. All of these are available through the Oak Contig Browser http://genotoul-contigbrowser.toulouse.inra.fr:9092/Quercus_robur/index.html. CONCLUSIONS: This genomic resource provides a unique tool to discover genes of interest, study the oak transcriptome, and develop new markers to investigate functional diversity in natural populations.

Transcriptional changes in two types of pre-mycorrhizal roots and in ectomycorrhizas of oak microcuttings inoculated with Piloderma croceum.

The formation of the ectomycorrhiza implies an alteration in gene expression of both the plant and fungal partners, a process which starts before the formation of any symbiotic interface. However, little is known on the regulation pattern occurring in different parts of the root system. Our experimental system consisting of a micropropagated oak with a hierarchical root system was shown to exhibit symbiosis functional traits prior to any mycorrhizal tissue differentiation after the inoculation with the basidiomycete Piloderma croceum. Using a cDNA array, the plant gene regulation was analyzed in the pre-mycorrhizal phase. Seventy-five transcripts showed differential expression in pre-mycorrhizal lateral and principal roots, and both root types exhibited different sets of responsive genes. For transcripts selected according to a statistical analysis, the alteration in gene expression was confirmed by RT-PCR and quantitative real-time PCR. Genes regulated in pre-mycorrhizal lateral roots displayed an almost identical expression in mycorrhizas. In contrast, genes regulated in pre-mycorrhizal principal roots were often regulated differently in ectomycorrhizas. Down-regulation affected most of the regulated genes involved in metabolism, whereas most of the regulated genes related to cell rescue functions, water regulation and defence response were up-regulated. Regulation of such genes could explain the increase of global resistance observed in mycorrhizal plants.

Transcriptome analysis of bud burst in sessile oak (Quercus petraea).

Expression patterns of hundreds of transcripts in apical buds were monitored during bud flushing in sessile oak (Quercus petraea), in order to identify genes differentially expressed between the quiescent and active stage of bud development. Different transcriptomic techniques combining the construction of suppression subtractive hybridization (SSH) libraries and the monitoring of gene expression using macroarray and real-time reverse transcriptase polymerase chain reaction (RT-PCR) were performed to dissect bud burst, with a special emphasis on the onset of the process. We generated 801 expressed sequence tags (ESTs) derived from six developmental stages of bud burst. Macroarray experiment revealed a total of 233 unique transcripts exhibiting differential expression during the process, and a putative function was assigned to 65% of them. Cell rescue/defense-, metabolism-, protein synthesis-, cell cycle- and transcription-related transcripts were among the most regulated genes. Macroarray and real-time RT-PCR showed that several genes exhibited contrasted expressions between quiescent and swelling buds, such as a putative homologue of the transcription factor DAG2 (Dof Affecting Germination 2), previously reported to be involved in the control of seed germination in Arabidopsis thaliana. These differentially expressed genes constitute relevant candidates for signaling pathway of bud burst in trees.

Comparison of quantitative trait loci for adaptive traits between oak and chestnut based on an expressed sequence tag consensus map.

A comparative genetic and QTL mapping was performed between Quercus robur L. and Castanea sativa Mill., two major forest tree species belonging to the Fagaceae family. Oak EST-derived markers (STSs) were used to align the 12 linkage groups of the two species. Fifty-one and 45 STSs were mapped in oak and chestnut, respectively. These STSs, added to SSR markers previously mapped in both species, provided a total number of 55 orthologous molecular markers for comparative mapping within the Fagaceae family. Homeologous genomic regions identified between oak and chestnut allowed us to compare QTL positions for three important adaptive traits. Colocation of the QTL controlling the timing of bud burst was significant between the two species. However, conservation of QTL for height growth was not supported by statistical tests. No QTL for carbon isotope discrimination was conserved between the two species. Putative candidate genes for bud burst can be identified on the basis of colocations between EST-derived markers and QTL.