Tracking population structure and phenology through time using ancient genomes from waterlogged white oak wood.

Whole genome characterizations of crop plants based on ancient DNA have provided unique keys for a better understanding of the evolutionary origins of modern cultivars, the pace and mode of selection underlying their adaptation to new environments and the production of phenotypes of interest. Although forests are among the most biologically rich ecosystems on earth and represent a fundamental resource for human societies, no ancient genome sequences have been generated for trees. This contrasts with the generation of multiple ancient reference genomes for important crops. Here, we sequenced the first ancient tree genomes using two white oak wood remains from Germany dating to the Last Little Ice Age (15th century CE, 7.3× and 4.0×) and one from France dating to the Bronze Age (1700 BCE, 3.4×). We assessed the underlying species and identified one medieval remains as a hybrid between two common oak species (Quercus robur and Q. petraea) and the other two remains as Q. robur. We found that diversity at the global genome level had not changed over time. However, exploratory analyses suggested that a reduction of diversity took place at different time periods. Finally, we determined the timing of leaf unfolding for ancient trees for the first time. The study extends the application of ancient wood beyond the classical proxies of dendroclimatology, dendrochronology, dendroarchaeology and dendroecology, thereby enhancing resolution of inferences on the responses of forest ecosystems to past environmental changes, epidemics and silvicultural practices.

Genome-wide evolutionary response of European oaks during the Anthropocene.

The pace of tree microevolution during Anthropocene warming is largely unknown. We used a retrospective approach to monitor genomic changes in oak trees since the Little Ice Age (LIA). Allelic frequency changes were assessed from whole-genome pooled sequences for four age-structured cohorts of sessile oak (Quercus petraea) dating back to 1680, in each of three different oak forests in France. The genetic covariances of allelic frequency changes increased between successive time periods, highlighting genome-wide effects of linked selection. We found imprints of parallel linked selection in the three forests during the late LIA, and a shift of selection during more recent time periods of the Anthropocene. The changes in allelic covariances within and between forests mirrored the documented changes in the occurrence of extreme events (droughts and frosts) over the last 300 years. The genomic regions with the highest covariances were enriched in genes involved in plant responses to pathogens and abiotic stresses (temperature and drought). These responses are consistent with the reported sequence of frost (or drought) and disease damage ultimately leading to the oak dieback after extreme events. They provide support for adaptive evolution of long-lived species during recent climatic changes. Although we acknowledge that other sources (e.g., gene flow, generation overlap) may have contributed to temporal covariances of allelic frequency changes, the consistent and correlated response across the three forests lends support to the existence of a systematic driving force such as natural selection.

Gene expression predictions and networks in natural populations supports the omnigenic theory.

BACKGROUND: Recent literature on the differential role of genes within networks distinguishes core from peripheral genes. If previous works have shown contrasting features between them, whether such categorization matters for phenotype prediction remains to be studied. RESULTS: We measured 17 phenotypic traits for 241 cloned genotypes from a Populus nigra collection, covering growth, phenology, chemical and physical properties. We also sequenced RNA for each genotype and built co-expression networks to define core and peripheral genes. We found that cores were more differentiated between populations than peripherals while being less variable, suggesting that they have been constrained through potentially divergent selection. We also showed that while cores were overrepresented in a subset of genes statistically selected for their capacity to predict the phenotypes (by Boruta algorithm), they did not systematically predict better than peripherals or even random genes. CONCLUSION: Our work is the first attempt to assess the importance of co-expression network connectivity in phenotype prediction. While highly connected core genes appear to be important, they do not bear enough information to systematically predict better quantitative traits than other gene sets.

The pulse of the tree is under genetic control: eucalyptus as a case study.

