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.

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.

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.

Genetic linkage maps of Eucalyptus grandis and Eucalyptus urophylla using a pseudo-testcross: mapping strategy and RAPD markers

We have used a "two-way pseudo-testcross" mapping strategy in combination with the random amplified polymorphic DNA (RAPD) assay to construct two moderate density genetic linkage maps for species of Eucalyptus. In the cross between two heterozygous individuals many single-dose RAPD markers will be heterozygous in one parent, null in the other and therefore segregate 1:1 in their F1 progeny following a testcross configuration. Meiosis and gametic segregation in each individual can be directly and efficiently analyzed using RAPD markers. We screened 305 primers of arbitrary sequence, and selected 151 to amplify a total of 558 markers. These markers were grouped at LOD 5.0, theta = 0.25, resulting in the maternal Eucalyptus grandis map having a total of 240 markers into 14 linkage groups (1552 cM) and the paternal Eucalyptus urophylla map with 251 markers in 11 linkage groups (1101 cM) (n = 11 in Eucalyptus). Framework maps ordered with a likelihood support > or = 1000:1 were assembled covering 95% of the estimated genome size in both individuals. Characterization of genome complexity of a sample of 48 mapped random amplified polymorphic DNA (RAPD) markers indicate that 53% amplify from low copy regions. These are the first reported high coverage linkage maps for any species of Eucalyptus and among the first for any hardwood tree species. We propose the combined use of RAPD markers and the pseudo-testcross configuration as a general strategy for the construction of single individual genetic linkage maps in outbred forest trees as well as in any highly heterozygous sexually reproducing living organisms. A survey of the occurrence of RAPD markers in different individuals suggests that the pseudo-testcross/RAPD mapping strategy should also be efficient at the intraspecific level and increasingly so with crosses of genetically divergent individuals. The ability to quickly construct single-tree genetic linkage maps in any forest species opens the way for a shift from the paradigm of a species index map to the heterodox proposal of constructing several maps for individual trees of a population, therefore mitigating the problem of linkage equilibrium between marker and trait loci for the application of marker assisted strategies in tree breeding.

Nucleotide sequence and evolution of the 18S ribosomal RNA gene in maize mitochondria

The nucleotide sequence of the gene coding for the 18S ribosomal RNA of maize mitochondria has been determined and a model for the secondary structure is proposed. Dot matrix analysis has been used to compare the extent and distribution of sequence similarities of the entire maize mitochondrial 18S rRNA sequence with that of 15 other small subunit rRNA sequences. The mitochondrial gene shows great similarity to the eubacterial sequences and to the maize chloroplast, and less similarity to mitochondrial rRNA genes in animals and fungi. We propose that this similarity is due to a slow rate of nucleotide divergence in plant mtDNA compared to the mtDNA of animals. Sequence comparisons indicate that the evolution of the maize mitochondrial 18S, chloroplast 16S and nuclear 17S ribosomal genes have been essentially independent, in spite of evidence for DNA transfer between organelles and the nucleus.

The characterization of a mutant affecting DNA metabolism in the development of D. melanogaster

Using a "single-fly" nucleic acid hybridization method, we have surveyed a collection of D. melanogaster strains in search of variants which affect DNA complementary to the polypyrimidine sequence corresponding to one strand of the 1.705 satellite. Hybridization of labelled polypyrimidine probe to polypurine sequence in nucleic acid extracts of single flies, followed by thermal chromatography over hydroxyapatite led to the identification of one variant. The strain Cy/M(2)S2(10) produced excess hybrid, much of which had low thermal stability. A developmental analysis of the low-melt hybrid phenotype showed that certain tissues, in particular the ovaries, were affected. In addition to the biochemical phenotype, the break down of nurse cell nuclei in Cy/M(2)S2(10) ovaries during oocyte maturation was abnormal. A genetic analysis demonstrated that both the biochemical and cytological phenotypes were the consequences of a single recessive mutation in the DNase-1 gene on chromosome III. Studies with purified DNA demonstrated that the low-melt hybrid phenotype resulted from the accumulation of low molecular weight DNA complementary to the polypyrimidine probe.

