Background: The 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. Results: We 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. Conclusions: In 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.

Background: In today's age of genomic discovery, no attempt has been made to comprehensively sequence a gymnosperm genome. The largest genus in the coniferous family Pinaceae is Pinus, whose 110-120 species have extremely large genomes (c. 20-40 Gb, 2N = 24). The size and complexity of these genomes have prompted much speculation as to the feasibility of completing a conifer genome sequence. Conifer genomes are reputed to be highly repetitive, but there is little information available on the nature and identity of repetitive units in gymnosperms. The pines have extensive genetic resources, with approximately 329000 ESTs from eleven species and genetic maps in eight species, including a dense genetic map of the twelve linkage groups in Pinus taeda. Results: We present here the Sanger sequence and annotation of ten P. taeda BAC clones and Genome Analyzer II whole genome shotgun (WGS) sequences representing 7.5% of the genome. Computational annotation of ten BACs predicts three putative protein-coding genes and at least fifteen likely pseudogenes in nearly one megabase of sequence. We found three conifer-specific LTR retroelements in the BACs, and tentatively identified at least 15 others based on evidence from the distantly related angiosperms. Alignment of WGS sequences to the BACs indicates that 80% of BAC sequences have similar copies (≥ 75% nucleotide identity) elsewhere in the genome, but only 23% have identical copies (99% identity). The three most common repetitive elements in the genome were identified and, when combined, represent less than 5% of the genome. Conclusions: This study indicates that the majority of repeats in the P. taeda genome are 'novel' and will therefore require additional BAC or genomic sequencing for accurate characterization. The pine genome contains a very large number of diverged and probably defunct repetitive elements. This study also provides new evidence that sequencing a pine genome using a WGS approach is a feasible goal.

In developing xylem, gene expression levels vary in different genotypes, at different stages of development, throughout a growing season, and in response to stresses. Commercially important characteristics such as wood-specific gravity are known to differ with seed source. For example, when grown on a common site, the specific gravity of Arkansas loblolly pine (Pinus taeda L.) trees is greater than that of Louisiana loblolly pine, and Texas loblolly pines have a greater specific gravity than loblolly pines from the Atlantic coast. A microarray analysis was performed to examine variation in gene expression among trees from different geographical sources when grown on a common site, and seasonal variation in gene expression in each seed source. We used microarrays containing 2171 expressed sequence tags (ESTs) with putative functions of interest, selected from several loblolly pine xylem partial cDNA libraries and a shoot tip library. Genes with significant variation in expression for each factor were identified. Many genes preferentially expressed in latewood compared with earlywood were for proteins involved in cell wall biosynthesis. Variation in gene expression among trees from the two seed sources in each growing season suggests that there may be more differences between South Arkansas trees and South Louisiana trees in latewood than in earlywood. Variation in gene expression among trees from different regions may reflect adaptation to different environments.

Gene expression analyses using native populations can contribute to the understanding of plant development and adaptation in multiple ways. These include the identification of candidate genes and genetic polymorphisms affecting expression and phenotypic traits and characterization of transcriptional networks. We analyzed the expression of 111 genes with probable roles in xylem/wood development in a population of loblolly pine (Pinus taeda L.) covering much of the natural range. Loblolly pine is one of the most commercially important forest tree species in the United States, and the discovery of genes and alleles contributing to desirable wood properties would be valuable. Of the 111 genes analyzed using quantitative reverse transcription–polymerase chain reaction, there were significant differences in gene expression between clones for 106 genes. Genes encoding lignin biosynthetic enzymes and arabinogalactan proteins were more variable than those encoding cellulose synthases or those involved in signal transduction. Several groups of genes with related functions form clusters. A network analysis identified transcription factors that may be key regulators of xylem development in pine. Secondary wall-associated NAC domain protein 1 (SND1) in particular appears to be involved in the regulation of many other genes. The cluster analysis using clones did not result in discrete populations but did identify some expression differences between regions. In the future, the gene expression data will be used for association analyses and promoter studies to understand how these gene expression differences are associated with specific genetic polymorphisms in other genes and promoters.

Variation in the expression of genes with putative roles in wood development was associated with single-nucleotide polymorphisms (SNPs) using a population of loblolly pine (Pinus taeda L.) that included individuals from much of the native range. Association studies were performed using 3938 SNPs and expression data obtained using quantitative real-time polymerase chain reaction (PCR) (qRT-PCR) for 106 xylem development genes in 400 clonally replicated loblolly pine individuals. A general linear model (GLM) approach, which takes the underlying population structure into consideration, was used to discover significant associations. After adjustment for multiple testing using a false discovery rate correction, 88 statistically significant associations (Q<0.05) were observed for 80 SNPs with the expression data of 33 xylem development genes. Thirty SNPs caused nonsynonymous mutations, 18 resulted in synonymous mutations, 11 were in 3′ untranslated regions (UTRs), 1 was in a 5′ UTR and 20 were in introns. Using AraNet, we found that Arabidopsis genes with high similarity to the loblolly pine genes involved in 21 of the 88 statistically significant associations are connected in functional gene networks. Comparisons of gene expression values revealed that in most cases the average expression in plants homozygous for the rare SNP allele was lower than that of plants that were heterozygous or homozygous for the abundant allele. Although there are association studies of SNPs and expression profiles for humans, Arabidopsis and white spruce, to the best of our knowledge, this is the first example of such an association genetic study in pines. Functional validation of these associations will lead to a deeper understanding of the molecular basis of phenotypic differences in wood development among individuals in conifer populations.