Nathalie Pavy

Involvement of Pinus taeda MYB1 and MYB8 in phenylpropanoid metabolism and secondary cell wall biogenesis: a comparative in planta analysis

The involvement of two R2R3-MYB genes from Pinus taeda L., PtMYB1 and PtMYB8, in phenylpropanoid metabolism and secondary cell wall biogenesis was investigated in planta. These pine MYBs were constitutively overexpressed (OE) in Picea glauca (Moench) Voss, used as a heterologous conifer expression system. Morphological, histological, chemical (lignin and soluble phenols), and transcriptional analyses, i.e. microarray and reverse transcription quantitative PCR (RT-qPCR) were used for extensive phenotyping of MYB-overexpressing spruce plantlets. Upon germination of somatic embryos, root growth was reduced in both transgenics. Enhanced lignin deposition was also a common feature but ectopic secondary cell wall deposition was more strongly associated with PtMYB8-OE. Microarray and RT-qPCR data showed that overexpression of each MYB led to an overlapping up-regulation of many genes encoding phenylpropanoid enzymes involved in lignin monomer synthesis, while misregulation of several cell wall-related genes and other MYB transcription factors was specifically associated with PtMYB8-OE. Together, the results suggest that MYB1 and MYB8 may be part of a conserved transcriptional network involved in secondary cell wall deposition in conifers.

Enhancing genetic mapping of complex genomes through the design of highly-multiplexed SNP arrays: application to the large and unsequenced genomes of white spruce and black spruce

BackgroundTo explore the potential value of high-throughput genotyping assays in the analysis of large and complex genomes, we designed two highly multiplexed Illumina bead arrays using the GoldenGate SNP assay for gene mapping in white spruce (Picea glauca [Moench] Voss) and black spruce (Picea mariana [Mill.] B.S.P.).ResultsEach array included 768 SNPs, identified by resequencing genomic DNA from parents of each mapping population. For white spruce and black spruce, respectively, 69.2% and 77.1% of genotyped SNPs had valid GoldenGate assay scores and segregated in the mapping populations. For each of these successful SNPs, on average, valid genotyping scores were obtained for over 99% of progeny. SNP data were integrated to pre-existing ALFP, ESTP, and SSR markers to construct two individual linkage maps and a composite map for white spruce and black spruce genomes. The white spruce composite map contained 821 markers including 348 gene loci. Also, 835 markers including 328 gene loci were positioned on the black spruce composite map. In total, 215 anchor markers (mostly gene markers) were shared between the two species. Considering lineage divergence at least 10 Myr ago between the two spruces, interspecific comparison of homoeologous linkage groups revealed remarkable synteny and marker colinearity.ConclusionThe design of customized highly multiplexed Illumina SNP arrays appears as an efficient procedure to enhance the mapping of expressed genes and make linkage maps more informative and powerful in such species with poorly known genomes. This genotyping approach will open new avenues for co-localizing candidate genes and QTLs, partial genome sequencing, and comparative mapping across conifers.

EgMYB1, an R2R3 MYB transcription factor from eucalyptus negatively regulates secondary cell wall formation in Arabidopsis and poplar

• The eucalyptus R2R3 transcription factor, EgMYB1 contains an active repressor motif in the regulatory domain of the predicted protein. It is preferentially expressed in differentiating xylem and is capable of repressing the transcription of two key lignin genes in vivo. • In order to investigate in planta the role of this putative transcriptional repressor of the lignin biosynthetic pathway, we overexpressed the EgMYB1 gene in Arabidopsis and poplar. • Expression of EgMYB1 produced similar phenotypes in both species, with stronger effects in transgenic Arabidopsis plants than in poplar. Vascular development was altered in overexpressors showing fewer lignified fibres (in phloem and interfascicular zones in poplar and Arabidopsis, respectively) and reduced secondary wall thickening. Klason lignin content was moderately but significantly reduced in both species. Decreased transcript accumulation was observed for genes involved in the biosynthesis of lignins, cellulose and xylan, the three main polymers of secondary cell walls. Transcriptomic profiles of transgenic poplars were reminiscent of those reported when lignin biosynthetic genes are disrupted. • Together, these results strongly suggest that EgMYB1 is a repressor of secondary wall formation and provide new opportunities to dissect the transcriptional regulation of secondary wall biosynthesis.

Identification of conserved core xylem gene sets: conifer cDNA microarray development, transcript profiling and computational analyses

