Genomics, domestication, and evolution of forest trees.

The forests of the world continue to be threatened by climate change, population growth, and loss to agriculture. Our ability to conserve natural forests and to meet the increasing demand for fuel, biomass, wood, and paper depends on our fundamental understanding of tree growth and adaptation (FAO 2001; Fenning and Gershenzon 2002; Campbell et al. 2003; Gray et al. 2006). Our knowledge of the unique biology of trees will be greatly advanced through the application of genomics. The purpose of this chapter is to describe this emergent genomic paradigm as it is being applied to trees.

Gene expression during formation of earlywood and latewood in loblolly pine: expression profiles of 350 genes.

The natural variability of wood formation in trees affords opportunities to correlate transcript profiles with the resulting wood properties. We have used cDNA microarrays to study transcript abundance in developing secondary xylem of loblolly pine (Pinus taeda) over a growing season. The cDNAs were selected from a collection of 75 000 ESTs that have been sequenced and annotated (http://web.ahc.umn.edu/biodata/nsfpine/). Cell wall thickness and climatic data were related to earlywood and latewood formation at different time points during the growing season. Seventy-one ESTs showed preferential expression in earlywood or latewood, including 23 genes with no significant similarity to genes in GenBank. Seven genes involved in lignin synthesis were preferentially expressed in latewood. The studies have provided initial insights into the variation of expression patterns of some of the genes related to the wood formation process.

Comparative genetic linkage maps of Eucalyptus grandis, Eucalyptus globulus and their F1 hybrid based on a double pseudo-backcross mapping approach.

Comparative genetic mapping in interspecific pedigrees presents a powerful approach to study genetic differentiation, genome evolution and reproductive isolation in diverging species. We used this approach for genetic analysis of an F(1) hybrid of two Eucalyptus tree species, Eucalyptus grandis (W. Hill ex Maiden.) and Eucalyptus globulus (Labill.). This wide interspecific cross is characterized by hybrid inviability and hybrid abnormality. Approximately 20% of loci in the genome of the F(1) hybrid are expected to be hemizygous due to a difference in genome size between E. grandis (640 Mbp) and E. globulus (530 Mbp). We investigated the extent of colinearity between the two genomes and the distribution of hemizygous loci in the F(1) hybrid using high-throughput, semi-automated AFLP marker analysis. Two pseudo-backcross families (backcrosses of an F(1) individual to non-parental individuals of the parental species) were each genotyped with more than 800 AFLP markers. This allowed construction of de novo comparative genetic linkage maps of the F(1) hybrid and the two backcross parents. All shared AFLP marker loci in the three single-tree parental maps were found to be colinear and little evidence was found for gross chromosomal rearrangements. Our results suggest that hemizygous AFLP loci are dispersed throughout the E. grandis chromosomes of the F(1) hybrid.

Conservation and synteny of SSR loci and QTLs for vegetative propagation in four Eucalyptus species.

Conservation of microsatellite loci, heterozygous in Eucalyptus grandis, Eucalyptus urophylla, Eucalyptus tereticornis and Eucalyptus globulus, allowed us to propose homeologies among genetic linkage groups in these species, supported by at least three SSR loci in two different linkage groups. Marker-trait associations for sprouting and adventitious rooting ability were also compared in the four species. Putative quantitative trait loci (QTLs) influencing vegetative propagation traits were located on homeologous linkage groups. Our findings indicate high transferability of microsatellite markers between Eucalyptus species of the Symphyomyrtus subgenus and establish foundations for the use of synteny in the genetic analysis of this genus. Microsatellite markers should help integrate eucalypt genetic linkage maps from various sources. The availability of comparative linkage maps provides a basis of more-efficient use of genetic information for molecular breeding and evolutionary studies in Eucalyptus.

Regeneration of transgenic loblolly pine (Pinus taeda L.) from zygotic embryos transformed with Agrobacterium tumefaciens.

Embryos of 24 open-pollinated families of loblolly pine (Pinus teade L.) were used as explants to conduct in vitro regeneration. Then, Agrobacterium tumefaciens strain GV3101 harboring the plasmid pPCV6NFHygGUSINT was used to transform mature zygotic embryos of seven families of loblolly pine. The frequency of transformation varied among families infected with A. tumefaciens. The highest frequency (100%) of transient beta-glucuronidase (GUS)-expressing embryos was obtained from family 11-1029 with over 300 blue spots per embryo. Expression of the GUS reporter gene was observed in cotyledons, hypocotyls, and radicles of co-cultivated mature zygotic embryos, as well as in callus and shoots derived from co-cultivated mature zygotic embryos. Ninety transgenic plants were regenerated from hygromycin-resistant callus derived from families W03. 8-1082 and 11-1029. and 19 transgenic plantlets were established in soil. The presence of the GUS gene in the plant genome was confirmed by polymerase chain reaction. Southern blot, and plant DNA/T-DNA junction analysis. These results suggest that an efficient A. tumefaciens-mediated transformation protocol for stable integration of foreign genes into loblolly pine has been developed and that this transformation system could be useful for future studies on transferring economically important genes to loblolly pine.

Functional genomics and cell wall biosynthesis in loblolly pine.

Loblolly pine (Pinus taeda L.) is the most widely planted tree species in the USA and an important tree in commercial forestry world-wide. The large genome size and long generation time of this species present obstacles to both breeding and molecular genetic analysis. Gene discovery by partial DNA sequence determination of cDNA clones is an effective means of building a knowledge base for molecular investigations of mechanisms governing aspects of pine growth and development, including the commercially relevant properties of secondary cell walls in wood. Microarray experiments utilizing pine cDNA clones can be used to gain additional information about the potential roles of expressed genes in wood formation. Different methods have been used to analyze data from first-generation pine microarrays, with differing degrees of success. Disparities in predictions of differential gene expression between cDNA sequencing experiments and microarray experiments arise from differences in the nature of the respective analyses, but both approaches provide lists of candidate genes which should be further investigated for potential roles in cell wall formation in differentiating pine secondary xylem. Some of these genes seem to be specific to pine, while others also occur in model plants such as Arabidopsis, where they could be more efficiently investigated.

Genotypic control of high-frequency adventitious shoot regeneration via somatic organogenesis in loblolly pine.

Mature zygotic embryos of 24 genotypes of loblolly pine (Pinus taeda L.) were used as explants to establish an adventitious shoot regeneration system through somatic organogenesis. Callus formation frequencies of 18.2 (genotype 11-1103) -77.7% (genotype 7-100) have been induced from mature zygotic embryos of all genotypes tested on callus induction medium (basal salts) containing 2,4-dichlorophenoxyacetic acid (2,4-D) or alpha-naphthaleneacetic acid (NAA), 6-benzyladenine (BA), and kinetin. Adventitious shoot regeneration via organogenesis with the frequency of 5.4 (genotype 11-1103 and 7-2) -77.2% (genotype 8-1082) was obtained from callus and tissue cultures derived from mature zygotic embryos of 24 genotypes of loblolly pine. The highest mean number of 18 adventitious buds per piece of callus 0.5x0.5 cm(2) in size was obtained from genotype 8-1082. Elongation of adventitious buds was achieved on TE medium supplemented with 0.5 mg/l indole-3-butyric acid (IBA) and 1 mg/l BA. After rooting, regenerated plantlets were established in soil. These results suggested that adventitious shoot regeneration via somatic organogenesis was influenced by the genotypes. The in vitro regeneration procedure established in this investigation could be used for clonal micropropagation of some genotypes of loblolly pine, as well as for establishing a transformation system in coniferous species.

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.

High-throughput AFLP analysis using infrared dye-labeled primers and an automated DNA sequencer.

Amplified fragment length polymorphism (AFLP) analysis is currently the most powerful and efficient technique for the generation of large numbers of anonymous DNA markers in plant and animal genomes. We have developed a protocol for high-throughput AFLP analysis that allows up to 70,000 polymorphic marker genotype determinations per week on a single automated DNA sequencer. This throughput is based on multiplexed PCR amplification of AFLP fragments using two different infrared dyelabeled primer combinations. The multiplexed AFLPs are resolved on a two-dye, model 4200 LI-COR automated DNA sequencer, and the digital images are scored using semi-automated scoring software specifically designed for complex AFLP banding patterns (AFLP-Quantar). Throughput is enhanced by using high-quality genomic DNA templates obtained by a 96-well DNA isolation procedure.

Lignin structure in a mutant pine deficient in cinnamyl alcohol dehydrogenase.

Cinnamyl alcohol dehydrogenase (CAD) activity is deficient in loblolly pine (Pinus taeda L.) harboring a mutated allele of the cad gene (cad-n1). We compared lignin structure of CAD-deficient and wild-type pines, both types segregating within full-sib families obtained by controlled crosses. The type and frequency of lignin building units and distribution of interunit bonds were determined from the GC-MS analysis of thioacidolysis monomers and dimers. While the lignin content was only slightly reduced, the lignin structure was dramatically modified by the mutation in both mature and juvenile trees. Lignins from CAD-deficient pine displayed unusually high levels of coniferaldehyde and dihydroconiferyl alcohol. In addition, biphenyl and biphenyl ether bonds were in large excess in these abnormal lignins. These results suggest that the CAD-deficient pines efficiently compensate for the shortage in normal lignin precursors by utilizing nontraditional wall phenolics to construct unusual lignins particularly enriched in resistant interunit bonds.

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.

High-efficiency Agrobacterium-mediated transformation of Norway spruce (Picea abies) and loblolly pine (Pinus taeda).

Agrobacterium-mediated gene transfer is the method of choice for many plant biotechnology laboratories; however, large-scale use of this organism in conifer transformation has been limited by difficult propagation of explant material, selection efficiencies and low transformation frequency. We have analyzed co-cultivation conditions and different disarmed strains of Agrobacterium to improve transformation. Additional copies of virulence genes were added to three common disarmed strains. These extra virulence genes included either a constitutively active virG or extra copies of virG and virB, both from pTiBo542. In experiments with Norway spruce, we increased transformation efficiencies 1000-fold from initial experiments where little or no transient expression was detected. Over 100 transformed lines expressing the marker gene beta-glucuronidase (GUS) were generated from rapidly dividing embryogenic suspension-cultured cells co-cultivated with Agrobacterium. GUS activity was used to monitor transient expression and to further test lines selected on kanamycin-containing medium. In loblolly pine, transient expression increased 10-fold utilizing modified Agrobacterium strains. Agrobacterium-mediated gene transfer is a useful technique for large-scale generation of transgenic Norway spruce and may prove useful for other conifer species.

Analysis of xylem formation in pine by cDNA sequencing.

Secondary xylem (wood) formation is likely to involve some genes expressed rarely or not at all in herbaceous plants. Moreover, environmental and developmental stimuli influence secondary xylem differentiation, producing morphological and chemical changes in wood. To increase our understanding of xylem formation, and to provide material for comparative analysis of gymnosperm and angiosperm sequences, ESTs were obtained from immature xylem of loblolly pine (Pinus taeda L.). A total of 1,097 single-pass sequences were obtained from 5' ends of cDNAs made from gravistimulated tissue from bent trees. Cluster analysis detected 107 groups of similar sequences, ranging in size from 2 to 20 sequences. A total of 361 sequences fell into these groups, whereas 736 sequences were unique. About 55% of the pine EST sequences show similarity to previously described sequences in public databases. About 10% of the recognized genes encode factors involved in cell wall formation. Sequences similar to cell wall proteins, most known lignin biosynthetic enzymes, and several enzymes of carbohydrate metabolism were found. A number of putative regulatory proteins also are represented. Expression patterns of several of these genes were studied in various tissues and organs of pine. Sequencing novel genes expressed during xylem formation will provide a powerful means of identifying mechanisms controlling this important differentiation pathway.

Abnormal lignin in a loblolly pine mutant.

Novel lignin is formed in a mutant loblolly pine (Pinus taeda L.) severely depleted in cinnamyl alcohol dehydrogenase (E.C. 1.1.1.195), which converts coniferaldehyde to coniferyl alcohol, the primary lignin precursor in pines. Dihydroconiferyl alcohol, a monomer not normally associated with the lignin biosynthetic pathway, is the major component of the mutant's lignin, accounting for approximately 30 percent (versus approximately 3 percent in normal pine) of the units. The level of aldehydes, including new 2-methoxybenzaldehydes, is also increased. The mutant pines grew normally indicating that, even within a species, extensive variations in lignin composition need not disrupt the essential functions of lignin.

Inheritance, gene expression, and lignin characterization in a mutant pine deficient in cinnamyl alcohol dehydrogenase.

We have discovered a mutant loblolly pine (Pinus taeda L.) in which expression of the gene encoding cinnamyl alcohol dehydrogenase (CAD; EC 1.1.1.195) is severely reduced. The products of CAD, cinnamyl alcohols, are the precursors of lignin, a major cell wall polymer of plant vascular tissues. Lignin composition in this mutant shows dramatic modifications, including increased incorporation of the substrate of CAD (coniferaldehyde), indicating that CAD may modulate lignin composition in pine. The recessive cad-n1 allele, which causes this phenotype, was discovered in a tree heterozygous for this mutant allele. It is inherited as a simple Mendelian locus that maps to the same genomic region as the cad locus. In mutant plants, CAD activity and abundance of cad RNA transcript are low, and free CAD substrate accumulates to a high level. The wood of the mutant is brown, whereas the wood in wild types is nearly white. The wood phenotype resembles that of brown midrib (bm) mutants and some transgenic plants in which xylem is red-brown due to a reduction in CAD activity. However, unlike transgenics with reduced CAD, the pine mutant has decreased lignin content. Wood in which the composition of lignin varies beyond previous expectations still provides vascular function and mechanical support.

Genetic mapping of quantitative trait loci controlling growth and wood quality traits in Eucalyptus grandis using a maternal half-sib family and RAPD markers.

Quantitative trait loci (QTL) mapping of forest productivity traits was performed using an open pollinated half-sib family of Eucalyptus grandis. For volume growth, a sequential QTL mapping approach was applied using bulk segregant analysis (BSA), selective genotyping (SG) and cosegregation analysis (CSA). Despite the low heritability of this trait and the heterogeneous genetic background employed for mapping, BSA detected one putative QTL and SG two out of the three later found by CSA. The three putative QTL for volume growth were found to control 13.7% of the phenotypic variation, corresponding to an estimated 43.7% of the genetic variation. For wood specific gravity five QTL were identified controlling 24.7% of the phenotypic variation corresponding to 49% of the genetic variation. Overlapping QTL for CBH, WSG and percentage dry weight of bark were observed. A significant case of digenic epistasis was found, involving unlinked QTL for volume. Our results demonstrate the applicability of the within half-sib design for QTL mapping in forest trees and indicate the existence of major genes involved in the expression of economically important traits related to forest productivity in Eucalyptus grandis. These findings have important implications for marker-assisted tree breeding.

Detection of a major gene for resistance to fusiform rust disease in loblolly pine by genomic mapping.

Genomic mapping has been used to identify a region of the host genome that determines resistance to fusiform rust disease in loblolly pine where no discrete, simply inherited resistance factors had been previously found by conventional genetic analysis over four decades. A resistance locus, behaving as a single dominant gene, was mapped by association with genetic markers, even though the disease phenotype deviated from the expected Mendelian ratio. The complexity of forest pathosystems and the limitations of genetic analysis, based solely on phenotype, had led to an assumption that effective long-term disease resistance in trees should be polygenic. However, our data show that effective long-term resistance can be obtained from a single qualitative resistance gene, despite the presence of virulence in the pathogen population. Therefore, disease resistance in this endemic coevolved forest pathosystem is not exclusively polygenic. Genomic mapping now provides a powerful tool for characterizing the genetic basis of host pathogen interactions in forest trees and other undomesticated, organisms, where conventional genetic analysis often is limited or not feasible.