The pulse of the tree (diurnal cycle of stem radius fluctuations) has been widely studied as a way of analyzing tree responses to the environment, including the phenotypic plasticity of tree-water relationships in particular. However, the genetic basis of this daily phenotype and its interplay with the environment remain largely unexplored. We characterized the genetic and environmental determinants of this response, by monitoring daily stem radius fluctuation (dSRF) on 210 trees from a Eucalyptus urophylla × E. grandis full-sib family over 2 years. The dSRF signal was broken down into hydraulic capacitance, assessed as the daily amplitude of shrinkage (DA), and net growth, estimated as the change in maximum radius between two consecutive days (ΔR). The environmental determinants of these two traits were clearly different: DA was positively correlated with atmospheric variables relating to water demand, while ΔR was associated with soil water content. The heritability for these two traits ranged from low to moderate over time, revealing a time-dependent or environment-dependent complex genetic determinism. We identified 686 and 384 daily quantitative trait loci (QTL) representing 32 and 31 QTL regions for DA and ΔR, respectively. The identification of gene networks underlying the 27 major genomics regions for both traits generated additional hypotheses concerning the biological mechanisms involved in response to water demand and supply. This study highlights that environmentally induced changes in daily stem radius fluctuation are genetically controlled in trees and suggests that these daily responses integrated over time shape the genetic architecture of mature traits.

Cytological, Biochemical and Molecular Events of the Embryogenic State in Douglas-fir (Pseudotsuga menziesii [Mirb.]).

Somatic embryogenesis techniques have been developed for most coniferous species, but only using very juvenile material. To extend the techniques' scope, better integrated understanding of the key biological, physiological and molecular characteristics of embryogenic state is required. Therefore, embryonal masses (EMs) and non-embryogenic calli (NECs) have been compared during proliferation at multiple levels. EMs and NECs originating from a single somatic embryo (isogenic lines) of each of three unrelated genotypes were used in the analyses, which included comparison of the lines' anatomy by transmission light microscopy, transcriptomes by RNAseq Illumina sequencing, proteomes by free-gel analysis, contents of endogenous phytohormones (indole-3-acetic acid, cytokinins and ABA) by LC-MS analysis, and soluble sugar contents by HPLC. EMs were characterized by upregulation (relative to levels in NECs) of transcripts, proteins, transcription factors and active cytokinins associated with cell differentiation accompanied by histological, carbohydrate content and genetic markers of cell division. In contrast, NECs were characterized by upregulation (relative to levels in EMs) of transcripts, proteins and products associated with responses to stimuli (ABA, degradation forms of cytokinins, phenols), oxidative stress (reactive oxygen species) and carbohydrate storage (starch). Sub-Network Enrichment Analyses that highlighted functions and interactions of transcripts and proteins that significantly differed between EMs and NECs corroborated these findings. The study shows the utility of a novel approach involving integrated multi-scale transcriptomic, proteomic, biochemical, histological and anatomical analyses to obtain insights into molecular events associated with embryogenesis and more specifically to the embryogenic state of cell in Douglas-fir.

Repetitive somatic embryogenesis induced cytological and proteomic changes in embryogenic lines of Pseudotsuga menziesii [Mirb.].

BACKGROUND: To explore poorly understood differences between primary and subsequent somatic embryogenic lines of plants, we induced secondary (2ry) and tertiary (3ry) lines from cotyledonary somatic embryos (SEs) of two Douglas-fir genotypes: SD4 and TD17. The 2ry lines exhibited significantly higher embryogenic potential (SE yields) than the 1ry lines initiated from zygotic embryos (SD4, 2155 vs 477; TD17, 240 vs 29 g- 1 f.w.). Moreover, we observed similar differences in yield between 2ry and 3ry lines of SD4 (2400 vs 3921 g- 1 f.w.). To elucidate reasons for differences in embryogenic potential induced by repetitive somatic embryogenesis we then compared 2ry vs 1ry and 2ry vs 3ry lines at histo-cytological (using LC-MS/MS) and proteomic levels. RESULTS: Repetitive somatic embryogenesis dramatically improved the proliferating lines' cellular organization (genotype SD4's most strongly). Frequencies of singulated, bipolar SEs and compact polyembryogenic centers with elongated suspensors and apparently cleavable embryonal heads increased in 2ry and (even more) 3ry lines. Among 2300-2500 identified proteins, 162 and 228 were classified significantly differentially expressed between 2ry vs 1ry and 3ry vs 2ry lines, respectively, with special emphasis on "Proteolysis" and "Catabolic process" Gene Ontology categories. Strikingly, most of the significant proteins (> 70%) were down-regulated in 2ry relative to 1ry lines, but up-regulated in 3ry relative to 2ry lines, revealing a down-up pattern of expression. GO category enrichment analyses highlighted the opposite adjustments of global protein patterns, particularly for processes involved in chitin catabolism, lignin and L-phenylalanine metabolism, phenylpropanoid biosynthesis, oxidation-reduction, and response to karrikin. Sub-Network Enrichment Analyses highlighted interactions between significant proteins and both plant growth regulators and secondary metabolites after first (especially jasmonic acid, flavonoids) and second (especially salicylic acid, abscisic acid, lignin) embryogenesis cycles. Protein networks established after each induction affected the same "Plant development" and "Defense response" biological processes, but most strongly after the third cycle, which could explain the top embryogenic performance of 3ry lines. CONCLUSIONS: This first report of cellular and molecular changes after repetitive somatic embryogenesis in conifers shows that each cycle enhanced the structure and singularization of EMs through modulation of growth regulator pathways, thereby improving the lines' embryogenic status.