One approach for investigating the molecular basis of wood formation is to integrate microarray profiling data sets and sequence analyses, comparing tree species with model plants such as Arabidopsis. Conifers may be included in comparative studies thanks to large-scale expressed sequence tag (EST) analyses, which enable the development of cDNA microarrays with very significant genome coverage. A microarray of 10,400 low-redundancy sequences was designed starting from white spruce (Picea glauca (Moench.) Voss) cDNAs. Computational procedures that were developed to ensure broad transcriptome coverage and efficient PCR amplification were used to select cDNA clones, which were re-sequenced in the microarray manufacture process. White spruce transcript profiling experiments that compared secondary xylem to phloem and needles identified 360 xylem-preferential gene sequences. The functional annotations of all differentially expressed sequences were highly consistent with the results of similar analyses carried out in angiosperm trees and herbaceous plants. Computational analyses comparing the spruce microarray sequences and core xylem gene sets from Arabidopsis identified 31 transcripts that were highly conserved in angiosperms and gymnosperms, in terms of both sequence and xylem expression. Several other spruce sequences have not previously been linked to xylem differentiation (including genes encoding TUBBY-like domain proteins (TLPs) and a gibberellin insensitive (gai) gene sequence) or were shown to encode proteins of unknown function encompassing diverse conserved domains of unknown function.

Towards decoding the conifer giga-genome.

Several new initiatives have been launched recently to sequence conifer genomes including pines, spruces and Douglas-fir. Owing to the very large genome sizes ranging from 18 to 35 gigabases, sequencing even a single conifer genome had been considered unattainable until the recent throughput increases and cost reductions afforded by next generation sequencers. The purpose of this review is to describe the context for these new initiatives. A knowledge foundation has been acquired in several conifers of commercial and ecological interest through large-scale cDNA analyses, construction of genetic maps and gene mapping studies aiming to link phenotype and genotype. Exploratory sequencing in pines and spruces have pointed out some of the unique properties of these giga-genomes and suggested strategies that may be needed to extract value from their sequencing. The hope is that recent and pending developments in sequencing technology will contribute to rapidly filling the knowledge vacuum surrounding their structure, contents and evolution. Researchers are also making plans to use comparative analyses that will help to turn the data into a valuable resource for enhancing and protecting the world’s conifer forests.

A spruce gene map infers ancient plant genome reshuffling and subsequent slow evolution in the gymnosperm lineage leading to extant conifers

BackgroundSeed plants are composed of angiosperms and gymnosperms, which diverged from each other around 300 million years ago. While much light has been shed on the mechanisms and rate of genome evolution in flowering plants, such knowledge remains conspicuously meagre for the gymnosperms. Conifers are key representatives of gymnosperms and the sheer size of their genomes represents a significant challenge for characterization, sequencing and assembling.ResultsTo gain insight into the macro-organisation and long-term evolution of the conifer genome, we developed a genetic map involving 1,801 spruce genes. We designed a statistical approach based on kernel density estimation to analyse gene density and identified seven gene-rich isochors. Groups of co-localizing genes were also found that were transcriptionally co-regulated, indicative of functional clusters. Phylogenetic analyses of 157 gene families for which at least two duplicates were mapped on the spruce genome indicated that ancient gene duplicates shared by angiosperms and gymnosperms outnumbered conifer-specific duplicates by a ratio of eight to one. Ancient duplicates were much more translocated within and among spruce chromosomes than conifer-specific duplicates, which were mostly organised in tandem arrays. Both high synteny and collinearity were also observed between the genomes of spruce and pine, two conifers that diverged more than 100 million years ago.ConclusionsTaken together, these results indicate that much genomic evolution has occurred in the seed plant lineage before the split between gymnosperms and angiosperms, and that the pace of evolution of the genome macro-structure has been much slower in the gymnosperm lineage leading to extent conifers than that seen for the same period of time in flowering plants. This trend is largely congruent with the contrasted rates of diversification and morphological evolution observed between these two groups of seed plants.

Development of high-density SNP genotyping arrays for white spruce (Picea glauca) and transferability to subtropical and nordic congeners.

High‐density SNP genotyping arrays can be designed for any species given sufficient sequence information of high quality. Two high‐density SNP arrays relying on the Infinium iSelect technology (Illumina) were designed for use in the conifer white spruce (Picea glauca). One array contained 7338 segregating SNPs representative of 2814 genes of various molecular functional classes for main uses in genetic association and population genetics studies. The other one contained 9559 segregating SNPs representative of 9543 genes for main uses in population genetics, linkage mapping of the genome and genomic prediction. The SNPs assayed were discovered from various sources of gene resequencing data. SNPs predicted from high‐quality sequences derived from genomic DNA reached a genotyping success rate of 64.7%. Nonsingleton in silico SNPs (i.e. a sequence polymorphism present in at least two reads) predicted from expressed sequenced tags obtained with the Roche 454 technology and Illumina GAII analyser resulted in a similar genotyping success rate of 71.6% when the deepest alignment was used and the most favourable SNP probe per gene was selected. A variable proportion of these SNPs was shared by other nordic and subtropical spruce species from North America and Europe. The number of shared SNPs was inversely proportional to phylogenetic divergence and standing genetic variation in the recipient species, but positively related to allele frequency in P. glauca natural populations. These validated SNP resources should open up new avenues for population genetics and comparative genetic mapping at a genomic scale in spruce species.

