Search

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.

Understanding the genomic basis of local adaptation is crucial to determine the potential of long-lived woody species to withstand changes in their natural environment. In the past, efforts to dissect the genomic architecture in gymnosperms species have been limited due to the absence of reference genomes. Recently, the genomes of some commercially important conifers, such as loblolly pine, have become available, allowing whole-genome studies of these species. In this study, we test for associations between 87k SNPs, obtained from whole-genome resequencing of loblolly pine individuals, and 270 environmental variables and combinations of them. We determine the geographic location of significant loci and identify their genomic location using our newly constructed ultradense 26k SNP linkage map. We found that water availability is the main climatic variable shaping local adaptation of the species, and found 821 SNPs showing significant associations with climatic variables or combinations of them based on the consistent results of three different genotype–environment association methods. Our results suggest that adaptation to climate in the species might have occurred by many changes in the frequency of alleles with moderate to small effect sizes, and by the smaller contribution of large effect alleles in genes related to moisture deficit, temperature and precipitation. Genomic regions of low recombination and high population differentiation harbored SNPs associated with groups of environmental variables, suggesting climate adaptation might have evolved as a result of different selection pressures acting on groups of genes associated with an aspect of climate rather than on individual environmental variables. [ABSTRACT FROM AUTHOR]

Summary: Dissecting the genetic and genomic architecture of complex traits is essential to understand the forces maintaining the variation in phenotypic traits of ecological and economical importance.Whole‐genome resequencing data were used to generate high‐resolution polymorphic single nucleotide polymorphism (SNP) markers and genotype individuals from common gardens across the loblolly pine (Pinus taeda) natural range. Genome‐wide associations were tested with a large phenotypic dataset comprising 409 variables including morphological traits (height, diameter, carbon isotope discrimination, pitch canker resistance), and molecular traits such as metabolites and expression of xylem development genes.Our study identified 2335 new SNP × trait associations for the species, with many SNPs located in physical clusters in the genome of the species; and the genomic location of hotspots for metabolic × genotype associations.We found a highly polygenic basis of quantitative inheritance, with significant differences in number, effects size, genomic location and frequency of alleles contributing to variation in phenotypes in the different traits. While mutation‐selection balance might be shaping the genetic variation in metabolic traits, balancing selection is more likely to shape the variation in expression of xylem development genes. Our work contributes to the study of complex traits in nonmodel plant species by identifying associations at a whole‐genome level. See also the Commentary on this article by Casola, 221: 1669–1671. [ABSTRACT FROM AUTHOR]

The 22-gigabase genome of loblolly pine (Pinus taeda) is one of the largest ever sequenced. The draft assembly published in 2014 was built entirely from short Illumina reads, with lengths ranging from 100 to 250 base pairs (bp). The assembly was quite fragmented, containing over 11 million contigs whose weighted average (N50) size was 8206 bp. To improve this result, we generated approximately 12-fold coverage in long reads using the Single Molecule Real Time sequencing technology developed at Pacific Biosciences. We assembled the long and short reads together using the MaSuRCA mega-reads assembly algorithm, which produced a substantially better assembly, P. taeda version 2.0. The new assembly has an N50 contig size of 25 361, more than three times as large as achieved in the original assembly, and an N50 scaffold size of 107 821, 61% larger than the previous assembly. [ABSTRACT FROM AUTHOR]

Background: The use of wood as an industrial raw material has led to development of plantation forestry, in which trees are planted, managed, and harvested as crops. The productivity of such plantations often exceeds that of lessintensively- managed forests, and land managers have the option of choosing specific planting stock to produce specific types of wood for industrial use. Stem forking, or division of the stem into two or more stems of roughly equal size, is a character trait important in determining the quality of the stem for production of solid wood products. This trait typically has very low individual-tree heritability, but can be more accurately assessed in clonally-replicated plantings where each genotype is represented by several individual trees. We report results from a quantitative trait mapping experiment in a clonally-replicated full-sibling family of loblolly pine (Pinus taeda L.). Results: Quantitative trait loci influencing forking defects were identified in an outbred full-sibling family of loblolly pine, using single-nucleotide polymorphism markers. Genetic markers in this family segregated either in 1:2:1 (F2 intercross-like segregation) or 1:1 ratio (backcross-like segregation). An integrated linkage map combining markers with different segregation ratios was assembled for this full-sib family, and a total of 409 SNP markers were mapped on 12 linkage groups, covering 1622 cM. Two and three trait loci were identified for forking and ramicorn branch traits, respectively, using the interval mapping method. Three trait loci were detected for both traits using multiple-trait analysis. Conclusions: The detection of three loci for forking and ramicorn branching in a multiple-trait analysis could mean that there are genes with pleiotropic effects on both traits, or that separate genes affecting different traits are clustered together. The detection of genetic loci associated with variation in stem quality traits in this study supports the hypothesis that marker-assisted selection can be used to decrease the rate of stem defects in breeding populations of loblolly pine. [ABSTRACT FROM AUTHOR]