Oak genome reveals facets of long lifespan.

Oaks are an important part of our natural and cultural heritage. Not only are they ubiquitous in our most common landscapes1 but they have also supplied human societies with invaluable services, including food and shelter, since prehistoric times2. With 450 species spread throughout Asia, Europe and America3, oaks constitute a critical global renewable resource. The longevity of oaks (several hundred years) probably underlies their emblematic cultural and historical importance. Such long-lived sessile organisms must persist in the face of a wide range of abiotic and biotic threats over their lifespans. We investigated the genomic features associated with such a long lifespan by sequencing, assembling and annotating the oak genome. We then used the growing number of whole-genome sequences for plants (including tree and herbaceous species) to investigate the parallel evolution of genomic characteristics potentially underpinning tree longevity. A further consequence of the long lifespan of trees is their accumulation of somatic mutations during mitotic divisions of stem cells present in the shoot apical meristems. Empirical4 and modelling5 approaches have shown that intra-organismal genetic heterogeneity can be selected for6 and provides direct fitness benefits in the arms race with short-lived pests and pathogens through a patchwork of intra-organismal phenotypes7. However, there is no clear proof that large-statured trees consist of a genetic mosaic of clonally distinct cell lineages within and between branches. Through this case study of oak, we demonstrate the accumulation and transmission of somatic mutations and the expansion of disease-resistance gene families in trees.

Non-cellulosic polysaccharide distribution during G-layer formation in poplar tension wood fibers: abundance of rhamnogalacturonan I and arabinogalactan proteins but no evidence of xyloglucan.

MAIN CONCLUSION: RG-I and AGP, but not XG, are associated to the building of the peculiar mechanical properties of tension wood. Hardwood trees produce tension wood (TW) with specific mechanical properties to cope with environmental cues. Poplar TW fibers have an additional cell wall layer, the G-layer responsible for TW mechanical properties. We investigated, in two poplar hybrid species, the molecules potentially involved in the building of TW mechanical properties. First, we evaluated the distribution of the different classes of non-cellulosic polysaccharides during xylem fiber differentiation, using immunolocalization. In parallel, G-layers were isolated and their polysaccharide composition determined. These complementary approaches provided information on the occurrence of non-cellulosic polysaccharides during G-fiber differentiation. We found no evidence of the presence of xyloglucan (XG) in poplar G-layers, whereas arabinogalactan proteins (AGP) and rhamnogalacturonan type I pectins (RG-I) were abundant, with an apparent progressive loss of RG-I side chains during G-layer maturation. Similarly, the intensity of immunolabeling signals specific for glucomannans and glucuronoxylans varies during G-layer maturation. RG-I and AGP are best candidate matrix components to be responsible for TW mechanical properties.

Implication of the suberin pathway in adaptation to waterlogging and hypertrophied lenticels formation in pedunculate oak (Quercus robur L.).

