ABSTRACT
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.
Subject(s)
Quercus , Wood , Humans , Quercus/genetics , Ecosystem , Forests , Trees/geneticsABSTRACT
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.
Subject(s)
Eucalyptus/physiology , Plant Stems/physiology , Trees/physiology , Circadian Rhythm/physiology , Environment , Eucalyptus/anatomy & histology , Eucalyptus/genetics , Plant Stems/anatomy & histology , Plant Stems/genetics , Quantitative Trait Loci/genetics , Trees/anatomy & histology , Trees/genetics , Water/metabolismABSTRACT
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.
Subject(s)
Computational Biology/methods , Gene Expression Profiling/methods , Gene Regulatory Networks , Populus/growth & development , Gene Expression Regulation, Developmental , Gene Expression Regulation, Plant , Genotype , Machine Learning , Phenotype , Plant Proteins/genetics , Populus/genetics , Quantitative Trait Loci , Sequence Analysis, RNAABSTRACT
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.
Subject(s)
Plant Somatic Embryogenesis Techniques/methods , Pseudotsuga/embryology , Seeds/growth & development , Gene Regulatory Networks , Mass Spectrometry , Plant Proteins/metabolism , Plant Proteins/physiology , Proteomics , Pseudotsuga/growth & development , Pseudotsuga/metabolism , Seeds/metabolismABSTRACT
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.
Subject(s)
Mucoproteins/analysis , Pectins/analysis , Polysaccharides/analysis , Populus/metabolism , Cell Wall/metabolism , Mannans/analysis , Mannans/metabolism , Mucoproteins/metabolism , Pectins/metabolism , Plant Proteins/analysis , Plant Proteins/metabolism , Polysaccharides/metabolism , Populus/cytology , Populus/growth & development , Trees , Wood/cytology , Wood/genetics , Wood/metabolism , Xylem/cytology , Xylem/growth & development , Xylem/metabolismABSTRACT
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.
Subject(s)
Aldehyde Oxidoreductases/deficiency , Biofuels , Ethanol/metabolism , Lignin/metabolism , Populus/metabolism , Belgium , Biomass , Fermentation , France , Plants, Genetically Modified , Populus/geneticsABSTRACT
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.
Subject(s)
Gene Expression Regulation, Plant , Phloem/growth & development , Plant Proteins/metabolism , Populus/growth & development , Cell Wall/metabolism , Chimera , Molecular Sequence Data , Phenotype , Phloem/genetics , Phloem/metabolism , Plant Proteins/genetics , Plant Stems/genetics , Plant Stems/metabolism , Plants, Genetically Modified , Populus/geneticsABSTRACT
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.
Subject(s)
Phosphoproteins/metabolism , Plant Proteins/metabolism , Populus/metabolism , Proteome/metabolism , Proteomics/methods , Wood/metabolism , Amino Acid Sequence , Cluster Analysis , Gene Expression Regulation, Plant , Gene Ontology , Gene Regulatory Networks , Gravitation , Gravitropism , Mass Spectrometry , Molecular Sequence Data , Peptides/genetics , Peptides/metabolism , Phosphopeptides/genetics , Phosphopeptides/metabolism , Phosphoproteins/genetics , Plant Proteins/classification , Plant Proteins/genetics , Populus/genetics , Proteome/classification , Proteome/genetics , Signal Transduction/genetics , Wood/genetics , Xylem/genetics , Xylem/metabolismABSTRACT
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.
Subject(s)
Adaptation, Physiological/genetics , Gene Expression Regulation, Plant , Plant Dormancy/genetics , Transcriptome/genetics , Base Sequence , Chromosome Mapping , Genetic Speciation , Genome, Plant , Quercus/genetics , Quercus/growth & development , Sequence Analysis, RNAABSTRACT
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.
Subject(s)
Gene Expression Regulation, Plant , Plant Physiological Phenomena/physiology , Quercus/genetics , Quercus/physiology , Transcriptome , Abscisic Acid/physiology , Base Sequence , Cold Temperature , Databases, Genetic , Gene Expression Profiling , Gene Library , Genes, Plant , Gibberellins/genetics , Globulins/genetics , Glycosyltransferases/genetics , Quercus/growth & development , RNA, Plant/genetics , Sequence Analysis, RNA , Stress, Physiological/genetics , Up-RegulationABSTRACT
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.
Subject(s)
Eucalyptus/genetics , Gene Expression Profiling/standards , Genes, Plant/genetics , Reverse Transcriptase Polymerase Chain Reaction/standards , Algorithms , Cold Temperature , DNA Primers/genetics , Droughts , Environment , Eucalyptus/growth & development , Eucalyptus/physiology , Fertilization , Gene Expression , Gene Expression Regulation, Plant , Organ Specificity , Reference Standards , Stress, Physiological , Xylem/genetics , Xylem/growth & development , Xylem/physiologyABSTRACT
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.
Subject(s)
Genetic Association Studies , Genome, Plant/genetics , Molecular Sequence Annotation , Populus/growth & development , Populus/genetics , Quantitative Trait Loci/genetics , Water/metabolism , Chromosome Mapping , Confidence Intervals , Crosses, Genetic , Genes, Plant/genetics , Genetic Linkage , Genetic Markers , Plant Leaves/genetics , Populus/anatomy & histology , Quantitative Trait, HeritableABSTRACT
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.
ABSTRACT
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.
Subject(s)
Gene Expression/genetics , Populus/enzymology , Wood/enzymology , Xylosidases/chemistry , Xylosidases/metabolism , Amino Acid Sequence , Cloning, Molecular , DNA, Plant/genetics , Gene Expression Regulation, Enzymologic , Mass Spectrometry , Molecular Sequence Data , Molecular Weight , Phylogeny , Populus/genetics , Wood/genetics , Xylosidases/geneticsABSTRACT
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.
Subject(s)
Gravitropism/physiology , Plant Proteins/metabolism , Plant Stems/metabolism , Populus/physiology , Proteomics , Electrophoresis, Gel, Two-Dimensional , Gene Expression Regulation, Plant , Plant Stems/physiology , Populus/metabolism , Spectrometry, Mass, Matrix-Assisted Laser Desorption-Ionization , Tandem Mass SpectrometryABSTRACT
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.
ABSTRACT
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.
Subject(s)
Genome, Plant/genetics , Quercus/genetics , Biological Evolution , DNA, Plant/genetics , Genetic Variation/genetics , Longevity/genetics , Mutation , Phylogeny , Sequence Analysis, DNAABSTRACT
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.
Subject(s)
Arabidopsis/genetics , Genes, Plant/genetics , Genome, Plant/genetics , Multigene Family/genetics , Oryza/genetics , Populus/genetics , Amino Acid Sequence , Consensus Sequence , DNA, Complementary/isolation & purification , Expressed Sequence Tags , Gene Library , Models, Genetic , Molecular Sequence Data , Phylogeny , Sequence Analysis, Protein , Structure-Activity RelationshipABSTRACT
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.
Subject(s)
Lignin/genetics , Lignin/metabolism , Plants, Genetically Modified , Populus/genetics , Populus/metabolism , Wood , Alcohol Oxidoreductases/genetics , Alcohol Oxidoreductases/metabolism , Antisense Elements (Genetics) , Ecosystem , Environmental Pollution/prevention & control , France , Gene Expression , Industrial Waste/prevention & control , Methyltransferases/genetics , Methyltransferases/metabolism , Models, Chemical , Paper , Plant Roots/metabolism , Reproducibility of Results , Sensitivity and Specificity , Soil Microbiology , Transgenes , United KingdomABSTRACT
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.