W. Walter Lorenz

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.

SAGE Analysis of Transcriptome Responses in Arabidopsis Roots Exposed to 2,4,6-Trinitrotoluene

Serial analysis of gene expression was used to profile transcript levels in Arabidopsis roots and assess their responses to 2,4,6-trinitrotoluene (TNT) exposure. SAGE libraries representing control and TNT-exposed seedling root transcripts were constructed, and each was sequenced to a depth of roughly 32,000 tags. More than 19,000 unique tags were identified overall. The second most highly induced tag (27-fold increase) represented a glutathione S-transferase. Cytochrome P450 enzymes, as well as an ABC transporter and a probable nitroreductase, were highly induced by TNT exposure. Analyses also revealed an oxidative stress response upon TNT exposure. Although some increases were anticipated in light of current models for xenobiotic metabolism in plants, evidence for unsuspected conjugation pathways was also noted. Identifying transcriptome-level responses to TNT exposure will better define the metabolic pathways plants use to detoxify this xenobiotic compound, which should help improve phytoremediation strategies directed at TNT and other nitroaromatic compounds.

Microarray analysis and scale-free gene networks identify candidate regulators in drought-stressed roots of loblolly pine (P. taeda L.).

BackgroundGlobal transcriptional analysis of loblolly pine (Pinus taeda L.) is challenging due to limited molecular tools. PtGen2, a 26,496 feature cDNA microarray, was fabricated and used to assess drought-induced gene expression in loblolly pine propagule roots. Statistical analysis of differential expression and weighted gene correlation network analysis were used to identify drought-responsive genes and further characterize the molecular basis of drought tolerance in loblolly pine.ResultsMicroarrays were used to interrogate root cDNA populations obtained from 12 genotype × treatment combinations (four genotypes, three watering regimes). Comparison of drought-stressed roots with roots from the control treatment identified 2445 genes displaying at least a 1.5-fold expression difference (false discovery rate = 0.01). Genes commonly associated with drought response in pine and other plant species, as well as a number of abiotic and biotic stress-related genes, were up-regulated in drought-stressed roots. Only 76 genes were identified as differentially expressed in drought-recovered roots, indicating that the transcript population can return to the pre-drought state within 48 hours. Gene correlation analysis predicts a scale-free network topology and identifies eleven co-expression modules that ranged in size from 34 to 938 members. Network topological parameters identified a number of central nodes (hubs) including those with significant homology (E-values ≤ 2 × 10-30) to 9-cis-epoxycarotenoid dioxygenase, zeatin O-glucosyltransferase, and ABA-responsive protein. Identified hubs also include genes that have been associated previously with osmotic stress, phytohormones, enzymes that detoxify reactive oxygen species, and several genes of unknown function.ConclusionPtGen2 was used to evaluate transcriptome responses in loblolly pine and was leveraged to identify 2445 differentially expressed genes responding to severe drought stress in roots. Many of the genes identified are known to be up-regulated in response to osmotic stress in pine and other plant species and encode proteins involved in both signal transduction and stress tolerance. Gene expression levels returned to control values within a 48-hour recovery period in all but 76 transcripts. Correlation network analysis indicates a scale-free network topology for the pine root transcriptome and identifies central nodes that may serve as drivers of drought-responsive transcriptome dynamics in the roots of loblolly pine.

A SNP resource for Douglas-fir: de novo transcriptome assembly and SNP detection and validation