Genetically improving constitutive resin canal development in Pinus stems may enhance the capacity to synthesize terpenes for bark beetle resistance, chemical feedstocks, and biofuels., To discover genes that potentially regulate axial resin canal number ( RCN), single nucleotide polymorphisms ( SNPs) in 4027 genes were tested for association with RCN in two growth rings and three environments in a complex pedigree of 520 Pinus taeda individuals ( CCLONES). The map locations of associated genes were compared with RCN quantitative trait loci ( QTLs) in a ( P. taeda × Pinus elliottii) × P. elliottii pseudo-backcross of 345 full-sibs ( BC1)., Resin canal number was heritable ( h2 ˜ 0.12-0.21) and positively genetically correlated with xylem growth ( rg ˜ 0.32-0.72) and oleoresin flow ( rg ˜ 0.15-0.51). Sixteen well-supported candidate regulators of RCN were discovered in CCLONES, including genes associated across sites and ages, unidirectionally associated with oleoresin flow and xylem growth, and mapped to RCN QTLs in BC1., Breeding is predicted to increase RCN 11% in one generation and could be accelerated with genomic selection at accuracies of 0.45-0.52 across environments. There is significant genetic variation for RCN in loblolly pine, which can be exploited in breeding for elevated terpene content. [ABSTRACT FROM AUTHOR]

Quantitative genetic analyses aim to estimate genetic parameters and breeding values to select superior parents, families, and individuals. For these estimates a relationship matrix derived from the pedigree typically is used in a mixed model framework. However, breeding is a complex, multistep process and errors in the pedigree are common. Because errors reduce the accuracy of genetic parameter estimates and affect genetic gain, it is important to correct these errors. Here we show that a realized relationship matrix (RRM) derived from single nucleotide polymorphism markers based on the normality of the relationship coefficients can be used to correct pedigree errors. For a loblolly pine (Pinus taeda L.) breeding population, errors in the pedigree were detected and corrected with the RRM. With the corrected pedigree, best linear unbiased predictor (BLUP) models fit the data significantly better for 14 out of 15 traits evaluated, and the predictive ability of the genomic selection models using ridge regression BLUP increased for 13 traits. The corrected pedigree based on the normality of the relationship coefficients improves accuracy of traditional estimations of heritability and breeding values as well as genomic selection predictions. As more breeding programs begin to use genomic selection, we recommend first using the dense panel of markers to correct pedigree errors and then using the improved information to develop genomic selection prediction models. [ABSTRACT FROM AUTHOR]

Fusiform rust resistance can involve gene-for-gene interactions where resistance (Fr) genes in the host interact with corresponding avirulence genes in the pathogen, Cronartium quercuum f.sp. fusiforme (Cqf). Here, we identify trees with Fr genes in a loblolly pine population derived from a complex mating design challenged with two Cqf inocula (one gall and 10 gall mixtures). We used single nucleotide polymorphism (SNP) genotypes at sufficient density to ensure linkage between segregating markers and Fr genes identifying SNPs that explained high proportions of variance in disease incidence using BayesCπ, that also were significant using Bayesian Association with Missing Data (BAMD) software. Two SNPs mapped near Fr1 and generated significant LOD scores in single marker regression analyses for Fr1/fr1 parent 17 as well as four other parents. One SNP mapped near Fr8 and was significant for parent 28. Two SNPs mapped to linkage groups not previously shown to contain Fr genes and were significant for three parents. Parent 2 showed evidence of Fr gene stacking. Our results suggest that it is feasible to identify trees segregating for Fr genes, and to map Fr genes, based on parental analysis of SNPs that cosegregate with disease incidence in designed resistance screening trials. [ABSTRACT FROM AUTHOR]