Waterlogging causes stressful conditions for perennial species. The temporary overabundance of water in waterlogged soil can induce hypoxia in the rhizosphere, leading to root death, tree decline and even dieback. Two closely related members of the European white oak complex, pedunculate (Quercus robur L.) and sessile (Quercus petraea Matt. Liebl.) oaks, have different ecological characteristics, especially regarding their adaptation to soil waterlogging. The tolerance of waterlogging observed in pedunculate oak is driven principally by its ability to produce adaptive structures, hypertrophied lenticels and adventitious roots, and to switch rapidly its metabolism to the fermentative pathway. This study had two objectives: (i) to identify genes important for adaptation to waterlogging and (ii) to gain insight into the molecular mechanisms involved in hypertrophied lenticel formation in pedunculate oak. We subjected seedlings of the two species to hypoxia by maintaining the water level 2 cm above the collar. The immersed part of the stem (i.e., containing hypertrophied lenticels in pedunculate oak) was sampled after 9 days of waterlogging stress and its gene expression was investigated by RNA-seq. Genes displaying differential expression between the two species were identified with the DESeq R package and a false discovery rate of 0.001. We found that 3705 contigs were differentially regulated between the two species. Twenty-two differentially expressed genes were validated by real-time quantitative polymerase chain reaction. The suberin biosynthesis pathway was found to be upregulated in pedunculate oak, consistent with molecular mechanisms analogous to those operating in the radial oxygen loss barrier in waterlogging-tolerant species.

Genome Sequences of Populus tremula Chloroplast and Mitochondrion: Implications for Holistic Poplar Breeding.

Complete Populus genome sequences are available for the nucleus (P. trichocarpa; section Tacamahaca) and for chloroplasts (seven species), but not for mitochondria. Here, we provide the complete genome sequences of the chloroplast and the mitochondrion for the clones P. tremula W52 and P. tremula x P. alba 717-1B4 (section Populus). The organization of the chloroplast genomes of both Populus clones is described. A phylogenetic tree constructed from all available complete chloroplast DNA sequences of Populus was not congruent with the assignment of the related species to different Populus sections. In total, 3,024 variable nucleotide positions were identified among all compared Populus chloroplast DNA sequences. The 5-prime part of the LSC from trnH to atpA showed the highest frequency of variations. The variable positions included 163 positions with SNPs allowing for differentiating the two clones with P. tremula chloroplast genomes (W52, 717-1B4) from the other seven Populus individuals. These potential P. tremula-specific SNPs were displayed as a whole-plastome barcode on the P. tremula W52 chloroplast DNA sequence. Three of these SNPs and one InDel in the trnH-psbA linker were successfully validated by Sanger sequencing in an extended set of Populus individuals. The complete mitochondrial genome sequence of P. tremula is the first in the family of Salicaceae. The mitochondrial genomes of the two clones are 783,442 bp (W52) and 783,513 bp (717-1B4) in size, structurally very similar and organized as single circles. DNA sequence regions with high similarity to the W52 chloroplast sequence account for about 2% of the W52 mitochondrial genome. The mean SNP frequency was found to be nearly six fold higher in the chloroplast than in the mitochondrial genome when comparing 717-1B4 with W52. The availability of the genomic information of all three DNA-containing cell organelles will allow a holistic approach in poplar molecular breeding in the future.

Decoding the oak genome: public release of sequence data, assembly, annotation and publication strategies.

The 1.5 Gbp/2C genome of pedunculate oak (Quercus robur) has been sequenced. A strategy was established for dealing with the challenges imposed by the sequencing of such a large, complex and highly heterozygous genome by a whole-genome shotgun (WGS) approach, without the use of costly and time-consuming methods, such as fosmid or BAC clone-based hierarchical sequencing methods. The sequencing strategy combined short and long reads. Over 49 million reads provided by Roche 454 GS-FLX technology were assembled into contigs and combined with shorter Illumina sequence reads from paired-end and mate-pair libraries of different insert sizes, to build scaffolds. Errors were corrected and gaps filled with Illumina paired-end reads and contaminants detected, resulting in a total of 17,910 scaffolds (>2 kb) corresponding to 1.34 Gb. Fifty per cent of the assembly was accounted for by 1468 scaffolds (N50 of 260 kb). Initial comparison with the phylogenetically related Prunus persica gene model indicated that genes for 84.6% of the proteins present in peach (mean protein coverage of 90.5%) were present in our assembly. The second and third steps in this project are genome annotation and the assignment of scaffolds to the oak genetic linkage map. In accordance with the Bermuda and Fort Lauderdale agreements and the more recent Toronto Statement, the oak genome data have been released into public sequence repositories in advance of publication. In this presubmission paper, the oak genome consortium describes its principal lines of work and future directions for analyses of the nature, function and evolution of the oak genome.