BackgroundDouglas-fir (Pseudotsuga menziesii), one of the most economically and ecologically important tree species in the world, also has one of the largest tree breeding programs. Although the coastal and interior varieties of Douglas-fir (vars. menziesii and glauca) are native to North America, the coastal variety is also widely planted for timber production in Europe, New Zealand, Australia, and Chile. Our main goal was to develop a SNP resource large enough to facilitate genomic selection in Douglas-fir breeding programs. To accomplish this, we developed a 454-based reference transcriptome for coastal Douglas-fir, annotated and evaluated the quality of the reference, identified putative SNPs, and then validated a sample of those SNPs using the Illumina Infinium genotyping platform.ResultsWe assembled a reference transcriptome consisting of 25,002 isogroups (unique gene models) and 102,623 singletons from 2.76 million 454 and Sanger cDNA sequences from coastal Douglas-fir. We identified 278,979 unique SNPs by mapping the 454 and Sanger sequences to the reference, and by mapping four datasets of Illumina cDNA sequences from multiple seed sources, genotypes, and tissues. The Illumina datasets represented coastal Douglas-fir (64.00 and 13.41 million reads), interior Douglas-fir (80.45 million reads), and a Yakima population similar to interior Douglas-fir (8.99 million reads). We assayed 8067 SNPs on 260 trees using an Illumina Infinium SNP genotyping array. Of these SNPs, 5847 (72.5%) were called successfully and were polymorphic.ConclusionsBased on our validation efficiency, our SNP database may contain as many as ~200,000 true SNPs, and as many as ~69,000 SNPs that could be genotyped at ~20,000 gene loci using an Infinium II array—more SNPs than are needed to use genomic selection in tree breeding programs. Ultimately, these genomic resources will enhance Douglas-fir breeding and allow us to better understand landscape-scale patterns of genetic variation and potential responses to climate change.

The phenylalanine ammonia lyase (PAL) gene family shows a gymnosperm-specific lineage

BackgroundPhenylalanine ammonia lyase (PAL) is a key enzyme of the phenylpropanoid pathway that catalyzes the deamination of phenylalanine to trans-cinnamic acid, a precursor for the lignin and flavonoid biosynthetic pathways. To date, PAL genes have been less extensively studied in gymnosperms than in angiosperms. Our interest in PAL genes stems from their potential role in the defense responses of Pinus taeda, especially with respect to lignification and production of low molecular weight phenolic compounds under various biotic and abiotic stimuli. In contrast to all angiosperms for which reference genome sequences are available, P. taeda has previously been characterized as having only a single PAL gene. Our objective was to re-evaluate this finding, assess the evolutionary history of PAL genes across major angiosperm and gymnosperm lineages, and characterize PAL gene expression patterns in Pinus taeda.MethodsWe compiled a large set of PAL genes from the largest transcript dataset available for P. taeda and other conifers. The transcript assemblies for P. taeda were validated through sequencing of PCR products amplified using gene-specific primers based on the putative PAL gene assemblies. Verified PAL gene sequences were aligned and a gene tree was estimated. The resulting gene tree was reconciled with a known species tree and the time points for gene duplication events were inferred relative to the divergence of major plant lineages.ResultsIn contrast to angiosperms, gymnosperms have retained a diverse set of PAL genes distributed among three major clades that arose from gene duplication events predating the divergence of these two seed plant lineages. Whereas multiple PAL genes have been identified in sequenced angiosperm genomes, all characterized angiosperm PAL genes form a single clade in the gene PAL tree, suggesting they are derived from a single gene in an ancestral angiosperm genome. The five distinct PAL genes detected and verified in P. taeda were derived from a combination of duplication events predating and postdating the divergence of angiosperms and gymnosperms.ConclusionsGymnosperms have a more phylogenetically diverse set of PAL genes than angiosperms. This inference has contrasting implications for the evolution of PAL gene function in gymnosperms and angiosperms.

Transcriptomic analysis highlights epigenetic and transcriptional regulation during zygotic embryo development of Pinus pinaster