Large-scale statistical analysis of secondary xylem ESTs in pine.

A computational analysis of pine transcripts was conducted to contribute to the functional annotation of conifer sequences. A statistical analysis of expressed sequential tags(ESTs) belonging the 7732 contigs in the TIGR Pinus Gene Index (PGI1.0) identified 260 differentially represented gene sequences across six cDNA libraries from loblolly pine secondary xylem. Cluster analysis of this subset of contigs resulted in five groups representing genes preferentially represented in one of the xylem samples (compression wood, plannings, root xylem, latewood) and one group containing mostly genes simultaneously present in compression and side wood libraries. To complement the sequence annotation, 27 cDNA clones representing selected transcripts were completely sequenced. Several genes were identified that could represent putative markers for xylem from different organs, at different stages of development. Several sequences encoding regulatory proteins were over-represented in root xylem as opposed to the other xylem samples. Some of them belonged to known families of plant transcription factors, but two genes were previously uncharacterized in plants. One transcript was homologous to the gene encoding the Smad4 interacting factor, a key co-activator in TGFβ (transforming growth factor) signalling in animals. Thus, the digital analysis of pine ESTs highlighted a putative gene function of potentially broad interest but that has yet to be investigated in plants. More generally, this study showed that the application of numerical approaches to EST databases should be helpful in establishing priorities among genes to consider for targeted functional studies. Thus, we illustrated the potential of extracting information from conifer sequences already accessible through well-structured public databases.

The heterogeneous levels of linkage disequilibrium in white spruce genes and comparative analysis with other conifers

In plants, knowledge about linkage disequilibrium (LD) is relevant for the design of efficient single-nucleotide polymorphism arrays in relation to their use in population and association genomics studies. Previous studies of conifer genes have shown LD to decay rapidly within gene limits, but exceptions have been reported. To evaluate the extent of heterogeneity of LD among conifer genes and its potential causes, we examined LD in 105 genes of white spruce (Picea glauca) by sequencing a panel of 48 haploid megagametophytes from natural populations and further compared it with LD in other conifer species. The average pairwise r2 value was 0.19 (s.d.=0.19), and LD dropped quickly with a half-decay being reached at a distance of 65 nucleotides between sites. However, LD was significantly heterogeneous among genes. A first group of 29 genes had stronger LD (mean r2=0.28), and a second group of 38 genes had weaker LD (mean r2=0.12). While a strong relationship was found with the recombination rate, there was no obvious relationship between LD and functional classification. The level of nucleotide diversity, which was highly heterogeneous across genes, was also not significantly correlated with LD. A search for selection signatures highlighted significant deviations from the standard neutral model, which could be mostly attributed to recent demographic changes. Little evidence was seen for hitchhiking and clear relationships with LD. When compared among conifer species, on average, levels of LD were similar in genes from white spruce, Norway spruce and Scots pine, whereas loblolly pine and Douglas fir genes exhibited a significantly higher LD.

High nitrogen fertilization and stem leaning have overlapping effects on wood formation in poplar but invoke largely distinct molecular pathways

Previous studies indicated that high nitrogen fertilization may impact secondary xylem development and alter fibre anatomy and composition. The resulting wood shares some resemblance with tension wood, which has much thicker cell walls than normal wood due to the deposition of an additional layer known as the G-layer. This report compares the short-term effects of high nitrogen fertilization and tree leaning to induce tension wood, either alone or in combination, upon wood formation in young trees of Populus trichocarpa (Torr. & Gray) × P. deltoides Bartr. ex Marsh. Fibre anatomy, chemical composition and transcript profiles were examined in newly formed secondary xylem. Each of the treatments resulted in thicker cell walls relative to the controls. High nitrogen and tree leaning had overlapping effects on chemical composition based on Fourier transform infrared analysis, specifically indicating that secondary cell wall composition was shifted in favour of cellulose and hemicelluloses relative to lignin content. In contrast, the high-nitrogen trees had shorter fibres, whilst the leaning trees had longer fibres that the controls. Microarray transcript profiling carried out after 28 days of treatment identified 180 transcripts that accumulated differentially in one or more treatments. Only 10% of differentially expressed transcripts were affected in all treatments relative to the controls. Several of the affected transcripts were related to carbohydrate metabolism, secondary cell wall formation, nitrogen metabolism and osmotic stress. RT-qPCR analyses at 1, 7 and 28 days showed that several transcripts followed very different accumulation profiles in terms of rate and level of accumulation, depending on the treatment. Our findings suggest that high nitrogen fertilization and tension wood induction elicit largely distinct and molecular pathways with partial overlap. When combined, the two types of environmental cue yielded additive effects.