Quantitative Proteomic and Phosphoproteomic Approaches for Deciphering the Signaling Pathway for Tension Wood Formation in Poplar.

Trees adjust their growth following forced changes in orientation to re-establish a vertical position. In angiosperms, this adjustment involves the differential regulation of vascular cambial activity between the lower (opposite wood) and upper (tension wood) sides of the leaning stem. We investigated the molecular mechanisms leading to the formation of differential wood types through a quantitative proteomic and phosphoproteomic analysis on poplar subjected to a gravitropic stimulus. We identified and quantified 675 phosphopeptides, corresponding to 468 phosphoproteins, and 3 763 nonphosphorylated peptides, corresponding to 1 155 proteins, in the differentiating xylem of straight-growing trees (control) and trees subjected to a gravitational stimulus during 8 weeks. About 1% of the peptides were specific to a wood type (straight, opposite, or tension wood). Proteins quantified in more than one type of wood were more numerous: a mixed linear model showed 389 phosphopeptides and 556 proteins to differ in abundance between tension wood and opposite wood. Twenty-one percent of the phosphoproteins identified here were described in their phosphorylated form for the first time. Our analyses revealed remarkable developmental molecular plasticity, with wood type-specific phosphorylation events, and highlighted the involvement of different proteins in the biosynthesis of cell wall components during the formation of the three types of wood.

PtaRHE1, a Populus tremula × Populus alba RING-H2 protein of the ATL family, has a regulatory role in secondary phloem fibre development.

REALLY INTERESTING NEW GENE (RING) proteins play important roles in the regulation of many processes by recognizing target proteins for ubiquitination. Previously, we have shown that the expression of PtaRHE1, encoding a Populus tremula × Populus alba RING-H2 protein with E3 ubiquitin ligase activity, is associated with tissues undergoing secondary growth. To further elucidate the role of PtaRHE1 in vascular tissues, we have undertaken a reverse genetic analysis in poplar. Within stem secondary vascular tissues, PtaRHE1 and its corresponding protein are expressed predominantly in the phloem. The downregulation of PtaRHE1 in poplar by artificial miRNA triggers alterations in phloem fibre patterning, characterized by an increased portion of secondary phloem fibres that have a reduced cell wall thickness and a change in lignin composition, with lower levels of syringyl units as compared with wild-type plants. Following an RNA-seq analysis, a biological network involving hormone stress signalling, as well as developmental processes, could be delineated. Several candidate genes possibly associated with the altered phloem fibre phenotype observed in amiRPtaRHE1 poplar were identified. Altogether, our data suggest a regulatory role for PtaRHE1 in secondary phloem fibre development.

The oak gene expression atlas: insights into Fagaceae genome evolution and the discovery of genes regulated during bud dormancy release.

BACKGROUND: Many northern-hemisphere forests are dominated by oaks. These species extend over diverse environmental conditions and are thus interesting models for studies of plant adaptation and speciation. The genomic toolbox is an important asset for exploring the functional variation associated with natural selection. RESULTS: The assembly of previously available and newly developed long and short sequence reads for two sympatric oak species, Quercus robur and Quercus petraea, generated a comprehensive catalog of transcripts for oak. The functional annotation of 91 k contigs demonstrated the presence of a large proportion of plant genes in this unigene set. Comparisons with SwissProt accessions and five plant gene models revealed orthologous relationships, making it possible to decipher the evolution of the oak genome. In particular, it was possible to align 9.5 thousand oak coding sequences with the equivalent sequences on peach chromosomes. Finally, RNA-seq data shed new light on the gene networks underlying vegetative bud dormancy release, a key stage in development allowing plants to adapt their phenology to the environment. CONCLUSION: In addition to providing a vast array of expressed genes, this study generated essential information about oak genome evolution and the regulation of genes associated with vegetative bud phenology, an important adaptive traits in trees. This resource contributes to the annotation of the oak genome sequence and will provide support for forward genetics approaches aiming to link genotypes with adaptive phenotypes.