BackgroundIt is during embryogenesis that the plant body plan is established and the meristems responsible for all post-embryonic growth are specified. The molecular mechanisms governing conifer embryogenesis are still largely unknown. Their elucidation may contribute valuable information to clarify if the distinct features of embryo development in angiosperms and gymnosperms result from differential gene regulation. To address this issue, we have performed the first transcriptomic analysis of zygotic embryo development in a conifer species (Pinus pinaster) focusing our study in particular on regulatory genes playing important roles during plant embryo development, namely epigenetic regulators and transcription factors.ResultsMicroarray analysis of P. pinaster zygotic embryogenesis was performed at five periods of embryo development from early developing to mature embryos. Our results show that most changes in transcript levels occurred in the first and the last embryo stage-to-stage transitions, namely early to pre-cotyledonary embryo and cotyledonary to mature embryo. An analysis of functional categories for genes that were differentially expressed through embryogenesis highlighted several epigenetic regulation mechanisms. While putative orthologs of transcripts associated with mechanisms that target transposable elements and repetitive sequences were strongly expressed in early embryogenesis, PRC2-mediated repression of genes seemed more relevant during late embryogenesis. On the other hand, functions related to sRNA pathways appeared differentially regulated across all stages of embryo development with a prevalence of miRNA functions in mid to late embryogenesis. Identification of putative transcription factor genes differentially regulated between consecutive embryo stages was strongly suggestive of the relevance of auxin responses and regulation of auxin carriers during early embryogenesis. Such responses could be involved in establishing embryo patterning. Later in development, transcripts with homology to genes acting on modulation of auxin flow and determination of adaxial-abaxial polarity were up-regulated, as were putative orthologs of genes required for meristem formation and function as well as establishment of organ boundaries. Comparative analysis with A. thaliana embryogenesis also highlighted genes involved in auxin-mediated responses, as well as epigenetic regulation, indicating highly correlated transcript profiles between the two species.ConclusionsThis is the first report of a time-course transcriptomic analysis of zygotic embryogenesis in a conifer. Taken together our results show that epigenetic regulation and transcriptional control related to auxin transport and response are critical during early to mid stages of pine embryogenesis and that important events during embryogenesis seem to be coordinated by putative orthologs of major developmental regulators in angiosperms.

Noctilisin, a Venom Glycopeptide of Sirex noctilio (Hymenoptera: Siricidae), Causes Needle Wilt and Defense Gene Responses in Pines.

ABSTRACT During oviposition, female Sirex noctilio (F.) (Siricidae) woodwasps inject their conifer hosts with a venom gland secretion. The secretion induces a variety of host physiological changes that facilitate subsequent lethal infection by a symbiotic fungus. A heat-stable factor that can migrate from the site of oviposition in the trunk through the xylem to needles in the crown of attacked pines was purified by size-fractionation and reversed-phase—high-performance liquid chromatography using activity assays based on defense gene induction as well as the needle wilt response in pine shoot explants. An 11-amino acid, posttranslationally modified peptide (SEGPROGTKRP) encoded by the most abundant transcript recovered from S. noctilio venom gland tissue comprised the backbone of the 1,850 Da active factor. Posttranslational modifications included hydroxylation of a Pro residue at position 6 as well as O-glycosylation of Ser and Thr residues at positions 1 and 8, respectively. The O-linked sugars were identical α-linked N-acetylgalactosamine residues modified at the C6 position by addition of phosphoethanolamine. In contrast to the native peptide, a synthetic version of the hydroxylated peptide backbone lacking the glycosyl side chains failed to induce pine defense genes or cause needle wilt in excised shoots. This peptide, hereafter called noctilisin, is related to the O-glycosylated short-chain proline-rich antimicrobial peptides exemplified by drosocin. The noctilisin structure contains motifs which may explain how it avoids detection by pine defense systems.

An Improved Method of RNA Isolation from Loblolly Pine (P. taeda L.) and Other Conifer Species

