Search

Until very recently, complete characterization of the megagenomes of conifers has remained elusive. The diploid genome of sugar pine (Pinus lambertiana Dougl.) has a highly repetitive, 31 billion bp genome. It is the largest genome sequenced and assembled to date, and the first from the subgenus Strobus, or white pines, a group that is notable for having the largest genomes among the pines. The genome represents a unique opportunity to investigate genome "obesity" in conifers and white pines. Comparative analysis of P. lambertiana and P. taeda L. reveals new insights on the conservation, age, and diversity of the highly abundant transposable elements, the primary factor determining genome size. Like most North American white pines, the principal pathogen of P. lambertiana is white pine blister rust (Cronartium ribicola J.C. Fischer ex Raben.). Identification of candidate genes for resistance to this pathogen is of great ecological importance. The genome sequence afforded us the opportunity to make substantial progress on locating the major dominant gene for simple resistance hypersensitive response, Cr1 We describe new markers and gene annotation that are both tightly linked to Cr1 in a mapping population, and associated with Cr1 in unrelated sugar pine individuals sampled throughout the species' range, creating a solid foundation for future mapping. This genomic variation and annotated candidate genes characterized in our study of the Cr1 region are resources for future marker-assisted breeding efforts as well as for investigations of fundamental mechanisms of invasive disease and evolutionary response.

Conifers are the predominant gymnosperm. The size and complexity of their genomes has presented formidable technical challenges for whole-genome shotgun sequencing and assembly. We employed novel strategies that allowed us to determine the loblolly pine (Pinus taeda) reference genome sequence, the largest genome assembled to date. Most of the sequence data were derived from whole-genome shotgun sequencing of a single megagametophyte, the haploid tissue of a single pine seed. Although that constrained the quantity of available DNA, the resulting haploid sequence data were well-suited for assembly. The haploid sequence was augmented with multiple linking long-fragment mate pair libraries from the parental diploid DNA. For the longest fragments, we used novel fosmid DiTag libraries. Sequences from the linking libraries that did not match the megagametophyte were identified and removed. Assembly of the sequence data were aided by condensing the enormous number of paired-end reads into a much smaller set of longer "super-reads," rendering subsequent assembly with an overlap-based assembly algorithm computationally feasible. To further improve the contiguity and biological utility of the genome sequence, additional scaffolding methods utilizing independent genome and transcriptome assemblies were implemented. The combination of these strategies resulted in a draft genome sequence of 20.15 billion bases, with an N50 scaffold size of 66.9 kbp.

The largest genus in the conifer family Pinaceae is Pinus, with over 100 species. The size and complexity of their genomes (∼20-40 Gb, 2n = 24) have delayed the arrival of a well-annotated reference sequence. In this study, we present the annotation of the first whole-genome shotgun assembly of loblolly pine (Pinus taeda L.), which comprises 20.1 Gb of sequence. The MAKER-P annotation pipeline combined evidence-based alignments and ab initio predictions to generate 50,172 gene models, of which 15,653 are classified as high confidence. Clustering these gene models with 13 other plant species resulted in 20,646 gene families, of which 1554 are predicted to be unique to conifers. Among the conifer gene families, 159 are composed exclusively of loblolly pine members. The gene models for loblolly pine have the highest median and mean intron lengths of 24 fully sequenced plant genomes. Conifer genomes are full of repetitive DNA, with the most significant contributions from long-terminal-repeat retrotransposons. In depth analysis of the tandem and interspersed repetitive content yielded a combined estimate of 82%.

A primary goal of evolutionary genetics is to discover and explain the genetic basis of fitness-related traits and how this genetic basis evolves within natural populations. Unprecedented technological advances have fueled the discovery of genetic variants associated with ecologically relevant phenotypes in many different life forms, as well as the ability to scan genomes for deviations from selectively neutral models of evolution. Theoretically, the degree of overlap between lists of genomic regions identified using each approach is related to the genetic architecture of fitness-related traits and the strength and type of natural selection molding variation at these traits within natural populations. Here we address for the first time in a plant the degree of overlap between these lists, using patterns of nucleotide diversity and divergence for >7000 unique amplicons described from the extensive expressed sequence tag libraries generated for loblolly pine (Pinus taeda L.) in combination with the >1000 published genetic associations. We show that loci associated with phenotypic traits are distinct with regard to neutral expectations. Phenotypes measured at the whole plant level (e.g., disease resistance) exhibit an approximately twofold increase in the proportion of adaptive nonsynonymous substitutions over the genome-wide average. As expected for polygenic traits, these signals were apparent only when loci were considered at the level of functional sets. The ramifications of this result are discussed in light of the continued efforts to dissect the genetic basis of quantitative traits.

