Ronald R. Sederoff

High-throughput gene and SNP discovery in Eucalyptus grandis, an uncharacterized genome.

BackgroundBenefits from high-throughput sequencing using 454 pyrosequencing technology may be most apparent for species with high societal or economic value but few genomic resources. Rapid means of gene sequence and SNP discovery using this novel sequencing technology provide a set of baseline tools for genome-level research. However, it is questionable how effective the sequencing of large numbers of short reads for species with essentially no prior gene sequence information will support contig assemblies and sequence annotation.ResultsWith the purpose of generating the first broad survey of gene sequences in Eucalyptus grandis, the most widely planted hardwood tree species, we used 454 technology to sequence and assemble 148 Mbp of expressed sequences (EST). EST sequences were generated from a normalized cDNA pool comprised of multiple tissues and genotypes, promoting discovery of homologues to almost half of Arabidopsis genes, and a comprehensive survey of allelic variation in the transcriptome. By aligning the sequencing reads from multiple genotypes we detected 23,742 SNPs, 83% of which were validated in a sample. Genome-wide nucleotide diversity was estimated for 2,392 contigs using a modified theta (θ) parameter, adapted for measuring genetic diversity from polymorphisms detected by randomly sequencing a multi-genotype cDNA pool. Diversity estimates in non-synonymous nucleotides were on average 4x smaller than in synonymous, suggesting purifying selection. Non-synonymous to synonymous substitutions (Ka/Ks) among 2,001 contigs averaged 0.30 and was skewed to the right, further supporting that most genes are under purifying selection. Comparison of these estimates among contigs identified major functional classes of genes under purifying and diversifying selection in agreement with previous researches.ConclusionIn providing an abundance of foundational transcript sequences where limited prior genomic information existed, this work created part of the foundation for the annotation of the E. grandis genome that is being sequenced by the US Department of Energy. In addition we demonstrated that SNPs sampled in large-scale with 454 pyrosequencing can be used to detect evolutionary signatures among genes, providing one of the first genome-wide assessments of nucleotide diversity and Ka/Ks for a non-model plant species.

Variation in Lignin Content and Composition (Mechanisms of Control and Implications for the Genetic Improvement of Plants)

Lignin, a complex phenolic polymer, is important for mechanical support, water transport, and defense in vaseular plants. Compressive strength and hydrophobicity of xylem cell walls are imparted by the lignin polymer, which is deposited during the terminal differentiation of tracheids and other cell types. The resistance of xylem to compressive stresses imposed by water transport and by the mass of the plants is important to growth and development. In addition, the insolubility and complexity of the lignin polymer makes it resistant to degradation by most microorganisms. Therefore, lignin serves an important function in plant defense. Variation in lignin content, composition, and location is likely to affect these essential processes. The constraints on the amount, composition, and localization of lignin for normal xylem function and plant defense are not known. Lignin composition, quantity, and distribution also affect the agroindustrial uses of plant material. Digestibility and dietary conversion of herbaceous crops are affected by differences in lignin content and composition (Akin et al., 1986, 1991). Lignin is an undesirable component in the conversion of wood into pulp and paper; remova1 of lignin is a major step in the paper making process. Furthermore, the resistance of lignin to microbial degradation enhances its persistence in soils. Lignin is, therefore, a significant component in the global carbon cycle. The mechanisms of control of lignin composition and quantity have wide implications regarding the adaptation and evolution of land plants and provide a basis for improved genetic manipulation of lignin for agroindustrial end uses. In this Update, we will focus on the levels of control of lignin variation, including (a) metabolic control, (b) regulation of individual enzymes in the biosynthetic pathway, and (c) regulation of gene expression. These levels of regulation affect variation in lignin content, quality, and distribution. Finally, the implications of these regulatory mechanisms for the genetic improvement of lignin for

A laccase associated with lignification in loblolly pine xylem

Peroxidase has been thought to be the only enzyme that oxidizes monolignol precursors to initiate lignin formation in plants. A laccase was purified from cell walls of differentiating xylem of loblolly pine and shown to coincide in time and place with lignin formation and to oxidize monolignols to dehydrogenation products in vitro. These results suggest that laccase participates in lignin biosynthesis and therefore could be an important target for genetic engineering to modify wood properties or to improve the digestibility of forage crops.