Tissues isolated from conifer species, particularly those belonging to the Pinaceae family, such as loblolly pine (Pinus taeda L.), contain high concentrations of phenolic compounds and polysaccharides that interfere with RNA purification. Isolation of high-quality RNA from these species requires rigorous tissue collection procedures in the field and the employment of an RNA isolation protocol comprised of multiple organic extraction steps in order to isolate RNA of sufficient quality for microarray and other genomic analyses. The isolation of high-quality RNA from field-collected loblolly pine samples can be challenging, but several modifications to standard tissue and RNA isolation procedures greatly improve results. The extent of general RNA degradation increases if samples are not properly collected and transported from the field, especially during large-scale harvests. Total RNA yields can be increased significantly by pulverizing samples in a liquid nitrogen freezer mill prior to RNA isolation, especially when samples come from woody tissues. This is primarily due to the presence of oxidizing agents, such as phenolic compounds, and polysaccharides that are both present at high levels in extracts from the woody tissues of most conifer species. If not removed, these contaminants can carry over leading to problems, such as RNA degradation, that result in low yields and a poor quality RNA sample. Carryover of phenolic compounds, as well as polysaccharides, can also reduce or even completely eliminate the activity of reverse transcriptase or other polymerases commonly used for cDNA synthesis. In particular, RNA destined to be used as template for double-stranded cDNA synthesis in the generation of cDNA libraries, single-stranded cDNA synthesis for PCR or qPCRs, or for the synthesis of microarray target materials must be of the highest quality if researchers expect to obtain optimal results. RNA isolation techniques commonly employed for many other plant species are often insufficient in their ability to remove these contaminants from conifer samples and thus do not yield total RNA samples suitable for downstream manipulations. In this video we demonstrate methods for field collection of conifer tissues, beginning with the felling of a forty year-old tree, to the harvesting of phloem, secondary xylem, and reaction wood xylem. We also demonstrate an RNA isolation protocol that has consistently yielded high-quality RNA for subsequent enzymatic manipulations.

Biomarker genes highlight intraspecific and interspecific variations in the responses of Pinus taeda L. and Pinus radiata D. Don to Sirex noctilio F. acid gland secretions

Sirex noctilio F., a Eurasian horntail woodwasp recently introduced into North America, oviposits in pines and other conifers and in the process spreads a phytopathogenic fungus that serves as a food source for its larvae. During oviposition the woodwasp also deposits mucus produced in its acid (venom) gland that alters pine defense responses and facilitates infection by the fungus. A 26,496-feature loblolly pine cDNA microarray was used to survey gene expression of pine tissue responding to S. noctilio venom. Six genes were selected for further assessment by quantitative real-time polymerase chain reaction (qRT-PCR), including one that encoded an apparent PR-4 protein and another that encoded a thaumatin-like protein. Expression of both was strongly induced in response to venom, while expression of an apparent actin gene (ACT1) was stable in response to the venom. The pattern of gene response was similar in Pinus taeda L. and Pinus radiata D. Don, but the magnitude of response in P. radiata was significantly stronger for each of the induced genes. The magnitude of the biomarker gene response to venom also varied according to genotype within these two species. The qRT-PCR assay was used to demonstrate that the primary bioactive component in S. noctilio venom is a polypeptide.

Processing the Loblolly Pine PtGen2 cDNA Microarray

PtGen2 is a 26,496 feature cDNA microarray containing amplified loblolly pine ESTs. The array is produced in our laboratory for use by researchers studying gene expression in pine and other conifer species. PtGen2 was developed as a result of our gene discovery efforts in loblolly pine, and is comprised of sequences identified primarily from root tissues, but also from needle and stem.1,2 PtGen2 has been tested by hybridizing different Cy-dye labeled conifer target cDNAs, using both amplified and non-amplified indirect labeling methods, and also tested with a number of hybridization and washing conditions. This video focuses on the handling and processing of slides before and after pre-hybridization, as well as after hybridization, using some modifications to procedures developed previously.3,4 Also included, in text form only, are the protocols used for the generation, labeling and clean up of target cDNA s, as well as information on software used for downstream data processing. PtGen2 is printed with a proprietary print buffer that contains high concentrations of salt that can be difficult to remove completely. The slides are washed first in a warm SDS solution prior to pre-hybridization. After pre-hybridization, the slides are washed vigorously in several changes of water to complete removal of remaining salts. LifterSlips™ are then cleaned and positioned on the slides and labeled cDNA is carefully loaded onto the microarray by way of capillary action which provides for even distribution of the sample across the slide, and reduces the chance of bubble incorporation. Hybridization of targets to the array is done at 48°C in high humidity conditions. After hybridization, a series of standard washes are done at 53°C and room temperature for extended times. Processing PtGen2 slides using this technique reduces salt and SDS-derived artifacts often seen when the array is processed less rigorously. Hybridizing targets derived from several different conifer RNA sources, this processing protocol yielded fewer artifacts, reduced background, and provided better consistency among different experimental groups of arrays.