A unigene set for European beech (Fagus sylvatica L.) and its use to decipher the molecular mechanisms involved in dormancy regulation.

Systematic sequencing is the method of choice for generating genomic resources for molecular marker development and candidate gene identification in nonmodel species. We generated 47,357 Sanger ESTs and 2.2M Roche-454 reads from five cDNA libraries for European beech (Fagus sylvatica L.). This tree species of high ecological and economic value in Europe is among the most representative trees of deciduous broadleaf forests. The sequences generated were assembled into 21,057 contigs with MIRA software. Functional annotations were obtained for 85% of these contigs, from the proteomes of four plant species, Swissprot accessions and the Gene Ontology database. We were able to identify 28,079 in silico SNPs for future marker development. Moreover, RNAseq and qPCR approaches identified genes and gene networks regulated differentially between two critical phenological stages preceding vegetative bud burst (the quiescent and swelling buds stages). According to climatic model-based projection, some European beech populations may be endangered, particularly at the southern and eastern edges of the European distribution range, which are strongly affected by current climate change. This first genomic resource for the genus Fagus should facilitate the identification of key genes for beech adaptation and management strategies for preserving beech adaptability.

Improved saccharification and ethanol yield from field-grown transgenic poplar deficient in cinnamoyl-CoA reductase.

Lignin is one of the main factors determining recalcitrance to enzymatic processing of lignocellulosic biomass. Poplars (Populus tremula x Populus alba) down-regulated for cinnamoyl-CoA reductase (CCR), the enzyme catalyzing the first step in the monolignol-specific branch of the lignin biosynthetic pathway, were grown in field trials in Belgium and France under short-rotation coppice culture. Wood samples were classified according to the intensity of the red xylem coloration typically associated with CCR down-regulation. Saccharification assays under different pretreatment conditions (none, two alkaline, and one acid pretreatment) and simultaneous saccharification and fermentation assays showed that wood from the most affected transgenic trees had up to 161% increased ethanol yield. Fermentations of combined material from the complete set of 20-mo-old CCR-down-regulated trees, including bark and less efficiently down-regulated trees, still yielded ∼ 20% more ethanol on a weight basis. However, strong down-regulation of CCR also affected biomass yield. We conclude that CCR down-regulation may become a successful strategy to improve biomass processing if the variability in down-regulation and the yield penalty can be overcome.

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.

Reference genes for high-throughput quantitative reverse transcription-PCR analysis of gene expression in organs and tissues of Eucalyptus grown in various environmental conditions.

Interest in the genomics of Eucalyptus has skyrocketed thanks to the recent sequencing of the genome of Eucalyptus grandis and to a growing number of large-scale transcriptomic studies. Quantitative reverse transcription-PCR (RT-PCR) is the method of choice for gene expression analysis and can now also be used as a high-throughput method. The selection of appropriate internal controls is becoming of utmost importance to ensure accurate expression results in Eucalyptus. To this end, we selected 21 candidate reference genes and used high-throughput microfluidic dynamic arrays to assess their expression among a large panel of developmental and environmental conditions with a special focus on wood-forming tissues. We analyzed the expression stability of these genes by using three distinct statistical algorithms (geNorm, NormFinder and ΔCt), and used principal component analysis to compare methods and rankings. We showed that the most stable genes identified depended not only on the panel of biological samples considered but also on the statistical method used. We then developed a comprehensive integration of the rankings generated by the three methods and identified the optimal reference genes for 17 distinct experimental sets covering 13 organs and tissues, as well as various developmental and environmental conditions. The expression patterns of Eucalyptus master genes EgMYB1 and EgMYB2 experimentally validated our selection. Our findings provide an important resource for the selection of appropriate reference genes for accurate and reliable normalization of gene expression data in the organs and tissues of Eucalyptus trees grown in a range of conditions including abiotic stresses.