Genetic association is a powerful method for dissecting complex adaptive traits due to (i) fine-scale mapping resulting from historical recombination, (ii) wide coverage of phenotypic and genotypic variation within a single experiment, and (iii) the simultaneous discovery of loci and alleles. In this article, genetic association among single nucleotide polymorphisms (58 SNPs) from 20 wood- and drought-related candidate genes and an array of wood property traits with evolutionary and commercial importance, namely, earlywood and latewood specific gravity, percentage of latewood, earlywood microfibril angle, and wood chemistry (lignin and cellulose content), was tested using mixed linear models (MLMs) that account for relatedness among individuals by using a pairwise kinship matrix. Population structure, a common systematic bias in association studies, was assessed using 22 nuclear microsatellites. Different phenotype: genotype associations were found, some of them confirming previous evidence from collocation of QTL and genes in linkage maps (for example, 4cl and percentage of latewood) and two that involve nonsynonymous polymorphisms (cad SNP M28 with earlywood specific gravity and 4cl SNP M7 with percentage of latewood). The strongest genetic association found in this study was between allelic variation in [alpha]-tubulin, a gene involved in the formation of cortical microtubules, and earlywood microfibril angle. Intragenic LD decays rapidly in conifers; thus SNPs showing genetic association are likely to be located in close proximity to the causative polymorphisms. This first multigene association genetic study in forest trees has shown the feasibility of candidate gene strategies for dissecting complex adaptive traits, provided that genes belonging to key pathways and appropriate statistical tools are used. This approach is of particular utility in species such as conifers, where genomewide strategies are limited by their large genomes.

Genetic association studies are rapidly becoming the experimental approach of choice to dissect complex traits, including tolerance to drought stress, which is the most common cause of mortality and yield losses in forest trees. Optimization of association mapping requires knowledge of the patterns of nucleotide diversity and linkage disequilibrium and the selection of suitable polymorphisms for genotyping. Moreover, standard neutrality tests applied to DNA sequence variation data can be used to select candidate genes or amino acid sites that are putatively under selection for association mapping. In this article, we study the pattern of polymorphism of 18 candidate genes for drought-stress response in Pinus taeda L., an important tree crop. Data analyses based on a set of 21 putatively neutral nuclear microsatellites did not show population genetic structure or genomewide departures from neutrality. Candidate genes had moderate average nucleotide diversity at silent sites ([[pi].sub.sil] = 0.00853), varying 100-fold among single genes. The level of within-gene LD was low, with an average pairwise [r.sup.2] of 0.30, decaying rapidly from ~0.50 to ~0.20 at 800 bp. No apparent LD among genes was found. A selective sweep may have occurred at the early-response-to-drought-3 (erd3) gene, although population expansion can also explain our results and evidence for selection was not conclusive. One other gene, ccoaomt-1, a methylating enzyme involved in lignification, showed dimorphism (i.e., two highly divergent haplotype lineages at equal frequency), which is commonly associated with the long-term action of balancing selection. Finally, a set of haplotype-tagging SNPs (htSNPs) was selected. Using htSNPs, a reduction of genotyping effort of ~30-40%, while sampling most common allelic variants, can be gained in our ongoing association studies for drought tolerance in pine.

Nuclear sequence variation and linkage disequilibrium (LD) were studied in 15 cold-hardiness- and 3 wood quality-related candidate genes in Douglas fir [Pseudotsuga menziesii (Mirb.) Franco]. This set of genes was selected on the basis of its function in other plants and collocation with cold-hardiness-related quantitative trait loci (QTL). The single-nucleotide polymorphism (SNP) discovery panel represented 24 different trees from six regions in Washington and Oregon plus parents of a segregating population used in the QTL study. The frequency of SNPs was one SNP per 46 bp across coding and noncoding regions on average. Haplotype and nucleotide diversities were also moderately high with [H.sub.d] = 0.827 [+ or -] 0.043 and [pi] = 0.00655 [+ or -] 0.00082 on average, respectively. The nonsynonymous (replacement) nucleotide substitutions were almost five times less frequent than synonymous ones and substitutions in noncoding regions. LD decayed relatively slowly but steadily within genes. Haploblock analysis was used to define haplotype tag SNPs (htSNPs). These data will help to select SNPs for association mapping, which is already in progress.