Genetic maps have been successfully applied to assist in the dissection of complex traits, provide insight on genome structure, and estimate recombination in conjunction with physical maps. Despite an extensive list of genetic maps developed for loblolly pine ( Pinus taeda L.) over the past two decades, a high-density consensus map has not yet been constructed. In this study, we used two reference three-generation outbred pedigrees, base and qtl, obtained from the North Carolina State University Cooperative Tree Improvement Program, to obtain a high-density genetic consensus map. Both populations were genotyped with ≈ 7,000 different markers (restriction fragment length polymorphisms, expressed sequence tag polymorphisms, simple sequence repeats, SNPs). The grouping, ordering, and spacing of the markers on each linkage group were performed with JoinMap® 4.1, which implements the multipoint maximum likelihood algorithm for outbred populations. The final consensus map contains 2,466 markers, with a total length of 1,476 centimorgans (cM). The average marker density across the 12 linkage groups was 0.62 cM/marker. This high-density map provides an important resource for breeders and geneticists and will enable comparative studies across species, as well as improve the loblolly pine genome sequence assembly. [ABSTRACT FROM AUTHOR]

Despite their prevalence and importance, the genome sequences of loblolly pine, Norway spruce, and white spruce, three ecologically and economically important conifer species, are just becoming available to the research community. Following the completion of these large assemblies, annotation efforts will be undertaken to characterize the reference sequences. Accurate annotation of these ancient genomes would be aided by a comprehensive repeat library; however, few studies have generated enough sequence to fully evaluate and catalog their non-genic content. In this paper, two sets of loblolly pine genomic sequence, 103 previously assembled BACs and 90,954 newly sequenced and assembled fosmid scaffolds, were analyzed. Together, this sequence represents 280 Mbp (roughly 1% of the loblolly pine genome) and one of the most comprehensive studies of repetitive elements and genes in a gymnosperm species. A combination of homology and de novo methodologies were applied to identify both conserved and novel repeats. Similarity analysis estimated a repetitive content of 27% that included both full and partial elements. When combined with the de novo investigation, the estimate increased to almost 86%. Over 60% of the repetitive sequence consists of full or partial LTR (long terminal repeat) retrotransposons. Through de novo approaches, 6,270 novel, full-length transposable element families and 9,415 sub-families were identified. Among those 6,270 families, 82% were annotated as single-copy. Several of the novel, high-copy families are described here, with the largest, PtPiedmont, comprising 133 full-length copies. In addition to repeats, analysis of the coding region reported 23 full-length eukaryotic orthologous proteins (KOGS) and another 29 novel or orthologous genes. These discoveries, along with other genomic resources, will be used to annotate conifer genomes and address long-standing questions about gymnosperm evolution. [ABSTRACT FROM AUTHOR]

Rapidly enhancing oleoresin production in conifer stems through genomic selection and genetic engineering may increase resistance to bark beetles and terpenoid yield for liquid biofuels., We integrated association genetic and genomic prediction analyses of oleoresin flow (g 24 h−1) using 4854 single nucleotide polymorphisms ( SNPs) in expressed genes within a pedigreed population of loblolly pine ( Pinus taeda) that was clonally replicated at three sites in the southeastern United States., Additive genetic variation in oleoresin flow ( h2 ≈ 0.12-0.30) was strongly correlated between years in which precipitation varied ( ra ≈ 0.95), while the genetic correlation between sites declined from 0.8 to 0.37 with increasing differences in soil and climate among sites. A total of 231 SNPs were significantly associated with oleoresin flow, of which 81% were specific to individual sites. SNPs in sequences similar to ethylene signaling proteins, ABC transporters, and diterpenoid hydroxylases were associated with oleoresin flow across sites., Despite this complex genetic architecture, we developed a genomic prediction model to accelerate breeding for enhanced oleoresin flow that is robust to environmental variation. Results imply that breeding could increase oleoresin flow 1.5- to 2.4-fold in one generation. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Summary: •The metabolome of a plant comprises all small molecule metabolites, which are produced during cellular processes. The genetic basis for metabolites in nonmodel plants is unknown, despite frequently observed correlations between metabolite concentrations and stress responses. A quantitative genetic analysis of metabolites in a nonmodel plant species is thus warranted.•Here, we use standard association genetic methods to correlate 3563 single nucleotide polymorphisms (SNPs) to concentrations of 292 metabolites measured in a single loblolly pine (Pinus taeda) association population.•A total of 28 single locus associations were detected, representing 24 and 20 unique SNPs and metabolites, respectively. Multilocus Bayesian mixed linear models identified 2998 additional associations for a total of 1617 unique SNPs associated to 255 metabolites. These SNPs explained sizeable fractions of metabolite heritabilities when considered jointly (56.6% on average) and had lower minor allele frequencies and magnitudes of population structure as compared with random SNPs.•Modest sets of SNPs (n = 1–23) explained sizeable portions of genetic effects for many metabolites, thus highlighting the importance of multi‐SNP models to association mapping, and exhibited patterns of polymorphism consistent with being linked to targets of natural selection. The implications for association mapping in forest trees are discussed. [ABSTRACT FROM AUTHOR]