Integrating genome annotation and QTL position to identify candidate genes for productivity, architecture and water-use efficiency in Populus spp.

BACKGROUND: Hybrid poplars species are candidates for biomass production but breeding efforts are needed to combine productivity and water use efficiency in improved cultivars. The understanding of the genetic architecture of growth in poplar by a Quantitative Trait Loci (QTL) approach can help us to elucidate the molecular basis of such integrative traits but identifying candidate genes underlying these QTLs remains difficult. Nevertheless, the increase of genomic information together with the accessibility to a reference genome sequence (Populus trichocarpa Nisqually-1) allow to bridge QTL information on genetic maps and physical location of candidate genes on the genome. The objective of the study is to identify QTLs controlling productivity, architecture and leaf traits in a P. deltoides x P. trichocarpa F1 progeny and to identify candidate genes underlying QTLs based on the anchoring of genetic maps on the genome and the gene ontology information linked to genome annotation. The strategy to explore genome annotation was to use Gene Ontology enrichment tools to test if some functional categories are statistically over-represented in QTL regions. RESULTS: Four leaf traits and 7 growth traits were measured on 330 F1 P. deltoides x P. trichocarpa progeny. A total of 77 QTLs controlling 11 traits were identified explaining from 1.8 to 17.2% of the variation of traits. For 58 QTLs, confidence intervals could be projected on the genome. An extended functional annotation was built based on data retrieved from the plant genome database Phytozome and from an inference of function using homology between Populus and the model plant Arabidopsis. Genes located within QTL confidence intervals were retrieved and enrichments in gene ontology (GO) terms were determined using different methods. Significant enrichments were found for all traits. Particularly relevant biological processes GO terms were identified for QTLs controlling number of sylleptic branches: intervals were enriched in GO terms of biological process like 'ripening' and 'adventitious roots development'. CONCLUSION: Beyond the simple identification of QTLs, this study is the first to use a global approach of GO terms enrichment analysis to fully explore gene function under QTLs confidence intervals in plants. This global approach may lead to identification of new candidate genes for traits of interest.

Expression analysis of LIM gene family in poplar, toward an updated phylogenetic classification.

BACKGROUND: Plant LIM domain proteins may act as transcriptional activators of lignin biosynthesis and/or as actin binding and bundling proteins. Plant LIM genes have evolved in phylogenetic subgroups differing in their expression profiles: in the whole plant or specifically in pollen. However, several poplar PtLIM genes belong to uncharacterized monophyletic subgroups and the expression patterns of the LIM gene family in a woody plant have not been studied. FINDINGS: In this work, the expression pattern of the twelve duplicated poplar PtLIM genes has been investigated by semi quantitative RT-PCR in different vegetative and reproductive tissues. As in other plant species, poplar PtLIM genes were widely expressed in the tree or in particular tissues. Especially, PtXLIM1a, PtXLIM1b and PtWLIM1b genes were preferentially expressed in the secondary xylem, suggesting a specific function in wood formation. Moreover, the expression of these genes and of the PtPLIM2a gene was increased in tension wood. Western-blot analysis confirmed the preferential expression of PtXLIM1a protein during xylem differentiation and tension wood formation. Genes classified within the pollen specific PLIM2 and PLIM2-like subgroups were all strongly expressed in pollen but also in cottony hairs. Interestingly, pairs of duplicated PtLIM genes exhibited different expression patterns indicating subfunctionalisations in specific tissues. CONCLUSIONS: The strong expression of several LIM genes in cottony hairs and germinating pollen, as well as in xylem fibers suggests an involvement of plant LIM domain proteins in the control of cell expansion. Comparisons of expression profiles of poplar LIM genes with the published functions of closely related plant LIM genes suggest conserved functions in the areas of lignin biosynthesis, pollen tube growth and mechanical stress response. Based on these results, we propose a novel nomenclature of poplar LIM domain proteins.