Improved estimation of DNA fragment lengths from agarose gels

An improved method for estimation of DNA fragment lengths from mobility on agarose gels is introduced. Commonly, this estimation is done by graphing the logarithm of length against the mobility of standards, although this procedure produces a noticeably curved plot. The method presented here is based on the relationship of E. M. Southern (1979, Anal. Biochem., 100, 319–323), who showed that the reciprocal of mobility plotted against fragment length is linear. A least-squares analysis has been developed to use this relationship together with any number of standards to estimate fragment length. The method also tests internal consistency of the assigned lengths of a set of fragments. This analysis has been applied to sets of restriction fragments of λ DNA, indicating errors or inconsistencies in previously assigned lengths. Using a more consistent set of λ values, the lengths of pBR322 fragments of known sequence were predicted with better than 1% accuracy. No evidence of sequence-dependent migration was seen. A computer program to do this analysis is included.

Decoding the massive genome of loblolly pine using haploid DNA and novel assembly strategies

BackgroundThe size and complexity of conifer genomes has, until now, prevented full genome sequencing and assembly. The large research community and economic importance of loblolly pine, Pinus taeda L., made it an early candidate for reference sequence determination.ResultsWe develop a novel strategy to sequence the genome of loblolly pine that combines unique aspects of pine reproductive biology and genome assembly methodology. We use a whole genome shotgun approach relying primarily on next generation sequence generated from a single haploid seed megagametophyte from a loblolly pine tree, 20-1010, that has been used in industrial forest tree breeding. The resulting sequence and assembly was used to generate a draft genome spanning 23.2 Gbp and containing 20.1 Gbp with an N50 scaffold size of 66.9 kbp, making it a significant improvement over available conifer genomes. The long scaffold lengths allow the annotation of 50,172 gene models with intron lengths averaging over 2.7 kbp and sometimes exceeding 100 kbp in length. Analysis of orthologous gene sets identifies gene families that may be unique to conifers. We further characterize and expand the existing repeat library based on the de novo analysis of the repetitive content, estimated to encompass 82% of the genome.ConclusionsIn addition to its value as a resource for researchers and breeders, the loblolly pine genome sequence and assembly reported here demonstrates a novel approach to sequencing the large and complex genomes of this important group of plants that can now be widely applied.

Coordinated genetic regulation of growth and lignin revealed by quantitative trait locus analysis of cDNA microarray data in an interspecific backcross of eucalyptus.

Phenotypic, genotypic, and transcript level (microarray) data from an interspecific backcross population of Eucalyptus grandis and Eucalyptus globulus were integrated to dissect the genetic and metabolic network underlying growth variation. Transcript abundance, measured for 2,608 genes in the differentiating xylem of a 91 (E. grandis × E. globulus) × E. grandis backcross progeny was correlated with diameter variation, revealing coordinated down-regulation of genes encoding enzymes of the lignin biosynthesis and associated methylation pathways in fast growing individuals. Lignin analysis of wood samples confirmed the content and quality predicted by the transcript levels measured on the microarrays. Quantitative trait locus (QTL) analysis of transcript levels of lignin-related genes showed that their mRNA abundance is regulated by two genetic loci, demonstrating coordinated genetic control over lignin biosynthesis. These two loci colocalize with QTLs for growth, suggesting that the same genomic regions are regulating growth, and lignin content and composition in the progeny. Genetic mapping of the lignin genes revealed that most of the key biosynthetic genes do not colocalize with growth and transcript level QTLs, with the exception of the locus encoding the enzyme S-adenosylmethionine synthase. This study illustrates the power of integrating quantitative analysis of gene expression data and genetic map information to discover genetic and metabolic networks regulating complex biological traits.

Towards a Systems Approach for Lignin Biosynthesis in Populus trichocarpa: Transcript Abundance and Specificity of the Monolignol Biosynthetic Genes