Background: Host-pathogen interactions that may lead to a competitive co-evolution of virulence and resistance mechanisms present an attractive system to study molecular evolution because strong, recent (or even current) selective pressure is expected at many genomic loci. However, it is unclear whether these selective forces would act to preserve existing diversity, promote novel diversity, or reduce linked neutral diversity during rapid fixation of advantageous alleles. In plants, the lack of adaptive immunity places a larger burden on genetic diversity to ensure survival of plant populations. This burden is even greater if the generation time of the plant is much longer than the generation time of the pathogen. Methodology/Principal Findings: Here, we present nucleotide polymorphism and substitution data for 41 candidate genes from the long-lived forest tree loblolly pine, selected primarily for their prospective influences on host-pathogen interactions. This dataset is analyzed together with 15 drought-tolerance and 13 wood-quality genes from previous studies. A wide range of neutrality tests were performed and tested against expectations from realistic demographic models. Conclusions/Significance: Collectively, our analyses found that axr (auxin response factor), caf1 (chromatin assembly factor) and gatabp1 (gata binding protein 1) candidate genes carry patterns consistent with directional selection and erd3 (early response to drought 3) displays patterns suggestive of a selective sweep, both of which are consistent with the arm-race model of disease response evolution. Furthermore, we have identified patterns consistent with diversifying selection at erf1-like (ethylene responsive factor 1), ccoaoemt (caffeoyl-CoA-O-methyltransferase), cyp450-like (cytochrome p450-like) and pr4.3 (pathogen response 4.3), expected under the trench-warfare evolution model. Finally, a drought-tolerance candidate related to the plant cell wall, lp5, displayed patterns consistent with balancing selection. In conclusion, both arms-race and trench-warfare models seem compatible with patterns of polymorphism found in different disease-response candidate genes, indicating a mixed strategy of disease tolerance evolution for loblolly pine, a major tree crop in southeastern United States. [ABSTRACT FROM AUTHOR]

Genetic variation is often arrayed in latitudinal or altitudinal clines, reflecting either adaptation along environmental gradients, migratory routes, or both. For forest trees, climate is one of the most important drivers of adaptive phenotypic traits. Correlations of single and multilocus genotypes with environmental gradients have been identified for a variety of forest trees. These correlations are interpreted normally as evidence of natural selection. Here, we use a genome-wide dataset of single nucleotide polymorphisms (SNPs) typed from 1730 loci in 682 loblolly pine ( Pinus taeda L.) trees sampled from 54 local populations covering the full-range of the species to examine allelic correlations to five multivariate measures of climate. Applications of a Bayesian generalized linear mixed model, where the climate variable was a fixed effect and an estimated variance–covariance matrix controlled random effects due to shared population history, identified several well-supported SNPs associating to principal components corresponding to geography, temperature, growing degree-days, precipitation and aridity. Functional annotation of those genes with putative orthologs in Arabidopsis revealed a diverse set of abiotic stress response genes ranging from transmembrane proteins to proteins involved in sugar metabolism. Many of these SNPs also had large allele frequency differences among populations ( FST = 0.10–0.35). These results illustrate a first step towards a ecosystem perspective of population genomics for non-model organisms, but also highlight the need for further integration of the methodologies employed in spatial statistics, population genetics and climate modeling during scans for signatures of natural selection from genomic data. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Summary Evidence for the molecular basis of a null allele of cinnamyl alcohol dehydrogenase (CAD) has been discovered in the loblolly pine (Pinus taeda L.) clone 7-56. The mutation is a two-base pair adenosine insertion located in exon 5 that causes a frame-shift which is predicted to result in premature termination of the protein. For routine detection of the mutation, a diagnostic assay was developed utilizing Template-directed Dye-terminator Incorporation and Fluorescence Polarization detection (FP-TDI). Loblolly pine is the most important commercial tree species in the USA, being harvested for pulp and solid wood products. Chemical pulping could be increased in efficiency by selecting for trees having a two-base pair adenosine insertion, by use of the rapid diagnostic assay developed in this study. [ABSTRACT FROM AUTHOR]

