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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Cloning and expression analysis of a wood-associated xylosidase gene (PtaBXL1) in poplar tension wood.

In stems of woody angiosperms responding to mechanical stress, imposed for instance by tilting the stem or formation of a branch, tension wood (TW) forms above the affected part, while anatomically distinct opposite wood (OW) forms below it. In poplar TW the S3 layer of the secondary walls is substituted by a "gelatinous layer" that is almost entirely composed of cellulose and has much lower hemicellulose contents than unstressed wood. However, changes in xylan contents (the predominant hemicelluloses), their interactions with other wall components and the mechanisms involved in TW formation have been little studied. Therefore, in the study reported here we determined the structure and distribution of xylans, cloned the genes encoding the xylan remodeling enzymes beta-xylosidases (PtaBXLi), and examined their expression patterns during tension wood, normal wood and opposite wood xylogenesis in poplar. We confirm that poplar wood xylans are substituted solely by 4-O-methylglucuronic acid in both TW and OW. However, although glucuronoxylans are strongly represented in both primary and secondary layers of OW, no 4-O-methylGlcA xylan was found in G-layers of TW. Four full-length BXL cDNAs encoding putative beta-xylosidases were cloned. One, PtaBXL1, for which xylosidase activity was confirmed by heterologous expression in Escherichia coli, exhibited a wood-specific expression pattern in TW. In conclusion, xylan as PtaBXL1, encoding beta4-xylosidase activity, are down-regulated in TW.

Proteome analysis of apical and basal regions of poplar stems under gravitropic stimulation.

Gravity is a constant force guiding the direction of plant growth. In young poplar stem, reorientation of the apical region is mainly obtained by differential growth of elongating primary tissues. At the base, where elongation is achieved but where the cambium is active, reorientation is due to asymmetrical formation of reaction wood. After 45 min of gravistimulation, the stem showed no reorientation, but 1 week later, reaction wood was observed at the base of the stem. To determine the molecular mechanisms taking place at the top and base of the stem, after 45 min or 1 week of inclination, the changes induced in protein accumulation were studied by two-dimensional polyacrylamide gel electrophoresis and quantitatively analyzed using image analysis software. Around 300 protein spots were reproducibly detected and analyzed. Forty percent of these proteins showed significant changes after inclination. Mass spectrometry analysis of 135 spots led to the identification of 60 proteins involved in a wide range of activities and metabolisms. Very different patterns of protein expression were obtained according to conditions tested, highlighting the complexity of gravitropic responses. Our results suggest that primary and secondary tissues present specific mechanisms to sense reorientation and to respond to inclination. Some selected proteins are discussed.

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.

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.

Genome-wide analysis of LIM gene family in Populus trichocarpa, Arabidopsis thaliana, and Oryza sativa.

In Eukaryotes, LIM proteins act as developmental regulators in basic cellular processes such as regulating the transcription or organizing the cytoskeleton. The LIM domain protein family in plants has mainly been studied in sunflower and tobacco plants, where several of its members exhibit a specific pattern of expression in pollen. In this paper, we finely characterized in poplar six transcripts encoding these proteins. In Populus trichocarpa genome, the 12 LIM gene models identified all appear to be duplicated genes. In addition, we describe several new LIM domain proteins deduced from Arabidopsis and rice genomes, raising the number of LIM gene models to six for both species. Plant LIM genes have a core structure of four introns with highly conserved coding regions. We also identified new LIM domain proteins in several other species, and a phylogenetic analysis of plant LIM proteins reveals that they have undergone one or several duplication events during the evolution. We gathered several LIM protein members within new monophyletic groups. We propose to classify the plant LIM proteins into four groups: alphaLIM1, betaLIM1, gammaLIM2, and deltaLIM2, subdivided according to their specificity to a taxonomic class and/or to their tissue-specific expression. Our investigation of the structure of the LIM domain proteins revealed that they contain many conserved motifs potentially involved in their function.

Field and pulping performances of transgenic trees with altered lignification.

The agronomic and pulping performance of transgenic trees with altered lignin has been evaluated in duplicated, long-term field trials. Poplars expressing cinnamyl alcohol dehydrogenase (CAD) or caffeate/5-hydroxy-ferulate O-methyltransferase (COMT) antisense transgenes were grown for four years at two sites, in France and England. The trees remained healthy throughout the trial. Growth indicators and interactions with insects were normal. No changes in soil microbial communities were detected beneath the transgenic trees. The expected modifications to lignin were maintained in the transgenics over four years, at both sites. Kraft pulping of tree trunks showed that the reduced-CAD lines had improved characteristics, allowing easier delignification, using smaller amounts of chemicals, while yielding more high-quality pulp. This work highlights the potential of engineering wood quality for more environmentally benign papermaking without interfering with tree growth or fitness.

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.