A comparative genetic map was constructed between two important genera of the family Pinaceae. Ten homologous linkage groups in loblolly pine (Pinus taeda L.) and Douglas fir (Pseudotsuga menziesii [Mirb.] Franco) were identified using orthologous expressed sequence tag polymorphism (ESTP) and restriction fragment length polymorphism (RFLP) markers. The comparative mapping revealed extensive synteny and colinearity between genomes of the Pinaceae, consistent with the hypothesis of conservative chromosomal evolution in this important plant family. This study reports the first comparative map in forest trees at the family taxonomic level and establishes a framework for comparative genomics in Pinaceae.

Quantitative trait loci (QTL) were mapped in the woody perennial Douglas fir (Pseudotsuga menziesii var. menziesii [Mirb.] Franco) for complex traits controlling the timing of growth initiation and growth cessation. QTL were estimated under controlled environmental conditions to identify QTL interactions with photoperiod, moisture stress, winter chilling, and spring temperatures. A three-generation mapping population of 460 cloned progeny was used for genetic mapping and phenotypic evaluations. An all-marker interval mapping method was used for scanning the genome for the presence of QTL and single-factor ANOVA was used for estimating QTL-by-environment interactions. A modest number of QTL were detected per trait, with individual QTL explaining up to 9.5% of the phenotypic variation. Two QTL-by-treatment interactions were found for growth initiation, whereas several QTL-by-treatment interactions were detected among growth cessation traits. This is the first report of QTL interactions with specific environmental signals in forest trees and will assist in the identification of candidate genes controlling these important adaptive traits in perennial plants.

A long-term series of experiments to map QTL influencing wood property traits in loblolly pine has been completed. These experiments were designed to identify and subsequently verify QTL in multiple genetic backgrounds, environments, and growing seasons. Verification of QTL is necessary to substantiate a biological basis for observed marker-trait associations, to provide precise estimates of the magnitude of QTL effects, and to predict QTL expression at a given age or in a particular environment. Verification was based on the repeated detection of QTL among populations, as well as among multiple growing seasons for each population. Temporal stability of QTL was moderate, with approximately half being detected in multiple seasons. Fewer QTL were common to different populations, but the results are nonetheless encouraging for restricted applications of marker-assisted selection. QTL from larger populations accounted for less phenotypic variation than QTL detected in smaller populations, emphasizing the need for experiments employing much larger families. Additionally, 18 candidate genes related to lignin biosynthesis and cell wall structure were mapped genetically. Several candidate genes colocated with wood property QTL; however, these relationships must be verified in future experiments.

Anchored reference loci provide a framework for comparative mapping. They are landmarks to denote conserved chromosomal segments, allowing the synthesis of genetic maps from multiple sources. We evaluated 90 expressed sequence tag polymorphisms (ESTPs) from loblolly pine (Pinus taeda L.) for this function. Primer sets were assayed for amplification and polymorphism in six pedigrees, representing two subgenera of Pinus and a distant member of the Pinaceae, Douglas-fir (Pseudotsuga menziesii [Mirb.] Franco). On average, 89% of primer sets amplified in four species of subgenus Pinus, 49% in one species of subgenus Strobus, and 22% in Douglas-fir. Polymorphisms were detected for 37-61% of the ESTPs within each pedigree. Comparative mapping in loblolly and slash pine (P. elliottii Englm.) revealed that ESTPs mapped to the same location. Disrupted synteny or significant disruptions in colinearity were not detected. Thirty-five ESTPs met criteria established for anchor loci. The majority of those that did not meet these criteria were excluded when map location was known in only a single species. Anchor loci provide a unifying tool for the community, facilitating the creation of a 'generic' pine map and serving as a foundation for studies on genome organization and evolution.

A consensus map for loblolly pine (Pinus taeda L.) was constructed from the integration of linkage data from two unrelated three-generation outbred pedigrees, The progeny segregation data from restriction fragment length polymorphism, random amplified polymorphic DNA, and isozyme genetic markers from each pedigree were recoded to reflect the two independent populations of parental meioses, and genetic maps were constructed to represent each parent. The rate of meiotic recombination was significantly greater for males than females, as was the average estimate of genome length for males {1983.7 cM [Kosambi mapping function (K)]} and females [1339.5 cM (K)]. The integration of individual maps allows for the synthesis of genetic information from independent sources onto a single consensus map and facilitates the consolidation of linkage groups to represent the chromosomes (n = 12) of loblolly pine. The resulting consensus map consists of 357 unique molecular markers and covers [approximately equal to] 1300 cM(K).