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.     

Genetic mapping of expressed sequence tag polymorphism (ESTP) markers in loblolly pine (Pinus taeda L.)

The development and mapping of genetic markers based upon expressed sequence tag polymorphisms (ESTPs) in loblolly pine (Pinus taeda L.) are reported. The new markers were generated by PCR-amplification of loblolly pine genomic DNAs with primers designed from sequenced cDNAs. The cDNA libraries were constructed from RNAs expressed in the needles of loblolly pine seedlings or in the xylem from young trees. DNA polymorphisms were identified by analyzing the amplified products for differences in fragment size or restriction sites, or by examining mobility differences using denaturing gradient gel electrophoresis (DGGE). DGGE revealed more DNA polymorphisms than the other two methods. Fifty six ESTPs were mapped using either of two mapping populations and positioned onto a loblolly pine consensus genetic map. Unlike many other markers commonly used in forestry, ESTPs can be used as orthologous markers for comparative mapping, to map genes of known function, or to identify candidate genes affecting important traits in loblolly pine. Received: 10 April 2000 / Accepted: 13 July 2000     

Three loblolly pine CesA genes expressed in developing xylem are orthologous to secondary cell wall CesA genes of angiosperms

Specific plant cellulose synthases (CesA), encoded by a multigene family, are necessary for secondary wall synthesis in vascular tissues and are critical to wood production. We obtained full-length clones for the three CesAs that are highly expressed in developing xylem and examined their phylogenetic relationships and expression patterns in loblolly pine tissues. Full-length CesA clones were isolated from cDNA of developing loblolly pine (Pinus taeda) xylem and phylogenetic inferences made from plant CesA protein sequences. Expression of the three genes was examined by Northern blot analysis and semiquantitative RT-PCR. Each of three PtCesA genes is orthologous to one of the three angiosperm secondary cell wall CesAs. The PtCesAs are coexpressed in tissues of loblolly pine with tissues undergoing secondary cell wall biosynthesis showing the highest levels of expression. Phylogenetic and expression analyses suggest that functional roles for these loblolly pine CesAs are analogous to those of orthologs in angiosperm taxa. Based upon evidence from this and other studies, we suggest division of seed plant CesA genes into six major paralogous groups, each containing orthologs from various taxa. Available evidence suggests that paralogous CesA genes and their distinct functional roles evolved before the divergence of gymnosperm and angiosperm lineages.     

Cloning and characterization of cellulose synthase-like gene, PtrCSLD2 from developing xylem of aspen trees

Genetic improvement of cell wall polymer synthesis in forest trees is one of the major goals of forest biotechnology that could possibly impact their end product utilization. Identification of genes involved in cell wall polymer biogenesis is essential for achieving this goal. Among various candidate cell wall-related genes, cellulose synthase-like D (CSLD) genes are intriguing due to their hitherto unknown functions in cell wall polymer synthesis but strong structural similarity with cellulose synthases (CesAs) involved in cellulose deposition. Little is known about CSLD genes from trees. In the present article PtrCSLD2, a first CSLD gene from an economically important tree, aspen (Populus tremuloides) is reported. PtrCSLD2 cDNA was isolated from an aspen xylem cDNA library and encodes a protein that shares 90% similarity with Arabidopsis AtCSLD3 protein involved in root hair tip growth. It is possible that xylem fibers that also grow by intrusive tip growth may need expression of PtrCSLD2 for controlling the length of xylem fibers, a wood quality trait of great economical importance. PtrCSLD2 protein has a N-terminal cysteine-rich putative zinc-binding domain; eight transmembrane domains; alternating conserved and hypervariable domains; and a processive glycosyltransferases signature, D, D, D, QXXRW; all similar to aspen CesA proteins. However, PtrCSLD2 shares only 43-48% overall identity with the known aspen CesAs suggesting its distinct functional role in cell wall polymer synthesis perhaps other than cellulose biosynthesis. Based on Southern analysis, the aspen CSLD gene family consists of at least three genes and this gene copy estimate is supported by phylogenetic analysis of available CSLDs from plants. Moreover, gene expression studies using RT-PCR and in situ mRNA hybridization showed that PtrCSLD2 is expressed at a low level in all aspen tissues examined with a slightly higher expression level in secondary cell wall-enriched aspen xylem as compared to primary cell wall enriched tissues. Together, these observations suggest that PtrCSLD2 gene may be involved in the synthesis of matrix polysaccharides that are dominant in secondary cell walls of poplar xylem. Future molecular genetic analyses will clarify the functional significance of CSLD genes in the development of woody trees.     

Carbohydrate-active enzymes involved in the secondary cell wall biogenesis in hybrid aspen

Wood formation is a fundamental biological process with significant economic interest. While lignin biosynthesis is currently relatively well understood, the pathways leading to the synthesis of the key structural carbohydrates in wood fibers remain obscure. We have used a functional genomics approach to identify enzymes involved in carbohydrate biosynthesis and remodeling during xylem development in the hybrid aspen Populus tremula x tremuloides. Microarrays containing cDNA clones from different tissue-specific libraries were hybridized with probes obtained from narrow tissue sections prepared by cryosectioning of the developing xylem. Bioinformatic analyses using the sensitive tools developed for carbohydrate-active enzymes allowed the identification of 25 xylem-specific glycosyltransferases belonging to the Carbohydrate-Active EnZYme families GT2, GT8, GT14, GT31, GT43, GT47, and GT61 and nine glycosidases (or transglycosidases) belonging to the Carbohydrate-Active EnZYme families GH9, GH10, GH16, GH17, GH19, GH28, GH35, and GH51. While no genes encoding either polysaccharide lyases or carbohydrate esterases were found among the secondary wall-specific genes, one putative O-acetyltransferase was identified. These wood-specific enzyme genes constitute a valuable resource for future development of engineered fibers with improved performance in different applications.     