As a step toward a comprehensive description of lignin biosynthesis in Populus trichocarpa, we identified from the genome sequence 95 phenylpropanoid gene models in 10 protein families encoding enzymes for monolignol biosynthesis. Transcript abundance was determined for all 95 genes in xylem, leaf, shoot and phloem using quantitative real-time PCR (qRT-PCR). We identified 23 genes that most probably encode monolignol biosynthesis enzymes during wood formation. Transcripts for 18 of the 23 are abundant and specific to differentiating xylem. We found evidence suggesting functional redundancy at the transcript level for phenylalanine ammonia-lyase (PAL), cinnamate 4-hydroxylase (C4H), 4-coumarate:CoA ligase (4CL), p-hydroxycinnamoyl-CoA:quinate shikimate p-hydroxycinnamoyltransferase (HCT), caffeoyl-CoA O-methyltransferase (CCoAOMT) and coniferyl aldehyde 5-hydroxylase (CAld5H). We carried out an enumeration-based motif identification and discriminant analysis on the promoters of all 95 genes. Five core motifs correctly discriminate the 18 xylem-specific genes from the 77 non-xylem genes. These motifs are similar to promoter elements known to regulate phenylpropanoid gene expression. This work suggests that genes in monolignol biosynthesis are regulated by multiple motifs, often related in sequence.

Unexpected variation in lignin

Recent studies on mutant and transgenic plants indicate that lignification may be far more flexible than previously realized. Pines with a mutation affecting the biosynthesis of the major lignin precursor, coniferyl alcohol, show a high level of an unusual subunit, dihydroconiferyl alcohol. These results argue in favor of an increased potential for genetic modification of lignin and indicate that our knowledge of the biosynthesis of lignin is far from complete.

Elucidation of new structures in lignins of CAD- and COMT-deficient plants by NMR.

Studying lignin-biosynthetic-pathway mutants and transgenics provides insights into plant responses to perturbations of the lignification system, and enhances our understanding of normal lignification. When enzymes late in the pathway are downregulated, significant changes in the composition and structure of lignin may result. NMR spectroscopy provides powerful diagnostic tools for elucidating structures in the difficult lignin polymer, hinting at the chemical and biochemical changes that have occurred. COMT (caffeic acid O-methyl transferase) downregulation in poplar results in the incorporation of 5-hydroxyconiferyl alcohol into lignins via typical radical coupling reactions, but post-coupling quinone methide internal trapping reactions produce novel benzodioxane units in the lignin. CAD (cinnamyl alcohol dehydrogenase) downregulation results in the incorporation of the hydroxycinnamyl aldehyde monolignol precursors intimately into the polymer. Sinapyl aldehyde cross-couples 8-O-4 with both guaiacyl and syringyl units in the growing polymer, whereas coniferyl aldehyde cross-couples 8-O-4 only with syringyl units, reflecting simple chemical cross-coupling propensities. The incorporation of hydroxycinnamyl aldehyde and 5-hydroxyconiferyl alcohol monomers indicates that these monolignol intermediates are secreted to the cell wall for lignification. The recognition that novel units can incorporate into lignins portends significantly expanded opportunities for engineering the composition and consequent properties of lignin for improved utilization of valuable plant resources.

Apparent homology of expressed genes from wood-forming tissues of loblolly pine (Pinus taeda L.) with Arabidopsis thaliana

Pinus taeda L. (loblolly pine) and Arabidopsis thaliana differ greatly in form, ecological niche, evolutionary history, and genome size. Arabidopsis is a small, herbaceous, annual dicotyledon, whereas pines are large, long-lived, coniferous forest trees. Such diverse plants might be expected to differ in a large number of functional genes. We have obtained and analyzed 59,797 expressed sequence tags (ESTs) from wood-forming tissues of loblolly pine and compared them to the gene sequences inferred from the complete sequence of the Arabidopsis genome. Approximately 50% of pine ESTs have no apparent homologs in Arabidopsis or any other angiosperm in public databases. When evaluated by using contigs containing long, high-quality sequences, we find a higher level of apparent homology between the inferred genes of these two species. For those contigs 1,100 bp or longer, ≈90% have an apparent Arabidopsis homolog (E value < 10-10). Pines and Arabidopsis last shared a common ancestor ≈300 million years ago. Few genes would be expected to retain high sequence similarity for this time if they did not have essential functions. These observations suggest substantial conservation of gene sequence in seed plants.