A consensus genetic map for Pinus taeda (loblolly pine) and Pinus elliottii (slash pine) was constructed by merging three previously published P. taeda maps with a map from a pseudo-backcross between P. elliottii and P. taeda. The consensus map positioned 3856 markers via genotyping of 1251 individuals from four pedigrees. It is the densest linkage map for a conifer to date. Average marker spacing was 0.6 cM and total map length was 2305 cM. Functional predictions of mapped genes were improved by aligning expressed sequence tags used for marker discovery to full-length P. taeda transcripts. Alignments to the P. taeda genome mapped 3305 scaffold sequences onto 12 linkage groups. The consensus genetic map was used to compare the genome-wide linkage disequilibrium in a population of distantly related P. taeda individuals (ADEPT2) used for association genetic studies and a multiplefamily pedigree used for genomic selection (CCLONES). The prevalence and extent of LD was greater in CCLONES as compared to ADEPT2; however, extended LD with LGs or between LGs was rare in both populations. The average squared correlations, r², between SNP alleles less than 1 cM apart were less than 0.05 in both populations and r² did not decay substantially with genetic distance. The consensus map and analysis of linkage disequilibrium establish a foundation for comparative association mapping and genomic selection in P. taeda and P. elliottii. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

Natural populations of forest trees exhibit striking phenotypic adaptations to diverse environmental gradients, thereby making them appealing subjects for the study of genes underlying ecologically relevant phenotypes. Here, we use a genome-wide data set of single nucleotide polymorphisms genotyped across 3059 functional genes to study patterns of population structure and identify loci associated with aridity across the natural range of loblolly pine (Pinus taeda L.). Overall patterns of population structure, as inferred using principal components and Bayesian cluster analyses, were consistent with three genetic clusters likely resulting from expansions out of Pleistocene refugia located in Mexico and Florida. A novel application of association analysis, which removes the confounding effects of shared ancestry on correlations between genetic and environmental variation, identified five loci correlated with aridity. These loci were primarily involved with abiotic stress response to temperature and drought. A unique set of 24 loci was identified as FST outliers on the basis of the genetic clusters identified previously and after accounting for expansions out of Pleistocene refugia. These loci were involved with a diversity of physiological processes. Identification of nonoverlapping sets of loci highlights the fundamental differences implicit in the use of either method and suggests a pluralistic, yet complementary, approach to the identification of genes underlying ecologically relevant phenotypes. [ABSTRACT FROM AUTHOR]

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 92 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 α-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. [ABSTRACT FROM AUTHOR]

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 (πsil = 0.00853), varying 100-fold among single genes. The level of within-gene LD was low, with an average pairwise r2 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. [ABSTRACT FROM AUTHOR]

Discusses a long-term series of experiments to map Quantitative Trait Loci (QTL) influencing wood properties traits in loblolly pine. Identification and subsequent verification of QTL in multiple genetic backgrounds, environments, and growing seasons; Importance of verification of QTL; Biological basis for observed marker-trait associations; Precise estimates of the magnitude of QTL effects; Temporal stability of QTL.

Evaluates the importance of anchored reference loci for integrating gene analysis of Loblolly pine. Examination of DNA sequences of the plant; Selection of anchored references loci; Amplification and detection of polymorphisms among conifers.