Evolution of polypyrimidines in Drosophila

We surveyed 101 different Drosophila species for the presence of a particular highly repetitive DNA sequence containing long tracts of polypyrimidine/polypurine DNA, first found in D. melanogaster. Out of 55 tested species in the melanogaster group, only the sibling species D. simulans and D. mauritiana, as well as one distant relative in the ananassae subgroup, D. varians, contained the same sequence. All four of these species have long pyrimidine tracts as shown by acid hydrolysis of labelled DNA. All four species have the same sequence, bu the amount of this polypyrimidine/polypurine DNA varies greatly. Four other species in the hydei subgroup were found to contain a polypyrimidine/polpurine sequence, with an oligonucleotide composition different from that of D. melanogaster. This polypyrimidine DNA varies from as much as 10% of the total DNA in D. nigrohydei, to as little as 0.4% in D. neohydei. The long pyrimidine tracts in the hydei subgroup are often more than a thousand nucleotides in length, representing exceedingly homogeneous repetitious sequences.--These results show a rapid but discontinuous pattern of evolution for polypyrimidine/polypurine DNA . These sequences are not species specific, yet closely related species have greatly different amounts of polypyrimidines. Drastic changes occur in the amounts of these satellite type DNA sequences, as if the sequence had no continuous selective advantage in evolution. The implications of these results with regard to the general function and evolution of satellite DNA are discussed.

Acid treatment of Drosophila deoxyribonucleic acid

In cytologic preparations of chromosomes, acid-treated deoxyribonucleic acid (DNA) is found largely in the native state. However, acid treatments widely used for chromosome preparations produce significant amounts of depurination in DNA. DNA is similarly sensitive to depurination in intact cells or as purified DNA. If treated with alkali, these apurinic gaps can be converted to single strand breaks. Acid treatment has widely different effects on specific fractions of DNA. In Drosophila melanogaster 3% of the DNA is composed of very long tracts of pyrimidines (polypyrimidines) which are resistant to acid hydrolysis. The implications of these results for molecular cytogenetics are discussed.

Polypyrimidine segments in Drosophila melanogaster DNA: II. Chromosome location and nucleotide sequence

Long pyrimidine tracts, purified from Drosophila melanogaster DNA after treatment with formic acid-diphenylamine, were used as template for E. coli RNA polymerase to produce a polynucleotide containing only purines. This polypurine RNA hybridized specifically to D. melanogaster DNA with high efficiency at low Cot values. The resulting hybrid showed high thermal stability. When polypurine RNA was subjected to complete hydrolysis with ribonuclease T1, over 90% of the nucleotide products were ApGp and ApApGp. Partial hydrolysis yielded a distinct additional component, ApApGpApGp + ApGpApApGp. We conclude that the major sequence in the polypurine transcript is (ApGpApApGp)n. In situ hybridization to salivary gland polytene chromosomes and to metaphase chromosomes from neural ganglia indicated that polypyrimidines complementary to polypurine RNA are located in heterochromatin. In femal cells, the predominant labeling was on centromeric heterochromatin of the 2nd chromosome. We have verified the location of polypyrimidines in neural ganglion cells, by using a cytological marker of chromosomes 2. In male cells, hybrid was also found on the Y chromosome.

RNA synthesis by exogenous RNA polymerase on cytological preparations of chromosomes

Cytological preparations of Drosophila polytene chromosomes serve as templates for RNA synthesis carried out by exogenous RNA polymerase (Escherichia coli). Incorporation of labeled ribonucleoside triphosphates into RNA may be observed directly by autoradiography. Because of the effects of rifampicin, actinomycin D, ribonuclease, high salt, and the requirement for all four nucleoside triphosphates, we conclude that the labeling observed over chromosomes is due to DNA-dependent RNA polymerase activity. Using this method, one can observe RNA synthesis in vitro on specific chromosome regions due to the activity of exogenous RNA polymerase. We find that much of the RNA synthesis in this system occurs on DNA sequences which appear to be in a nondenatured state.