Analysis by molecular cloning of the human class II genes

The HLA class II genes control immune responsiveness to defined antigens; they encode cell surface heterodimers composed of alpha and beta glycopeptides. Recently, cDNA and genomic clones encoding these chains have been isolated, which allows molecular analysis of the class II genes. cDNA clones encoding the alpha chain of the HLA-DR antigen as well as that of another HLA class II antigen have been identified and characterized by nucleotide sequence analysis. These clones have been used as probes to isolate additional class II alpha cDNA clones in cDNA libraries and to identify polymorphisms in genomic DNA. Polymorphic restriction sites have been localized within the HLA-DR alpha gene and used as genetic markers in the analysis of families and of disease (insulin-dependent diabetes mellitus) and control populations. In addition, cDNA clones encoding the DR beta and DC beta chains were used as hybridization probes to identify DNA polymorphism. cDNA clones encoding the DR gamma (Ii) chain have also been identified; unlike the DR alpha and DR beta loci, the DR gamma gene is located on some chromosome other than chromosome 6. The genetic complexity of the human class II alpha and beta loci, as revealed by analysis with cDNA and genomic clones, is greater than that of the murine class II genes. The extent of that complexity will be defined by future work in this area.     

Functional mapping in pea, as an aid to the candidate gene selection and for investigating synteny with the model legume Medicago truncatula

The identification of the molecular polymorphisms giving rise to phenotypic trait variability—both quantitative and qualitative—is a major goal of the present agronomic research. Various approaches such as positional cloning or transposon tagging, as well as the candidate gene strategy have been used to discover the genes underlying this variation in plants. The construction of functional maps, i.e. composed of genes of known function, is an important component of the candidate gene approach. In the present paper we report the development of 63 single nucleotide polymorphism markers and 15 single-stranded conformation polymorphism markers for genes encoding enzymes mainly involved in primary metabolism, and their genetic mapping on a composite map using two pea recombinant inbred line populations. The complete genetic map covers 1,458 cM and comprises 363 loci, including a total of 111 gene-anchored markers: 77 gene-anchored markers described in this study, 7 microsatellites located in gene sequences, 16 flowering time genes, the Tri gene, 5 morphological markers, and 5 other genes. The mean spacing between adjacent markers is 4 cM and 90% of the markers are closer than 10 cM to their neighbours. We also report the genetic mapping of 21 of these genes in Medicago truncatula and add 41 new links between the pea and M. truncatula maps. We discuss the use of this new composite functional map for future candidate gene approaches in pea.     

Chloroplast DNA Diversity among Trees, Populations and Species in the California Closed-Cone Pines (Pinus Radiata, Pinus Muricata and Pinus Attenuata)

The amount, distribution and mutational nature of chloroplast DNA polymorphisms were studied via analysis of restriction fragment length polymorphisms in three closely related species of conifers, the California closed-cone pines-knobcone pine: Pinus attenuata Lemm.; bishop pine: Pinus muricata D. Don; and Monterey pine: Pinus radiata D. Don. Genomic DNA from 384 trees representing 19 populations were digested with 9-20 restriction enzymes and probed with cloned cpDNA fragments from Douglas-fir [Pseudotsuga menziesii (Mirb.) Franco] that comprise 82% of the chloroplast genome. Up to 313 restriction sites were surveyed, and 25 of these were observed to be polymorphic among or within species. Differences among species accounted for the majority of genetic (haplotypic) diversity observed [G(st) = 84(+/-13)%]; nucleotide diversity among species was estimated to be 0.3(+/-0.1)%. Knobcone pine and Monterey pine displayed almost no genetic variation within or among populations. Bishop pine also showed little variability within populations, but did display strong population differences [G(st) = 87(+/-8)%] that were a result of three distinct geographic groups. Mean nucleotide diversity within populations was 0.003(+/-0.002)%; intrapopulation polymorphisms were found in only five populations. This pattern of genetic variation contrasts strongly with findings from study of nuclear genes (allozymes) in the group, where most genetic diversity resides within populations rather than among populations or species. Regions of the genome subject to frequent length mutations were identified; estimates of subdivision based on length variant frequencies in one region differed strikingly from those based on site mutations or allozymes. Two trees were identified with a major chloroplast DNA inversion that closely resembled one documented between Pinus and Pseudotsuga.     

Genes induced by WDS are differentially expressed in two populations of aleppo pine (<Emphasis Type="Italic">Pinus halepensis</Emphasis>)

Gene expression in response to water-deficit stress (WDS) is a critical factor determining the survival and growth of pine seedlings. To understand how gene expression varies with different water stress levels, we differentially screened a cDNA library from roots of aleppo pine (Pinus halepensis) seedlings subjected to WDS. We found 156 clones of genes up-regulated and 56 down-regulated in response to WDS. Of the 14 clones selected for further characterization, 8 resemble WDS-responsive genes previously identified in angiosperms. The genes showing similarities to known proteins included an aldehyde dehydrogenase, a late embryogenesis abundant (LEA) protein, a chitinase, a cyclophilin, an MYB factor, an inorganic pyrophosphatase, a sucrose synthase, and a 4-coumarate ligase. Six of the clones did not have any similarity to previously identified proteins of known function. Quantitative polymerase chain reaction (qPCR) was used to compare the differential expression of these genes during control (no stress), moderate, and severe WDS treatments of seedling progeny from two different geographic origins within Israel, Yirka, and Beit Jann. The difference in expression between the treatments for various genes ranged from 1.9 to 8.0 cycle threshold. Most of the genes were expressed at similar levels in trees from the two populations or at higher levels in trees from Beit Jann, Israel. Northern blot analyses done for four highly expressed genes verify the results of the qPCR.